_RL psixh ! latent & sensib. stability function
psixh = -psim_fac*htol*stable + (exf_one-stable)*
rh = rhn/(exf_one + rhn*(ztln-psixh)/karman)
re = ren/(exf_one + ren*(ztln-psixh)/karman)
psixh = -psim_fac*htol*stable
rh(i,j)= rhn/( exf_one + rhn*(ztln-psixh)/karman )
re(i,j)= ren/( exf_one + ren*(ztln-psixh)/karman )
psixh = -5. _d 0*huol*stable + (1. _d 0-stable)*
rh = rhn/(1. _d 0 + rhn*(aln-psixh)/xkar)
psixh = -5. _d 0*htol*stable
rh = rhn/(1. _d 0 + rhn*(ztln-psixh)/xkar)
re = ren/(1. _d 0 + ren*(ztln-psixh)/xkar)
_RL huol,stable,xsq,x,psimh,psixh