Correted a OMP bug

NOAA-GFDL · Jul 31, 2024 · 4d55159 · 4d55159
1 parent 8f2c547
commit 4d55159
Showing 1 changed file with 6 additions and 2 deletions.
diff --git a/src/core/MOM_isopycnal_slopes.F90 b/src/core/MOM_isopycnal_slopes.F90
@@ -250,7 +250,9 @@ subroutine calc_isoneutral_slopes(G, GV, US, h, e, tv, dt_kappa_smooth, use_stan
   !$OMP                                  haB,haL,haR,dzaL,dzaR,wtA,wtB,wtL,wtR,drdz,  &
   !$OMP                                  drdx,mag_grad2,slope,l_seg)
   if (present_N2_u .or. (present(dzSxN))) then
-    GxSpV_u(:) = G_Rho0
+    do I=is-1,ie
+      GxSpV_u(I) = G_Rho0
+    enddo
   endif
   do j=js,je ; do K=nz,2,-1
     if (.not.(use_EOS)) then
@@ -395,7 +397,9 @@ subroutine calc_isoneutral_slopes(G, GV, US, h, e, tv, dt_kappa_smooth, use_stan
   !$OMP                                  haB,haL,haR,dzaL,dzaR,wtA,wtB,wtL,wtR,drdz,  &
   !$OMP                                  drdy,mag_grad2,slope,l_seg)
   if ((present_N2_v) .or. (present(dzSyN))) then
-    GxSpV_v(:) = G_Rho0
+    do i=is,ie
+      GxSpV_v(i) = G_Rho0
+    enddo
   endif
   do J=js-1,je ; do K=nz,2,-1
     if (.not.(use_EOS)) then