_RL GM_Visbeck_depth
& GM_Visbeck_depth,
deltaH=abs(GM_Visbeck_depth)-abs(rF(k))
deltaH=deltaH/GM_Visbeck_depth
GM_Visbeck_depth = 1000.D0