List of Parameters by Number

Number	Name	Format	Default	Unit	Synopsis
1	dtnew	float	1.	sec	Initial or current timestep.
2	timesec	float	0.	sec	Initial or current time.
3	extrap	float	1.	-	Extrapolation Parameter used in guessing new , , and values.
4	q2fac	float	4.	-	Quadratic artificial viscosity () factor.
5	maxit	integer	40	-	Maximum number of times subroutine regess is called to iterate R, T, and L estimates before the timestep is reduced by dtcut (p 53) and the step repeated (termed a backup).
6	dtcr	float	0.05	-	Maximum desired fractional change in radius per step.
7	dtct	float	0.05	-	Maximum desired fractional change in temperature per step.
8	dtcd	float	0.1	-	Maximum desired fractional change in density per step.
9	dtcq	float	0.05	-	Maximum desired fractional linear contraction per step.
10	dtcdt	float	0.99	-	Maximum fractional change in the timestep per step.
11	fcrmax	float	1.e-6	-	Maximum allowed relative convergence error in radius.
12	fctmax	float	1.e-6	-	Maximum allowed relative convergence error in temperature.
13	q1fac	float	0.1	-	Linear artificial viscosity factor.
14	nstop	integer	100000	-	Maximum number of cycles.
15	tstop	float	1.e+30	sec	Stop time.
16	nedit	integer	20	-	Number of cycles between ASCII edits.
17	dtedit	float	0.	sec	Time between ASCII edits (0. = ).
18	ndump	integer	10	-	Number of cycles between restart dumps.
19	xmlen	float	1.	-	Ratio of the convective mixing length to the pressure scale height.
20	fudgc	float	0.01	-	Fudge factor for convection (about 0.01).
21	difim	float	1.	-	Multiplier for the rate of convective mixing.
22	fcrbu	float	0.01	-	Reduce timestep by dtcut (p 53) and redo step if convergence error in radius is still greater than fcrbu (p 22) after maxit (p 5) iterations.
23	fctbu	float	0.01	-	Reduce timestep by dtcut (p 53) and redo step if convergence error in temperature is still greater than fctbu (p 23) after maxit (p 5) iterations.
24	drmult	float	0.1	-	Semi-convective mixing will be slower than thermal transport by at least drmult (p 24) (about 0.1) in zones with mean atomic weight, , is abarsemi (p 324).
25	dtmax	float	1.e+99	sec	Maximum timestep allowed.
26	fcrext	float	10.	-	Factor used to reduce noise in radius extrapolation.
27	fctext	float	10.	-	Factor used to reduce noise in temperature extrapolation.
28	ipup	integer	1	-	Abundance update parameter. Update abundances only if ipup (p 28) :math:’ne’ 0.
29	xkimt	float	1.	-	Multiplier on iben1 opacity.
30	xk2mt	float	1.	-	Multiplier on iben2 opacity.
31	xk3mt	float	1.	-	Multiplier on Christy opacity.
32	xk4mt	float	1.	-	Multiplier on Compton opacity.
33	rxkcmt	float	1.	-	Multiplier on conductive opacity.
34	pimult	float	1.	-	Ion energy and pressure multiplier.
35	prmult	float	1.	-	Radiation energy and pressure multiplier.
36	pemult	float	1.	-	Electron energy and pressure multiplier.
37	transm	float	1.	-	Multiplier on diffusive heat transport.
38	timex0	float	-1.e+99	sec	Time at which any ni56 present in low temperature regions (as determined by dypmin (p 89) and tnucmin (p 65)) of the star is assumed to have been produced and begun to decay, first to and then to .
39	dyemult	float	200.	-	If iytsflag (p 67) 1, increase the timestep sensitivity to changes in ye by a factor of dyemult (p 39).
40	dyqmult	float	2.5	-	If iytsflag (p 67) 1, increase the timestep sensitivity to changes in yq by a factor of dyqmult (p 40).
41	yemin	float	0.02	-	No further changes in ye are allowed for ye yemin (p 41) :unit: mol/g.
42	iwinsize	integer	10240750	-	Size of graphics window to be created in the form xxxxyyyy.
43	hstatm	float	1.	-	Multiplier on the inertial terms in the momentum balance equation.
44	lenqmax	integer	3000000	bytes	Maximum length of a post-processor dump file.
45	iwinloc	integer	0	-	Location of the top-left-hand corner of the graphics window relative to top-left corner of screen, given in the form xxxxyyyy.
46	dtcp	float	0.1	-	Maximum desired fractional change in abundances per step.
47	yfloorx	float	0.001	-	Minimum elemental mass fraction that effects the timestep.
48	cenu	float	1.	-	Velocity centring parameter (Range: 0.5 - 1.).
49	t7peek	float	1.e+50	cgs	Opacity will be no larger than xkmin (p 50) + t7peek (p 49) * .
50	xkmin	float	1.e-10	cm**2/g	Least upper opacity bound.
51	thickfac	float	0.001	-	Determine whether zones are considered in thermal equilibrium depending on time step and optical thickness.
52	maxbak	integer	5	-	Maximum number of times a given step is redone before the code quits.
53	dtcut	float	0.1	-	Fractional timestep reduction when a step is redone.
54	tfcrbu	float	2.	-	If the maximum fractional change in radius during a timestep exceeds tfcrbu (p 54) dtcr (p 6) then redo step.
55	tfctbu	float	2.	-	If the maximum fractional change in temperature during a timestep exceeds tfctbu (p 55) dtct (p 7) then redo step.
56	xipot	float	13.6053	eV	Ionization potential.
57	dzero	float	0.1	g/cc	Characteristic density for pressure ionization.
58	npflag	integer	1	-	Pairs are included in EOS calculation only if npflag (p 58) 0.
59	fracneut	float	0.05	-	If the semi convective test parameter, , is 0. but greater than - fracneut (p 59) * abs(log(T1/T0)), then the zonal interface is flagged convectively neutral (“NEUT” or “,”).
60	radius0	float	0.	cm	Radius of inner boundary.
61	summ0	float	0.	g	Mass inside inner boundary.
62	xlum0	float	0.	erg/sec	Luminosity emerging from inner boundary.
63	irnet	integer	0	-	Network switch parameter.
64	npixedit	integer	50	-	Graphics edits to the monitor are made every npixedit (p 64) KEPLER cycles (if itvstart (p 127) 0).
65	tnucmin	float	0.	K	Don’t calculate nuclear burning in APPROX if the temperature is less than tnucmin (p 65).
66	nsetparm	integer	170001	-	Initialize new parameters in subroutine restart and re-arrange/add/remove/initialize internal arrays and edit quantities if nsetparm (p 66) current version.
67	iytsflag	integer	5	-	Consider only specific changes in ye, yq, yf, o16 and other select species in ISE zones when computing the timestep.
68	tbound	float	0.	K	Temperature at outer boundary.
69	pbound	float	0.	erg/cc	Pressure at outer boundary.
70	fclmax	float	1.	-	Maximum allowed relative convergence error in luminosity.
71	fclext	float	1.e-5	-	Factor used to reduce noise in luminosity extrapolation.
72	dtcl	float	1.e+99	-	Maximum desired fractional change in luminosity per step.
73	fclbu	float	1.	-	Reduce timestep by dtcut (p 53) and redo step if convergence error in luminosity is still greater than fclbu (p 73) after maxit (p 5) iterations.
74	tfclbu	float	1.e+99	-	If the maximum fractional change in luminosity during a timestep exceeds tfclbu (p 74) dtcl (p 72) then redo step.
75	dtsmult	float	1.e+99	-	The fractional amount of semi convective mixing that can occur in one timestep is limited to approximately dtsmult (p 75).
76	rnratmax	float	0.25	-	Maximum fractional radius change allowed between zones before adzoning.
77	rnratmin	float	0.1	-	Minimum fractional radius change allowed between zones before dezoning.
78	tnratmax	float	0.25	-	Maximum fractional temperature change allowed between zones before adzoning.
79	tnratmin	float	0.1	-	Minimum fractional temperature change allowed between zones before dezoning.
80	dnratmax	float	0.25	-	Maximum fractional density change allowed between zones before adzoning.
81	dnratmin	float	0.1	-	Minimum fractional density change allowed between zones before dezoning.
82	rnmin	float	1.e+99	cm	Minimum radius for which adzoning is considered.
83	tnmin	float	1.	K	Minimum temperature for which adzoning is considered.
84	dnmin	float	1.e-99	g/cc	Minimum density for which adzoning is considered.
85	jmmin	integer	10	-	Minimum number of zones allowed after dezoning.
86	izonef	integer	1	-	General rezoning flag.
87	idzonef	integer	1	-	Dezoning Flags.
88	rnmax	float	1.e+99	cm	Maximum radius for which rezoning is considered.
89	dypmin	float	1.e-10	mol/g	Minimum amount of potential hydrogen burn in the present timestep needed to trigger the use of any reaction network.
90	dymgmin	float	1.e-10	mol/g	Minimum amount of Mg24 production required in the last timestep to trigger the use of the full reaction network.
91	zbound	float	0.1	-	Mass fraction of heavy elements above which iben1 opacities are used.
92	etaconv	float	1.e-5	-	Relative convergence required in calculating the electron Fermi degeneracy parameter, .
93	jshell0	integer	1	-	Innermost zone in which there is neutrino deposition.
94	jshell1	integer	1000	-	Outermost zone in which there is neutrino deposition.
95	eexplode	float	0.	erg	Total neutrino deposition energy.
96	texplode	float	1.e+99	sec	Time of neutrino deposition.
97	tauexp	float	0.01	sec	Time scale for neutrino deposition.
98	p0core	float	0.	cgs	Artificial neutrino core pressure parameter.
99	d0core	float	3.e+11	g/cc	Core density cutoff.
100	d0power	float	1.	-	Core pressure density dependence exponent.
101	sneutmt	float	1.	-	Neutrino energy loss rate multiplier (APPROX only).
102	snucmt	float	1.	-	Non-neutrino nuclear energy generation rate multiplier (APPROX only).
103	dsnum	float	1.e-6	-	Fractional density and temperature perturbations used to get nuclear energy generation rate derivatives.
104	tcorefac	float	1.	-	Increase the rezoner’s sensitivity to temperature gradients by a factor of tcorefac (p 104) in the region where fracm, the cumulative mass fraction measured from the center of the star, is less than or equal fracrz1 (p 138).
105	tqsemin	float	2.e+9	K	Floor on the temperature used in the ISE calculation.
106	jqse	integer	0	-	Change all APPROX network zones with j jqse (p 106) to the ISE network.
107	artv1	float	1.	-	Multiply q1fac (p 13) by artv1 (p 107) in zone j if jartv1 (p 110) j artv2 (p 108).
108	artv2	float	1.	-	Multiply q1fac (p 13) by artv2 (p 108) in zone j if jartv2 (p 111) j jm - jartv3 (p 112).
109	artv3	float	1.	-	Multiply q1fac (p 13) by artv3 (p 109) in zone j if j max(jartv1 (p 110), jartv2 (p 111), jm - jartv3 (p 112)).
110	jartv1	float	1000.	-	See artv1 (p 107) - artv3 (p 109).
111	jartv2	integer	1000	-	See artv1 (p 107) - artv3 (p 109).
112	jartv3	integer	1000	-	See artv1 (p 107) - artv3 (p 109).
113	ipixtype	integer	31	-	Graphics picture-type control parameter:.
114	nplotsym	integer	8	-	Label each curve with nplotsym (p 114) character symbols.
115	tscalem	float	1.e+7	K	Characteristic temperature for choosing the location of the temperature plotting grid in Thermodynamics (TD) graph.
116	igridt	integer	370	-	Grid choice for the radial coordinate (irtype (p 132) 1).
117	igridm	integer	498	-	Obsolete.
118	igridv	integer	496	-	Obsolete.
119	jp0	integer	1	-	Innermost zone to plot.
120	jp1	integer	0	-	Outermost zone to plot is jm (outermost zone) - jp1 (p 120).
121	rpmino	float	1.e+7	cm	Minimum radius plotted when jp0 (p 119) 0.
122	rpmaxo	float	1.e+14	cm	Maximum radius plotted when jp1 (p 120) 0.
123	ymintd	float	1.	g/cc	Lower bound of thermodynamics graph ordinate shall be no greater than ymintd (p 123).
124	ymaxtd	float	1.e+8	g/cc	Upper bound of thermodynamics graph ordinate shall be no less than ymaxtd (p 124).
125	sscalem	float	0.001	-	Energy generation rate graph scale multiplier in the TD graph.
126	pscalem	float	0.1	-	Pressure graph scale multiplier in the TD graph.
127	itvstart	integer	0	-	Graphics device control parameter.
128	abunlim	float	0.001	-	Least elemental mass fraction plotted or listed in a terminal ion edit.
129	ivplot	integer	1	-	Obsolete.
130	vscalem	float	1.e-5	-	Velocity graph scale multiplier in the TD graph.
131	rscalem	float	1.e-5	-	Radius graph scale multiplier in the TD graph.
132	irtype	integer	3	-	-axix type for plots.
133	convord	float	3.	g/cc	Ordinate value at which to plot convection sentinels in the TD graph.
134	yplotmin	float	0.	-	Minimum fractional mass coordinate () plotted (fraction of total mass).
135	yplotmax	float	1.	-	Maximum fractional mass coordinate () plotted (fraction of total mass).
136	etacut	float	10.	-	Electron degeneracy parameter, , above which to use temperature interpolation in subroutine adzone.
137	nadapb	integer	0	-	Enable Adaptive BURN network adjustment.
138	fracrz1	float	1.	-	Multiplier in the effective values of the density, temperature, and radius gradients used to determine the necessity for adzoning or dezoning.
139	fracrz2	float	1.	-	See fracrz1 (p 138) and fmax0 (p 150).
140	rzmult0	float	1.	-	See fracrz1 (p 138).
141	rzmult1	float	1.	-	See fracrz1 (p 138).
142	rzmult2	float	1.	-	See fracrz1 (p 138).
143	abarrat0	float	2.	-	Limiting ratio for for rezoning.
144	abarratm	float	1.3	-	See abarrat0 (p 143).
145	fracmlim	float	0.01	-	See abarrat0 (p 143).
146	frcsound	float	0.1	-	Don’t do convection if the absolute value of the zone velocity exceeds frcsound (p 146) times the local sound speed.
147	convlim	float	1.	-	Limit the convective velocity to a fraction convlim (p 147) of the local sound speed.
148	woversht	float	0.01	-	The semiconvective test parameter, , is taken to be woversht (p 148) * abs ( log (T1/T0)) for the special overshoot semiconvective zones where would otherwise be 0 and when abarsemi (p 324).
149	xmimult	float	1.	-	The mass used in calculating ion degeneracy is xmimult (p 149) times the mass of a neutron.
150	fmax0	float	1.	-	Adzone mass fraction parameter.
151	fmax1	float	1.	-	See fmax0 (p 150).
152	fmax2	float	1.	-	See fmax0 (p 150).
153	ddsfrac	float	0.01	-	Maximum fractional change in density since the last iteration for which changes and energy generation rates can be extrapolated from their previous values (APPROX only).
154	dtsfrac	float	0.001	-	Maximum fractional change in temperature since the last iteration for which ABAR changes and energy generation rates can be extrapolated from their previous values (APPROX only).
155	iudflag	integer	0	-	Extrapolate energy generation rates and changes in cases of small temperature and density change only if iudflag (p 155) is also 0.
156	nsdump	integer	10	-	Save every nsdump (p 156) restart dumps.
157	niterbar	integer	1	-	Total number of cycles used to compute iterbar, the average number of iterations per cycle. (Internally set.).
158	iterbarm	integer	100	-	Maximum allowed value of niterbar without terminating the problem.
159	iautoout	integer	0	-	Autumated ASCII and graphical output.
160	iflgabar	integer	0	-	The mean atomic weight, , calculated in subroutine sdot is implicitly coupled to the ion equation of state provided iflgabar (p 160) bit 0 (value 1) is set and the ISE or NSE network is being used; or iflgabar (p 160) bit 1 is set and the APPROX network is being used.
161	fraccore	float	1.1	-	The effective number of zones used in regess (jmcalc (p 162)) shall be the minimum needed to contain the inner fraccore (p 161) mass-fraction of the initial mass of the star.
162	jmcalc	integer	jm	-	Effective number of zones used in subroutine regess.
163	iupdflag	integer	1	-	Don’t calculate derivatives in subroutine sdot when called from subroutine update if iupdflag (p 163) 0.
164	niondump	integer	40	-	Number of dump cycles between forced dumps of all dump-grid points of all ‘ions’ specified as dump variables.
165	nisodump	integer	400	-	Number of dump cycles between forced dumps of all dump-grid points of all BURN ‘isotopes’ specified as dump variables.
166	nzondump	integer	40	-	Number of dump cycles between forced dumps of all dump-grid points of all arrays specified as dump variables except ‘ions’ and ‘isotopes’.
167	iflagye	integer	1	-	The electron EOS is implicitly coupled to nuclear burning-induced changes in the electron abundance, , provided iflagye (p 167) bit 0 (value 1) is set and the ISE or NSE network is being used; or iflagye (p 167) bit 1 is set and the APPROX network is being used.
168	njeditq	integer	5	-	Make a ISE edit every njeditq (p 168) zones.
169	dtqnum	float	0.0001	-	Relative temperature change used in calculating numerical derivatives in sdotq.
170	ddqnum	float	0.0001	-	Relative density change used in calculating numerical derivatives in sdotq.
171	iexciteh	integer	1	-	Include excited states in ISE calculations only if iexciteh (p 171) 0.
172	iterqmh	integer	1000	-	Maximum number of iterations allowed in ISE calculation.
173	ypconvh	float	1.e-8	-	Allowed relative convergence error in the proton abundance in the ISE calculation.
174	ynconvh	float	1.e-8	-	Allowed relative convergence error in the neutron abundance in the ISE calculation.
175	ysiconvh	float	1.e-8	-	Allowed relative convergence error in the Si28 abundance in the ISE calculation.
176	cnseh	float	1.	-	Increment the proton, neutron, and abundances by a fraction cnseh (p 176) of that calculated by the Newton-Raphson method for an ISE iteration.
177	fyph	float	0.05	-	Maximum allowed relative change in proton abundance during an ISE iteration.
178	fynh	float	0.05	-	Maximum allo~ed relative change in neutron abundance during an ISE iteration.
179	fysih	float	0.15	-	Maximum allowed relative change in Si28 abundance during an ISE iteration.
180	neditq	integer	20	-	Make a detailed edit of the ISE zones every neditq (p 180) general numerical edits.
181	xiter1qe	float	0.5	-	If the relative change of the proton, neutron, or abundance is opposite in sign and more than xiter1qe (p 181) in magnitude with respect to the corresponding change during the previous ISE iteration cycle, then cut the current step size in half if more than iter1qe (p 182) iterations have been done.
182	iter1qe	integer	10	-	See xiter1qe (p 181).
183	xthres	float	0.0001	-	Edit only those ISE isotopes with mass fractions exceeding xthres (p 183).
184	tqselim	float	1.5e+9	K	A sufficient condition to change a zone from the APPROX to ISE network is if its temperature exceeds tqselim (p 184), its o16 mass fraction is less than o16lim (p 185), its iron peak mass fraction exceeds qn56lim (p 186), and its density exceeds dqselim (p 203).
185	o16lim	float	0.04	-	See tqselim (p 184).
186	qn56lim	float	0.	-	See tqselim (p 184).
187	snuwmult	float	1.	-	Multiplier on neutrino energy losses from weak processes on nuclei and nucleons.
188	jnse	integer	0	-	Change all ISE zones with j jnse (p 188) to the NSE network.
189	siqselim	float	0.001	-	A sufficient condition to change a zone from the ISE to the NSE network is for the sum of the silicon and sulphur “group” elemental mass fractions to be less than or equal to siqselim (p 189).
190	jlcalc	integer	0	-	Remove the inner jlcalc (p 190) zones from the problem and reset the inner boundary conditions on radius, radius0 (p 60), and mass, summ0 (p 61), but don’t change the central luminosity, xlum0 (p 62).
191	vlimset	float	-1.	-	The ordinate bounds of the velocity graph are vlimset (p 191) if vlimset (p 191) 0. Otherwise the ordinate bounds are rescaled automatically. :unit: cm/s.
192	jpause	integer	-1	-	Pause the code at the end of the calculation for zone j = jpause (p 192) in regess (”pause 1”) and in update (”pause 2”). Do not pause if jpause (p 192) < 0. Type a carriage return to continue.
193	fracrz0	float	-1.	-	See fracrz1 (p 138) and fmax0 (p 150).
194	rzmultm	float	1.	-	See fracrz1 (p 138).
195	fmaxm	float	1.	-	See fmax0 (p 150).
196	eionmult	float	1.	-	Multiplier on the contribution of ionization potential energy to the equation of state. See subroutine ES.
197	nnewoutf	integer	2000	-	Start a new labelled ASCII output file every nnewoutf (p 197) cycles.
198	neditq1	integer	5	-	Make an ISE edit for the central zone every neditq1 (p 198) nedit (p 16) cycles.
199	wilsonmt	float	-1.	-	Multiplier on the Wilson-based nuclear EOS (except for thee thermal ion component) if it is 0. Otherwise, the old non-relativistic, partial degeneracy model for the ion EOS is used.
200	t11cut	float	1.	1.e+11 K	Upper bound on the temperature used to calculate the energy in nuclear excited states and nuclear partition functions.
201	y56gessm	float	0.0001	mol/g	Initial guess made in sdot for the y56 abundance when initializing a new ISE zone.
202	isi30brn	integer	1	-	Neutron-rich-silicon-burning flows are included in the ISE network if isi30brn (p 202) 1.
203	dqselim	float	100000.	g/cc	Nominal minimum density for changing a zone from the APPROX to the ISE network.
204	abunminx	float	-1.e-5	-	Force a backup after a call to subroutine BURN from subroutine sdot if a calculated elemental mass fraction change would cause the resulting mass fraction to be less than abunminx (p 204) or to change by a fractional amount greater than tfcybu (p 205).
205	tfcybu	float	0.5	-	Maximim allowed fractional change in abundance before an abundance backup is made. (See abunminx (p 204)).
206	yfloorbx	float	0.001	-	Elemental mass fraction floor for making abundance backups. (See abunminx (p 204)).
207	iqerrflg	integer	1	-	Fatal errors in qburn encountered during a normal regess iteration will cause KEPLER to back-up if iqerrflg (p 207) 1, or to terminate if iqerrflg (p 207) 1.
208	c12agmlt	float	1.	-	Multiplier on the overall C12(a,g) rate (APPROX only).
209	xltaucon	float	0.	sec	Limit how fast convective velocity can increase.
210	noqsecon	integer	0	-	Flags regulating transition to QSE for convectively coupled zones.
211	accrate	float	0.	Msun/yr	Rate at which mass in the form of new zones is added to the surface of the star.
212	xmacrete	float	0.	g	Mass of phantom outer zone used to mediate mass accretion. See accrate (p 211).
213	denconv	float	1.e-7	-	Allowed fractional convergence error in density when calculating a hydrostatic initial stellar configuration in the generator using the dstat option.
214	flamerad	float	0.	cm	Characteristic flame radius for carbon deflagration studies.
215	coulmult	float	1.	-	Multiplier on the Wigner-Seitz Coulomb corrections to the ion energy and pressure (See [Cla], p. 152, and subroutine es).
216	irzopt	integer	0	-	Rezoning Option Flag.
217	nsurfz	integer	-10	-	If nsurfz (p 217) 0 and izonef (p 86) 0, maintain a logarithmic ramp in zonal mass of the NSURFZ zones at the surface of the star, where fracsz0 (p 218) is twice the desired mass-fraction of the surface zone and fracsz1 (p 219) is twice the desired mass-fraction of the nsurfz (p 217)th zone from the surface.
218	fracsz0	float	1.	-	Inner surface zoning parameter.
219	fracsz1	float	1.	-	Outer surface zoning parameter.
220	xmlossm	float	0.	-	Multiplier on the de Jager mass-loss rate.
221	xmloss0	float	0.	Msun/yr	Nominal mass loss rate from the surface of the star.
222	totm0	float	0.	g	The original total mass of the star.
223	fracdez	float	0.5	-	Dezone the mass-losing zone if its mass drops below fracdez (p 223) times the average mass of the zones on either side of it.
224	xmratbak	float	0.2	-	Redo the timestep (“backup”) if the fractional change of mass in the mass-losing zone exceeds xmratbak (p 224).
225	xfracml	float	0.01	-	Subtract any “surface” mass loss specified by xmlossm (p 220) or xmloss0 (p 221) from the zone closest to the surface that still has more than fracs grams of material overlying it.
226	mlcompf	integer	0	-	Obsolete.
227	e1mltc12	float	1.	-	Multiplier on the E1 part of the C12(a,g) cross-section.
228	e2mltc12	float	1.	-	Multiplier on the E2 part of the C12(a,g) cross-section.
229	nupdate	integer	100	-	Do a forced update of BURN isotopic abundances in every zone every nupdate (p 229) cycles.
230	dtfrac	float	0.01	-	When BURN co-processing is initiated, or when the zerotime or reset commands are given, (re)set the zonal timesteps used by the BURN coprocessor to dtfrac (p 230) * dtnew (p 1).
231	bmassmin	float	-1.	g	BURN co-processing is skipped if a zone’s exterior mass coordinate is less than bmassmin (p 231).
232	bmassmax	float	1.e+99	g	BURN co-processing is skipped if a zone’s exterior mass coordinate is greater than bmassmax (p 232).
233	btempmin	float	1.e+6	K	BURN co-processing is skipped if a zone’s temperature is less than btempmin (p 233).
234	btempmax	float	1.e+99	K	BURN co-processing is skipped if a zone’s temperature is greater than btempmax (p 234).
235	snucmin	float	1.e-99	erg/g/sec	BURN co-processing is skipped if the absolute value of a zone’s normal nuclear energy generation rate (snn) is less than snucmin (p 235).
236	bdenmin	float	1.e-99	g/cc	BURN co-processing is skipped if a zone’s density is less than bdenmin (p 236).
237	bdenmax	float	1.e+99	g/cc	BURN co-processing is skipped if a zone’s density is greater than bdenmax (p 237).
238	tchange	float	0.02	-	Force a zone to be updated by the BURN coprocessor if its fractional temperature change since its last BURN processing exceeds tchange (p 238).
239	dchange	float	0.05	-	Force a zone to be updated by the BURN coprocessor if its fractional density change since its last BURN processing exceeds dchange (p 239).
240	netmax	integer	1	-	BURN co-processing is skipped if a zone’s network number, netnum, is greater than netmax (p 240).
241	neditb	integer	5	-	Make a BURN isotopic abundance edit every neditb (p 241) nedit (p 16) KEPLER cycles.
242	edmassl	float	0.0001	-	Obsolete.
243	nedita	integer	10	-	Make an elemental abundance edit every nedita (p 243) nedit (p 16) cycles.
244	jmeditb	integer	1	-	Edit only zones whose zone number is jmeditb (p 244) during normal BURN isotopic abundance edits.
245	neditall	integer	100	-	Make a BURN isotopic abundance edit for all zones every neditall (p 245) nedit (p 16) cycles regardless of the value of neditb (p 241).
246	chimin	float	1.e-6	mol/g	Minimum isotopic abundance that affects the calculation of the zonal timestep in the BURN coprocessor.
247	delchi	float	0.15	-	Maximum desired fractional change of an isotopic abundance used in determining the zonal timestep used in the BURN coprocessor.
248	fdtn	float	2.	-	Maximum factor by which the zonal timestep in the BURN coprocessor can be increased in one zonal cycle.
249	dtbkup	float	0.5	-	Back up to the previous zonal cycle in the BURN coprocessor if the new zonal timestep is less than dtbkup (p 249) times the previous timestep.
250	ncomp	integer	999999	-	Make a complete ASCII edit of the nuclear processes in all zones considered by the BURN coprocessor every ncomp (p 250) cycles.
251	ncent	integer	999999	-	Make an ASCII edit of the nuclear processes in the central zone considered by the BURN coprocessor every ncent (p 251) cycles.
252	nedt	integer	999999	-	Make an ASCII edit of the nuclear processes involved every nedt (p 252) BURN coprocessor matrix inversions (debugging only).
253	nzro	integer	1	-	Set newly calculated negative BURN isotope abundances to 0. if nzro (p 253) 0.
254	amaglim	float	0.	-	Minimum absolute magnitude of a BURN matrix element for which Gauss-Jordan elimination is carried out in solving for changes in isotopic abundances.
255	ninv	integer	0	-	Total number of BURN-coprocessor matrix inversions so far (internally incremented).
256	nbkup	integer	0	-	Total number of BURN-coprocessor backups so far (internally incremented).
257	ntty	integer	6	-	FORTRAN I/O unit number for zburn messages sent to the user’s terminal.
258	ileqs	integer	1	-	Flag to switch on sparse matrix solver.
259	negbkup	integer	0	-	Cumulative number of BURN coprocessor backups due to encountering negative isotopic abundances (incremented internally).
260	bkupdiv	float	10.	-	Factor by which to reduce the BURN coprocessor zonal timestep when a negative isotopic abundance backup occurs.
261	bkupmass	float	1.e-13	mol/g	Minimum absolute magnitude ofa negative isotopic abundance (except of protons) that can cause a zonal backup in the BURN coprocessor.
262	bkupmp	float	1.e-19	mol/g	Minimum absolute magnitude of a negative proton abundance that can cause a zonal backup in the BURN coprocessor.
263	nbkupmax	integer	10	-	Maximum number of consecutive negative isotopic abundance backups allowed in the BURN coprocessor before quitting.
264	nsubcycm	integer	1	-	Number of coupled BURN co-processing / BURN isotope convection subcycles per KEPLER cycle.
265	mazful	integer	1	-	Use Fuller et. al.’s weak rates in the BURN coprocessor if mazful (p 265) 1, otherwise use the old rates of Mazurek and Hansen.
266	al26mult	float	100.	-	Decrease the abundance threshold for Al26 to affect the BURN coprocessor zonal timestep from chimin (p 246) to chimin (p 246)/ al26mult (p 266).
267	inburn	integer	0	-	BURN coprocessing and related edits are done only if inburn (p 267) 0.
268	npage	integer	53	-	Number of lines printed per “page” of ASCII output.
269	nsavez	integer	43	-	Total number of non-BURN zonal arrays to save in restart dumps.
270	nsaveb	integer	12	-	Total number of BURN zonal arrays to save in restart dumps.
271	vloss	float	1.e+99	cm/sec	Remove the outer zone if its velocity exceeds vloss (p 271) and its radius is larger than rlossmin (p 436), but do not change the previous values of pbound (p 69) and tbound (p 68).
272	abunlimb	float	1.e-5	-	Edit only those BURN isotope mass-fractions whose values exceed abunlimb (p 272) in making terminal edits (including isotopic mass-fraction sums over zones).
273	scalem	float	1.9892e+33	g	Mass unit used for the mass coordinate employed in making ASCII and terminal edits.
274	neditz1	integer	10	-	Make an ascii edit of primary zonal quantities every neditz1 (p 274) nedit (p 16) cycles.
275	neditz2	integer	20	-	Make an ascii edit of secondary zonal quantities every neditz2 (p 275) nedit (p 16) cycles.
276	medit	integer	0	-	Flag determining the minimum amount of information printed in an ASCII cycle edit regardless of the settings of other edit parameters:.
277	meditfin	integer	7	-	Effective value of medit (p 276) used in determining the scope of the final edit made when the problem is finished.
278	neditp	integer	50	-	Make an ascii edit of all changeable parameters every neditp (p 278) nedit (p 16) cycles.
279	c12flame	float	0.	-	Carbon mass fraction ahead of flame.
280	o16flame	float	0.	-	Oxygen mass fraction for ne-ox problems.
281	timenew	float	0.	sec	When the flame became RT unstable.
282	sharp1	float	0.	-	Sharp1 parameter.
283	sharpr	float	1.	cm/sec	Minimum speed from convection.
284	sharp2	float	0.	-	Sharp2 parameter.
285	newflam	integer	0	-	(Edit only.).
286	taunu	float	0.	sec	Time scale for the neutrino pulse considered by the BURN coprocessor in calculating neutrino-induced nucleosynthesis after core collapse.
287	enu53	float	3.	1.e+53 erg	Total energy for the neutrino pulse considered by the BURN coprocessor in calculating neutrino- induced nucleosynthesis after core collapse.
288	tmunu	float	8.	MeV	Temperature of the and neutrinos in the neutrino pulse considered by the BURN coprocessor in calculating neutrino- induced nucleosynthesis after core collapse.
289	tenu	float	4.	MeV	Temperature of the electron neutrinos in the neutrino pulse considered by the BURN coprocessor in calculating neutrino-induced nucleosynthesis after core collapse.
290	noiland	integer	3	-	Prevent separated islands of zones employing a new network from developing by not letting zone j go to the new network unless zone j - 1 has already done so.
291	charsizg	float	1.	-	Graphics character size for grid labels (relative to MONGO’s default character size).
292	charsizc	float	0.85	-	Graphics character size for curve labels (relative to MONGO’s default character size).
293	charsizz	float	0.85	-	Graphics character size for zone sentinels (relative to MONGO’s default character size).
294	charsizh	float	1.	-	Graphics character size for header info (relative to MONGO’s default character size).
295	widthtd	float	60.	-	Obsolete.
296	ibackgnd	integer	-1	-	Graphics window background color is white if ibackgnd (p 296) 0 and black if ibackgnd (p 296) 0.
297	dscalem	float	1.9892e+33	g	Mass units used for the internal mass coordinate, ymass, used for dumping purpos~s.
298	ngridmin	integer	10	-	Minimum number of points allowed in a dump grid.
299	ncycqq	integer	5	-	Number of KEPLER cycles between post-processor dump cycles, i.e., calls to dumpq.
300	lentrack	integer	16384	bytes	Length of the track(s) assigned to each dump variable in the post-processor dump(s).
301	idtmaxl	integer	40	-	Maximum number of post-processor dump cycles beyond the currently-specified dump cycle (LTIME) that readq will search to try to get values for each point in the ‘advanced’ grid formed from dump variable values at or beyond ltime.
302	idtlook	integer	10	-	Default number of post-processor dump cycles between LOOK plots or prints or reconstructed TIMEMAP grids.
303	backfacq	float	0.5	-	If the fractional change in a dump variable since the last dump cycle exceeds backfacq (p 303) ratzdump(IDAT), then also dump the old value of that variable at the previous (dump cycle) time point, where backfacq (p 303) ratzdump(idat) is the (previously specified) maximum allowed fractional change between dumps of this zonal dump variable (indexed by idat).
304	tempstop	float	1.e+99	K	Terminate the problem when the central temperature reaches TEMPSTOP.
305	denstop	float	1.e+99	g/cc	Terminate the problem when the central density reaches denstop (p 305).
306	vinstop	float	1.e+99	cm/sec	Terminate the problem when the infall velocity below vinstopm (p 462) exceeds vinstop (p 306).
307	o16stop	float	-1.	-	Terminate the problem when the o16 mass-fraction drops below o16stop (p 307) provided the central temperature exceeds tqselim (p 184).
308	timezms	float	1.e+99	sec	Time at which to make the zero-age-main-sequence (ZAMS) parameter changes and restart dump (typically 1.e+12).
309	izonezms	integer	1	-	Reset the value of izonef (p 86) to izonezms (p 309) at the time specified by timezms (p 308).
310	q1faczms	float	0.1	-	Reset the value of q1fac (p 13) to q1faczms (p 310) at the time specified by timezms (p 308).
311	tempcig	float	1.e+99	K	Central temperature at which to make the pre-carbon-ignition parameter changes and restart dump (typically 5.e+8).
312	yflrxcig	float	0.003	-	Reset the value of yfloorx (p 47) to yflrxcig (p 312) when the central temperature specified by tempcig (p 311) is reached.
313	fmaxmcig	float	1.	-	Reset the value of fmaxm (p 195) to fmaxmcig (p 313) when the central temperature specified by tempcig (p 311) is reached if fmaxmcig (p 313) 0.
314	fmax0cig	float	1.	-	Reset the value of fmax0 (p 150) to fmax0cig (p 314) when the central temperature specified by tempcig (p 311) is reached if fmax0cig (p 314) 0.
315	toffset	float	0.	sec	Cumulative amount of time by which the problem time has been offset by zerotime commands.
316	abunminb	float	0.0001	-	Lower mass-fraction limit of the isotopic abundance plot.
317	abunmaxb	float	1.	-	Upper mass-fraction limit of the isotopic abundance plot.
318	numiso	integer	0	-	Number of BURN isotopes to be plotted, starting from the first one listed by the most recent setiso command.
319	timeref	float	-1.e+99	sec	Reference time used in calculating the time coordinate in timeplots and timemaps.
320	tosetref	float	0.	sec	Reference offset time used in calculating the time-coordinate for timeplots and timemaps.
321	timecmin	float	0.	-	Minimum value of the time-coordinate to be plotted in timeplots and timemaps (time coordinate units, see maptime (p 327)).
322	timecmax	float	0.	-	Maximum value of the time-coordinate to be plotted in timeplots and timemaps (time coordinate units, see maptime (p 327)).
323	yemax	float	0.5	mol/g	Maximum value of allowed when initializing a new ISE zone.
324	abarsemi	float	4.	g/mol	Value of the zonal mean atomic weight, , used to divide the star into two regions with separately specifiable values of the semiconvective mixing rate and the overshoot mixing coefficient.
325	drmultlo	float	-1.	-	Semiconvective mixing will be slower than thermal transport by at least drmultlo (p 325) (about 0.1) in zones where the mean atomic weight, , is below abarsemi (p 324).
326	woverslo	float	-1.	-	The semiconvective test parameter, , is taken to be = woverslo (p 326) * abs ( log(T1/T0)) for the special overshoot semiconvective zones where would otherwise be less than 0 and when < abarsemi (p 324).
327	maptime	integer	2	-	Flag indicating the desired time coordinate in timeplots and timemaps.
328	vminmap	float	1.e+99	-	Minimum value of the time map variable to be mapped.
329	vmaxmap	float	-1.e+99	-	Maximum value of the timemap variable to be mapped.
330	vratmap	float	1.e-99	-	Minimum ratio of the minimum timemap variable limit to the maximum timemap variable limit in the case when the actual minimum value of the current timemap variable would otherwise be used as the minimum timemap limit.
331	tempcdep	float	1.e+99	K	If the central temperature is tempcdep (p 331), then make a restart dump labeled #cdep, execute the alias-defined cdep command, and reset tempcdep (p 331) to 1.e+99.
332	o16odep	float	-1.e+99	-	If the central oxygen abundance is o16odep (p 332) and the central temperature is tqselim (p 184), then make a restart dump labelled #odep, execute the alias-defined odep command, and reset o16odep (p 332) to -1.e+99.
333	tempchar	float	1.e+99	K	If the central temperature is tempchar (p 333), then make a restart dump labeled with #tn and the current cycle number, reset tempchar (p 333) to 1.e+99, and finally execute the alias-defined tnchar command.
334	denchar	float	1.e+99	g/cc	If the central density is denchar (p 334), then make a restart dump labeled with #dn and the current cycle number, reset denchar (p 334) to 1.e+99, and finally execute the alias-defined dnchar command.
335	abarchar	float	1.e+99	-	If the central mean atomic weight (abar) is abarchar (p 335), then make a restart dump labelled with #ab and the current cycle number, reset abarchar (p 335) to 1.e+99, and finally execute the alias-defined abchar command.
336	zonemmin	float	1.e-99	g	Minimum mass that a pair of zone may have and still be allowed to be adzoned.
337	icalcne	integer	0	-	Calculate more accurate electron densities in partially ionized regions if icalcne (p 337) 0 using Lisa Ensman’s multiple-ion, Saba equilibrium subroutine, calcne.
338	xneconv	float	1.e-5	-	Maximum allowable fractional convergence error in the electron density calculated by subroutine calcne.
339	ionflag	integer	0	-	Determine calculation of ionization.
340	xnemin	float	1.e-5	-	Minimum mass fraction for which an element is included in the calculation of Saha ionization equilibrium done in subroutine calcne.
341	xkapgam	float	0.054	cm**2/g	Assumed effective opacity for the deposition of gamma ray energy from the radioactive decay of and , as controlled by timex0 (p 38).
342	egamp	float	0.	-	Dimensionless correction factor used in calculating the escape of gamma rays from the radioactive decay of Ni56 and , as controlled by timex0 (p 38).
343	tshock	float	1.e+99	sec	If the problem time is tshock (p 343), then make a restart dump labelled #shock, reset tshock (p 343) to 1.e+99, and finally execute the alias-defined tshock command.
344	tnucleo	float	1.e+99	-	If the problem time is tnucleo (p 344), then make a restart dump labelled #nucleo, reset tnucleo (p 344) to 1.e+99, and finally execute the alias-defined tnucleo command.
345	tenvel	float	1.e+99	-	If the problem time is tenvel (p 345), then make a restart dump labeled #envel , reset tenvel (p 345) to 1.e+99, and finally execute the alias-defined tenvel command.
346	nfirstq	integer	0	-	Default value of the first cycle to be read or plotted in making post-processor edits, time plots, or timemaps.
347	xlmxmult	float	1.	-	Obsolete.
348	binm10	float	:p:`totm0`	Msun	The initial mass of the star being evolved (primary).
349	binm20	float	0.	Msun	The initial mass of the binary companion star.
350	binalp	float	1.	-	[PJH92]‘s alpha parameter, related to the angular momentum of the mass lost in binary transfer.
351	binbet	float	1.	-	[PJH92]‘s beta parameter.
352	bina0	float	2.	AU	Initial binary separation.
353	binmdt	float	0.001	Msun/yr	Mass loss rate due to binary mass transfer assumed when the primary stars exceeds its Roche radius.
354	rocher	float	1.e+99	cm	Current Roche radius.
355	xlanger1	float	0.	?	Langer mass loss rate parameter 1.
356	xlanger2	float	0.	?	Langer mass loss rate parameter 2.
357	iburnye	integer	0	-	Initialization of ISE zones with taken from the BURN network.
358	relmult	float	0.	-	Multiplier on GR corrections.
359	geemult	float	1.	-	Multiplier on gravitational constant.
360	grbparm	float	0.	?	Energy deposition for GRB modelling.
361	swmult	float	0.	-	Obsolete.
362	tsharp	float	0.	?	A parameter for Type Ia SNe simulations.
363	xmlossn	float	0.	-	Multiplier on Niewenhuijzen & de Jager mass loss rate.
364	nangmix	integer	0	-	Enable rotationally-induced mixing.
365	angfmu	float	1.	-	Molecular weight sensitivity of rotational mixing processes.
366	angfc	float	1.	-	Efficiency of chemical mixing by rotational instabilities.
367	angfjc	float	1.	-	Efficiency of angular momentum transport by (semi)convection.
368	angrcrit	float	2500.	-	Critical Reynolds number (affects secular shear instability).
369	angric	float	0.25	-	Critical Richardson number (do not change).
370	angfjdsi	float	1.	-	General efficiency multiplier for dynamical shaer instability.
371	angfjshi	float	1.	-	General efficiency multiplier for Solberg-Hoiland instability.
372	angfjssi	float	1.	-	General efficiency multiplier for secular shear instability.
373	angfjez	float	1.	-	General efficiency multiplier for Eddington-Sweet circulation.
374	angfjgsf	float	1.	-	General efficiency multiplier for Goldreich-Schubert-Fricke instability.
375	fmin	float	1.	-	Under-relaxation factor on the corrections taken each iteration in the Henyey-solver.
376	ncnvout	integer	0	-	Write out convection plot file data (*.cnv) every ncnvout (p 376) cycles. Off when 0.
377	kaptab	integer	1	-	Select opacity table.
378	fkapz	float	1.	-	Multiplier on metallicity used in OPAL opacities.
379	zfakexp	float	0.5	-	Metallicity-dependence of the mass loss.
380	angsmt	float	1.	-	Secular rotational mixing processes may not change by more than that per time-step.
381	nangsmg	integer	2	-	Smooth some gradients used for computation of the rotational instabilities over that much grid points on either side. A Gaussian smoothing profile is used.
382	angsml	float	0.001	-	Secular rotational mixing processes may change by at least that fraction for the local zonal diffusion time-scale.
383	angsmm	float	0.001	-	Secular rotational mixing processes may change by at least that fraction of the total diffusion coefficient.
384	ipapsize	integer	6120792	-	Paper dimensions in points (inch/72).
385	hstatxm	float	1.e+99	g	Interior mass coordinate above which hydrostatic stratification is assumed.
386	hstatym	float	-1.e+99	g	Surface mass coordinate above which hydrostatic stratification is assumed.
387	xmlossw	float	0.	-	Multiplier on WR mass loss rates.
388	iold	integer	0	-	Set to value other than 0 to use old physics - mostly fix that energy generation in APPROX did not include neutrino losses and mass excess but only considered differences in binding energy. These are used as flags.
389	rhotrans	float	1.e+7	-	Some SNIa stuff.
390	nwndout	integer	0	-	Write out wind data to wind file nameprob.wnd every nwndout (p 390) cycles. Off when 0.
391	kapverb	integer	0	-	Verbostiy of opacity subroutine. Zero gives no messages.
392	xl0limf	float	4.	-	Multiplier on limiting flux in radiation flux limiter.
393	xl0limk	float	0.	-	Limiting flux multiplied by.
394	llimout	integer	1	-	Limit to radiative flux of outer zone only if set to 1 (old behaviour), otherwise limit to flux of the hotter zone (upper or lower).
395	nenout	integer	0	-	Write out total energies in file nameprob.ent every nenout cycles. Set to zero to turn off.
396	ipromin	integer	1	-	Minimum zone for production factor/yield plot.
397	ipromax	integer	99999	-	Maximum zone for production factor/yield plot.
398	iprownd	integer	1	-	Take into account wind when computing production factors/yields.
399	iproyld	integer	0	-	Determines the BURN abundance plot type (plot 9).
400	minapro	integer	-1000	-	Minimum mass number for production factor/yield plot; automatic determination if set .
401	maxapro	integer	-1000	-	Maximum mass number for production factor/yield plot; automatic determination if set LE -99.
402	proymin	float	1.e-10	-	Minimum value for yield in BURN plot types 1, 2, 6, 8, 11, and 13 (iproyld (p 399)).
403	proymax	float	1000.	-	Maximum value for yield in BURN plot types 1, 2, 6, 8, 11, and 13 (iproyld (p 399)).
404	xkdmin	float	1.e-20	-	Minimum value for density used in subroutine kappa.
405	h1hdep	float	0.01	-	Central hydrogen abundance at which the #hdep dump is made.
406	he4hedep	float	0.01	-	Central helium abundance at which the #hedep dump is made.
407	nenuout	integer	0	-	Cycle frequency used to write *.enu data file.
408	optconv	float	2./3.	-	Minimum optical depth below (i.e., outside of) which zone are not allowed to become convective.
409	rloss	float	1.e+99	-	Maximum radius beyond which zones are removed from star, similar to vloss (p 271).
410	tloss	float	-1.e+99	-	Minimum temperature below which zones are removed from star, similar to vloss (p 271) and rloss (p 409).
411	tapprox	float	-1.e+99	-	Change zones back to APPROX network if temperature drops below tapprox (p 411).
412	semilan	float	0.	-	efficiency parameter for semiconvection according to [Lan83].
413	profmin	float	0.0001	-	Minimum value for production factor in BURN plot types 0, 5, and 10 (iproyld (p 399)).
414	profmax	float	1000.	-	Maximum value for production factor in BURN plot types:n:0, 5, and 10 (iproyld (p 399)).
415	proamin	float	1.e-12	-	Minimum value for the mass fraction BURN plot types 3, 4, 7, 9, 12, and 14 (iproyld (p 399)).
416	proamax	float	1000.	-	Maximum value for the mass fraction BURN plot types 3, 4 7, 9, 12, and 14 (iproyld (p 399)).
417	yelimb	float	0.497	-	Minimum value of yeburn below which no BURN co-processing is considered.
418	irecb	integer	0	-	Record maximum abundances in BURN network (and save in the restart dumps if set to 1.
419	bmasslow	float	-1.e+99	-	Mass coordinate below which no BURN co-processing is considered.
420	si28dep	float	0.0001	-	Central si28 mass fraction at which the #sidep dump is made, if the mass fractions of o16 and he4 are below 0.01.
421	ivspecl	integer	0	-	Version of the special rate set to use.
422	ivrate	integer	1	-	Version of the APPROX network rate subroutine to use.
423	magnet	integer	0	-	Switch to enable magnetic fields according to [Spr02].
424	nosht	integer	2	-	Convective zones are bound by overshooting layer only of they are at least nosht (p 424) zones thick.
425	alpth	float	1.	-	Efficiency factor for thermohaline convection.
426	abarstep	float	1.e+99	-	No rezoning is performed when the region under consideration contains a step in of more than abarstep (p 426).
427	zbarstep	float	1.e+99	-	No rezoning is performed when the region under consideration contains a step in of more than zbarstep (p 427).
428	xmustep	float	1.e+99	-	No rezoning is performed when the region under consideration contains a step in the mean molecular weight, of more than xmustep (p 428).
429	netmin	integer	1	-	Minimum APPROX network number until which BURN coprocessing is followed.
430	awcorotz	float	0.	-	-component of co-rotation angular velocity.
431	ymcorot	float	-1.	-	Make the outermost ymcorot (p 431) rotate with angular velocity awcorotx (p 581), awcoroty (p 582), awcorotz (p 430). Off if 0..
432	nstrout	integer	0	-	Write out structure data file nameprob.str every nstrout (p 432) time steps.
433	mixcycl	integer	0	-	Determine when during cycle “”mixing” is being done.
434	lburn	integer	0	-	Substitute BURN network for APPROX network (including energy generation, , , ) when set to 1.
435	lbbkup	integer	2	-	This parameter regulates the behaviour of lburn (p 434).
436	rlossmin	float	1.	-	Minimum radius for which to apply vloss (p 271).
437	lcout	integer	10	-	Number of outer layers to be written in light curve output file, .lc.
438	xmagfmu	float	1.	-	Multiply -gradient considered in magnetic field model by this factor.
439	xmagft	float	1.	-	Multiply -gradient considered in magnetic field model by this factor.
440	xmagfnu	float	1.	-	Multiply eddy viscosity by considered in magnetic field model this factor.
441	xmagfdif	float	1.	-	Multiply eddy diffusivity in magnetic field model by this factor.
442	dxncbkup	float	1.e-7	-	Backup if abundance change vector dy returned from burn implies a mass non-conservation larger than dxncbkup.
443	iplotb	integer	0	-	Control isotope/network used in abundance plots, Plot 3.
444	minzone	integer	1	-	Do not rezone the innermost minzone (p 444) zones.
445	zonemmax	float	1.e+99	-	Do not dezone zones bigger than zonemmax (p 445).
446	tenubar	float	-1.	MeV	Electron anti-neutrino temperature for core collapse neutrino flux.
447	levcnv	integer	1	-	Number of levels per dex for .cnv output file.
448	mingain	integer	-1	-	log of minimum energy generation (nuclear + neutrinos) for .cnv output file.
449	minloss	integer	-1	-	log of minimum energy loss (nuclear + neutrinos) for .cnv output file.
450	minnucl	integer	-1	-	log of minimum nuclear energy loss for .cnv output file.
451	ddmin	float	-1.	-	Minimum density for dezoning.
452	iazonef	integer	1	-	Flag to switch off adzoning.
453	dynfac	float	1.	-	Multiplier on dynamic time-scale used to determine whether to treat problem (zones) as “dynamic” or “static” in subroutine update.
454	h1hburn	float	0.4	-	Hydrogen mass fraction at which to make the #hburn dump.
455	c12heign	float	0.01	-	Carbon mass fraction at which to make the #heign dump.
456	he4hebrn	float	0.5	-	Helium mass fraction at which to make the #heburn dump.
457	zonermin	float	-1.	-	Do not adzone a pair of zones with thickness below zonermin (p 457).
458	zonermax	float	1.e+99	-	Do not dezone zones with thickness greater than zonermax (p 458).
459	xmixnova	float	0.	-	Extent (in mass) of a linear composition gradient between substrate and newly accreted material.
460	accmass	float	0.	-	Determine accretion mass.
461	accmassf	float	1.4142	-	The mass of newly accreted zone is may not differ by more than a factor (or fraction) accmassf (p 461) from the current outermost zone.
462	vinstopm	float	1.e+99	-
463	lowamul	integer	0	-	Minimum mass number for which weak rates are modified.
464	pdmult	float	1.	-	Multiplier on positron decay/electron capture rate for lowamul (p 463) ihwamul (p 466).
465	edmult	float	1.	-	Multiplier on electron decay rate for lowamul (p 463) ihwamul (p 466).
466	ihwamul	integer	0	-	Maximum mass number for which weak rates are modified.
467	kapburn	integer	0	-	Use of BURN abundances for opacity.
468	fackap	float	1.	-	Multiply opacity and its derivatives by fackap (p 468).
469	awwkloss	float	1.e+99	-	Remove all outer shells that have an angular velocity angwwk -times bigger than Keplerian velocity.
470	lossrot	integer	0	-	Use Langer’s (1998) formula for rotationally enhanced mass loss.
471	ymjkep	float	-1.e+99	-	Reduce angular momentum to Keplerian angular momentum in the surface layes down to an exterior mass of ymjkep (p 471) if it exceeds Keplerian rotation.
472	maxzone	integer	0	-	Do not rezone outermost maxzone (p 472) zones.
473	cfakexp	float	0.	-	Metallicity dependent mass loss scaling with C abundance for cool Pop III stars.
474	minnucg	integer	-1	-	log of minimum nuclear energy generation for .cnv output file.
475	mingaind	integer	21	-	log of minimum energy generation (nuclear + neutrinos) per cm for .cnv output file.
476	minlossd	integer	21	-
477	minnucgd	integer	21	-	log of minimum nuclear energy generation per cm for .cnv output file.
478	minnucld	integer	21	-	log of minimum nuclear energy loss per cm for .cnv output file. output file.
479	tweakmin	float	0.5	1.e+9 K	Minimum temperature for weak rates.
480	centmult	float	0.	-	Multiplier on simplistic centrifugal force.
481	mixout	integer	0	-	Write out mixing file nameprob.mix every mixout (p 481) cycles.
482	irprox	integer	1	-	Determine use of RPROX network.
483	n14pg	integer	3	-	Version of N14(p,g) rate to use.
484	r3amult	float	1.	-	Multiplier on triple-alpha reaction rate.
485	ibwarn	integer	1	-	Set BURN date warning messages.
486	ifallbk	integer	0	-	Activate fallback treatment.
487	xnumu12	float	0.	-	Neutrino magnetic moment in units of .
488	nzsave	integer	0	-	Number of older z files to save.
489	axion	float	0.	-	Axion mass eV.
490	zmhiconv	float	1.e+99	g	Maximum mass for convection.
491	rnhiconv	float	1.e+99	cm	Maximum radius for convection.
492	icutbin	integer	1	-	Truncate binary output files on restart/generation.
493	nconvers	integer	<most	-	Recent version> Version of the current convection output file.
494	nwndvers	integer	<most	-	Recent version> Version of the current wind output file.
495	h1hign	float	0.01	-	Mass fraction of h1 burnt to write out the #hign dump.
496	wimp	float	0.	GeV	WIMP mass.
497	ipdtmin	integer	0	-	Minimum time between plot outputs in seconds.
498	minneug	integer	-1	-	log of minimum neutrino energy generation (deposition?) for .cnv output file.
499	minneul	integer	-1	-	log of minimum neutrino energy loss (deposition?) for .cnv output file.
500	minneugd	integer	21	-	log of minimum neutrino energy generation (deposition?) per cm for .cnv output file.
501	minneuld	integer	21	-	log of minimum neutrino energy loss per cm for .cnv output file.
502	wimpsip	float	1.e-43	cm**2	Spin-independent cross section of WIMPs on protons.
503	wimpsin	float	1.e-43	cm**2	Spin-independent cross section of WIMPs on protons.
504	wimpsdp	float	1.e-38	cm**2	Spin-dependent cross section of WIMPs on protons.
505	wimpsdn	float	1.e-38	cm**2	Spin-dependent cross section of WIMPs on protons.
506	wimprho0	float	1.e+13	GeV/cc	WIMP density.
507	wimpv0	float	1.e+6	cm/sec	WIMP velocity dispersion.
508	wimpvelo	float	0.	cm/sec	Velocity of star relative to WIMP dark matter halo.
509	iwimpb	integer	1	-	Network to use for WIMP cross section calculations.
510	angw0z	float	0.	rad/sec	-component of angular velocity of the inner angw0m (p 511) mass.
511	angw0m	float	-1.	g	Mass coordinate (not including summ0 (p 61)) for which inner angular velocities angw0x (p 579), angw0y (p 580), and angw0z (p 510) are set.
512	angjaccz	float	0.	cm/sec**2	-component of specific angular momentum of newly accreted material.
513	zoneymax	float	1.e+99	-	Do not dezone zones with xm (i)/ym (i) zoneymax (p 513).
514	accdepth	float	0.	-	Depth where mass is to be accreted.
515	pulse051	float	0.	B	Initial pulsar rotational energy.
516	pulsb15	float	0.	1.e+15 G	Asymptotic pulsar magnetic field at late times.
517	nzpuls	integer	10	-	Number of zones over which to distribute pulsar energy.
518	fracadz	float	2.	-	Maximum ratio of mass accreting zone before it is forced to adzone.
519	losseadv	integer	1	-	Include energy term from advection in mass loss.
520	iacceadv	integer	1	-	Include energy term from advection in accretion.
521	iaccadv	integer	1	-	Advection of composition during accretion.
522	tnumin	float	1.e+7	-	Minimum temperature for neutrino losses if APPROX and BURN are not active.
523	isurf	integer	0	-	Choose atmosphere model for added boundary pressure 0.
524	nlogout	integer	1	-	Output log file if set to 1.
525	ipnuc	integer	1	-	Do nuclear burning/energy generation.
526	ipnu	integer	1	-	Do neutrino losses.
527	amasslow	float	-1.e+99	-	Minimum mass coordinate for APPROX network.
528	umasslow	float	-1.e+99	-	Minimum mass for neutrino losses.
529	idecmode	integer	1	-	Mode for “decretion” model.
530	decrate	float	0.	-	Rate of mass decretion from inner zone.
531	fracdec	float	1.	-	Fraction of mass of inner zone for dezoning if in decretion mode.
532	jmdec	integer	1	-	Zone from which to remove mass in decretion.
533	lumfile	integer	0	-	Use nameprob.lum file for time-dependent base luminosity.
534	acctimef	float	1.	-	Multiply accretion time and time scale by this factor for accretion data from a file.
535	xl0timef	float	1.	-	Multiply base luminosity time and time scale by this factor for accretion data from a file.
536	nsekout	integer	0	-	Write out data file for NuGrid every nsekout (p 536) cycles.
537	iadapv	integer	1	-	Verbosity of adapent output.
538	ittyv	integer	-1	-	Verbosity of tty output.
539	ihe4cc	integer	1	-	Switch on charged current neutrino reactions on he4 if set to 1.
540	inuenc	integer	1	-	Switch on neutral current due to electron neutrinos if set to 1.
541	inuebnc	integer	1	-	Switch on neutral current due to electron anti-neutrinos if set to 1.
542	ibdatov	integer	0	-	Overwriting of hard-coded and bdat rates.
543	h1hm2	float	0.02	-	Mass fraction of h1 burnt to write out the #hm2 dump.
544	h1hm5	float	0.05	-	Mass fraction of h1 burnt to write out the #hm5 dump.
545	h1hm10	float	0.1	-	Mass fraction of h1 burnt to write out the #hm10 dump.
546	nentout	integer	0	-	Frequency to write out *.ent file.
547	nentlev	integer	0	-	Number of levels to write out in *.ent file.
548	entdm	float	0.5	Msun	Mass of level interval to write out into *.ent file.
549	xl0ratef	float	1.	-	Multiplier on luminosity from base from time-dependent input file.
550	accratef	float	1.	-	Multiplier on accretion rate from time-dependent input file.
551	nangdis	integer	0	-	Rotational energy source from dissipation.
552	xmagfbr	float	0.0001	-	Strength of equipartition radial magnetic field used for simple dynamo with magnet (p 423) 8.
553	xmagfbt	float	0.01	-	Strength of equipartition toroidal magnetic field used for simple dynamo with magnet (p 423) 8.
554	accdnfac	float	0.5	-	Density of newly accreted zone relative to outermost zone.
555	taucorot	float	0.	sec	Time scale on which co-rotation is established.
556	iccrate	integer	0	-	Select BR for rate.
557	nidecay	integer	1	-	Switch for decay energy depositon in supernova after explosion.
558	nnuout	integer	0	-	Write out data file for neutrino information every nnuout (p 558) cycles.
559	pulsb150	float	1.	1.e+15 G	Initial magnetic field strength.
560	pulstau	float	100.	1.e+15 G	Decay time of initial magnetic field toward.
561	xheatl	float	0.	erg/sec	Rate of energy deposition from extra heating.
562	xheatym	float	-1.e+99	g	Center of mass where energy is deposited; mass is measured relative to surface.
563	xheatdm	float	0.	g	Width () of the Gaussian in mass where energy is deposited.
564	rmaxeadv	float	1.e+99	-	Maximum radius for which energy deposition for advection from accretion is considered.
565	xmlosse	float	0.	-	Mass loss for super-eddignton luminosity.
566	jslosse	integer	1	-	Determine zone where values for Eddington luminosity mass loss are taken.
567	loczone	integer	0	-	Flag to switch on local criterion for new mass grid on adzoning and dezoning.
568	iaccunit	integer	0	-	Unit in which to measure accretion depth.
569	convmass	float	1.e+99	-	Don’t do convcetion outside convmass (p 569).
570	taumin	float	0.	-	Use optical depth-based zoning where taumin (p 570).
571	taurat	float	0.2	-	Optical depth ratio used in rezoning.
572	taulim	float	0.1	-	Optical depth limit used in rezoning.
573	ir1212	integer	0	-	rate to use.
574	frsurf	float	0.	-	Multiplier on centrifugal force used for atmosphere boundary pressure.
575	psurfmlt	float	1.	-	Multiplier atmosphere boundary pressure.
576	nang3d	integer	0	-	Switch on 3D angular momentum.
577	angjaccx	float	0.	cm/sec**2	-component of specific angular momentum of newly accreted material.
578	angjaccy	float	0.	cm/sec**2	-component of specific angular momentum of newly accreted material.
579	angw0x	float	0.	rad/sec	-component of angular velocity of the inner angw0m (p 511) mass.
580	angw0y	float	0.	rad/sec	-component of angular velocity of the inner angw0m (p 511) mass.
581	awcorotx	float	0.	-	-component of co-rotation angular velocity.
582	awcoroty	float	0.	-	-component of co-rotation angular velocity.
583	nrotout	integer	0	-	Write out rotation data file.
584	dweakmin	float	100000.	g/cc	Minimum density for weak rates.
585	angfjdyn	float	0.	-	General efficiency multiplier for generalised dynamical instability.
586	angfsh	float	1.	-	Relative weighing of Solberg-Hoiland instability contribution for generalised dynamical instability (angfjdyn (p 585)).
587	impoparm	integer	1	-	When set to 1, optional parameters can be defined implicitly.
588	rossbyf	float	0.	-	Multiplier on Rossby reference rotation profile.
589	rossbym	float	0.333333	-	Maximum value of gradient in Rossby reference rotation profile.
590	rossbyp	float	0.6	-	Power of Rossby number (Ro) in gradient of Rossby reference rotation profile.
591	ibintype	integer	4	-	Binary interaction type.
592	xmlossrf	float	-0.2	-	Maximum allowed relative change in mass loss rate.
593	cbfmult	float	0.	-	Multiplier for convectively bound flame for O burning flame.
594	cbfkapm	float	0.	-
595	ncbfmult	integer	15	-
596	lburndtc	integer	1	-	If set to 1 track BURN species for time step rather than APPROX species.
597	fnue	float	0.167	-	Fraction of neutrino energy in electron neutrinos.
598	fnueb	float	0.167	-	Fraction of neutrino energy in electron anti-neutrinos.
599	nserec	integer	1	-	Record and plot NSE/QSE isotope abundances.
600	icnet	integer	0	-	Switch on CNET.
601	cnet1	string	Placeholder	-	For future CNET use.
602	cnet2	string	Placeholder	-	For future CNET use.
603	cnet3	string	Placeholder	-	For future CNET use.
604	cnet4	string	Placeholder	-	For future CNET use.
605	cnet5	string	Placeholder	-	For future CNET use.
606	cnet6	string	Placeholder	-	For future CNET use.
607	cnet7	string	Placeholder	-	For future CNET use.
608	cnet8	string	Placeholder	-	For future CNET use.
609	cnet9	string	Placeholder	-	For future CNET use.
610	cnet0	string	Placeholder	-	For future CNET use.
611	dtcrc	float	1.	-	First zone multiplier on change in radius per step.
612	dtctc	float	1.	-	First zone multiplier on change in temp per step.
613	dtclc	float	1.	-	First zone multiplier on change in lumin per step.
614	dtcdc	float	1.	-	First zone multiplier on change in density per step.
615	dtcqc	float	1.	-	First zone multiplier on linear contraction per step.
616	mlosswr	integer	3	-	WR mass loss rate formula to use.
617	xmlossr0	float	0.	Msun/yr	Mass loss bass for limiting relative change in mass loss rates.
618	temp0	float	0.	K	Temperature of inner boundary.
619	gridconv	float	1.e-10	-	Grid generation convergence limit.
620	xksmin	float	0.01	-	Minimum opacity used for atmosphere model.
621	xksmax	float	100.	-	Maximum opacity used for atmosphere model.
622	aentrain	float	0.	-	Efficiency factor for convective entrainment.
623	eentrain	float	1.	-	Exponent for convective entrainment.
624	cnebmax	float	0.03	-
625	nifcout	integer	1	-	Write out interface data to file nameprob.ifc every nifcout (p 625) cycles. Off when 0.
626	rifcout	float	0.	-	Write out interface data for radius rifcout (p 626). Off when 0.
627	nifcvers	integer	<most	-	Recent version> Version of the current interface output file.
628	mrun	integer	0	-	Flags for run termination.
629	dtmin	float	0.	-	Minimum time step for run.
630	nlcout	integer	0	-	Write out slight curve data file nameprob.lc every nlcout (p 630) time steps.
631	nlcvers	integer	<most	-	Recent version> Version of the current lightcurve output file.
632	ioldfree	integer	0	-	Free memory of old variable storage after convergence. May make code use slightly less memory at the expense of overhead for re-allocation of memory at the next cycle.