Merge branch 'development' of github.com:amrex-astro/Microphysics int…

…o development
AMReX-Astro · Sep 10, 2023 · 3c75538 · 3c75538
2 parents a7d5990 + 743d388
commit 3c75538
Show file tree

Hide file tree

Showing 3 changed files with 10 additions and 2 deletions.
diff --git a/integration/nse_update_simplified_sdc.H b/integration/nse_update_simplified_sdc.H
@@ -73,7 +73,7 @@ void sdc_nse_burn(BurnT& state, const Real dt) {
     // reaction terms and advection to dt
 
     eos_re_t eos_state;
-    eos_state.T = 1.e6_rt; // initial guess
+    eos_state.T = T_in; // initial guess
 
     // initial aux_sources
     Real aux_source[NumAux] = {0.0_rt};
@@ -238,7 +238,7 @@ void sdc_nse_burn(BurnT& state, const Real dt) {
     // updates with advection to dt
 
     BurnT burn_state;
-    burn_state.T = 1.e6_rt; // initial guess
+    burn_state.T = T_in; // initial guess
     burn_state.mu_p = state.mu_p;
     burn_state.mu_n = state.mu_n;
 

diff --git a/unit_test/burn_cell_sdc/_parameters b/unit_test/burn_cell_sdc/_parameters
@@ -69,3 +69,5 @@ Adv_Aux3          real       0.0e0
 
 
 
+mu_p              real       -5.0
+mu_n              real      -12.0
diff --git a/unit_test/burn_cell_sdc/burn_cell.H b/unit_test/burn_cell_sdc/burn_cell.H
@@ -293,6 +293,12 @@ void burn_cell_c()
 
     burn_state.T_fixed = -1.0;
 
+#ifdef NSE_NET
+    burn_state.mu_p = unit_test_rp::mu_p;
+    burn_state.mu_n = unit_test_rp::mu_n;
+    burn_state.y_e = -1.0;
+#endif
+
     std::ofstream state_over_time("state_over_time.txt");
 
     // we will divide the total integration time into nsteps that are
Original file line number	Diff line number	Diff line change
Expand Up		@@ -69,3 +69,5 @@ Adv_Aux3 real 0.0e0



		mu_p real -5.0
		mu_n real -12.0