Merge pull request #39 from ofuhrer/patch-1

Fix mesh generation in init_grid()
NOAA-GFDL · May 4, 2020 · eee6555 · eee6555
2 parents 75b6ba3 + 022c121
commit eee6555
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