61 real :: u1(2),v1(2),w1(2),uh(2),vh(2),wh(2)
62 real :: usl,vsl,wsl,usq,vsq,wsq,xaux
63 integer :: i,m,n,indexh,indzh
64 real,
parameter :: eps=1.0e-30
86 dt1=
real(itime-memtime(1))
87 dt2=
real(memtime(2)-itime)
94 if (height(i).gt.zt)
then
105 dz=1./(height(indz+1)-height(indz))
106 dz1=(zt-height(indz))*dz
107 dz2=(height(indz+1)-zt)*dz
129 u1(n)=p1*uupol(ix ,jy ,indzh,indexh) &
130 +p2*uupol(ixp,jy ,indzh,indexh) &
131 +p3*uupol(ix ,jyp,indzh,indexh) &
132 +p4*uupol(ixp,jyp,indzh,indexh)
133 v1(n)=p1*vvpol(ix ,jy ,indzh,indexh) &
134 +p2*vvpol(ixp,jy ,indzh,indexh) &
135 +p3*vvpol(ix ,jyp,indzh,indexh) &
136 +p4*vvpol(ixp,jyp,indzh,indexh)
137 usl=usl+uupol(ix ,jy ,indzh,indexh)+ &
138 uupol(ixp,jy ,indzh,indexh) &
139 +uupol(ix ,jyp,indzh,indexh)+uupol(ixp,jyp,indzh,indexh)
140 vsl=vsl+vvpol(ix ,jy ,indzh,indexh)+ &
141 vvpol(ixp,jy ,indzh,indexh) &
142 +vvpol(ix ,jyp,indzh,indexh)+vvpol(ixp,jyp,indzh,indexh)
144 usq=usq+uupol(ix ,jy ,indzh,indexh)* &
145 uupol(ix ,jy ,indzh,indexh)+ &
146 uupol(ixp,jy ,indzh,indexh)*uupol(ixp,jy ,indzh,indexh)+ &
147 uupol(ix ,jyp,indzh,indexh)*uupol(ix ,jyp,indzh,indexh)+ &
148 uupol(ixp,jyp,indzh,indexh)*uupol(ixp,jyp,indzh,indexh)
149 vsq=vsq+vvpol(ix ,jy ,indzh,indexh)* &
150 vvpol(ix ,jy ,indzh,indexh)+ &
151 vvpol(ixp,jy ,indzh,indexh)*vvpol(ixp,jy ,indzh,indexh)+ &
152 vvpol(ix ,jyp,indzh,indexh)*vvpol(ix ,jyp,indzh,indexh)+ &
153 vvpol(ixp,jyp,indzh,indexh)*vvpol(ixp,jyp,indzh,indexh)
155 u1(n)=p1*uu(ix ,jy ,indzh,indexh) &
156 +p2*uu(ixp,jy ,indzh,indexh) &
157 +p3*uu(ix ,jyp,indzh,indexh) &
158 +p4*uu(ixp,jyp,indzh,indexh)
159 v1(n)=p1*vv(ix ,jy ,indzh,indexh) &
160 +p2*vv(ixp,jy ,indzh,indexh) &
161 +p3*vv(ix ,jyp,indzh,indexh) &
162 +p4*vv(ixp,jyp,indzh,indexh)
163 usl=usl+uu(ix ,jy ,indzh,indexh)+uu(ixp,jy ,indzh,indexh) &
164 +uu(ix ,jyp,indzh,indexh)+uu(ixp,jyp,indzh,indexh)
165 vsl=vsl+vv(ix ,jy ,indzh,indexh)+vv(ixp,jy ,indzh,indexh) &
166 +vv(ix ,jyp,indzh,indexh)+vv(ixp,jyp,indzh,indexh)
168 usq=usq+uu(ix ,jy ,indzh,indexh)*uu(ix ,jy ,indzh,indexh)+ &
169 uu(ixp,jy ,indzh,indexh)*uu(ixp,jy ,indzh,indexh)+ &
170 uu(ix ,jyp,indzh,indexh)*uu(ix ,jyp,indzh,indexh)+ &
171 uu(ixp,jyp,indzh,indexh)*uu(ixp,jyp,indzh,indexh)
172 vsq=vsq+vv(ix ,jy ,indzh,indexh)*vv(ix ,jy ,indzh,indexh)+ &
173 vv(ixp,jy ,indzh,indexh)*vv(ixp,jy ,indzh,indexh)+ &
174 vv(ix ,jyp,indzh,indexh)*vv(ix ,jyp,indzh,indexh)+ &
175 vv(ixp,jyp,indzh,indexh)*vv(ixp,jyp,indzh,indexh)
177 w1(n)=p1*ww(ix ,jy ,indzh,indexh) &
178 +p2*ww(ixp,jy ,indzh,indexh) &
179 +p3*ww(ix ,jyp,indzh,indexh) &
180 +p4*ww(ixp,jyp,indzh,indexh)
181 wsl=wsl+ww(ix ,jy ,indzh,indexh)+ww(ixp,jy ,indzh,indexh) &
182 +ww(ix ,jyp,indzh,indexh)+ww(ixp,jyp,indzh,indexh)
183 wsq=wsq+ww(ix ,jy ,indzh,indexh)*ww(ix ,jy ,indzh,indexh)+ &
184 ww(ixp,jy ,indzh,indexh)*ww(ixp,jy ,indzh,indexh)+ &
185 ww(ix ,jyp,indzh,indexh)*ww(ix ,jyp,indzh,indexh)+ &
186 ww(ixp,jyp,indzh,indexh)*ww(ixp,jyp,indzh,indexh)
194 uh(m)=dz2*u1(1)+dz1*u1(2)
195 vh(m)=dz2*v1(1)+dz1*v1(2)
196 wh(m)=dz2*w1(1)+dz1*w1(2)
204 u=(uh(1)*dt2+uh(2)*dt1)*dtt
205 v=(vh(1)*dt2+vh(2)*dt1)*dtt
206 w=(wh(1)*dt2+wh(2)*dt1)*dtt
213 if (xaux.lt.eps)
then
220 if (xaux.lt.eps)
then
228 if (xaux.lt.eps)
then
subroutine interpol_wind(itime, xt, yt, zt)