Properly normalize the initial composition in Detonation

yut23 · Mar 27, 2024 · 7fae7f6 · 7fae7f6
1 parent aa8d3b2
commit 7fae7f6
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/Exec/science/Detonation/problem_initialize.H b/Exec/science/Detonation/problem_initialize.H
@@ -55,10 +55,10 @@ void problem_initialize ()
 
     if (problem::in14 >= 0) {
        problem::xn[problem::in14] = amrex::max(problem::nfrac, problem::smallx);
-       problem::xn[problem::ihe4] = 1.0_rt - problem::cfrac - problem::nfrac - problem::ofrac - (NumSpec - 3) * problem::smallx;
+       problem::xn[problem::ihe4] = 1.0_rt - problem::xn[problem::ic12] - problem::xn[problem::in14] - problem::xn[problem::io16] - (NumSpec - 4) * problem::smallx;
     }
     else {
-       problem::xn[problem::ihe4] = 1.0_rt - problem::cfrac - problem::ofrac - (NumSpec - 2) * problem::smallx;
+       problem::xn[problem::ihe4] = 1.0_rt - problem::xn[problem::ic12] - problem::xn[problem::io16] - (NumSpec - 3) * problem::smallx;
     }
 
     // Set the ambient material