czol = hu*karman*gravity_mks
usn = ws/( exf_one + rdn*(zwln-psimh)/karman )
rd = rdn/( exf_one + rdn*(zwln-psimh)/karman )
rh = rhn/( exf_one + rhn*(ztln-psixh)/karman )
re = ren/( exf_one + ren*(ztln-psixh)/karman )
_RL karman
PARAMETER ( karman = 0.4 _d 0 )
dzTmp = (zwln-psimh)/karman
usn = sh(i,j,bi,bj)/(exf_one - rdn(i,j)/karman*psimh)
rd(i,j) = rdn(i,j)/(exf_one - rdn(i,j)/karman*psimh)
rh = rhn/(exf_one + rhn*(ztln-psixh)/karman)
re = ren/(exf_one + ren*(ztln-psixh)/karman)