rdn=xkar/(log(zref/zice))
rhn=rdn
ren=rdn
ustar=rdn*usm
rd = rdn/(1. _d 0 + rdn*(aln-psimh)/xkar)
_RL rdn ,ren, rhn
_RL rdn
rdn = ustar/us
rdn = 0. _d 0
rdn = sqrt(tmpbulk)
ustar = rdn*sh
rd = rdn/(exf_one - rdn/karman*psimh )
shn = sh*rd/rdn
rd = rdn/(exf_one - rdn/karman*psimh)