60 if (h/abs(ol).lt.1.)
then
62 sigw=1.3*exp(-2.e-4*z/ust)
65 tlw=0.5*z/sigw/(1.+1.5e-3*z/ust)
72 else if (ol.lt.0.)
then
78 sigw=sqrt(1.2*wst**2*(1.-.9*zeta)*zeta**0.66666+ &
79 (1.8-1.4*zeta)*ust**2)+1.e-2
80 dsigwdz=0.5/sigw/h*(-1.4*ust**2+wst**2* &
81 (0.8*max(zeta,1.e-3)**(-.33333)-1.8*zeta**0.66666))
87 if (z.lt.abs(ol))
then
88 tlw=0.1*z/(sigw*(0.55-0.38*abs(z/ol)))
89 else if (zeta.lt.0.1)
then
92 tlw=0.15*h/sigw*(1.-exp(-5*zeta))
101 sigw=1.e-2+1.3*ust*(1.-zeta)
103 tlw=0.1*h/sigw*zeta**0.8
110 if (dsigwdz.eq.0.) dsigwdz=1.e-10
subroutine hanna_short(z)