Fix mesh generation in init_grid()

Tile index to `grid_global` was not correct inside a loop over tiles.
NOAA-GFDL · Apr 30, 2020 · 022c121 · 022c121
1 parent 75b6ba3
commit 022c121
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/tools/fv_grid_tools.F90 b/tools/fv_grid_tools.F90
@@ -617,8 +617,8 @@ subroutine init_grid(Atm, grid_name, grid_file, npx, npy, npz, ndims, nregions,
 !----------------------------------------------------------------------------------------------------
                       if ( grid_global(i,j,1,n) < 0. )              &
                            grid_global(i,j,1,n) = grid_global(i,j,1,n) + 2.*pi
-                      if (ABS(grid_global(i,j,1,1)) < 1.d-10) grid_global(i,j,1,1) = 0.0
-                      if (ABS(grid_global(i,j,2,1)) < 1.d-10) grid_global(i,j,2,1) = 0.0
+                      if (ABS(grid_global(i,j,1,n)) < 1.d-10) grid_global(i,j,1,n) = 0.0
+                      if (ABS(grid_global(i,j,2,n)) < 1.d-10) grid_global(i,j,2,n) = 0.0
                    enddo
                    enddo
                    enddo