23 maxnests,ngrid,nxn,nyn,memind,xt,yt,level,itime1,itime2,itime, &
75 integer :: maxnests,ngrid
76 integer :: nxn(maxnests),nyn(maxnests),nxmaxn,nymaxn,nzmax,memind(2)
77 integer :: m,ix,jy,ixp,jyp,itime,itime1,itime2,level,indexh
78 real :: yy1(0:nxmaxn-1,0:nymaxn-1,nzmax,2,maxnests)
79 real :: yy2(0:nxmaxn-1,0:nymaxn-1,nzmax,2,maxnests)
80 real :: yy3(0:nxmaxn-1,0:nymaxn-1,nzmax,2,maxnests)
81 real :: ddx,ddy,rddx,rddy,dt1,dt2,dt,y1(2),y2(2),y3(2)
82 real :: xt,yt,yint1,yint2,yint3,p1,p2,p3,p4
89 if (xt.ge.(
real(nxn(ngrid)-1)-0.0001)) &
90 xt=
real(nxn(ngrid)-1)-0.0001
91 if (yt.ge.(
real(nyn(ngrid)-1)-0.0001)) &
92 yt=
real(nyn(ngrid)-1)-0.0001
124 y1(m)=p1*yy1(ix ,jy ,level,indexh,ngrid) &
125 + p2*yy1(ixp,jy ,level,indexh,ngrid) &
126 + p3*yy1(ix ,jyp,level,indexh,ngrid) &
127 + p4*yy1(ixp,jyp,level,indexh,ngrid)
128 y2(m)=p1*yy2(ix ,jy ,level,indexh,ngrid) &
129 + p2*yy2(ixp,jy ,level,indexh,ngrid) &
130 + p3*yy2(ix ,jyp,level,indexh,ngrid) &
131 + p4*yy2(ixp,jyp,level,indexh,ngrid)
132 y3(m)=p1*yy3(ix ,jy ,level,indexh,ngrid) &
133 + p2*yy3(ixp,jy ,level,indexh,ngrid) &
134 + p3*yy3(ix ,jyp,level,indexh,ngrid) &
135 + p4*yy3(ixp,jyp,level,indexh,ngrid)
143 dt1=
real(itime-itime1)
144 dt2=
real(itime2-itime)
147 yint1=(y1(1)*dt2+y1(2)*dt1)/dt
148 yint2=(y2(1)*dt2+y2(2)*dt1)/dt
149 yint3=(y3(1)*dt2+y3(2)*dt1)/dt
subroutine interpol_rain_nests(yy1, yy2, yy3, nxmaxn, nymaxn, nzmax, maxnests, ngrid, nxn, nyn, memind, xt, yt, level, itime1, itime2, itime, yint1, yint2, yint3)