_RL Ts2
Ts2 = Tsf*Tsf
dEvdT(i,j) = ssq(i,j)*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]