60 if (h/abs(ol).lt.1.)
then
63 sigu=2.0*ust*exp(-3.e-4*z/ust)
65 sigv=1.3*ust*exp(-2.e-4*z/ust)
68 dsigw2dz=-6.76e-4*ust*exp(-4.e-4*z/ust)
69 tlu=0.5*z/sigw/(1.+1.5e-3*z/ust)
78 else if (ol.lt.0.)
then
84 sigu=ust*(12-0.5*h/ol)**0.33333
88 if (zeta.lt.0.03)
then
89 sigw=0.96*wst*(3*zeta-ol/h)**0.33333
90 dsigw2dz=1.8432*wst*wst/h*(3*zeta-ol/h)**(-0.33333)
91 else if (zeta.lt.0.4)
then
92 s1=0.96*(3*zeta-ol/h)**0.33333
96 dsigw2dz=1.8432*wst*wst/h*(3*zeta-ol/h)**(-0.33333)
99 dsigw2dz=0.203759*wst*wst/h*zeta**(-0.65)
101 else if (zeta.lt.0.96)
then
102 sigw=0.722*wst*(1-zeta)**0.207
103 dsigw2dz=-.215812*wst*wst/h*(1-zeta)**(-0.586)
104 else if (zeta.lt.1.00)
then
116 if (z.lt.abs(ol))
then
117 tlw=0.1*z/(sigw*(0.55-0.38*abs(z/ol)))
118 else if (zeta.lt.0.1)
then
121 tlw=0.15*h/sigw*(1.-exp(-5*zeta))
130 sigu=2.*ust*(1.-zeta)
131 sigv=1.3*ust*(1.-zeta)
135 dsigw2dz=3.38*ust*ust*(zeta-1.)/h
136 tlu=0.15*h/sigu*(sqrt(zeta))
138 tlw=0.1*h/sigw*zeta**0.8