_RL qsat_fac, qsat_exp
qsat_fac = cvapor_fac
qsat_fac = cvapor_fac_ice
tmpbulk= qsat_fac*EXP(-qsat_exp/Tsf)