BackwardEuler

parameter	description	default value
`maximum_timestep_change_factor`	Maximum amount any quantity can change by in a timestep	1.01
`max_iter`	Maximum number of iterations for the Newton solve	25
`tol`	tolerance for the Newton solve	1.e-10

ForwardEuler

parameter	description	default value
`maximum_timestep_change_factor`	Maximum amount any quantity can change by in a timestep	1.001

QSS

parameter	description	default value
`predictor_corrector_tolerance`	Allowable difference between the predictor and corrector	0.01
`tolerance_safety_factor`	Threshold factor on the error criterion used in timestep selection	10.0
`num_timestep_iters`	Maximum number of iterations on the timestep constraint loop	10
`num_corrector_iters`	Maximum number of iterations on the corrector loop	1
`dt_max_change_factor`	Maximum factor that dt is allowed to grow per timestep	1.05
`dt_cut_factor`	Multiply the timestep by this factor when we outright reject it	0.5
`species_tolerance`	Reject a species update if X < -tol or X > 1.0 + tol	0.01
`dt_init_fraction`	Scale factor for initial timestep	0.01

RKC

parameter	description	default value
`scale_system`	do we scale the ODE system we integrate to make it O(1)? for Strang, this simply means scaling e by the initial energy?	1
`use_circle_theorem`	use the Gershgorin circle theorem to estimate the spectral radius? note: requires integrator.scale_system = 1	1

VODE

parameter	description	default value
`X_reject_buffer`	for the step rejection logic on mass fractions, we only consider species that are > X_reject_buffer * atol_spec	1.0

ase

parameter	description	default value
`disable_p_C12_to_N13`		0
`disable_He4_N13_to_p_O16`		0

breakout

parameter	description	default value
`eos_gamma`		0.e0

burn_cell

parameter	description	default value
`run_prefix`		“”
`small_temp`		1.e5
`small_dens`		1.e5
`tmax`	the final time to integrate to	1.e-2
`tfirst`	first output time – we will output in nsteps logarithmically spaced steps between [tfirst, tmax]	0.0
`nsteps`	number of steps (logarithmically spaced)	100
`density`		1.e7
`temperature`		3.e9
`skip_initial_normalization`		0

burn_cell_primordial_chem

parameter	description	default value
`run_prefix`		“burn_cell_primordial_chem”
`small_temp`	floor values of temperature and density	1.e1
`small_dens`		1.e-30
`tmax`	the final time to integrate to	1.e20
`tff_reduc`	tff_reduc reduces the calculated freefall time to accordingly increase the density during the single zone burn	1.e-1
`tfirst`	first output time – we will output in nsteps logarithmically spaced steps between [tfirst, tmax]	0.0
`nsteps`	number of steps for the single zone burn	1000
`temperature`	initial temperature	1e2
`primary_species_1`	list of species and their number densities used in the network (14 if including deuterium)	1.0e0
`primary_species_2`		0.0e0
`primary_species_3`		0.0e0
`primary_species_4`		0.0e0
`primary_species_5`		0.0e0
`primary_species_6`		0.0e0
`primary_species_7`		0.0e0
`primary_species_8`		0.0e0
`primary_species_9`		0.0e0
`primary_species_10`		0.0e0
`primary_species_11`		0.0e0
`primary_species_12`		0.0e0
`primary_species_13`		0.0e0
`primary_species_14`		0.0e0

burn_cell_sdc

parameter	description	default value
`small_temp`		1.e5
`small_dens`		1.e5
`tmax`	the final time to integrate to	1.e-2
`tfirst`	first output time – we will output in nsteps logarithmically spaced steps between [tfirst, tmax]	0.0
`nsteps`	number of steps (logarithmically spaced)	100
`recompute_aux`	do we recompute the aux quantities? or do we take them as given in the inputs?	0
`density`		1.e7
`temperature`		3.e9
`rhoe`		-1.e0
`Aux1`		0.0e0
`Aux2`		0.0e0
`Aux3`		0.0e0
`Adv_rho`		0.0e0
`Adv_rhoe`		0.0e0
`Adv_X1`		0.0e0
`Adv_X2`		0.0e0
`Adv_X3`		0.0e0
`Adv_X4`		0.0e0
`Adv_X5`		0.0e0
`Adv_X6`		0.0e0
`Adv_X7`		0.0e0
`Adv_X8`		0.0e0
`Adv_X9`		0.0e0
`Adv_X10`		0.0e0
`Adv_X11`		0.0e0
`Adv_X12`		0.0e0
`Adv_X13`		0.0e0
`Adv_X14`		0.0e0
`Adv_X15`		0.0e0
`Adv_X16`		0.0e0
`Adv_X17`		0.0e0
`Adv_X18`		0.0e0
`Adv_X19`		0.0e0
`Adv_X20`		0.0e0
`Adv_X21`		0.0e0
`Adv_X22`		0.0e0
`Adv_X23`		0.0e0
`Adv_X24`		0.0e0
`Adv_X25`		0.0e0
`Adv_X26`		0.0e0
`Adv_X27`		0.0e0
`Adv_X28`		0.0e0
`Adv_X29`		0.0e0
`Adv_X30`		0.0e0
`Adv_X31`		0.0e0
`Adv_X32`		0.0e0
`Adv_X33`		0.0e0
`Adv_X34`		0.0e0
`Adv_X35`		0.0e0
`Adv_Aux1`		0.0e0
`Adv_Aux2`		0.0e0
`Adv_Aux3`		0.0e0
`mu_p`		-5.0
`mu_n`		-12.0

cj_detonation

parameter	description	default value
`smallx`		1.e-10
`small_temp`		1.e5
`small_dens`		1.e5

constant

parameter	description	default value
`const_conductivity`		1.0e0

constant_opacity

parameter	description	default value
`const_opacity`		7.0e-2

eos_cell

parameter	description	default value
`small_temp`		1.e5
`small_dens`		1.e5
`density`		1.e7
`temperature`		3.e9

gamma_law

parameter	description	default value
`eos_gamma`		5.e0/3.e0
`eos_assume_neutral`		1

helmholtz

parameter	description	default value
`use_eos_coulomb`	use the Coulomb corrections	1
`eos_input_is_constant`	Force the EOS output quantities to match input	1
`eos_ttol`	Tolerance for iterations with respect to temperature	1.0e-8
`eos_dtol`	Tolerance for iterations with respect to density	1.0e-8
`prad_limiter_rho_c`	Density target for radiation pressure smoothing (negative means smoothing is disabled)	-1.0e0
`prad_limiter_delta_rho`	Density gradient for radiation pressure smoothing (negative means smoothing is disabled)	-1.0e0

integration

parameter	description	default value
`call_eos_in_rhs`	Normally we update the temperature during a burn to be consistent with the current internal energy. This is done with an EOS call, which can be turned off if desired. This will freeze the temperature and specific heat to the values at the beginning of the burn, which is inaccurate but cheaper.	1
`integrate_energy`	Allow the energy integration to be disabled by setting the RHS to zero.	1
`jacobian`	Whether to use an analytical or numerical Jacobian. 1 == Analytical 2 == Numerical	1
`burner_verbose`	Should we print out diagnostic output after the solve?	0
`rtol_spec`	Tolerances for the solver (relative and absolute), for the species and energy equations.	1.e-12
`rtol_enuc`		1.e-6
`atol_spec`		1.e-8
`atol_enuc`		1.e-6
`renormalize_abundances`	Whether to renormalize the mass fractions at each step in the evolution so that they sum to unity.	0
`SMALL_X_SAFE`	The absolute cutoff for species – note that this might be larger than `small_x`, but the issue is that we need to prevent underflow issues and keep mass fractions positive in the integrator. You may have to increase the floor to, e.g. 1.e-20 if your rates are large.	1.0e-30
`MAX_TEMP`	The maximum temperature for reactions in the integration.	1.0e11
`react_boost`	boost the reaction rates by a factor > 1	-1.e0
`ode_max_steps`	maximum number of timesteps for the integrator	150000
`ode_max_dt`	maximum timestep for the integrator	1.e30
`use_jacobian_caching`	Whether to use Jacobian caching in VODE	1
`nonaka_i`	Inputs for generating a Nonaka Plot (TM)	0
`nonaka_j`		0
`nonaka_k`		0
`nonaka_level`		0
`nonaka_file`		“nonaka_plot.dat”
`use_burn_retry`	do we retry a failed burn with different parameters?	0
`retry_swap_jacobian`	do we swap the Jacobian (from analytic to numerical or vice versa) on a retry?	1
`retry_rtol_spec`	Tolerances for the solver (relative and absolute), for the species and energy equations.	1.e-12
`retry_rtol_enuc`		1.e-6
`retry_atol_spec`		1.e-8
`retry_atol_enuc`		1.e-6
`do_species_clip`	in the clean_state process, do we clip the species such that they are in [0, 1]?	1
`use_number_densities`	flag for turning on the use of number densities for all species	0
`subtract_internal_energy`	flag for tuning on the subtraction of internal energy	1
`sdc_burn_tol_factor`	SDC iteration tolerance adjustment factor	1.e0
`scale_system`	do we scale the ODE system we integrate to make it O(1)? for Strang, this simply means scaling e by the initial energy?	0
`nse_iters`	for the NSE update predictor-corrector, how many iterations do we take to get the new time NSE state	3
`nse_deriv_dt_factor`	for SDC+NSE, when estimating the derivatives of the NSE table quantities, what fraction of dt do we use for the finite-difference estimate	0.05
`nse_include_enu_weak`	for NSE update, do we include the weak rate neutrino losses?	1
`linalg_do_pivoting`	for the linear algebra, do we allow pivoting?	1

make_table

parameter	description	default value
`rho_min`		1.e6
`rho_max`		1.e10
`nrho`		5
`T_min`		1.e9
`T_max`		1.e10
`nT`		5
`Ye_min`		0.4
`Ye_max`		0.7
`nye`		7
`small_dens`		1.e5
`small_temp`		1.e5
`rho_min`		1.e6
`rho_max`		1.e10
`nrho`		5
`T_min`		4.0e9
`T_max`		1.e10
`nT`		10
`Ye_min`		0.4
`Ye_max`		0.7
`nye`		7

multigamma

parameter	description	default value
`eos_gamma_default`		1.4
`species_a_name`		“”
`species_a_gamma`		1.4
`species_b_name`		“”
`species_b_gamma`		1.4
`species_c_name`		“”
`species_c_gamma`		1.4

networks

parameter	description	default value
`small_x`	cutoff for species mass fractions	1.e-30
`use_tables`	Should we use rate tables if they are present in the network?	0
`use_c12ag_deboer17`	Should we use Deboer + 2017 rate for c12(a,g)o16?	0

nse_compatibility

parameter	description	default value
`run_prefix`		“”
`tmax`	the final time to integrate to	1.e3
`tfirst`	first output time – we will output in nsteps logarithmically spaced steps between [tfirst, tmax]	0.0
`nsteps`	number of steps (logarithmically spaced)	100
`rho_min`		1.e7
`rho_max`		1.e9
`nrho`		4
`T_min`		6.e9
`T_max`		8.e9
`nT`		4

nse_solver

parameter	description	default value
`max_nse_iters`		500
`use_hybrid_solver`		1
`ase_tol`		0.1
`nse_abs_tol`		0.005
`nse_rel_tol`		0.2
`nse_dx_independent`		0
`nse_molar_independent`		0
`nse_skip_molar`		0
`T_nse_net`		-1.0

nse_table_cell

parameter	description	default value
`small_temp`		1.e5
`small_dens`		1.e5
`density`		1.23e9
`temperature`		5.18e9
`ye`		0.472

nse_tabular

parameter	description	default value
`rho_nse`		3.e8
`T_nse`		3.e9
`C_nse`		0.01
`O_nse`		0.01
`Si_nse`		0.01
`He_Fe_nse`		0.88
`nse_relax_factor`	factor (< 1) by which to relax the criteria for entering NSE. This is only applied after a failed burn.	1.0
`nse_table_interp_linear`	do we do tri-linear or tri-cubic interpolation on the table?	0

polytrope

parameter	description	default value
`polytrope_type`		0
`polytrope_gamma`		0.0e0
`polytrope_K`		0.0e0
`polytrope_mu_e`		2.0e0

powerlaw

parameter	description	default value
`cond_coeff`		1.0
`cond_exponent`		1.0
`rtilde`	the coefficient for the reaction rate	1.e0
`nu`	exponent for the temperature	4.e0
`specific_q_burn`	reaction specific q-value (in erg/g)	10.e0
`T_burn_ref`	reaction thresholds (for the power law)	1.0e0
`rho_burn_ref`		1.0e0
`f_act`		1.0e0

primordial_chem

parameter	description	default value
`small_x`	cutoff for species mass fractions	1.e-100
`redshift`	assumed redshift for primordial chem (Pop III star formation)	30e0
`eos_gamma_default`		1.4
`species_a_name`	define the specie names, and their masses and gammas	“elec”
`species_a_gamma`		5./3.
`species_a_mass`		9.10938188e-28
`species_b_name`		“hp”
`species_b_gamma`		5./3.
`species_b_mass`		1.67262158e-24
`species_c_name`		“h”
`species_c_gamma`		5./3.
`species_c_mass`		1.67353251819e-24
`species_d_name`		“hm”
`species_d_gamma`		5./3.
`species_d_mass`		1.67444345638e-24
`species_e_name`		“dp”
`species_e_gamma`		5./3.
`species_e_mass`		3.34512158e-24
`species_f_name`		“d”
`species_f_gamma`		5./3.
`species_f_mass`		3.34603251819e-24
`species_g_name`		“h2p”
`species_g_gamma`		1.4
`species_g_mass`		3.34615409819e-24
`species_h_name`		“dm”
`species_h_gamma`		5./3.
`species_h_mass`		3.34694345638e-24
`species_i_name`		“h2”
`species_i_gamma`		1.4
`species_i_mass`		3.34706503638e-24
`species_j_name`		“hdp”
`species_j_gamma`		1.4
`species_j_mass`		5.01865409819e-24
`species_k_name`		“hd”
`species_k_gamma`		1.4
`species_k_mass`		5.01956503638e-24
`species_l_name`		“hepp”
`species_l_gamma`		5./3.
`species_l_mass`		6.69024316e-24
`species_m_name`		“hep”
`species_m_gamma`		5./3.
`species_m_mass`		6.69115409819e-24
`species_n_name`		“he”
`species_n_gamma`		5./3.
`species_n_mass`		6.69206503638e-24

rad_power_law

parameter	description	default value
`const_kappa_p`	Opacity constant (Planck)	-1.0e0
`kappa_p_exp_m`	Density exponent (Planck)	0.0e0
`kappa_p_exp_n`	Temperature exponent (Planck)	0.0e0
`kappa_p_exp_p`	Frequency exponent (Planck)	0.0e0
`const_kappa_r`	Opacity constant (Rosseland)	-1.0e0
`kappa_r_exp_m`	Density exponent (Rosseland)	0.0e0
`kappa_r_exp_n`	Temperature exponent (Rosseland)	0.0e0
`kappa_r_exp_p`	Frequency exponent (Rosseland)	0.0e0
`const_scatter`	Opacity constant (scattering)	0.0e0
`scatter_exp_m`	Density exponent (scattering)	0.0e0
`scatter_exp_n`	Temperature exponent (scattering)	0.0e0
`scatter_exp_p`	Frequency exponent (scattering)	0.0e0
`kappa_floor`	Opacity floor	0.0e0
`rad_temp_floor`	Temperature floor	1.e-10
`eos_const_c_v`		-1.e0
`eos_c_v_exp_m`		0.e0
`eos_c_v_exp_n`		0.e0

rprox

parameter	description	default value
`atol_spec`	override the default tolerances for backwards compatibility	1.0e-11
`rtol_spec`		1.0e-12
`atol_enuc`		1.0e-8
`rtol_enuc`		1.0e-8
`jacobian`	override so that the default is an analytical Jacobian	1

screening

parameter	description	default value
`enable_chabrier1998_quantum_corr`		0

subch_approx

parameter	description	default value
`disable_p_C12_to_N13`		0
`disable_He4_N13_to_p_O16`		0

subch_base

parameter	description	default value
`disable_p_C12_to_N13`		0
`disable_He4_N13_to_p_O16`		0

subch_full

parameter	description	default value
`disable_p_C12_to_N13`		0
`disable_He4_N13_to_p_O16`		0

subch_simple

parameter	description	default value
`disable_p_C12_to_N13`		0
`disable_He4_N13_to_p_O16`		0

test_aprox_rates

parameter	description	default value
`dens_min`		1.e6
`dens_max`		1.e9
`temp_min`		1.e6
`temp_max`		1.e12
`metalicity_max`		0.1e0
`small_temp`		1.e4
`small_dens`		1.e-4

test_ase

parameter	description	default value
`run_prefix`		“”
`small_temp`		1.e5
`small_dens`		1.e5
`density`		1.e7
`temperature`		3.e9
`ye`		0.75
`mu_p`		-3.0
`mu_n`		-12.0
`X1`		1.0e0
`X2`		0.0e0
`X3`		0.0e0
`X4`		0.0e0
`X5`		0.0e0
`X6`		0.0e0
`X7`		0.0e0
`X8`		0.0e0
`X9`		0.0e0
`X10`		0.0e0
`X11`		0.0e0
`X12`		0.0e0
`X13`		0.0e0
`X14`		0.0e0
`X15`		0.0e0
`X16`		0.0e0
`X17`		0.0e0
`X18`		0.0e0
`X19`		0.0e0
`X20`		0.0e0
`X21`		0.0e0

test_conductivity

parameter	description	default value
`dens_min`		1.e6
`dens_max`		1.e9
`temp_min`		1.e6
`temp_max`		1.e12
`metalicity_max`		0.1e0
`small_temp`		1.e4
`small_dens`		1.e-4

test_eos

parameter	description	default value
`dens_min`		1.e6
`dens_max`		1.e9
`temp_min`		1.e6
`temp_max`		1.e12
`metalicity_max`		0.1e0
`small_temp`		1.e4
`small_dens`		1.e-4

test_jac

parameter	description	default value
`dens_min`		1.e6
`dens_max`		1.e9
`temp_min`		1.e6
`temp_max`		1.e15
`uniform_xn`		0
`tmax`		0.1e0
`small_temp`		1.e5
`small_dens`		1.e5
`do_acc`		1

test_linear_algebra

parameter	description	default value
`small_temp`		1.e5
`small_dens`		1.e5

test_neutrino_cooling

parameter	description	default value
`dens_min`		1.e6
`dens_max`		1.e9
`temp_min`		1.e6
`temp_max`		1.e12
`metalicity_max`		0.1e0
`small_temp`		1.e4
`small_dens`		1.e-4

test_nse

parameter	description	default value
`run_prefix`		“”
`small_temp`		1.e5
`small_dens`		1.e5
`density`		1.e7
`temperature`		3.e9
`mu_p`		-3.0
`mu_n`		-12.0
`y_e`		0.5

test_nse_interp

parameter	description	default value
`small_temp`		1.e5
`small_dens`		1.e5
`density`		1.23e9
`temperature`		5.18e9
`ye`		0.472

test_part_func

parameter	description	default value
`small_temp`		1.e5
`small_dens`		1.e5
`density`		1.e9
`temperature`		5.e9

test_react

parameter	description	default value
`dens_min`		1.e6
`dens_max`		1.e9
`temp_min`		1.e6
`temp_max`		1.e15
`uniform_xn`		0
`tmax`		0.1e0
`small_temp`		1.e5
`small_dens`		1.e5
`do_acc`		1

test_rhs

parameter	description	default value
`dens_min`		1.e6
`dens_max`		1.e9
`temp_min`		1.e6
`temp_max`		1.e15
`uniform_xn`		0
`small_temp`		1.e5
`small_dens`		1.e5

test_screening

parameter	description	default value
`dens_min`		1.e6
`dens_max`		1.e9
`temp_min`		1.e6
`temp_max`		1.e12
`metalicity_max`		0.1e0
`small_temp`		1.e4
`small_dens`		1.e-4

test_screening_templated

parameter	description	default value
`dens_min`		1.e6
`dens_max`		1.e9
`temp_min`		1.e6
`temp_max`		1.e12
`metalicity_max`		0.1e0
`small_temp`		1.e4
`small_dens`		1.e-4

test_sdc

parameter	description	default value
`dens_min`		1.e6
`dens_max`		1.e9
`temp_min`		1.e6
`temp_max`		1.e15
`uniform_xn`		0
`tmax`		0.1e0
`small_temp`		1.e5
`small_dens`		1.e3
`do_acc`		1

test_sdc_vode_rhs

parameter	description	default value
`small_temp`		1.e5
`small_dens`		1.e5
`density`		1.e7
`temperature`		3.e9
`tmax`		1.e-7
`Aux1`		0.0e0
`Aux2`		0.0e0
`Aux3`		0.0e0
`Adv_rho`		0.0e0
`Adv_rhoe`		0.0e0
`Adv_X1`		0.0e0
`Adv_X2`		0.0e0
`Adv_X3`		0.0e0
`Adv_X4`		0.0e0
`Adv_X5`		0.0e0
`Adv_X6`		0.0e0
`Adv_X7`		0.0e0
`Adv_X8`		0.0e0
`Adv_X9`		0.0e0
`Adv_X10`		0.0e0
`Adv_X11`		0.0e0
`Adv_X12`		0.0e0
`Adv_X13`		0.0e0
`Adv_X14`		0.0e0
`Adv_X15`		0.0e0
`Adv_X16`		0.0e0
`Adv_X17`		0.0e0
`Adv_X18`		0.0e0
`Adv_X19`		0.0e0
`Adv_X20`		0.0e0
`Adv_X21`		0.0e0
`Adv_X22`		0.0e0
`Adv_X23`		0.0e0
`Adv_X24`		0.0e0
`Adv_X25`		0.0e0
`Adv_X26`		0.0e0
`Adv_X27`		0.0e0
`Adv_X28`		0.0e0
`Adv_X29`		0.0e0
`Adv_X30`		0.0e0
`Adv_Aux1`		0.0e0
`Adv_Aux2`		0.0e0
`Adv_Aux3`		0.0e0

tillotson

parameter	description	default value
`eos_la`		0.5
`eos_lb`		1.3
`eos_e_0`		1.6e11
`eos_rho_0`		2.7
`eos_A`		1.8e11
`eos_B`		1.8e11
`eos_e_s`		3.5e10
`eos_e_s_prime`		1.8e11
`eos_alpha`		5.0
`eos_beta`		5.0
`eos_c_v`		7.9e6

triple_alpha_plus_cago

parameter	description	default value
`atol_spec`	override the default tolerances for backwards compatibility	1.0e-12
`rtol_spec`		1.0e-12
`atol_enuc`		1.0e-8
`rtol_enuc`		1.0e-6
`jacobian`	override so that the default is an analytical Jacobian	1

unit_test

parameter	description	default value
`primary_species_1`		“”
`primary_species_2`		“”
`primary_species_3`		“”
`X1`		1.0e0
`X2`		0.0e0
`X3`		0.0e0
`X4`		0.0e0
`X5`		0.0e0
`X6`		0.0e0
`X7`		0.0e0
`X8`		0.0e0
`X9`		0.0e0
`X10`		0.0e0
`X11`		0.0e0
`X12`		0.0e0
`X13`		0.0e0
`X14`		0.0e0
`X15`		0.0e0
`X16`		0.0e0
`X17`		0.0e0
`X18`		0.0e0
`X19`		0.0e0
`X20`		0.0e0
`X21`		0.0e0
`X22`		0.0e0
`X23`		0.0e0
`X24`		0.0e0
`X25`		0.0e0
`X26`		0.0e0
`X27`		0.0e0
`X28`		0.0e0
`X29`		0.0e0
`X30`		0.0e0
`X31`		0.0e0
`X32`		0.0e0
`X33`		0.0e0
`X34`		0.0e0
`X35`		0.0e0