huol = czol/ustar**2 *(tstar/t0 +
huol = sign( min(abs(huol),10. _d 0), huol)
stable = 5. _d -1 + sign(5. _d -1 , huol)
xsq = max(sqrt(abs(1. _d 0 - 16. _d 0*huol)),1. _d 0)
psimh = -5. _d 0*huol*stable + (1. _d 0-stable)*
psixh = -5. _d 0*huol*stable + (1. _d 0-stable)*
_RL huol
huol = czol*(tstar/t0 +
huol = max(huol,zolmin)
stable = exf_half + sign(exf_half, huol)
htol = huol*ht/hu
hqol = huol*hq/hu
xsq = max(sqrt(abs(exf_one - 16.*huol)),exf_one)
psimh = -psim_fac*huol*stable +