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