63 real :: ust1(2),wst1(2),oli1(2),oliaux
64 real :: y1(2),y2(2),y3(2),rho1(2),rhograd1(2)
65 real :: usl,vsl,wsl,usq,vsq,wsq,xaux
66 integer :: i,m,n,indexh
67 real,
parameter :: eps=1.0e-30
89 dt1=
real(itime-memtime(1))
90 dt2=
real(memtime(2)-itime)
103 ust1(m)=p1*ustarn(ix ,jy ,1,indexh,ngrid) &
104 + p2*ustarn(ixp,jy ,1,indexh,ngrid) &
105 + p3*ustarn(ix ,jyp,1,indexh,ngrid) &
106 + p4*ustarn(ixp,jyp,1,indexh,ngrid)
107 wst1(m)=p1*wstarn(ix ,jy ,1,indexh,ngrid) &
108 + p2*wstarn(ixp,jy ,1,indexh,ngrid) &
109 + p3*wstarn(ix ,jyp,1,indexh,ngrid) &
110 + p4*wstarn(ixp,jyp,1,indexh,ngrid)
111 oli1(m)=p1*olin(ix ,jy ,1,indexh,ngrid) &
112 + p2*olin(ixp,jy ,1,indexh,ngrid) &
113 + p3*olin(ix ,jyp,1,indexh,ngrid) &
114 + p4*olin(ixp,jyp,1,indexh,ngrid)
119 ust=(ust1(1)*dt2+ust1(2)*dt1)*dtt
120 wst=(wst1(1)*dt2+wst1(2)*dt1)*dtt
121 oliaux=(oli1(1)*dt2+oli1(2)*dt1)*dtt
123 if (oliaux.ne.0.)
then
139 if (height(i).gt.zt)
then
164 y1(m)=p1*uun(ix ,jy ,n,indexh,ngrid) &
165 +p2*uun(ixp,jy ,n,indexh,ngrid) &
166 +p3*uun(ix ,jyp,n,indexh,ngrid) &
167 +p4*uun(ixp,jyp,n,indexh,ngrid)
168 y2(m)=p1*vvn(ix ,jy ,n,indexh,ngrid) &
169 +p2*vvn(ixp,jy ,n,indexh,ngrid) &
170 +p3*vvn(ix ,jyp,n,indexh,ngrid) &
171 +p4*vvn(ixp,jyp,n,indexh,ngrid)
172 y3(m)=p1*wwn(ix ,jy ,n,indexh,ngrid) &
173 +p2*wwn(ixp,jy ,n,indexh,ngrid) &
174 +p3*wwn(ix ,jyp,n,indexh,ngrid) &
175 +p4*wwn(ixp,jyp,n,indexh,ngrid)
176 rhograd1(m)=p1*drhodzn(ix ,jy ,n,indexh,ngrid) &
177 +p2*drhodzn(ixp,jy ,n,indexh,ngrid) &
178 +p3*drhodzn(ix ,jyp,n,indexh,ngrid) &
179 +p4*drhodzn(ixp,jyp,n,indexh,ngrid)
180 rho1(m)=p1*rhon(ix ,jy ,n,indexh,ngrid) &
181 +p2*rhon(ixp,jy ,n,indexh,ngrid) &
182 +p3*rhon(ix ,jyp,n,indexh,ngrid) &
183 +p4*rhon(ixp,jyp,n,indexh,ngrid)
185 usl=usl+uun(ix ,jy ,n,indexh,ngrid)+uun(ixp,jy ,n,indexh,ngrid) &
186 +uun(ix ,jyp,n,indexh,ngrid)+uun(ixp,jyp,n,indexh,ngrid)
187 vsl=vsl+vvn(ix ,jy ,n,indexh,ngrid)+vvn(ixp,jy ,n,indexh,ngrid) &
188 +vvn(ix ,jyp,n,indexh,ngrid)+vvn(ixp,jyp,n,indexh,ngrid)
189 wsl=wsl+wwn(ix ,jy ,n,indexh,ngrid)+wwn(ixp,jy ,n,indexh,ngrid) &
190 +wwn(ix ,jyp,n,indexh,ngrid)+wwn(ixp,jyp,n,indexh,ngrid)
192 usq=usq+uun(ix ,jy ,n,indexh,ngrid)*uun(ix ,jy ,n,indexh,ngrid)+ &
193 uun(ixp,jy ,n,indexh,ngrid)*uun(ixp,jy ,n,indexh,ngrid)+ &
194 uun(ix ,jyp,n,indexh,ngrid)*uun(ix ,jyp,n,indexh,ngrid)+ &
195 uun(ixp,jyp,n,indexh,ngrid)*uun(ixp,jyp,n,indexh,ngrid)
196 vsq=vsq+vvn(ix ,jy ,n,indexh,ngrid)*vvn(ix ,jy ,n,indexh,ngrid)+ &
197 vvn(ixp,jy ,n,indexh,ngrid)*vvn(ixp,jy ,n,indexh,ngrid)+ &
198 vvn(ix ,jyp,n,indexh,ngrid)*vvn(ix ,jyp,n,indexh,ngrid)+ &
199 vvn(ixp,jyp,n,indexh,ngrid)*vvn(ixp,jyp,n,indexh,ngrid)
200 wsq=wsq+wwn(ix ,jy ,n,indexh,ngrid)*wwn(ix ,jy ,n,indexh,ngrid)+ &
201 wwn(ixp,jy ,n,indexh,ngrid)*wwn(ixp,jy ,n,indexh,ngrid)+ &
202 wwn(ix ,jyp,n,indexh,ngrid)*wwn(ix ,jyp,n,indexh,ngrid)+ &
203 wwn(ixp,jyp,n,indexh,ngrid)*wwn(ixp,jyp,n,indexh,ngrid)
205 uprof(n)=(y1(1)*dt2+y1(2)*dt1)*dtt
206 vprof(n)=(y2(1)*dt2+y2(2)*dt1)*dtt
207 wprof(n)=(y3(1)*dt2+y3(2)*dt1)*dtt
208 rhoprof(n)=(rho1(1)*dt2+rho1(2)*dt1)*dtt
209 rhogradprof(n)=(rhograd1(1)*dt2+rhograd1(2)*dt1)*dtt
210 indzindicator(n)=.false.
216 if (xaux.lt.eps)
then
219 usigprof(n)=sqrt(xaux/7.)
223 if (xaux.lt.eps)
then
226 vsigprof(n)=sqrt(xaux/7.)
231 if (xaux.lt.eps)
then
234 wsigprof(n)=sqrt(xaux/7.)
subroutine interpol_all_nests(itime, xt, yt, zt)