_RL rqh ! rho * q * h for a layer
rqh = rq * hs
if (etop .lt. rqh) then
etop = etop - rqh
rqh = rq * hnew(k)
if (ebot .lt. rqh) then
ebot = ebot - rqh
rqh = rhoi*qicen(1)*hnew(1) + rhos*qsnow*dhs
qicen(1) = rqh / (rhoi*hnew(1))
rqh = rq * hi + rs * hs
if (extend .lt. rqh) then
compact=(1. _d 0-extend/rqh)*compact
fresh=fresh+extend/rqh*freshe
fsalt=fsalt+extend/rqh*salte
qleft=qleft+(extend-rqh)/dt