update apply_source flag and add geometric source terms for 1d spherical

zhichen3 · Nov 13, 2024 · 6766e37 · 6766e37
1 parent b38a018
commit 6766e37
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 5 deletions.
diff --git a/Source/sources/Castro_geom.cpp b/Source/sources/Castro_geom.cpp
@@ -26,7 +26,7 @@ Castro::construct_old_geom_source(MultiFab& source, MultiFab& state_in, Real tim
       return;
   }
 
-#if AMREX_SPACEDIM == 2
+#if AMREX_SPACEDIM <= 2
   const Real strt_time = ParallelDescriptor::second();
 
   MultiFab geom_src(grids, dmap, source.nComp(), 0);
@@ -83,7 +83,7 @@ Castro::construct_new_geom_source(MultiFab& source, MultiFab& state_old,
       return;
   }
 
-#if AMREX_SPACEDIM == 2
+#if AMREX_SPACEDIM <= 2
   const Real strt_time = ParallelDescriptor::second();
 
   MultiFab geom_src(grids, dmap, source.nComp(), 0);

diff --git a/Source/sources/Castro_sources.cpp b/Source/sources/Castro_sources.cpp
@@ -61,10 +61,10 @@ Castro::source_flag(int src)
 #endif
 
     case geom_src:
-        if (geom.Coord() == 1) {
-          return true;
-        } else {
+        if (geom.Coord() == 0) {
           return false;
+        } else {
+          return true;
         }
 
 #ifdef DIFFUSION