49 real :: y1(2),y2(2),y3(2),rho1(2),rhograd1(2)
50 real :: usl,vsl,wsl,usq,vsq,wsq,xaux
52 real,
parameter :: eps=1.0e-30
76 y1(m)=p1*uun(ix ,jy ,n,indexh,ngrid) &
77 +p2*uun(ixp,jy ,n,indexh,ngrid) &
78 +p3*uun(ix ,jyp,n,indexh,ngrid) &
79 +p4*uun(ixp,jyp,n,indexh,ngrid)
80 y2(m)=p1*vvn(ix ,jy ,n,indexh,ngrid) &
81 +p2*vvn(ixp,jy ,n,indexh,ngrid) &
82 +p3*vvn(ix ,jyp,n,indexh,ngrid) &
83 +p4*vvn(ixp,jyp,n,indexh,ngrid)
84 y3(m)=p1*wwn(ix ,jy ,n,indexh,ngrid) &
85 +p2*wwn(ixp,jy ,n,indexh,ngrid) &
86 +p3*wwn(ix ,jyp,n,indexh,ngrid) &
87 +p4*wwn(ixp,jyp,n,indexh,ngrid)
88 rho1(m)=p1*rhon(ix ,jy ,n,indexh,ngrid) &
89 +p2*rhon(ixp,jy ,n,indexh,ngrid) &
90 +p3*rhon(ix ,jyp,n,indexh,ngrid) &
91 +p4*rhon(ixp,jyp,n,indexh,ngrid)
92 rhograd1(m)=p1*drhodzn(ix ,jy ,n,indexh,ngrid) &
93 +p2*drhodzn(ixp,jy ,n,indexh,ngrid) &
94 +p3*drhodzn(ix ,jyp,n,indexh,ngrid) &
95 +p4*drhodzn(ixp,jyp,n,indexh,ngrid)
97 usl=usl+uun(ix ,jy ,n,indexh,ngrid)+uun(ixp,jy ,n,indexh,ngrid) &
98 +uun(ix ,jyp,n,indexh,ngrid)+uun(ixp,jyp,n,indexh,ngrid)
99 vsl=vsl+vvn(ix ,jy ,n,indexh,ngrid)+vvn(ixp,jy ,n,indexh,ngrid) &
100 +vvn(ix ,jyp,n,indexh,ngrid)+vvn(ixp,jyp,n,indexh,ngrid)
101 wsl=wsl+wwn(ix ,jy ,n,indexh,ngrid)+wwn(ixp,jy ,n,indexh,ngrid) &
102 +wwn(ix ,jyp,n,indexh,ngrid)+wwn(ixp,jyp,n,indexh,ngrid)
104 usq=usq+uun(ix ,jy ,n,indexh,ngrid)*uun(ix ,jy ,n,indexh,ngrid)+ &
105 uun(ixp,jy ,n,indexh,ngrid)*uun(ixp,jy ,n,indexh,ngrid)+ &
106 uun(ix ,jyp,n,indexh,ngrid)*uun(ix ,jyp,n,indexh,ngrid)+ &
107 uun(ixp,jyp,n,indexh,ngrid)*uun(ixp,jyp,n,indexh,ngrid)
108 vsq=vsq+vvn(ix ,jy ,n,indexh,ngrid)*vvn(ix ,jy ,n,indexh,ngrid)+ &
109 vvn(ixp,jy ,n,indexh,ngrid)*vvn(ixp,jy ,n,indexh,ngrid)+ &
110 vvn(ix ,jyp,n,indexh,ngrid)*vvn(ix ,jyp,n,indexh,ngrid)+ &
111 vvn(ixp,jyp,n,indexh,ngrid)*vvn(ixp,jyp,n,indexh,ngrid)
112 wsq=wsq+wwn(ix ,jy ,n,indexh,ngrid)*wwn(ix ,jy ,n,indexh,ngrid)+ &
113 wwn(ixp,jy ,n,indexh,ngrid)*wwn(ixp,jy ,n,indexh,ngrid)+ &
114 wwn(ix ,jyp,n,indexh,ngrid)*wwn(ix ,jyp,n,indexh,ngrid)+ &
115 wwn(ixp,jyp,n,indexh,ngrid)*wwn(ixp,jyp,n,indexh,ngrid)
117 uprof(n)=(y1(1)*dt2+y1(2)*dt1)*dtt
118 vprof(n)=(y2(1)*dt2+y2(2)*dt1)*dtt
119 wprof(n)=(y3(1)*dt2+y3(2)*dt1)*dtt
120 rhoprof(n)=(rho1(1)*dt2+rho1(2)*dt1)*dtt
121 rhogradprof(n)=(rhograd1(1)*dt2+rhograd1(2)*dt1)*dtt
122 indzindicator(n)=.false.
128 if (xaux.lt.eps)
then
131 usigprof(n)=sqrt(xaux/7.)
135 if (xaux.lt.eps)
then
138 vsigprof(n)=sqrt(xaux/7.)
143 if (xaux.lt.eps)
then
146 wsigprof(n)=sqrt(xaux/7.)
subroutine interpol_misslev_nests(n)