czol = hu*karman*gravity_mks
& )*czol/(ustar*ustar)
_RL czol
czol = zref*xkar*gravity
huol = czol/ustar**2 *(tstar/t0 +
_RL czol ! = zref*Karman_cst*gravity
_RL aln,cdalton,czol,psim_fac
& )*czol/(ustar(i,j)*ustar(i,j))