update B/A in the SDC+NSE

AMReX-Astro · Nov 19, 2023 · c6d4d70 · c6d4d70
1 parent 2383aa2
commit c6d4d70
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/integration/nse_update_simplified_sdc.H b/integration/nse_update_simplified_sdc.H
@@ -138,6 +138,8 @@ void sdc_nse_burn(BurnT& state, const Real dt) {
         // update the new state for the next pass
 
         aux_source[iye] = 0.5_rt * rho_half * (dyedt_old + dyedt);
+        aux_source[ibea] = rho_dBEA / dt;
+
         rho_aux_new[iye] = state.y[SFX+iye] + dt * state.ydot_a[SFX+iye] + dt * aux_source[iye];
 
         rho_aux_new[iabar] = state.y[SRHO] * abar_out;