_RL dhi
dhi = -ebot / (qbot * rhoi)
qicen(nlyr) = (hnew(nlyr)*qicen(nlyr)+dhi*qbot) /
& (hnew(nlyr)+dhi)
hnew(nlyr) = hnew(nlyr) + dhi
dhi = evap/rhoi*dt
IF (dhi.GE.hnew(k)) THEN
hq = hnew(k)*qicen(k)-dhi*Lfresh
hnew(k) = hnew(k) - dhi
dhi = dhs * rhos / rhoi
hnew(1) = hnew(1) + dhi
hi = hi + dhi