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
dflhdT = -lath*dEvdT
_RL lath