55 real :: y1(2),y2(2),y3(2),rho1(2),rhograd1(2)
56 real :: usl,vsl,wsl,usq,vsq,wsq,xaux
58 real,
parameter :: eps=1.0e-30
83 y1(m)=p1*uupol(ix ,jy ,n,indexh) &
84 +p2*uupol(ixp,jy ,n,indexh) &
85 +p3*uupol(ix ,jyp,n,indexh) &
86 +p4*uupol(ixp,jyp,n,indexh)
87 y2(m)=p1*vvpol(ix ,jy ,n,indexh) &
88 +p2*vvpol(ixp,jy ,n,indexh) &
89 +p3*vvpol(ix ,jyp,n,indexh) &
90 +p4*vvpol(ixp,jyp,n,indexh)
91 usl=usl+uupol(ix ,jy ,n,indexh)+uupol(ixp,jy ,n,indexh) &
92 +uupol(ix ,jyp,n,indexh)+uupol(ixp,jyp,n,indexh)
93 vsl=vsl+vvpol(ix ,jy ,n,indexh)+vvpol(ixp,jy ,n,indexh) &
94 +vvpol(ix ,jyp,n,indexh)+vvpol(ixp,jyp,n,indexh)
96 usq=usq+uupol(ix ,jy ,n,indexh)*uupol(ix ,jy ,n,indexh)+ &
97 uupol(ixp,jy ,n,indexh)*uupol(ixp,jy ,n,indexh)+ &
98 uupol(ix ,jyp,n,indexh)*uupol(ix ,jyp,n,indexh)+ &
99 uupol(ixp,jyp,n,indexh)*uupol(ixp,jyp,n,indexh)
100 vsq=vsq+vvpol(ix ,jy ,n,indexh)*vvpol(ix ,jy ,n,indexh)+ &
101 vvpol(ixp,jy ,n,indexh)*vvpol(ixp,jy ,n,indexh)+ &
102 vvpol(ix ,jyp,n,indexh)*vvpol(ix ,jyp,n,indexh)+ &
103 vvpol(ixp,jyp,n,indexh)*vvpol(ixp,jyp,n,indexh)
105 y1(m)=p1*uu(ix ,jy ,n,indexh) &
106 +p2*uu(ixp,jy ,n,indexh) &
107 +p3*uu(ix ,jyp,n,indexh) &
108 +p4*uu(ixp,jyp,n,indexh)
109 y2(m)=p1*vv(ix ,jy ,n,indexh) &
110 +p2*vv(ixp,jy ,n,indexh) &
111 +p3*vv(ix ,jyp,n,indexh) &
112 +p4*vv(ixp,jyp,n,indexh)
113 usl=usl+uu(ix ,jy ,n,indexh)+uu(ixp,jy ,n,indexh) &
114 +uu(ix ,jyp,n,indexh)+uu(ixp,jyp,n,indexh)
115 vsl=vsl+vv(ix ,jy ,n,indexh)+vv(ixp,jy ,n,indexh) &
116 +vv(ix ,jyp,n,indexh)+vv(ixp,jyp,n,indexh)
118 usq=usq+uu(ix ,jy ,n,indexh)*uu(ix ,jy ,n,indexh)+ &
119 uu(ixp,jy ,n,indexh)*uu(ixp,jy ,n,indexh)+ &
120 uu(ix ,jyp,n,indexh)*uu(ix ,jyp,n,indexh)+ &
121 uu(ixp,jyp,n,indexh)*uu(ixp,jyp,n,indexh)
122 vsq=vsq+vv(ix ,jy ,n,indexh)*vv(ix ,jy ,n,indexh)+ &
123 vv(ixp,jy ,n,indexh)*vv(ixp,jy ,n,indexh)+ &
124 vv(ix ,jyp,n,indexh)*vv(ix ,jyp,n,indexh)+ &
125 vv(ixp,jyp,n,indexh)*vv(ixp,jyp,n,indexh)
127 y3(m)=p1*ww(ix ,jy ,n,indexh) &
128 +p2*ww(ixp,jy ,n,indexh) &
129 +p3*ww(ix ,jyp,n,indexh) &
130 +p4*ww(ixp,jyp,n,indexh)
131 rhograd1(m)=p1*drhodz(ix ,jy ,n,indexh) &
132 +p2*drhodz(ixp,jy ,n,indexh) &
133 +p3*drhodz(ix ,jyp,n,indexh) &
134 +p4*drhodz(ixp,jyp,n,indexh)
135 rho1(m)=p1*rho(ix ,jy ,n,indexh) &
136 +p2*rho(ixp,jy ,n,indexh) &
137 +p3*rho(ix ,jyp,n,indexh) &
138 +p4*rho(ixp,jyp,n,indexh)
139 wsl=wsl+ww(ix ,jy ,n,indexh)+ww(ixp,jy ,n,indexh) &
140 +ww(ix ,jyp,n,indexh)+ww(ixp,jyp,n,indexh)
141 wsq=wsq+ww(ix ,jy ,n,indexh)*ww(ix ,jy ,n,indexh)+ &
142 ww(ixp,jy ,n,indexh)*ww(ixp,jy ,n,indexh)+ &
143 ww(ix ,jyp,n,indexh)*ww(ix ,jyp,n,indexh)+ &
144 ww(ixp,jyp,n,indexh)*ww(ixp,jyp,n,indexh)
146 uprof(n)=(y1(1)*dt2+y1(2)*dt1)*dtt
147 vprof(n)=(y2(1)*dt2+y2(2)*dt1)*dtt
148 wprof(n)=(y3(1)*dt2+y3(2)*dt1)*dtt
149 rhoprof(n)=(rho1(1)*dt2+rho1(2)*dt1)*dtt
150 rhogradprof(n)=(rhograd1(1)*dt2+rhograd1(2)*dt1)*dtt
151 indzindicator(n)=.false.
158 if (xaux.lt.eps)
then
161 usigprof(n)=sqrt(xaux/7.)
165 if (xaux.lt.eps)
then
168 vsigprof(n)=sqrt(xaux/7.)
173 if (xaux.lt.eps)
then
176 wsigprof(n)=sqrt(xaux/7.)
subroutine interpol_misslev(n)