O flwupa, flha, fsha, df0dT,
flwupa=ocean_emissivity*stefan*tsf**4
flwupa=snow_emissivity*stefan*tsf**4
flwupa=ice_emissivity*stefan*tsf**4
_RL flwupa ! upward long wave radiation (>0 upward) [W/m2]
flwupa = ocean_emissivity*stefan*Ts2*Ts2
flwupa = snow_emissivity*stefan*Ts2*Ts2
flwupa = ice_emissivity*stefan*Ts2*Ts2