update the massive star to work with the new NSE table

massive_star/GNUmakefile

@@ -17,7 +17,7 @@ USE_CXX_REACTIONS = TRUE
 
 USE_FORT_MICROPHYSICS = FALSE
 
-MAX_ZONES := 16384
+MAX_ZONES := 32768
 
 DEFINES += -DNPTS_MODEL=$(MAX_ZONES)
 

massive_star/README.md

@@ -8,3 +8,14 @@ convert_21_to_19.py
 
 This setup uses Ye as the primary composition variable from the initial
 model in regions that are in NSE.
+
+Note: you should ensure that the NSE conditions in the inputs file match
+those of your simulation, so the model will be properly in HSE.
+
+Also note that when running with 32768 zones, you need to do:
+
+```
+ulimit -s 16384
+```
+
+Since the arrays are put on the stack.

massive_star/init_1d.H

@@ -77,11 +77,16 @@ set_aux(eos_t& eos_state) {
         // we are in NSE, so leave ye alone, but get abar and redo xn
 
         Real abar_pass;
-        Real dq_pass;
+        Real bea_pass;
         Real dyedt_pass;
+        Real dabardt_pass;
+        Real dbeadt_pass;
+        Real e_nu_pass;
 
         nse_interp(eos_state.T, eos_state.rho, eos_state.aux[AuxZero::iye],
-                   abar_pass, dq_pass, dyedt_pass, eos_state.xn);
+                   abar_pass, bea_pass,
+                   dyedt_pass, dabardt_pass, dbeadt_pass, e_nu_pass,
+                   eos_state.xn);
 
         eos_state.aux[AuxZero::iabar] = abar_pass;
 

massive_star/inputs

@@ -1,9 +1,11 @@
 problem.model_file = 15m_500_sec.aprox19.dat
 
 problem.xmin = 0.0
-problem.xmax = 2.4e10
+problem.xmax = 32768000000.0
 
-problem.nx = 12000
+problem.nx = 32768
 
-network.rho_nse = 2.e6
+network.rho_nse = 1.e7
 network.T_nse = 3.e9
+
+network.nse_table_interp_linear = 0