_RL esurp
esurp = 0. _d 0
esurp = sHeating * dt
esurp = esurp - rhos*qsnow*hs
esurp = esurp - rhoi*qicen(k)*hnew(k)
& 'ThSI_CALC_TH: -1 : esurp=',esurp
esurp = esurp - hnew(k)*rhoi*(qicen(k)-Lfresh)
& 'ThSI_CALC_TH: -2 : esurp,fresh=', esurp, fresh
qleft = qleft + (Fbot+(esurp+etop+ebot)/dt)
& 'ThSI_CALC_TH: [esurp,etop+ebot]/dt =',esurp/dt,etop/dt,ebot/dt
esurp = fzMlOc(i,j)
newIce = esurp*deltaTice/qicAv