62 real :: u1(2),v1(2),w1(2),uh(2),vh(2),wh(2)
63 real :: usl,vsl,wsl,usq,vsq,wsq,xaux
64 integer :: i,m,n,indexh,indzh
65 real,
parameter :: eps=1.0e-30
87 dt1=
real(itime-memtime(1))
88 dt2=
real(memtime(2)-itime)
95 if (height(i).gt.zt)
then
106 dz=1./(height(indz+1)-height(indz))
107 dz1=(zt-height(indz))*dz
108 dz2=(height(indz+1)-zt)*dz
130 u1(n)=p1*uun(ix ,jy ,indzh,indexh,ngrid) &
131 +p2*uun(ixp,jy ,indzh,indexh,ngrid) &
132 +p3*uun(ix ,jyp,indzh,indexh,ngrid) &
133 +p4*uun(ixp,jyp,indzh,indexh,ngrid)
134 v1(n)=p1*vvn(ix ,jy ,indzh,indexh,ngrid) &
135 +p2*vvn(ixp,jy ,indzh,indexh,ngrid) &
136 +p3*vvn(ix ,jyp,indzh,indexh,ngrid) &
137 +p4*vvn(ixp,jyp,indzh,indexh,ngrid)
138 w1(n)=p1*wwn(ix ,jy ,indzh,indexh,ngrid) &
139 +p2*wwn(ixp,jy ,indzh,indexh,ngrid) &
140 +p3*wwn(ix ,jyp,indzh,indexh,ngrid) &
141 +p4*wwn(ixp,jyp,indzh,indexh,ngrid)
143 usl=usl+uun(ix ,jy ,indzh,indexh,ngrid)+ &
144 uun(ixp,jy ,indzh,indexh,ngrid) &
145 +uun(ix ,jyp,indzh,indexh,ngrid)+ &
146 uun(ixp,jyp,indzh,indexh,ngrid)
147 vsl=vsl+vvn(ix ,jy ,indzh,indexh,ngrid)+ &
148 vvn(ixp,jy ,indzh,indexh,ngrid) &
149 +vvn(ix ,jyp,indzh,indexh,ngrid)+ &
150 vvn(ixp,jyp,indzh,indexh,ngrid)
151 wsl=wsl+wwn(ix ,jy ,indzh,indexh,ngrid)+ &
152 wwn(ixp,jy ,indzh,indexh,ngrid) &
153 +wwn(ix ,jyp,indzh,indexh,ngrid)+ &
154 wwn(ixp,jyp,indzh,indexh,ngrid)
156 usq=usq+uun(ix ,jy ,indzh,indexh,ngrid)* &
157 uun(ix ,jy ,indzh,indexh,ngrid)+ &
158 uun(ixp,jy ,indzh,indexh,ngrid)*uun(ixp,jy ,indzh,indexh,ngrid)+ &
159 uun(ix ,jyp,indzh,indexh,ngrid)*uun(ix ,jyp,indzh,indexh,ngrid)+ &
160 uun(ixp,jyp,indzh,indexh,ngrid)*uun(ixp,jyp,indzh,indexh,ngrid)
161 vsq=vsq+vvn(ix ,jy ,indzh,indexh,ngrid)* &
162 vvn(ix ,jy ,indzh,indexh,ngrid)+ &
163 vvn(ixp,jy ,indzh,indexh,ngrid)*vvn(ixp,jy ,indzh,indexh,ngrid)+ &
164 vvn(ix ,jyp,indzh,indexh,ngrid)*vvn(ix ,jyp,indzh,indexh,ngrid)+ &
165 vvn(ixp,jyp,indzh,indexh,ngrid)*vvn(ixp,jyp,indzh,indexh,ngrid)
166 wsq=wsq+wwn(ix ,jy ,indzh,indexh,ngrid)* &
167 wwn(ix ,jy ,indzh,indexh,ngrid)+ &
168 wwn(ixp,jy ,indzh,indexh,ngrid)*wwn(ixp,jy ,indzh,indexh,ngrid)+ &
169 wwn(ix ,jyp,indzh,indexh,ngrid)*wwn(ix ,jyp,indzh,indexh,ngrid)+ &
170 wwn(ixp,jyp,indzh,indexh,ngrid)*wwn(ixp,jyp,indzh,indexh,ngrid)
178 uh(m)=dz2*u1(1)+dz1*u1(2)
179 vh(m)=dz2*v1(1)+dz1*v1(2)
180 wh(m)=dz2*w1(1)+dz1*w1(2)
188 u=(uh(1)*dt2+uh(2)*dt1)*dtt
189 v=(vh(1)*dt2+vh(2)*dt1)*dtt
190 w=(wh(1)*dt2+wh(2)*dt1)*dtt
197 if (xaux.lt.eps)
then
204 if (xaux.lt.eps)
then
212 if (xaux.lt.eps)
then
subroutine interpol_wind_nests(itime, xt, yt, zt)