61 real :: ust1(2),wst1(2),oli1(2),oliaux
62 real :: y1(2),y2(2),y3(2),rho1(2),rhograd1(2)
63 real :: usl,vsl,wsl,usq,vsq,wsq,xaux
64 integer :: i,m,n,indexh
65 real,
parameter :: eps=1.0e-30
87 dt1=
real(itime-memtime(1))
88 dt2=
real(memtime(2)-itime)
101 ust1(m)=p1*ustar(ix ,jy ,1,indexh) &
102 + p2*ustar(ixp,jy ,1,indexh) &
103 + p3*ustar(ix ,jyp,1,indexh) &
104 + p4*ustar(ixp,jyp,1,indexh)
105 wst1(m)=p1*wstar(ix ,jy ,1,indexh) &
106 + p2*wstar(ixp,jy ,1,indexh) &
107 + p3*wstar(ix ,jyp,1,indexh) &
108 + p4*wstar(ixp,jyp,1,indexh)
109 oli1(m)=p1*oli(ix ,jy ,1,indexh) &
110 + p2*oli(ixp,jy ,1,indexh) &
111 + p3*oli(ix ,jyp,1,indexh) &
112 + p4*oli(ixp,jyp,1,indexh)
117 ust=(ust1(1)*dt2+ust1(2)*dt1)*dtt
118 wst=(wst1(1)*dt2+wst1(2)*dt1)*dtt
119 oliaux=(oli1(1)*dt2+oli1(2)*dt1)*dtt
121 if (oliaux.ne.0.)
then
137 if (height(i).gt.zt)
then
163 y1(m)=p1*uupol(ix ,jy ,n,indexh) &
164 +p2*uupol(ixp,jy ,n,indexh) &
165 +p3*uupol(ix ,jyp,n,indexh) &
166 +p4*uupol(ixp,jyp,n,indexh)
167 y2(m)=p1*vvpol(ix ,jy ,n,indexh) &
168 +p2*vvpol(ixp,jy ,n,indexh) &
169 +p3*vvpol(ix ,jyp,n,indexh) &
170 +p4*vvpol(ixp,jyp,n,indexh)
171 usl=usl+uupol(ix ,jy ,n,indexh)+uupol(ixp,jy ,n,indexh) &
172 +uupol(ix ,jyp,n,indexh)+uupol(ixp,jyp,n,indexh)
173 vsl=vsl+vvpol(ix ,jy ,n,indexh)+vvpol(ixp,jy ,n,indexh) &
174 +vvpol(ix ,jyp,n,indexh)+vvpol(ixp,jyp,n,indexh)
176 usq=usq+uupol(ix ,jy ,n,indexh)*uupol(ix ,jy ,n,indexh)+ &
177 uupol(ixp,jy ,n,indexh)*uupol(ixp,jy ,n,indexh)+ &
178 uupol(ix ,jyp,n,indexh)*uupol(ix ,jyp,n,indexh)+ &
179 uupol(ixp,jyp,n,indexh)*uupol(ixp,jyp,n,indexh)
180 vsq=vsq+vvpol(ix ,jy ,n,indexh)*vvpol(ix ,jy ,n,indexh)+ &
181 vvpol(ixp,jy ,n,indexh)*vvpol(ixp,jy ,n,indexh)+ &
182 vvpol(ix ,jyp,n,indexh)*vvpol(ix ,jyp,n,indexh)+ &
183 vvpol(ixp,jyp,n,indexh)*vvpol(ixp,jyp,n,indexh)
185 y1(m)=p1*uu(ix ,jy ,n,indexh) &
186 +p2*uu(ixp,jy ,n,indexh) &
187 +p3*uu(ix ,jyp,n,indexh) &
188 +p4*uu(ixp,jyp,n,indexh)
189 y2(m)=p1*vv(ix ,jy ,n,indexh) &
190 +p2*vv(ixp,jy ,n,indexh) &
191 +p3*vv(ix ,jyp,n,indexh) &
192 +p4*vv(ixp,jyp,n,indexh)
193 usl=usl+uu(ix ,jy ,n,indexh)+uu(ixp,jy ,n,indexh) &
194 +uu(ix ,jyp,n,indexh)+uu(ixp,jyp,n,indexh)
195 vsl=vsl+vv(ix ,jy ,n,indexh)+vv(ixp,jy ,n,indexh) &
196 +vv(ix ,jyp,n,indexh)+vv(ixp,jyp,n,indexh)
198 usq=usq+uu(ix ,jy ,n,indexh)*uu(ix ,jy ,n,indexh)+ &
199 uu(ixp,jy ,n,indexh)*uu(ixp,jy ,n,indexh)+ &
200 uu(ix ,jyp,n,indexh)*uu(ix ,jyp,n,indexh)+ &
201 uu(ixp,jyp,n,indexh)*uu(ixp,jyp,n,indexh)
202 vsq=vsq+vv(ix ,jy ,n,indexh)*vv(ix ,jy ,n,indexh)+ &
203 vv(ixp,jy ,n,indexh)*vv(ixp,jy ,n,indexh)+ &
204 vv(ix ,jyp,n,indexh)*vv(ix ,jyp,n,indexh)+ &
205 vv(ixp,jyp,n,indexh)*vv(ixp,jyp,n,indexh)
207 y3(m)=p1*ww(ix ,jy ,n,indexh) &
208 +p2*ww(ixp,jy ,n,indexh) &
209 +p3*ww(ix ,jyp,n,indexh) &
210 +p4*ww(ixp,jyp,n,indexh)
211 rhograd1(m)=p1*drhodz(ix ,jy ,n,indexh) &
212 +p2*drhodz(ixp,jy ,n,indexh) &
213 +p3*drhodz(ix ,jyp,n,indexh) &
214 +p4*drhodz(ixp,jyp,n,indexh)
215 rho1(m)=p1*rho(ix ,jy ,n,indexh) &
216 +p2*rho(ixp,jy ,n,indexh) &
217 +p3*rho(ix ,jyp,n,indexh) &
218 +p4*rho(ixp,jyp,n,indexh)
219 wsl=wsl+ww(ix ,jy ,n,indexh)+ww(ixp,jy ,n,indexh) &
220 +ww(ix ,jyp,n,indexh)+ww(ixp,jyp,n,indexh)
221 wsq=wsq+ww(ix ,jy ,n,indexh)*ww(ix ,jy ,n,indexh)+ &
222 ww(ixp,jy ,n,indexh)*ww(ixp,jy ,n,indexh)+ &
223 ww(ix ,jyp,n,indexh)*ww(ix ,jyp,n,indexh)+ &
224 ww(ixp,jyp,n,indexh)*ww(ixp,jyp,n,indexh)
226 uprof(n)=(y1(1)*dt2+y1(2)*dt1)*dtt
227 vprof(n)=(y2(1)*dt2+y2(2)*dt1)*dtt
228 wprof(n)=(y3(1)*dt2+y3(2)*dt1)*dtt
229 rhoprof(n)=(rho1(1)*dt2+rho1(2)*dt1)*dtt
230 rhogradprof(n)=(rhograd1(1)*dt2+rhograd1(2)*dt1)*dtt
231 indzindicator(n)=.false.
237 if (xaux.lt.eps)
then
240 usigprof(n)=sqrt(xaux/7.)
244 if (xaux.lt.eps)
then
247 vsigprof(n)=sqrt(xaux/7.)
252 if (xaux.lt.eps)
then
255 wsigprof(n)=sqrt(xaux/7.)
subroutine interpol_all(itime, xt, yt, zt)