O Tice, sHeating, flxCnB,
flxCnB = 4. _d 0*kice *(Tice(2)-Tf)/hi
& flx0,fct,flx0-fct,flxCnB
_RL flxCnB
I sHeating, flxCnB, evpAtm,
ebot = (flxCnB-Fbot) * dt