remove 2d spherical abort and add some asserts setup

AMReX-Astro · Oct 7, 2024 · fc38b7d · fc38b7d
1 parent f79e5a8
commit fc38b7d
Showing 1 changed file with 9 additions and 1 deletion.
diff --git a/Source/driver/Castro.cpp b/Source/driver/Castro.cpp
@@ -306,7 +306,15 @@ Castro::read_params ()
 #elif (AMREX_SPACEDIM == 2)
     if ( dgeom.IsSPHERICAL() )
       {
-        amrex::Abort("We don't support spherical coordinate systems in 2D");
+        if ( (dgeom.ProbLo(1) >= 0.0_rt) && (dgeom.ProbHi(1) <= 3.14159265359_rt) )
+        {
+          amrex::Abort("Theta must be within [0, Pi] for spherical coordinate system in 2D");
+        }
+
+        if ( dgeom.ProbLo(0) > static_cast<Real>(NUM_GROW)*dgeom.CellSize(0) )
+        {
+          amrex::Abort("R-min must be large enough so ghost cells doesn't extend to negative R");
+        }
       }
 #elif (AMREX_SPACEDIM == 3)
     if ( dgeom.IsRZ() )