_RL deltaH, integrDepth
integrDepth = -rC( kLowC(i,j,bi,bj) )
integrDepth = MIN( integrDepth, GM_Visbeck_depth )
integrDepth = MAX( integrDepth, GM_Visbeck_minDepth )
deltaH = integrDepth + rC(k-1)
deltaH = deltaH/( integrDepth+rC(1) )
integrDepth = drF(k)
deltaH=min(deltaH,integrDepth)