use balance charge function

AMReX-Astro · Aug 29, 2024 · 6e2141f · 6e2141f
1 parent 7056045
commit 6e2141f
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/unit_test/burn_cell_primordial_chem/burn_cell.H b/unit_test/burn_cell_primordial_chem/burn_cell.H
@@ -8,6 +8,7 @@
 #include <eos.H>
 #include <extern_parameters.H>
 #include <network.H>
+#include <actual_network.H>
 
 amrex::Real grav_constant = 6.674e-8;
 
@@ -209,8 +210,7 @@ auto burn_cell_c() -> int {
     }
 
     // update the number density of electrons due to charge conservation
-    state.xn[0] = -state.xn[3] - state.xn[7] + state.xn[1] + state.xn[12] +
-                  state.xn[6] + state.xn[4] + state.xn[9] + 2.0 * state.xn[11];
+    balance_charge(state);
 
     // reconserve mass fractions post charge conservation
     insum = 0;