_RL flx0 ! net surf heat flux, from Atmos. to sea-ice (W m-2)
flx0 = fswdn + flxExceptSw
flx0 = fswdn + flxExcSw(1)
& flx0,df0dT,k12,k12-df0dT
b1 = b10 - k12*(flx0-df0dT*Tsf) / (k12-df0dT)
dTsf = (flx0 + k12*(Tice(1)-Tsf)) / (k12-df0dT)
flx0 = fswdn + flxExcSw(0)
write (6,*) 'BB: thermw conv err: Tsf, flx0', Tsf,flx0
flx0 = flx0 + df0dT*dTsf
sHeating = flx0 - fct
& flx0,fct,flx0-fct,flxCnB