From fc7cfce6a08ce3d09bd2cd0f4169046976774458 Mon Sep 17 00:00:00 2001 From: zhi Date: Thu, 14 Nov 2024 19:43:19 -0500 Subject: [PATCH] fix ptheta reflux --- Source/driver/Castro.cpp | 19 +++++++------------ 1 file changed, 7 insertions(+), 12 deletions(-) diff --git a/Source/driver/Castro.cpp b/Source/driver/Castro.cpp index 184f930340..b55e38d9ec 100644 --- 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.