fix

AMReX-Astro · Sep 19, 2024 · 058febb · 058febb
1 parent 950b9c4
commit 058febb
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 6 deletions.
diff --git a/Source/driver/timestep.cpp b/Source/driver/timestep.cpp
@@ -82,7 +82,7 @@ Castro::estdt_cfl (int is_new)
       Real dt2;
 #if AMREX_SPACEDIM >= 2
       dt2 = dx[1]/(c + std::abs(uy));
-      if (geom.IsSPHERICAL) {
+      if (geom.IsSPHERICAL()) {
           // dx[1] in Spherical2D is just dtheta, need rdtheta for physical length
           // so just multiply by the smallest r
           dt2 *= geom.ProbLo(0) + 0.5_rt * dx[0];
@@ -211,7 +211,7 @@ Castro::estdt_mhd (int is_new)
       Real dt2;
 #if AMREX_SPACEDIM >= 2
       dt2 = dx[1]/(cy + std::abs(uy));
-      if (geom.IsSPHERICAL) {
+      if (geom.IsSPHERICAL()) {
           dt2 *= geom.ProbLo(0) + 0.5_rt * dx[0];
       }
 #else
@@ -295,7 +295,7 @@ Castro::estdt_temp_diffusion (int is_new)
           Real dt2;
 #if AMREX_SPACEDIM >= 2
           dt2 = 0.5_rt * dx[1]*dx[1] / D;
-          if (geom.IsSPHERICAL) {
+          if (geom.IsSPHERICAL()) {
               Real r = geom.ProbLo(0) + 0.5_rt * dx[0];
               dt2 *= r * r;
           }
@@ -382,7 +382,7 @@ Castro::estdt_burning (int is_new)
 #else
         Real r = 1.0_rt;
 #if AMREX_SPACEDIM >= 2
-        if (geom.IsSPHERICAL) {
+        if (geom.IsSPHERICAL()) {
             r = geom.ProbLo(0) + 0.5_rt * dx[0];
         }
 #endif
@@ -541,7 +541,7 @@ Castro::estdt_rad (int is_new)
             Real dt1 = dx[0] / (c + std::abs(ux));
 #if AMREX_SPACEDIM >= 2
             Real dt2 = dx[1] / (c + std::abs(uy));
-            if (geom.IsSPHERICAL) {
+            if (geom.IsSPHERICAL()) {
                 dt2 *= geom.ProbLo(0) + 0.5_rt * dx[0];
             }
 #else

diff --git a/Source/hydro/Castro_hydro.cpp b/Source/hydro/Castro_hydro.cpp
@@ -243,7 +243,7 @@ Castro::check_for_cfl_violation(const MultiFab& State, const Real dt)
     Real dtdy = 0.0_rt;
     if (AMREX_SPACEDIM >= 2) {
       dtdy = dt / dx[1];
-      if (geom.IsSPHERICAL) {
+      if (geom.IsSPHERICAL()) {
           // dx[1] in Spherical2D is just rdtheta, need rdtheta for physical length
           // Just choose to divide by the smallest r
           dtdy /= geom.ProbLo(0) + 0.5_rt * dx[0];