stop if invalid grid_type in topography_number_seas

COSIMA · Nov 3, 2024 · 6c49dd9 · 6c49dd9
1 parent 29cb4f5
commit 6c49dd9
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