_RL qsat_fac, qsat_exp
qsat_exp = cvapor_exp
qsat_exp = cvapor_exp_ice
tmpbulk= qsat_fac*EXP(-qsat_exp/Tsf)
dEvdT = (tau*re)*ssq*qsat_exp/Ts2
dEvdT = (atmrho*exf_iceCe*wsm)*(ssq*qsat_exp/Ts2)