List of Parameters by Number¶
Number 
Name 
Format 
Default 
Unit 
Synopsis 





Initial or current timestep. 





Initial or current time. 





Extrapolation Parameter used in guessing new , , and values. 





Quadratic artificial viscosity () factor. 





Maximum number of times subroutine 





Maximum desired fractional change in radius per step. 





Maximum desired fractional change in temperature per step. 





Maximum desired fractional change in density per step. 





Maximum desired fractional linear contraction per step. 





Maximum fractional change in the timestep per step. 





Maximum allowed relative convergence error in radius. 





Maximum allowed relative convergence error in temperature. 





Linear artificial viscosity factor. 





Maximum number of cycles. 





Stop time. 





Number of cycles between ASCII edits. 





Time between ASCII edits ( 





Number of cycles between restart dumps. 





Ratio of the convective mixing length to the pressure scale height. 





Fudge factor for convection (about 





Multiplier for the rate of convective mixing. 





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. 





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. 





Semiconvective mixing will be slower than thermal transport by
at least drmult (p 24) (about 





Maximum timestep allowed. 





Factor used to reduce noise in radius extrapolation. 





Factor used to reduce noise in temperature extrapolation. 





Abundance update parameter. Update abundances only if ipup (p 28)
:math:’ne’ 





Multiplier on 





Multiplier on 





Multiplier on Christy opacity. 





Multiplier on Compton opacity. 





Multiplier on conductive opacity. 





Ion energy and pressure multiplier. 





Radiation energy and pressure multiplier. 





Electron energy and pressure multiplier. 





Multiplier on diffusive heat transport. 





Time at which any 





If iytsflag (p 67) 





If iytsflag (p 67) 





No further changes in 





Size of graphics window to be created in the form 





Multiplier on the inertial terms in the momentum balance equation. 





Maximum length of a postprocessor dump file. 





Location of the toplefthand corner of the graphics window relative
to topleft corner of screen, given in the form 





Maximum desired fractional change in abundances per step. 





Minimum elemental mass fraction that effects the timestep. 





Velocity centring parameter (Range: 





Opacity will be no larger than xkmin (p 50) + t7peek (p 49) * . 





Least upper opacity bound. 





Determine whether zones are considered in thermal equilibrium depending on time step and optical thickness. 





Maximum number of times a given step is redone before the code quits. 





Fractional timestep reduction when a step is redone. 





If the maximum fractional change in radius during a timestep exceeds tfcrbu (p 54) dtcr (p 6) then redo step. 





If the maximum fractional change in temperature during a timestep exceeds tfctbu (p 55) dtct (p 7) then redo step. 





Ionization potential. 





Characteristic density for pressure ionization. 





Pairs are included in EOS calculation only if npflag (p 58) 0. 





If the semi convective test parameter, , is






Radius of inner boundary. 





Mass inside inner boundary. 





Luminosity emerging from inner surface. 





Network switch parameter. 





Graphics edits to the monitor are made every npixedit (p 64) KEPLER cycles (if itvstart (p 127) 0). 





Don’t calculate nuclear burning in APPROX if the temperature is less than tnucmin (p 65). 





Initialize new parameters in subroutine 





Consider only changes in 





Temperature at outer boundary. 





Pressure at outer boundary. 





Maximum allowed relative convergence error in luminosity. 





Factor used to reduce noise in luminosity extrapolation. 





Maximum desired fractional change in luminosity per step. 





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. 





If the maximum fractional change in luminosity during a timestep exceeds tfclbu (p 74) dtcl (p 72) then redo step. 





The fractional amount of semi convective mixing that can occur in one timestep is limited to approximately dtsmult (p 75). 





Maximum fractional radius change allowed betWeen zones before adzoning. 





Minimum fractional radius change allowed between zones before dezoning. 





Maximum fractional temperature change allowed between zones before adzoning. 





Minimum fractional temperature change allowed between zones before dezoning. 





Maximum fractional density change allowed between zones before adzoning. 





Minimum fractional density change allowed between zones before dezoning. 





Minimum radius for which adzoning is considered. 





Minimum temperature for which adzoning is considered. 





Minimum density for which adzoning is considered. 





Minimum number of zones allowed after dezoning. 





General rezoning flag. 





Dezoning Flags. 





Maximum radius for which rezoning is considered. 





Minimum amount of potential hydrogen burn in the present timestep needed to trigger the use of any reaction network. 





Minimum amount of 





Mass fraction of heavy elements above which 





Relative convergence required in calculating the electron Fermi degeneracy parameter, . 





Innermost zone in which there is neutrino deposition. 





Outermost zone in which there is neutrino deposition. 





Total neutrino deposition energy. 





Time of neutrino deposition. 





Time scale for neutrino deposition. 





Artificial neutrino core pressure parameter. 





Core density cutoff. 





Core pressure density dependence exponent. 





Neutrino energy loss rate multiplier (APPROX only). 





Nonneutrino nuclear energy generation rate multiplier (APPROX only). 





Fractional density and temperature perturbations used to get nuclear energy generation rate derivatives. 





Increase the rezoner’s sensitivity to temperature gradients by a
factor of tcorefac (p 104) in the region where 





Floor on the temperature used in the ISE calculation. 





Change all APPROX network zones with 





Multiply q1fac (p 13) by artv1 (p 107) in zone 





Multiply q1fac (p 13) by artv2 (p 108) in zone 





Multiply q1fac (p 13) by artv3 (p 109) in zone 





See artv1 (p 107)  artv3 (p 109). 





See artv1 (p 107)  artv3 (p 109). 





See artv1 (p 107)  artv3 (p 109). 





Graphics picturetype control parameter:. 





Label each curve with nplotsym (p 114) character symbols. 





Characteristic temperature for choosing the location of the temperature plotting grid in Thermodynamics (TD) graph. 





Grid choice for the radial coordinate (irtype (p 132) 





Obsolete. 





Obsolete. 





Innermost zone to plot. 





Outermost zone to plot is 





Minimum radius plotted when jp0 (p 119) 0. 





Maximum radius plotted when jp1 (p 120) 





Lower bound of thermodynamics graph ordinate shall be no greater than ymintd (p 123). 





Upper bound of thermodynamics graph ordinate shall be no less than ymaxtd (p 124). 





Energy generation rate graph scale multiplier in the TD graph. 





Pressure graph scale multiplier in the TD graph. 





Graphics device control parameter. 





Least elemental mass fraction plotted or listed in a terminal ion edit. 





Obsolete. 





Velocity graph scale multiplier in the TD graph. 





Radius graph scale multiplier in the TD graph. 





axix type for plots. 





Ordinate value at which to plot convection sentinels in the TD graph. 





Minimum fractional mass coordinate () plotted (fraction of total mass). 





Maximum fractional mass coordinate () plotted (fraction of total mass). 





Electron degeneracy parameter, , above which to use
temperature interpolation in subroutine 





Enable Adaptive BURN network adjustment. 





Multiplier in the effective values of the density, temperature, and radius gradients used to determine the necessity for adzoning or dezoning. 





See fracrz1 (p 138) and fmax0 (p 150). 





See fracrz1 (p 138). 





See fracrz1 (p 138). 





See fracrz1 (p 138). 





Limiting ratio for for rezoning. 





See abarrat0 (p 143). 





See abarrat0 (p 143). 





Don’t do convection if the absolute value of the zone velocity exceeds frcsound (p 146) times the local sound speed. 





Limit the convective velocity to a fraction convlim (p 147) of the local sound speed. 





The semiconvective test parameter, , is taken to be
woversht (p 148) * 





The mass used in calculating ion degeneracy is xmimult (p 149) times the mass of a neutron. 





Adzone mass fraction parameter. 





See fmax0 (p 150). 





See fmax0 (p 150). 





Maximum fractional change in density since the last iteration for which changes and energy generation rates can be extrapolated from their previous values (APPROX only). 





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





Extrapolate energy generation rates and changes in
cases of small temperature and density change only if
iudflag (p 155) is also 





Save every nsdump (p 156) restart dumps. 





Total number of cycles used to compute 





Maximum allowed value of 





Autumated ASCII and graphical output. 





The mean atomic weight, , calculated in subroutine






The effective number of zones used in 





Effective number of zones used in subroutine 





Don’t calculate derivatives in subroutine 





Number of dump cycles between forced dumps of all dumpgrid points of all ‘ions’ specified as dump variables. 





Number of dump cycles between forced dumps of all dumpgrid points of all BURN ‘isotopes’ specified as dump variables. 





Number of dump cycles between forced dumps of all dumpgrid points of all arrays specified as dump variables except ‘ions’ and ‘isotopes’. 





The electron EOS is implicitly coupled to
nuclear burninginduced changes in the electron abundance,
, provided iflagye (p 167) flag 





Make a ISE edit every njeditq (p 168) zones. 





Relative temperature change used in calculating numerical
derivatives in 





Relative density change used in calculating numerical
derivatives in 





Include excited states in ISE calculations only if iexciteh (p 171)






Maximum number of iterations allowed in ISE calculation. 





Allowed relative convergence error in the proton abundance in the ISE calculation. 





Allowed relative convergence error in the neutron abundance in the ISE calculation. 





Allowed relative convergence error in the Si28 abundance in the ISE calculation. 





Increment the proton, neutron, and abundances by a fraction cnseh (p 176) of that calculated by the NewtonRaphson method for an ISE iteration. 





Maximum allowed relative change in proton abundance during an ISE iteration. 





Maximum allo~ed relative change in neutron abundance during an ISE iteration. 





Maximum allowed relative change in Si28 abundance during an ISE iteration. 





Make a detailed edit of the ISE zones every neditq (p 180) general numerical edits. 





If the relative change of the proton, neutron, or abundance is opposite in sign and more than xiter1qe (p 181) in magnitude with respect to the corresponding change during the previous ISE iteration cycle, then cut the current step size in half if more than iter1qe (p 182) iterations have been done. 





See xiter1qe (p 181). 





Edit only those ISE isotopes with mass fractions exceeding xthres (p 183). 





A sufficient condition to change a zone from the APPROX to
ISE network is if its temperature exceeds tqselim (p 184),
its 





See tqselim (p 184). 





See tqselim (p 184). 





Multiplier on neutrino energy losses from weak processes on nuclei and nucleons. 





Change all ISE zones with 





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





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





The ordinate bounds of the velocity graph are
vlimset (p 191) if vlimset (p 191) 





Pause the code at the end of the calculation for zone






See fracrz1 (p 138) and fmax0 (p 150). 





See fracrz1 (p 138). 





See fmax0 (p 150). 





Multiplier on the contribution of ionization potential energy
to the equation of state. See subroutine 





Start a new labelled ASCII output file every nnewoutf (p 197) cycles. 





Make an ISE edit for the central zone every neditq1 (p 198) nedit (p 16) cycles. 





Multiplier on the Wilsonbased nuclear EOS (except for thee thermal
ion component) if it is 





Upper bound on the temperature used to calculate the energy in nuclear excited states and nuclear partition functions. 





Initial guess made in 





Neutronrichsiliconburning flows are included in the ISE
network if isi30brn (p 202) 





Nominal minimum density for changing a zone from the APPROX to the ISE network. 





Force a backup after a call to subroutine BURN from subroutine






Maximim allowed fractional change in abundance before an abundance backup is made. (See abunminx (p 204)). 





Elemental mass fraction floor for making abundance backups. (See abunminx (p 204)). 





Fatal errors in 





Multiplier on the overall 





Limit how fast convective velocity can increase. 





Flags regulating transition to QSE for convectively coupled zones. 





Rate at which mass in the form of new zones is added to the surface of the star. 





Mass of phantom outer zone used to mediate mass accretion. See accrate (p 211). 





Allowed fractional convergence error in density when calculating a hydrostatic initial stellar configuration in the generator using the dstat option. 





Characteristic flame radius for carbon deflagration studies. 





Multiplier on the WignerSeitz Coulomb corrections to the ion
energy and pressure (See [Cla], p. 152, and subroutine 





Rezoning Option Flag. 





If nsurfz (p 217) 





Inner surface zoning parameter. 





Outer surface zoning parameter. 





Multiplier on the de Jager massloss rate. 





Nominal mass loss rate from the surface of the star. 





The original total mass of the star. 





Dezone the masslosing zone if its mass drops below fracdez (p 223) times the average mass of the zones on either side of it. 





Redo the timestep (“backup”) if the fractional change of mass in the masslosing zone exceeds xmratbak (p 224). 





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 





Obsolete. 





Multiplier on the E1 part of the 





Multiplier on the E2 part of the 





Do a forced update of BURN isotopic abundances in every zone every nupdate (p 229) cycles. 





When BURN coprocessing 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). 





BURN coprocessing is skipped if a zone’s exterior mass coordinate is less than bmassmin (p 231). 





BURN coprocessing is skipped if a zone’s exterior mass coordinate is greater than bmassmax (p 232). 





BURN coprocessing is skipped if a zone’s temperature is less than btempmin (p 233). 





BURN coprocessing is skipped if a zone’s temperature is greater than btempmax (p 234). 





BURN coprocessing is skipped if the absolute value of a zone’s
normal nuclear energy generation rate ( 





BURN coprocessing is skipped if a zone’s density is less than bdenmin (p 236). 





BURN coprocessing is skipped if a zone’s density is greater than bdenmax (p 237). 





Force a zone to be updated by the BURN coprocessor if its fractional temperature change since its last BURN processing exceeds tchange (p 238). 





Force a zone to be updated by the BURN coprocessor if its fractional density change since its last BURN processing exceeds dchange (p 239). 





BURN coprocessing is skipped if a zone’s network number,






Make a BURN isotopic abundance edit every neditb (p 241) nedit (p 16) KEPLER cycles. 





Obsolete. 





Make an elemental abundance edit every nedita (p 243) nedit (p 16) cycles. 





Edit only zones whose zone number is jmeditb (p 244) during normal BURN isotopic abundance edits. 





Make a BURN isotopic abundance edit for all zones every neditall (p 245) nedit (p 16) cycles regardless of the value of neditb (p 241). 





Minimum isotopic abundance that affects the calculation of the zonal timestep in the BURN coprocessor. 





Maximum desired fractional change of an isotopic abundance used in determining the zonal timestep used in the BURN coprocessor. 





Maximum factor by which the zonal timestep in the BURN coprocessor can be increased in one zonal cycle. 





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. 





Make a complete ASCII edit of the nuclear processes in all zones considered by the BURN coprocessor every ncomp (p 250) cycles. 





Make an ASCII edit of the nuclear processes in the central zone considered by the BURN coprocessor every ncent (p 251) cycles. 





Make an ASCII edit of the nuclear processes involved every nedt (p 252) BURN coprocessor matrix inversions (debugging only). 





Set newly calculated negative BURN isotope abundances to 





Minimum absolute magnitude of a BURN matrix element for which GaussJordan elimination is carried out in solving for changes in isotopic abundances. 





Total number of BURNcoprocessor matrix inversions so far (internally incremented). 





Total number of BURNcoprocessor backups so far (internally incremented). 





FORTRAN I/O unit number for 





Flag to switch on sparse matrix solver. 





Cumulative number of BURN coprocessor backups due to encountering negative isotopic abundances (incremented internally). 





Factor by which to reduce the BURN coprocessor zonal timestep when a negative isotopic abundance backup occurs. 





Minimum absolute magnitude ofa negative isotopic abundance (except of protons) that can cause a zonal backup in the BURN coprocessor. 





Minimum absolute magnitude of a negative proton abundance that can cause a zonal backup in the BURN coprocessor. 





Maximum number of consecutive negative isotopic abundance backups allowed in the BURN coprocessor before quitting. 





Number of coupled BURN coprocessing / BURN isotope convection subcycles per KEPLER cycle. 





Use Fuller et. al.’s weak rates in the BURN coprocessor if
mazful (p 265) 





Decrease the abundance threshold for 





BURN coprocessing and related edits are done only if
inburn (p 267) 





Number of lines printed per “page” of ASCII output. 





Total number of nonBURN zonal arrays to save in restart dumps. 





Total number of BURN zonal arrays to save in restart dumps. 





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





Edit only those BURN isotope massfractions whose values exceed abunlimb (p 272) in making terminal edits (including isotopic massfraction sums over zones). 





Mass unit used for the mass coordinate employed in making ASCII and terminal edits. 





Make an ascii edit of primary zonal quantities every neditz1 (p 274) nedit (p 16) cycles. 





Make an ascii edit of secondary zonal quantities every neditz2 (p 275) nedit (p 16) cycles. 





Flag determining the minimum amount of information printed in an ASCII cycle edit regardless of the settings of other edit parameters:. 





Effective value of medit (p 276) used in determining the scope of the final edit made when the problem is finished. 





Make an ascii edit of all changeable parameters every neditp (p 278) nedit (p 16) cycles. 





Carbon mass fraction ahead of flame. 





Oxygen mass fraction for neox problems. 





When the flame became RT unstable. 





Sharp1 parameter. 





Minimum speed from convection. 





Sharp2 parameter. 





(Edit only.). 





Time scale for the neutrino pulse considered by the BURN coprocessor in calculating neutrinoinduced nucleosynthesis after core collapse. 





Total energy for the neutrino pulse considered by the BURN coprocessor in calculating neutrino induced nucleosynthesis after core collapse. 





Temperature of the and neutrinos in the neutrino pulse considered by the BURN coprocessor in calculating neutrino induced nucleosynthesis after core collapse. 





Temperature of the electron neutrinos in the neutrino pulse considered by the BURN coprocessor in calculating neutrinoinduced nucleosynthesis after core collapse. 





If noiland (p 290) 





Graphics character size for grid labels (relative to MONGO’s default character size). 





Graphics character size for curve labels (relative to MONGO’s default character size). 





Graphics character size for zone sentinels (relative to MONGO’s default character size). 





Graphics character size for header info (relative to MONGO’s default character size). 





Obsolete. 





Graphics window background color is white if ibackgnd (p 296)






Mass units used for the internal mass coordinate, 





Minimum number of points allowed in a dump grid. 





Number of KEPLER cycles between postprocessor dump
cycles, i.e., calls to 





Length of the track(s) assigned to each dump variable in the postprocessor dump(s). 





Maximum number of postprocessor dump cycles beyond the
currentlyspecified dump cycle ( 





Default number of postprocessor dump cycles between LOOK plots or prints or reconstructed TIMEMAP grids. 





If the fractional change in a dump variable since the last dump
cycle exceeds backfacq (p 303) 





Terminate the problem when the central temperature reaches TEMPSTOP. 





Terminate the problem when the central density reaches denstop (p 305). 





Terminate the problem when the infall velocity below vinstopm (p 462) exceeds vinstop (p 306). 





Terminate the problem when the 





Time at which to make the zeroagemainsequence (ZAMS) parameter
changes and restart dump (typically 





Reset the value of izonef (p 86) to izonezms (p 309) at the time specified by timezms (p 308). 





Reset the value of q1fac (p 13) to q1faczms (p 310) at the time specified by timezms (p 308). 





Central temperature at which to make the precarbonignition
parameter changes and restart dump (typically 





Reset the value of yfloorx (p 47) to yflrxcig (p 312) when the central temperature specified by tempcig (p 311) is reached. 





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) 





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) 





Cumulative amount of time by which the problem time has been offset by zerotime commands. 





Lower massfraction limit of the isotopic abundance plot. 





Upper massfraction limit of the isotopic abundance plot. 





Number of BURN isotopes to be plotted, starting from the first one listed by the most recent setiso command. 





Reference time used in calculating the time coordinate in timeplots and timemaps. 





Reference offset time used in calculating the timecoordinate for timeplots and timemaps. 





Minimum value of the timecoordinate to be plotted in timeplots and timemaps (time coordinate units, see maptime (p 327)). 





Maximum value of the timecoordinate to be plotted in timeplots and timemaps (time coordinate units, see maptime (p 327)). 





Maximum value of allowed when initializing a new ISE zone. 





Value of the zonal mean atomic weight, , used to divide the star into two regions with separately specifiable values of the semiconvective mixing rate and the overshoot mixing coefficient. 





Semiconvective mixing will be slower than thermal transport by
at least drmultlo (p 325) (about 





The semiconvective test parameter, , is taken to be
= woverslo (p 326) * 





Flag indicating the desired time coordinate in timeplots and timemaps. 





Minimum value of the time map variable to be mapped. 





Maximum value of the timemap variable to be mapped. 





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. 





If the central temperature is tempcdep (p 331), then make
a restart dump labeled #cdep, execute the aliasdefined
cdep command, and reset tempcdep (p 331) to 





If the central oxygen abundance is o16odep (p 332) and the central
temperature is tqselim (p 184), then make a restart dump
labelled #odep, execute the aliasdefined odep command, and
reset o16odep (p 332) to 





If the central temperature is tempchar (p 333), then make
a restart dump labeled with #tn and the current cycle
number, reset tempchar (p 333) to 





If the central density is denchar (p 334), then make a
restart dump labeled with #dn and the current cycle number,
reset denchar (p 334) to 





If the central mean atomic weight ( 





Minimum mass that a pair of zone may have and still be allowed to be adzoned. 





Calculate more accurate electron densities in partially ionized
regions if icalcne (p 337) 





Maximum allowable fractional convergence error in the electron
density calculated by subroutine 





Determine calculation of ionization. 





Minimum mass fraction for which an element is included in
the calculation of Saha ionization equilibrium done in
subroutine 





Assumed effective opacity for the deposition of gamma ray energy from the radioactive decay of and , as controlled by timex0 (p 38). 





Dimensionless correction factor used in calculating the escape of gamma rays from the radioactive decay of Ni56 and , as controlled by timex0 (p 38). 





If the problem time is tshock (p 343), then make a restart
dump labelled #shock, reset tshock (p 343) to 





If the problem time is tnucleo (p 344), then make a restart dump
labelled #nucleo, reset tnucleo (p 344) to 





If the problem time is tenvel (p 345), then make a restart
dump labeled #envel , reset tenvel (p 345) to 





Default value of the first cycle to be read or plotted in making postprocessor edits, time plots, or timemaps. 





Obsolete. 





The initial mass of the star being evolved (primary). 





The initial mass of the binary companion star. 





[PJH92]‘s alpha parameter, related to the angular momentum of the mass lost in binary transfer. 





[PJH92]‘s beta parameter. 





Initial binary separation. 





Mass loss rate due to binary mass transfer assumed when the primary stars exceeds its Roche radius. 





Current Roche radius. 





Langer mass loss rate parameter 1. 





Langer mass loss rate parameter 2. 





Initialization of ISE zones with taken from the BURN network. 





Multiplier on GR corrections. 





Multiplier on gravitational constant. 





Energy deposition for GRB modelling. 





Obsolete. 





A parameter for Type Ia SNe simulations. 





Multiplier on Niewenhuijzen & de Jager mass loss rate. 





Enable rotationallyinduced mixing. 





Molecular weight sensitivity of rotational mixing processes. 





Efficiency of chemical mixing by rotational instabilities. 





Efficiency of angular momentum transport by (semi)convection. 





Critical Reynolds number (affects secular shear instability). 





Critical Richardson number (do not change). 





General efficiency multiplier for dynamical shaer instability. 





General efficiency multiplier for SolbergHoiland instability. 





General efficiency multiplier for secular shear instability. 





General efficiency multiplier for EddingtonSweet circulation. 





General efficiency multiplier for GoldreichSchubertFricke instability. 





Underrelaxation factor on the corrections taken each iteration in the Henyeysolver. 





Write out convection plot file data ( 





Select opacity table. 





Multiplier on metallicity used in OPAL opacities. 





Metallicitydependence of the mass loss. 





Secular rotational mixing processes may not change by more than that per timestep. 





Smooth some gradients used for computation of the rotational instabilities over that much grid points on either side. A Gaussian smoothing profile is used. 





Secular rotational mixing processes may change by at least that fraction for the local zonal diffusion timescale. 





Secular rotational mixing processes may change by at least that fraction of the total diffusion coefficient. 





Paper dimensions in points (inch/72). 





Interior mass coordinate above which hydrostatic stratification is assumed. 





Surface mass coordinate above which hydrostatic stratification is assumed. 





Multiplier on WR mass loss rate. 





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. 





Some SNIa stuff. 





Write out wind data to wind file 





Verbostiy of opacity subroutine. Zero gives no messages. 





Multiplier on limiting flux in radiation flux limiter. 





Limiting flux multiplied by. 





Limit to radiative flux of outer zone only if set to 





Write out total energies in file 





Minimum zone for production factor/yield plot. 





Maximum zone for production factor/yield plot. 





Take into account wind when computing production factors/yields. 





Determines the BURN abundance plot type (plot 9). 





Minimum mass number for production factor/yield plot; automatic determination if set . 





Maximum mass number for production factor/yield plot; automatic determination if set LE 99. 





Minimum value for yield in BURN plot types 





Maximum value for yield in BURN plot types 





Minimum value for density used in subroutine 





Central hydrogen abundance at which the #hdep dump is made. 





Central helium abundance at which the #hedep dump is made. 





Cycle frequency used to write *.enu data file. 





Minimum optical depth below (i.e., outside of) which zone are not allowed to become convective. 





Maximum radius beyond which zones are removed from star, similar to vloss (p 271). 





Minimum temperature below which zones are removed from star, similar to vloss (p 271) and rloss (p 409). 





Change zones back to APPROX network if temperature drops below tapprox (p 411). 





efficiency parameter for semiconvection according to [Lan83]. 





Minimum value for production factor in BURN plot types 





Maximum value for production factor in BURN plot types:n:0,






Minimum value for the mass fraction BURN plot types 





Maximum value for the mass fraction BURN plot types 





Minimum value of 





Record maximum abundances in BURN network (and save in the restart dumps if set to 1. 





Mass coordinate below which no BURN coprocessing is considered. 





Central 





Version of the special rate set to use. 





Version of the APPROX network rate subroutine to use. 





Switch to enable magnetic fields according to [Spr02]. 





Convective zones are bound by overshooting layer only of they are at least nosht (p 424) zones thick. 





Efficiency factor for thermohaline convection. 





No rezoning is performed when the region under consideration contains a step in of more than abarstep (p 426). 





No rezoning is performed when the region under consideration contains a step in of more than zbarstep (p 427). 





No rezoning is performed when the region under consideration contains a step in the mean molecular weight, of more than xmustep (p 428). 





Minimum APPROX network number until which BURN coprocessing is followed. 





component of corotation angular velocity. 





Make the outermost ymcorot (p 431) rotate with angular velocity
awcorotx (p 581), awcoroty (p 582), awcorotz (p 430). Off
if 





Write out structure data file 





Determine when during cycle “”mixing” is being done. 





Substitute BURN network for APPROX network (including energy
generation, , , ) when
set to 





This parameter regulates the behaviour of lburn (p 434). 





Minimum radius for which to apply vloss (p 271). 





Number of outer layers to be written in light curve output file, .lc. 





Multiply gradient considered in magnetic field model by this factor. 





Multiply gradient considered in magnetic field model by this factor. 





Multiply eddy viscosity by considered in magnetic field model this factor. 





Multiply eddy diffusivity in magnetic field model by this factor. 





Backup if abundance change vector 





Control isotope/network used in abundance plots, Plot 





Do not rezone the innermost minzone (p 444) zones. 





Do not dezone zones bigger than zonemmax (p 445). 





Electron antineutrino temperature for core collapse neutrino flux. 





Number of levels per dex for .cnv output file. 























Minimum density for dezoning. 





Flag to switch off adzoning. 





Multiplier on dynamic timescale used to determine whether to treat
problem (zones) as “dynamic” or “static” in subroutine






Hydrogen mass fraction at which to make the #hburn dump. 





Carbon mass fraction at which to make the #heign dump. 





Helium mass fraction at which to make the #heburn dump. 





Do not adzone a pair of zones with thickness below zonermin (p 457). 





Do not dezone zones with thickness greater than zonermax (p 458). 





Extent (in mass) of a linear composition gradient between substrate and newly accreted material. 





Determine accretion mass. 





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. 










Minimum mass number for which weak rates are modified. 





Multiplier on positron decay/electron capture rate for lowamul (p 463) ihwamul (p 466). 





Multiplier on electron decay rate for lowamul (p 463) ihwamul (p 466). 





Maximum mass number for which weak rates are modified. 





Use of BURN abundances for opacity. 





Multiply opacity and its derivatives by fackap (p 468). 





Remove all outer shells that have an angular velocity






Use Langer’s (1998) formula for rotationally enhanced mass loss. 





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. 





Do not rezone outermost maxzone (p 472) zones. 





Metallicity dependent mass loss scaling with C abundance for cool Pop III stars. 


































Minimum temperature for weak rates. 





Multiplier on simplistic centrifugal force. 





Write out mixing file 





Determine use of RPROX network. 





Version of N14(p,g) rate to use. 





Multiplier on triplealpha reaction rate. 





Set BURN date warning messages. 





Activate fallback treatment. 





Neutrino magnetic moment in units of . 





Number of older z files to save. 





Axion mass






Maximum mass for convection. 





Maximum radius for convection. 





Truncate binary output files on restart/generation. 





Recent version> Version of the current convection output file. 





Recent version> Version of the current wind output file. 





Mass fraction of 





WIMP mass. 





Minimum time between plot outputs in seconds. 





























Spinindependent cross section of WIMPs on protons. 





Spinindependent cross section of WIMPs on protons. 





Spindependent cross section of WIMPs on protons. 





Spindependent cross section of WIMPs on protons. 





WIMP density. 





WIMP velocity dispersion. 





Velocity of star relative to WIMP dark matter halo. 





Network to use for WIMP cross section calculations. 





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





Mass coordinate (not including summ0 (p 61)) for which inner angular velocities angw0x (p 579), angw0y (p 580), and angw0z (p 510) are set. 





component of specific angular momentum of newly accreted material. 





Do not dezone zones with 





Depth where mass is to be accreted. 





Initial pulsar rotational energy. 





Asymptotic pulsar magnetic field at late times. 





Number of zones over which to distribute pulsar energy. 





Maximum ratio of mass accreting zone before it is forced to






Include energy term from advection in mass loss. 





Include energy term from advection in accretion. 





Advection of composition during accretion. 





Minimum temperature for neutrino losses if APPROX and BURN are not active. 





Choose atmosphere model for added boundary pressure






Output log file if set to 





Do nuclear burning/energy generation. 





Do neutrino losses. 





Minimum mass coordinate for APPROX network. 





Minimum mass for neutrino losses. 





Mode for “decretion” model. 





Rate of mass decretion from inner zone. 





Fraction of mass of inner zone for dezoning if in decretion mode. 





Zone from which to remove mass in decretion. 





Use 





Multiply accretion time and time scale by this factor for accretion data from a file. 





Multiply base luminosity time and time scale by this factor for accretion data from a file. 





Write out data file for NuGrid every nsekout (p 536) cycles. 





Verbosity of 





Verbosity of 





Switch on charged current neutrino reactions on 





Switch on neutral current due to electron neutrinos if set to 





Switch on neutral current due to electron antineutrinos if set to






Overwriting of hardcoded and 





Mass fraction of 





Mass fraction of 





Mass fraction of 





Frequency to write out *.ent file. 





Number of levels to write out in *.ent file. 





Mass of level interval to write out into *.ent file. 





Multiplier on luminosity from base from timedependent input file. 





Multiplier on accretion rate from timedependent input file. 





Rotational energy source from dissipation. 





Strength of equipartition radial magnetic field used for simple
dynamo with magnet (p 423) 





Strength of equipartition toroidal magnetic field used for simple
dynamo with magnet (p 423) 





Density of newly accreted zone relative to outermost zone. 





Time scale on which corotation is established. 





Select BR for rate. 





Switch for decay energy depositon in supernova after explosion. 





Write out data file for neutrino information every nnuout (p 558) cycles. 





Initial magnetic field strength. 





Decay time of initial magnetic field toward. 





Rate of energy deposition from extra heating. 





Center of mass where energy is deposited; mass is measured relative to surface. 





Width () of the Gaussian in mass where energy is deposited. 





Maximum radius for which energy deposition for advection from accretion is considered. 





Mass loss for supereddignton luminosity. 





Determine zone where values for Eddington luminosity mass loss are taken. 





Flag to switch on local criterion for new mass grid on adzoning and dezoning. 





Unit in which to measure accretion depth. 





Don’t do convcetion outside convmass (p 569). 





Use optical depthbased zoning where taumin (p 570). 





Optical depth ratio used in rezoning. 





Optical depth limit used in rezoning. 





rate to use. 





Multiplier on centrifugal force use for atmosphere boundary pressure. 





Multiplier atmosphere boundary pressure. 





Switch on 3D angular momentum. 





component of specific angular momentum of newly accreted material. 





component of specific angular momentum of newly accreted material. 





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





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





component of corotation angular velocity. 





component of corotation angular velocity. 





Write out rotation data file. 





Minimum density for weak rates. 





General efficiency multiplier for generalised dynamical instability. 





Relative weighing of SolbergHoiland instability contribution for generalised dynamical instability (angfjdyn (p 585)). 





When set to 





Multiplier on Rossby reference rotation profile. 





Maximum value of gradient in Rossby reference rotation profile. 





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





Binary interaction type. 





Maximum allowed relative change in mass loss rate. 





Multiplier for convectively bound flame for O burning flame. 














