flwupa=ocean_emissivity*stefan*tsf**4
dflwupdT=4. _d 0*ocean_emissivity*stefan*tsf**3
flwupa=snow_emissivity*stefan*tsf**4
dflwupdT=4. _d 0*snow_emissivity*stefan*tsf**3
flwupa=ice_emissivity*stefan*tsf**4
dflwupdT=4. _d 0*ice_emissivity*stefan*tsf**3
WRITE(iUnit,*) 'BlkF: stefan =',stefan
& stefan, xkar, Rvap,
stefan = 5.67 _d -8
_RL stefan