deltaH = integrDepth + rC(k-1)
deltaH = MIN( deltaH, drC(k) )
deltaH = deltaH/( integrDepth+rC(1) )
& +deltaH*GM_Visbeck_alpha
deltaH=abs(GM_Visbeck_depth)-abs(rF(k))
deltaH=min(deltaH,integrDepth)
deltaH=max(deltaH, 0. _d 0)
deltaH=deltaH/GM_Visbeck_depth
_RL deltaH, integrDepth