_RL rqh
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
iceFrac=(1. _d 0-extend/rqh)*iceFrac
IF ( extend.LT.rqh .AND. iceFrac.GE.iceMaskMin ) THEN
fresh=fresh+extend/rqh*freshe
fsalt(i,j)=fsalt(i,j)+extend/rqh*salte
qleft=qleft+(extend-rqh)/dt