don't include Flat dimensions in substep calculation

CliMA · Jun 4, 2023 · b395445 · b395445
1 parent a7ec07e
commit b395445
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/Models/HydrostaticFreeSurfaceModels/split_explicit_free_surface.jl b/src/Models/HydrostaticFreeSurfaceModels/split_explicit_free_surface.jl
@@ -270,9 +270,9 @@ function SplitExplicitSettings(FT::DataType=Float64;
             throw(ArgumentError("Need to specify the grid and max_Δt kwargs to calculate the barotropic substeps from the cfl"))
         end
 
-        Δx = minimum_xspacing(grid)
-        Δy = minimum_yspacing(grid)
-        Δs = sqrt(1 / (1 / Δx^2 + 1 / Δy^2))
+        Δx⁻² = topology(grid)[1] == Flat ? 0 : 1 / minimum_xspacing(grid)^2
+        Δy⁻² = topology(grid)[2] == Flat ? 0 : 1 / minimum_yspacing(grid)^2
+        Δs = sqrt(1 / (Δx⁻² + Δy⁻²))
 
         wave_speed = sqrt(gravitational_acceleration * grid.Lz)