58 real :: u1(2),v1(2),w1(2),uh(2),vh(2),wh(2)
59 integer :: i,m,n,indexh,indzh
78 dt1=
real(itime-memtime(1))
79 dt2=
real(memtime(2)-itime)
86 if (height(i).gt.zt)
then
97 dz=1./(height(indz+1)-height(indz))
98 dz1=(zt-height(indz))*dz
99 dz2=(height(indz+1)-zt)*dz
115 u1(n)=p1*uun(ix ,jy ,indzh,indexh,ngrid) &
116 +p2*uun(ixp,jy ,indzh,indexh,ngrid) &
117 +p3*uun(ix ,jyp,indzh,indexh,ngrid) &
118 +p4*uun(ixp,jyp,indzh,indexh,ngrid)
119 v1(n)=p1*vvn(ix ,jy ,indzh,indexh,ngrid) &
120 +p2*vvn(ixp,jy ,indzh,indexh,ngrid) &
121 +p3*vvn(ix ,jyp,indzh,indexh,ngrid) &
122 +p4*vvn(ixp,jyp,indzh,indexh,ngrid)
123 w1(n)=p1*wwn(ix ,jy ,indzh,indexh,ngrid) &
124 +p2*wwn(ixp,jy ,indzh,indexh,ngrid) &
125 +p3*wwn(ix ,jyp,indzh,indexh,ngrid) &
126 +p4*wwn(ixp,jyp,indzh,indexh,ngrid)
135 uh(m)=dz2*u1(1)+dz1*u1(2)
136 vh(m)=dz2*v1(1)+dz1*v1(2)
137 wh(m)=dz2*w1(1)+dz1*w1(2)
145 u=(uh(1)*dt2+uh(2)*dt1)*dtt
146 v=(vh(1)*dt2+vh(2)*dt1)*dtt
147 w=(wh(1)*dt2+wh(2)*dt1)*dtt
subroutine interpol_wind_short_nests(itime, xt, yt, zt)