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
116 u1(n)=p1*uupol(ix ,jy ,indzh,indexh) &
117 +p2*uupol(ixp,jy ,indzh,indexh) &
118 +p3*uupol(ix ,jyp,indzh,indexh) &
119 +p4*uupol(ixp,jyp,indzh,indexh)
120 v1(n)=p1*vvpol(ix ,jy ,indzh,indexh) &
121 +p2*vvpol(ixp,jy ,indzh,indexh) &
122 +p3*vvpol(ix ,jyp,indzh,indexh) &
123 +p4*vvpol(ixp,jyp,indzh,indexh)
125 u1(n)=p1*uu(ix ,jy ,indzh,indexh) &
126 +p2*uu(ixp,jy ,indzh,indexh) &
127 +p3*uu(ix ,jyp,indzh,indexh) &
128 +p4*uu(ixp,jyp,indzh,indexh)
129 v1(n)=p1*vv(ix ,jy ,indzh,indexh) &
130 +p2*vv(ixp,jy ,indzh,indexh) &
131 +p3*vv(ix ,jyp,indzh,indexh) &
132 +p4*vv(ixp,jyp,indzh,indexh)
134 w1(n)=p1*ww(ix ,jy ,indzh,indexh) &
135 +p2*ww(ixp,jy ,indzh,indexh) &
136 +p3*ww(ix ,jyp,indzh,indexh) &
137 +p4*ww(ixp,jyp,indzh,indexh)
145 uh(m)=dz2*u1(1)+dz1*u1(2)
146 vh(m)=dz2*v1(1)+dz1*v1(2)
147 wh(m)=dz2*w1(1)+dz1*w1(2)
156 u=(uh(1)*dt2+uh(2)*dt1)*dtt
157 v=(vh(1)*dt2+vh(2)*dt1)*dtt
158 w=(wh(1)*dt2+wh(2)*dt1)*dtt
subroutine interpol_wind_short(itime, xt, yt, zt)