Bug - indices are transposed

gustavo-marques · Jun 15, 2021 · e551526 · e551526
1 parent 445b865
commit e551526
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/core/MOM_stoch_eos.F90 b/src/core/MOM_stoch_eos.F90
@@ -158,8 +158,8 @@ subroutine MOM_calc_varT(G,Gv,h,tv,stoch_eos_CS)
   ! still a poor approximation in the interior when coordinates are strongly tilted.
   if (.not. associated(tv%varT)) call safe_alloc_ptr(tv%varT, G%isd, G%ied, G%jsd, G%jed, GV%ke)
   do k=1,G%ke
-     do j=G%isc,G%iec
-        do i=G%jsc,G%jec
+     do j=G%jsc,G%jec
+        do i=G%isc,G%iec
            hl(1) = h(i,j,k) * G%mask2dT(i,j)
            hl(2) = h(i-1,j,k) * G%mask2dCu(I-1,j)
            hl(3) = h(i+1,j,k) * G%mask2dCu(I,j)