O dTsf, flxAtm, evpAtm,
dTsf = (flx0 + k12*(Tice(1)-Tsf)) / (k12-df0dT)
Tsf = Tsf + dTsf
& k,Tsf,Tice(1),dTsf
dTsf = 0. _d 0
dTsf = 1000.
if (abs(dTsf).lt.Terrmax) then
write (6,*) 'BB: thermw conv err ',i,j,bi,bj,dTsf
flx0 = flx0 + df0dT*dTsf
evpAtm = evap + dEvdT*dTsf
flxAtm = flxSW + flxExceptSw + df0dT*dTsf
_RL dTsf