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