60 if (h/abs(ol).lt.1.)
then
63 sigu=1.e-2+2.0*ust*exp(-3.e-4*corr)
64 sigw=1.3*ust*exp(-2.e-4*corr)
68 tlu=0.5*z/sigw/(1.+1.5e-3*corr)
77 else if (ol.lt.0.)
then
83 sigu=1.e-2+ust*(12-0.5*h/ol)**0.33333
85 sigw=sqrt(1.2*wst**2*(1.-.9*zeta)*zeta**0.66666+ &
86 (1.8-1.4*zeta)*ust**2)+1.e-2
87 dsigwdz=0.5/sigw/h*(-1.4*ust**2+wst**2* &
88 (0.8*max(zeta,1.e-3)**(-.33333)-1.8*zeta**0.66666))
96 if (z.lt.abs(ol))
then
97 tlw=0.1*z/(sigw*(0.55-0.38*abs(z/ol)))
98 else if (zeta.lt.0.1)
then
101 tlw=0.15*h/sigw*(1.-exp(-5*zeta))
110 sigu=1.e-2+2.*ust*(1.-zeta)
111 sigv=1.e-2+1.3*ust*(1.-zeta)
114 tlu=0.15*h/sigu*(sqrt(zeta))
116 tlw=0.1*h/sigw*zeta**0.8
124 if (dsigwdz.eq.0.) dsigwdz=1.e-10