_RL GM_Visbeck_depth
& GM_Visbeck_depth,
& -rC(k-1).LT.GM_Visbeck_depth ) THEN
integrDepth = MIN( integrDepth, GM_Visbeck_depth )
deltaH=abs(GM_Visbeck_depth)-abs(rF(k))
deltaH=deltaH/GM_Visbeck_depth
GM_Visbeck_depth = 1000. _d 0