_RL lath
lath = flamb
lath = flamb+flami
hlLocal = -lath*evapLoc
dflhdT = -lath*dEvdT
lath=Lvap
if (iceornot.gt.0) lath=Lvap+Lfresh
ssq = ssq0*exp( lath*(ssq1-ssq2/Tsf) ) / p0
clha = rhoa*lath*us*re*rd
evp = -flha/lath
_RL lath ! latent heat of vaporization or sublimation
lath = Lvap+Lfresh
clha = lath*tau*re
_RL lath ! latent heat of vaporization or sublimation [J/kg]
lath=2.5d6
ssq = ssq0*exp( lath*(ssq1-ssq2/toa) ) / p0
clha = rhoa*lath*usm*re*rd
_RL lath,p0,deltap,delq
ssqa = ssq0*exp( lath*(ssq1-ssq2/local)) / p0
_RL ssq0,ssq1,ssq2,lath,p0,ssqa,q