_RL N2loc, SNloc
SNloc = Sloc*SQRT(N2loc)
SNloc = 0. _d 0
& *GM_Visbeck_length*GM_Visbeck_length*SNloc
SNloc = SQRT(Ssq(i,j)*N2loc )