Check if number of points in MOM5 vertical grid is odd, stop otherwise.

COSIMA · Sep 18, 2024 · e1374f5 · e1374f5
1 parent 8c7fd1d
commit e1374f5
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/vgrid.f90 b/src/vgrid.f90
@@ -74,6 +74,10 @@ type(vgrid_t) function vgrid_constructor(filename, type) result(vgrid)
     ! Handle the different types of grids
     select case (type)
     case ("mom5")
+      if (mod(zeta_len, 2) == 0) then
+        write(error_unit,'(a)') "ERROR: MOM5 vertical grid has an even number of points, which should never happen."
+        error stop
+      end if
       vgrid%nlevels = zeta_len/2
       allocate(vgrid%zeta_super(zeta_len))
       allocate(vgrid%zeta(0:vgrid%nlevels))