23 ny,memind,xt,yt,level,itime1,itime2,itime,yint1,yint2,yint3)
71 integer :: nx,ny,nxmax,nymax,nzmax,memind(2),m,ix,jy,ixp,jyp
72 integer :: itime,itime1,itime2,level,indexh
73 real :: yy1(0:nxmax-1,0:nymax-1,nzmax,2)
74 real :: yy2(0:nxmax-1,0:nymax-1,nzmax,2)
75 real :: yy3(0:nxmax-1,0:nymax-1,nzmax,2)
76 real :: ddx,ddy,rddx,rddy,dt1,dt2,dt,y1(2),y2(2),y3(2)
77 real :: xt,yt,yint1,yint2,yint3,p1,p2,p3,p4
84 if (xt.ge.
real(nx-1)) xt=
real(nx-1)-0.00001
85 if (yt.ge.
real(ny-1)) yt=
real(ny-1)-0.00001
117 y1(m)=p1*yy1(ix ,jy ,level,indexh) &
118 + p2*yy1(ixp,jy ,level,indexh) &
119 + p3*yy1(ix ,jyp,level,indexh) &
120 + p4*yy1(ixp,jyp,level,indexh)
121 y2(m)=p1*yy2(ix ,jy ,level,indexh) &
122 + p2*yy2(ixp,jy ,level,indexh) &
123 + p3*yy2(ix ,jyp,level,indexh) &
124 + p4*yy2(ixp,jyp,level,indexh)
125 y3(m)=p1*yy3(ix ,jy ,level,indexh) &
126 + p2*yy3(ixp,jy ,level,indexh) &
127 + p3*yy3(ix ,jyp,level,indexh) &
128 + p4*yy3(ixp,jyp,level,indexh)
136 dt1=
real(itime-itime1)
137 dt2=
real(itime2-itime)
140 yint1=(y1(1)*dt2+y1(2)*dt1)/dt
141 yint2=(y2(1)*dt2+y2(2)*dt1)/dt
142 yint3=(y3(1)*dt2+y3(2)*dt1)/dt
subroutine interpol_rain(yy1, yy2, yy3, nxmax, nymax, nzmax, nx, ny, memind, xt, yt, level, itime1, itime2, itime, yint1, yint2, yint3)