List of Parameters by Name¶

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