23 usigold,vsigold,wsigold,xt,yt,zt,icbt)
76 integer(kind=2) :: icbt
77 real :: zt,dz,dz1,dz2,up,vp,wp,usigold,vsigold,wsigold,
ran3
78 real(kind=dp) :: xt,yt
81 integer :: idummy = -7
85 nrand=int(
ran3(idummy)*
real(maxrand-1))+1
100 h=max(hmix(ix ,jy ,1,memind(1)), &
101 hmix(ixp,jy ,1,memind(1)), &
102 hmix(ix ,jyp,1,memind(1)), &
103 hmix(ixp,jyp,1,memind(1)), &
104 hmix(ix ,jy ,1,memind(2)), &
105 hmix(ixp,jy ,1,memind(2)), &
106 hmix(ix ,jyp,1,memind(2)), &
107 hmix(ixp,jyp,1,memind(2)))
133 u=(dz1*uprof(indzp)+dz2*uprof(indz))*dz
134 v=(dz1*vprof(indzp)+dz2*vprof(indz))*dz
135 w=(dz1*wprof(indzp)+dz2*wprof(indz))*dz
153 if (nrand+2.gt.maxrand) nrand=1
154 up=rannumb(nrand)*sigu
155 vp=rannumb(nrand+1)*sigv
157 if (.not.turbswitch) wp=wp*sigw
164 ldt=int(min(tlw,h/max(2.*abs(wp*sigw),1.e-5), &
165 0.5/abs(dsigwdz),600.)*ctl)
167 ldt=int(min(tlw,h/max(2.*abs(wp),1.e-5),600.)*ctl)
172 usig=(usigprof(indzp)+usigprof(indz))/2.
173 vsig=(vsigprof(indzp)+vsigprof(indz))/2.
174 wsig=(wsigprof(indzp)+wsigprof(indz))/2.
200 if (nrand+1.gt.maxrand) nrand=1
201 up=rannumb(nrand)*0.3
202 vp=rannumb(nrand+1)*0.3
221 if (nrand+2.gt.maxrand) nrand=1
222 usigold=rannumb(nrand)*usig
223 vsigold=rannumb(nrand+1)*vsig
224 wsigold=rannumb(nrand+2)*wsig
subroutine interpol_all(itime, xt, yt, zt)
subroutine initialize(itime, ldt, up, vp, wp, usigold, vsigold, wsigold, xt, yt, zt, icbt)
subroutine interpol_wind(itime, xt, yt, zt)