ebot = (flxCnB-Fbot) * dt
if (ebot.gt.0. _d 0) then
ebote = frace*ebot
ebot = ebot-ebote
& etop,etope,ebot,ebote
if (ebot .lt. 0. _d 0) then
dhi = -ebot / (qbot * rhoi)
ebot = 0. _d 0
if (ebot.gt.0. _d 0 .and. hnew(k).gt.0. _d 0) then
if (ebot .lt. rqh) then
hnew(k) = hnew(k) - ebot / rq
ebot = ebot - rqh
if (ebot.gt.0. _d 0 .and. hs.gt.0. _d 0) then
hs = hs - ebot / rq
& etop, ebot, hi, hs
qleft = qleft + (Fbot+(esurp+etop+ebot)/dt)
& ,esurp/dt,etop/dt,ebot/dt
_RL ebot ! energy for bottom melting (J m-2)