flwupa=ocean_emissivity*stefan*tsf**4
dflwupdT=4. _d 0*ocean_emissivity*stefan*tsf**3
WRITE(iUnit,*) 'BlkF: ocean_emissivity=',ocean_emissivity
& atm_emissivity, ocean_emissivity,
ocean_emissivity= .985 _d 0
_RL ocean_emissivity