ztln = LOG(ht/zref)
rh = rhn/( exf_one + rhn*(ztln-psixh)/karman )
re = ren/( exf_one + ren*(ztln-psixh)/karman )
_RL ztln ! = log(zth/zref)
ztln = LOG(zth/zref)
rh = rhn/(1. _d 0 + rhn*(ztln-psixh)/xkar)
re = ren/(1. _d 0 + ren*(ztln-psixh)/xkar)
rh = rhn/(exf_one + rhn*(ztln-psixh)/karman)
re = ren/(exf_one + ren*(ztln-psixh)/karman)