fix ptheta reflux

zhichen3 · Nov 15, 2024 · fc7cfce · fc7cfce
1 parent b11b943
commit fc7cfce
Showing 1 changed file with 7 additions and 12 deletions.
diff --git a/Source/driver/Castro.cpp b/Source/driver/Castro.cpp
@@ -2875,9 +2875,9 @@ Castro::reflux (int crse_level, int fine_level, bool in_post_timestep)
 
             reg = &getLevel(lev).pres_reg;
 
-            reg->ClearInternalBorders(crse_lev.geom);
+            Reg->ClearInternalBorders(crse_lev.geom);
 
-            reg->Reflux(crse_state, crse_lev.volume, 1.0, 0, UMX, 1, crse_lev.geom);
+            reg->Reflux(crse_state, crse_lev.volume, 0, 1.0, 0, UMX, 1, crse_lev.geom);
 
             if (update_sources_after_reflux || !in_post_timestep) {
 
@@ -2899,19 +2899,10 @@ Castro::reflux (int crse_level, int fine_level, bool in_post_timestep)
 
             }
 
-            reg->setVal(0.0);
-
-        }
-#endif
-
 #if (AMREX_SPACEDIM == 2)
         if (Geom().IsSPHERICAL()) {
 
-            reg = &getLevel(lev).pres_reg;
-
-            reg->ClearInternalBorders(crse_lev.geom);
-
-            reg->Reflux(crse_state, crse_lev.volume, 1.0, 0, UMY, 1, crse_lev.geom);
+            reg->Reflux(crse_state, crse_lev.volume, 1, 1.0, 0, UMY, 1, crse_lev.geom);
 
             if (update_sources_after_reflux || !in_post_timestep) {
 
@@ -2933,11 +2924,15 @@ Castro::reflux (int crse_level, int fine_level, bool in_post_timestep)
 
             }
 
+        }
+#endif
+
             reg->setVal(0.0);
 
         }
 #endif
 
+
 #ifdef RADIATION
 
         // This follows the same logic as the pure hydro fluxes; see above for details.