WIP: Minor fixes

Author: James471

integration/VODE/vode_dvstep.H

@@ -309,6 +309,7 @@ int dvstep (BurnT& state, DvodeT& vstate)
 
         }
 
+#ifdef PHOTOCHEMISTRY
         for (int i = net_ienuc + 1; i <= int_neqs; i+=MicrophysicsNumRadVarsPerGroup) {
 
 #ifdef STRANG
@@ -349,6 +350,7 @@ int dvstep (BurnT& state, DvodeT& vstate)
                 break;
             }
 #endif
+#endif
 
 #ifdef SDC
             // Abort the program. This is unimplemented

integration/_parameters

@@ -23,28 +23,28 @@ burner_verbose           bool     0
 # species and energy equations.
 
 # relative tolerance for species
-rtol_spec                real      1.e-12
+rtol_spec                real      1.e-10
 
 # relative tolerance for energy
-rtol_enuc                real      1.e-6
+rtol_enuc                real      1.e-10
 
 # absolute tolerance for species
-atol_spec                real      1.e-8
+atol_spec                real      1.e-10
 
 # absolute tolerance for energy
-atol_enuc                real      1.e-6
+atol_enuc                real      1.e-25
 
-# relative tolerance for radiation density
-rtol_rad_num             real      1.e-12
+# relative tolerance for radiation number density
+rtol_rad_num             real      1.e-10
 
-# relative teolerance for radiation flux
-rtol_rad_flux            real      1.e-6
+# relative tolerance for radiation flux
+rtol_rad_flux            real      1.e-10
 
-# absolute tolerance for radiation density
-atol_rad_num             real      1.e-12
+# absolute tolerance for radiation number density
+atol_rad_num             real      1.e-10
 
 # absolute tolerance for radiation flux
-atol_rad_flux            real      1.e-6
+atol_rad_flux            real      1.e-10
 
 # Whether to renormalize the mass fractions at each step in the evolution
 # so that they sum to unity.

integration/integrator_data.H

@@ -15,10 +15,12 @@ constexpr int INT_NEQS = NumRadEqs + NumSpec + 1;
 // abundance is unreasonably small or large (each X must satisfy
 // -failure_tolerance <= X <= 1.0 + failure_tolerance).
 constexpr amrex::Real species_failure_tolerance = 1.e-2_rt;
+#ifdef PHOTOCHEMISTRY
 // We will use this parameter to determine if 
 // a given radiation quantity (number density or flux) is 
 // unreasonably small (each quantity must satisfy X >= -radiation_failure_tolerance).
 constexpr amrex::Real radiation_failure_tolerance = 1.e-2_rt;
+#endif
 
 #ifdef NSE
 constexpr int MIN_NSE_BAILOUT_STEPS = 10;

integration/integrator_setup_strang.H

@@ -207,11 +207,13 @@ void integrator_cleanup (IntegratorT& int_state, BurnT& state,
             }
         }
     }
+#ifdef PHOTOCHEMISTRY
     for (int n = net_ienuc+1; n <= net_ienuc+NumRadEqs; n += MicrophysicsNumRadVarsPerGroup) {
         if (int_state.y(n) < -radiation_failure_tolerance) {
             state.success = false;
         }
     }
+#endif
 
 #ifndef AMREX_USE_GPU
     if (integrator_rp::burner_verbose) {