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.

angfjdyn

585

float

0.

-

General efficiency multiplier for generalised dynamical 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.

angfsh

586

float

1.

-

Relative weighing of Solberg-Hoiland instability contribution for generalised dynamical instability (angfjdyn (p 585)).

angjaccx

577

float

0.

cm/sec**2

x-component of specific angular momentum of newly accreted material.

angjaccy

578

float

0.

cm/sec**2

y-component of specific angular momentum of newly accreted material.

angjaccz

512

float

0.

cm/sec**2

z-component of 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.

angw0m

511

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.

angw0x

579

float

0.

rad/sec

x-component of angular velocity of the inner angw0m (p 511) mass.

angw0y

580

float

0.

rad/sec

y-component of angular velocity of the inner angw0m (p 511) mass.

angw0z

510

float

0.

rad/sec

z-component of angular velocity of the inner angw0m (p 511) mass.

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)).

awcorotx

581

float

0.

-

x-component of co-rotation angular velocity.

awcoroty

582

float

0.

-

y-component of co-rotation angular velocity.

awcorotz

430

float

0.

-

z-component of co-rotation angular velocity.

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.

cbfkapm

594

float

0.

-

cbfmult

593

float

0.

-

Multiplier for convectively bound flame for O burning flame.

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.

dnratmax

80

float

0.25

-

Maximum fractional density change allowed between zones before adzoning.

dnratmin

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

-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, 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).

dweakmin

584

float

100000.

g/cc

Minimum density for weak rates.

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.

ibintype

591

integer

4

-

Binary interaction type.

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.

impoparm

587

integer

1

-

When set to 1, optional parameters can be defined implicitly.

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) :math:’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.

lburndtc

596

integer

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.

nang3d

576

integer

0

-

Switch on 3D angular momentum.

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.

ncbfmult

595

integer

15

-

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.

nrotout

583

integer

0

-

Write out rotation data file.

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.

rnratmin

77

float

0.1

-

Minimum fractional radius change allowed between zones before dezoning.

rocher

354

float

1.e+99

cm

Current Roche radius.

rossbyf

588

float

0.

-

Multiplier on Rossby reference rotation profile.

rossbym

589

float

0.333333

-

Maximum value of gradient in Rossby reference rotation profile.

rossbyp

590

float

0.6

-

Power of Rossby number (Ro) in gradient of Rossby reference rotation profile.

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.

tnratmax

78

float

0.25

-

Maximum fractional temperature change allowed between zones before adzoning.

tnratmin

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

-1.

-

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.

xmlossrf

592

float

-0.2

-

Maximum allowed relative change in 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 (awcorotx (p 581), awcoroty (p 582), awcorotz (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).