_RL GM_K3D_minN2
& GM_maxK3D, GM_K3D_minCori, GM_K3D_minN2,
N2loc(i,j,k)=GM_K3D_minN2
N2W(i,j,k) = GM_K3D_minN2
N2S(i,j,k) = GM_K3D_minN2
N2(i,j,k) = MAX(N2(i,j,k),GM_K3D_minN2)*maskC(i,j,k,bi,bj)
maxDRhoDz = -rhoConst*GM_K3D_minN2/gravity
GM_K3D_minN2 = 1.0e-8
& GM_K3D_maxC, GM_K3D_minCori, GM_K3D_minN2,