Merge branch 'burn_cell_all_equal' into he-simple

AMReX-Astro · Nov 14, 2024 · 0c23f03 · 0c23f03
2 parents 2c12496 + a0a3983
commit 0c23f03
Show file tree

Hide file tree

Showing 3 changed files with 10 additions and 5 deletions.
diff --git a/unit_test/burn_cell/GNUmakefile b/unit_test/burn_cell/GNUmakefile
@@ -24,7 +24,7 @@ EOS_DIR     := helmholtz
 
 # This sets the network directory
 NETWORK_DIR := aprox13
-
+SCREEN_METHOD := screen5
 CONDUCTIVITY_DIR := stellar
 
 INTEGRATOR_DIR =  VODE

diff --git a/unit_test/burn_cell/_parameters b/unit_test/burn_cell/_parameters
@@ -19,3 +19,5 @@ density       real       1.e7
 temperature   real       3.e9
 
 skip_initial_normalization    bool    0
+
+init_species_all_equal        bool    0
diff --git a/unit_test/burn_cell/burn_cell.H b/unit_test/burn_cell/burn_cell.H
@@ -21,12 +21,15 @@ void burn_cell_c()
     // Make sure user set all the mass fractions to values in the interval [0, 1]
     for (int n = 1; n <= NumSpec; ++n) {
 
-        massfractions[n-1] = get_xn(n);
+        if (unit_test_rp::init_species_all_equal) {
+            massfractions[n-1] = 1.0_rt / static_cast<Real>(NumSpec);
+        } else {
+            massfractions[n-1] = get_xn(n);
 
-        if (massfractions[n-1] < 0 || massfractions[n-1] > 1) {
-            amrex::Error("mass fraction for " + short_spec_names_cxx[n-1] + " not initialized in the interval [0,1]!");
+            if (massfractions[n-1] < 0 || massfractions[n-1] > 1) {
+                amrex::Error("mass fraction for " + short_spec_names_cxx[n-1] + " not initialized in the interval [0,1]!");
+            }
         }
-
     }
 
     // Echo initial conditions at burn and fill burn state input
Original file line number	Diff line number	Diff line change
Expand Up		@@ -19,3 +19,5 @@ density real 1.e7
		temperature real 3.e9

		skip_initial_normalization bool 0

		init_species_all_equal bool 0