diff --git a/Source/MaestroBurner.cpp b/Source/MaestroBurner.cpp index 816199b32..df863a2c6 100644 --- a/Source/MaestroBurner.cpp +++ b/Source/MaestroBurner.cpp @@ -117,6 +117,7 @@ void Maestro::Burner(const Vector& s_in, Vector& s_out, // fraction is above the cutoff. if ((rho > burning_cutoff_density_lo && rho < burning_cutoff_density_hi) && + (T_in > burning_cutoff_temperature_lo) && (ispec_threshold < 0 || (ispec_threshold > 0 && x_test > burner_threshold_cutoff))) { diff --git a/Source/param/_cpp_parameters b/Source/param/_cpp_parameters index 04c765ae0..f56a65ea4 100644 --- a/Source/param/_cpp_parameters +++ b/Source/param/_cpp_parameters @@ -308,6 +308,9 @@ anelastic_cutoff_density Real -1.0 # Refer to Section \ref{Sec:Base Cutoff Density} base_cutoff_density Real -1.0 +# The temperature below which we disable burning +burning_cutoff_temperature_lo Real -1.0 + # The density below which we disable burning burning_cutoff_density_lo Real -1.0