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