_RL Ts2
Ts2 = Tsf*Tsf
dEvdT = (tau*re)*ssq*qsat_exp/Ts2
dEvdT = (atmrho*exf_iceCe*wsm)*(ssq*qsat_exp/Ts2)
flwup = emiss*stefanBoltzmann*Ts2*Ts2
dflwupdT = emiss*stefanBoltzmann*Ts2*Tsf * 4. _d 0
dEvdT = tau*re*ssq*dEvdT/Ts2
flwupa = ocean_emissivity*stefan*Ts2*Ts2
dflwupdT= ocean_emissivity*stefan*Ts2*Tsf*4. _d 0
flwupa = snow_emissivity*stefan*Ts2*Ts2
dflwupdT = snow_emissivity*stefan*Ts2*Tsf*4. _d 0
flwupa = ice_emissivity*stefan*Ts2*Ts2
dflwupdT = ice_emissivity*stefan*Ts2*Tsf*4. _d 0
_RL Ts2 ! surface temperature square [K^2]