stop if invalid grid_type in topography_number_seas

COSIMA · Nov 4, 2024 · 1b255f1 · 1b255f1
1 parent f7d11bd
commit 1b255f1
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/src/topography.f90 b/src/topography.f90
@@ -350,6 +350,9 @@ subroutine topography_number_seas(this, sea_number, number_of_seas, silent)
                 mask=[choke_west, choke_east, choke_south, choke_north]))
             case ('C')
               new_sea = min(sea(i, j), sea(im, j), sea(ip, j), sea(i, jm), sea(i, jp))
+            case default
+              write(error_unit, '(a)') "topogtools: grid_type must be B or C"
+              error stop
             end select
             if (sea(i, j) /= new_sea) then
               sea(i, j) = new_sea
@@ -397,6 +400,9 @@ subroutine topography_number_seas(this, sea_number, number_of_seas, silent)
                 mask=[choke_west, choke_east, choke_south, choke_north]))
             case ('C')
               new_sea = min(sea(i, j), sea(im, j), sea(ip, j), sea(i, jm), sea(i, jp))
+            case default
+              write(error_unit, '(a)') "topogtools: grid_type must be B or C"
+              error stop
             end select
             if (sea(i, j) /= new_sea) then
               sea(i, j) = new_sea