Edit Formats¶

Selected Ascii And Graphics Edit Formats¶

Terminal Edits¶

Below are listed examples of the results of issuing several of the terminal edit commands described on the Zonal Edit Varlables page. To save space and constant repetition, the output fonnat of these commands is fairly terse, but is explained in detail below. The internal KEPLER variables corresponding to the edited quantities are provided where they are not obvious. Each example starts with the $>$ sign that KEPLER uses to prompt the user for input, followed by the command issued. Subsequent line(s) show the results of the command.

Time Edit¶

>t
11981 5.125483408712E+05  dtd(  1) = 9.66E-03  iter=4  jm=505  zb=0

Here

value	result
”`11981`”	the current cycle number ncyc (q 4)
“`5.1254834087l2E+05`”	the current problem time time (p 2) in seconds (relative to the reference time offset toffset (p 315)),
”`dtd`”	the symbol of the ion or other quantity whose rate of change is determining the timestep, `idtcsym` $($ idtcon (q 13) $)$
”`1`”	the index of the zone where the rate of change of this quantity is most time-step-limiting, `jdtc` $($ idtcon (q 13) $)$
”`9.66E-03`”	the value of the current timestep, dt (q 1) (`sec`)
”`iter=4`”	gives the number of iterations that KEPLER took to converge during the latest cycle, iter (q 7).
”`jm=505`”	the number of zones currently in the problem (“`505`”) “`zb=0`” the number of zones processed by the BURN coprocessor during the last KEPLER timestep, nburnz (q 51), if BURN coprocessing is being done.

Time and Sucface Edit¶

>tq
11981  5.12S483408712E+05   dtd(  1) = 9.66E-03 iter=4  jm=505  zb=0
radius = 9.06016E+13  lum = 2.32789E+39  teff = 4.46663E+03

The first row of this edit is the same as for the time edit, as given above. The second row lists respectively,

value	result
”code:radius = 9.06016E+13”	the star’s photospheric radius, radius (q 47), in `cm`
”`lum = 2.32789E+39`”	the star’s surface (non-neutrino) luminosity, xlum (q 46), in `erg/sec`, and
”`teff = 4.46663E+03`”	the star’s effective temperature, teff (q 48) in `K`.

Energy Edit¶

>e
ke = 1.298S6E+49 ie = 9.S7900E+51 pe =-1.12123E+52 te =-1.62027E+51
kd = 5.9669SE+49 id = 2.89362E+51 pd =-3.00413E+51 sd =-5.01087E+49

The first row lists respectively, all in erg

value	result
”code:ke”	the star’s total kinetic energy, enk (q 32),
”code:ie”	total internal energy, eni (q 31),
”code:pe”	total potential energy enp (q 33), and
”code:te”	the sum of these three, ent (q 34).

The second row lists, respectively the current time derivatives of these quantities in erg/sec

value	result
”code:kd”	enkd (q 40)
”code:id”	enid (q 39)
”code:pd”	enpd (q 41)
”code:sd”	entd (q 42)

Ion Mass-Fraction Edit¶

>1 i
nt1    5.688E-03
he4    1.672E-02
ti44   3.982E-03
cr48   2.429E-0l
'fe'   7.303E-01
1-tot  2.883E-07

This edit makes a list of mass fractions of all the ions in the zone specified (here “1”) that exceed abunlim (p 128). The row labeled ‘1-tot gives the difference between unity and the sum of the mass fractions of all ions in the current network. The ion symbols are usually self-explanatory.

BURN Isotope Mass-Fraction Edit¶

This edit makes a list of mass fractions of all the BURN isotopes in the zone specified (here ‘1’) that exceed abunlimb (p 272). The row labeled 1-tot” gives the difference between unity and the sum of the mass fractions of all BURN isotopes in the current network. The isotope symbols are self-explanatory. A sample list of isotopes used in the BURN network can be found here.

Zone Edit¶

>1
   1 m= 0.00000E+00 z= 2.011E-03 t= 2.556E+09 d= 9.897E+07 v=-5.643E-01
conv r= 2.12892E+07 k= 1.020E-02 e= 5.413E+17 p= 2.357E+25 q= 0.000E+00
ise  eta= 6.195E+00 1= 1.282E+42 s= 3.425E+11 n= 3.580E+11  ye=0.49799
      p33(7390,6322) = 7.148E+01             nb= 1.971E+11 yeb=0.49710
     sig= 1.238E+00 x= 2.258E+03 a= 2.751E+01 w= 6.370E-04 j= 4.973E+09
      Bt= 0.000E+00 Br=0.000E+00

Here:

value	result
”`1`”	Index of the zone being edited, $j$
”`m`”	Mass interior to the inner boundary of the zone in solar mass units.
”`z`”	Mass of the zone in solar masses, xm $(j)$ (`Msun`).
”`t`”	Zonal temperature, tn $(j)$ (`K`)
”`d`”	Zonal density, dn $(j)$ (`g/cc`)
”`v`”	Velocity of the outer interface of the zone, un $(j)$ (`cm/sec`)
”`r`”	Radius of the outer interface of the zone, rn $(j)$ (`cm`)
”`k`”	Zonal opacity, xkn $(j)$ (`cm**2/g`)
”`e`”	Specific energy of the zone, en $(j)$ (`erg/g`)
”`p`”	Zonal pressure, pn $(j)$ (`erg/cc`)
”`q`”	Viscous stress (including artificial viscosity), zn $(j)$ (`erg`)
”`eta`”	Electron degeneracy parameter, etan $(j)$
”`l`”	Luminosity across the outer interface of the zone, xln $(j)$ (`erg/sec`)
”`s`”	Net energy generation rate (including neutrino losses), sn $(j)$ (`erg/g/sec`)
”`n`”	Net nuclear energy generation rate (excludes neutrinos losses), snn $(j)$ (`erg/g/sec`)
”`ye`”	Non-BURN electron abundance, zbar $(j)$ $/$ abar $(j)$ (`mol/g`)
”`nb`”	Nuclear energy generation rate last calculated by BURN, sburn $(j)$ (`erg/g/sec`)
”`ye`”	BURN electron abundance yeburn $(j)$ (`mol/g`)
”`sig`”	Total entropy per baryon, stot $(j)$ (`kb/baryon`)
”`x`”	Convective mixing timescale across the outer zonal interface, taumix $(j)$ (`sec`)
”`a`”	Zone’s average atomic weight, abar $(j)$ (`g/mol`)
”`w`”	Zonal angular velocity, angw $(j)$ (`rad/sec`)
”`j`”	Zonal specific angular momentum, angj $(j)$ (`cm**2/sec`)
”`Bt`”	Toroidal component of magnetic field on zone outer boundary, bfbt $(j)$ (`G`), only shown if magnet (p 423) $>$ `0`
”`Br`”	Radial component of magnetic field on zone outer boundary, bfbr $(j)$ (`G`) , only shown if magnet (p 423) $>$ `0`

The symbol “code:conv” exemplifies the convection sentinel, icon $($ j $)$ , that describes the status of convection between zones j and j $+$ 1. Possible values of the convection sentinel are describe with the description of icon.

The symbol “ise” exemplifies the network sentinel. Its value is

value	result
”`ise`”	if the zone is currently using the ISE network, netnum $(j)=$ `2`
”`nse`”	if the NSE network is being used, netnum $(j)=$ `3`
(blank)	if the APPROX network is being used, netnum $(j)=$ `3`

The construct “code:p33(7390,6322) = 7.148E+0l gives the current timestep being used by BURN for this zone, dtimen $(j)$ , and the symbol of the BURN isotope whose rate of change is most limiting this timestep (formed from limnuc $(j)$ ). The two numbers in parentheses are the current KEPLER cycle, ncyc (q 4), and the KEPLER cycle at which this zone was last processed by BURN (also formed from limnuc $(J)$ . Note that BURN processing is usually not done in ISE or NSE zones (see netmax (p 240) and bmasslow (p 419)). When BURN coprocessing has not been requested, the fourth line is omitted from this edit.

Changed In Version 15.2.0: Added angular velocity and specific angular momentum.

Zonal Variable Edit¶

>z dn tn sig sigi sige 1 10
      zone      mass         dn         tn        sig       sigi       sige
 thal    1 0.00000000  2.319E+06  8.374E+07  1.087E+00  9.878E-01  9.269E-02
 thal    2 0.00000043  2.319E+06  8.374E+07  1.087E+00  9.878E-01  9.270E-02
 thal    3 0.00000079  2.319E+06  8.374E+07  1.087E+00  9.878E-01  9.270E-02
 thal    4 0.00000134  2.319E+06  8.375E+07  1.087E+00  9.878E-01  9.271E-02
 thal    5 0.00000217  2.318E+06  8.375E+07  1.088E+00  9.878E-01  9.272E-02
 thal    6 0.00000452  2.318E+06  8.376E+07  1.088E+00  9.879E-01  9.274E-02
 thal    7 0.00000827  2.317E+06  8.377E+07  1.088E+00  9.879E-01  9.276E-02
 thal    8 0.00001514  2.316E+06  8.379E+07  1.088E+00  9.879E-01  9.280E-02
 thal    9 0.00002772  2.314E+06  8.382E+07  1.088E+00  9.880E-01  9.286E-02
 thal   10 0.00005411  2.312E+06  8.385E+07  1.088E+00  9.881E-01  9.294E-02

This edit lists the current values of each zonal quantity named in the header for the range of zones specified (Here, zones 1 through 10. The default is all current zones.) The first column, labeled “zone” contains the zones number and a convection sentinel, icon $(j)$ , that describes the status of convection between zones $j$ and $j+1$ . Possible values of the convection sentinel are describe with the description of icon.

The second column, labeled “mass” gives the mass interior to the inner boundary of each zone in units of scalem (p 273) (usually solar masses). This is converted from the internal KEPLER array ym $(j)$ , the mass exterior to the outer boundary of each zone in units of g.

The following columns list the values of the edit variables specified on the command line for these zones.

Graphics Edits¶

Graphics Edit Commands describes the graphics edits that the user can ask KEPLER to make, whereas Device-Control Graphics Parameters describes the graphics parameters that control many aspects of the resulting plots. Generally the resulting plots are self-explanatory, but sometimes terse formats or labels are used to save space. These cases are explained below.

Standard Plot Header¶

Standard plots produced by the plot command begin with a three line header that briefly describes the current status of the problem.

s15s7b2     2670   1.43666511135682E+17   c12( 185)= 1.0000E+14     Fri Jun 21 15:54:22 2013
 R = 6.5565E+10   Teff = 5.9352E+03  L = 3.8006E+33   Iter = 8   Zb = 3361   inv = 3361
 Dc = 1.4910E+02   Tc = 1.5516E+07  Ln = 1.0421E+32   Jm = 954   Etot = -3.091E+48

The first line lists the problem name nameprob, the current cycle number, ncyc (q 4), the current problem time, time (p 2) (sec), the current timestep, dt (q 1) (sec) along with the quantity and zone controlling it (see the description of the Time Edit for a more detailed discussion of these quantities), and the date and time when the plot was made.

The labeled quantities listed on the second and third lines have the meaning:

value	result
”`R`”	Photospheric radius, radius (q 47) (`cm`)
”`Teff`”	Photospheric temperature, teff (q 48) (`K`)
”`L`”	Non-neutrino surface luminosity, xlum (q 46) (`erg/sec`)
”`Iter`”	Number of KEPLER iterations it took the current cycle to converge, iter (q 7)
“`Dc`”	Density of the innermost zone, dn $(1)$ (`g/cc`)
”`Tc`”	Temperature of the innermost zone, tn $(1)$ (`K`)
”`Ln`”	Neutrino luminosity xlumn (q 25) (`erg/sec`)
”`Jm`”	Current number of zones jm (q 2)
“`Zb`”	Number of zones processed by BURN in the last time step, nburnz (q 51). If lburn (p 434) $>$ `1` this can be larger than the number of zones, jm (q 2), because BURN may be called multiple times for each zone, to obtain dreivatives and in case temperature and density of the zone change significantly during iterations. Only shown if BURN is used, inburn (p 267) $>$ `1`.
”`inv`”	Number of inversions done by BURN, ninvl (q 90). This should be $\ge$ nburnz (q 51). Only shown if BURN is used, inburn (p 267) $>$ `1`.

Note

Deprecated Since Version <15.0.0: The LOOK command generates plots with similar headers for previous cycles using infonnation stored in the “.lib” file(s).

Standard Zone Sentinels¶

The current zoning of problem and its convective state are indicated by a set of marks plotted at the position of the outer interface of each zone. Possible values of the convection sentinel are

value	result
(semicolon)	to indicate semi convection
”`\|`”	(vertical bar) to indicate full convection
(colon)	to indicate thermohaline regions (controlled by alpth (p 425))
”`'`”	(apostrophe) to indicate that overshoot mixing is being done in zones adjacent to convective regions
”`,`”	(comma) to indicate that while still radiative, the zonal interface is within a small amount (quantified and controlled by fracneut (p 59)) of being convective
”`.`”	(period) to indicate that the interface is well into the regime where only radiative transfer takes place across the interface.

Thermodvnamics Plot (Plot Type 1)¶

This plot plots temperature (“ $T$ ”, tn in K), the density (“ $\rho$ ”, dn in g/cc), and the standard zone sentinels described above. It also plots, in the units given, the pressure (“ $P$ ”, pn), net nuclear energy generation rate (if positive), “ $S$ ” (sgain, erg/g/sec), rate of energy loss including neutrino losses and net nuclear energy generation if negative, “ $L$ ” (sloss, erg/g/sec), and either the radius, “ $R$ ” (rn, cm, if the the radius is not already being used as the radial coordinate) or otherwise the velocity, un (cm/sec). Positive (outward) velocities are plotted as a curve labeled “ $U$ ”, and negative (inward) velocities are labeled “ $V$ ”.

Velocity Plot (Plot Type 2)¶

This plot should by self-explanatory. See above for an explanation of its header.

Ion Abundance Plot (Plot Tvpe 3)¶

This plot should be largely self-explanatory given the infonnation on headers and zone sentinels discussed above. Here “elemental” mass fraction means “ion” mass fraction. A detailed explanation of the ions corresponding to the symbols on the curves is given here.

Note

Pairs of adjacent points in different networks are not joined by lines.

Entropy Plot (Plot Type 4)¶

This plot graphs

value	result
unlabeled curve	total entropy sig
“ $E$ ”	the electron entropy, sige
“ $I$ ”	the ion entropy, sigi
“ $R$ ”	radiation entropy, sigr
“ $P$ ”	the positron entropy, sigp
“ $Z$ ”	the entropy due to ionization, sigion

Note

Currently these entropies are calculated assuming an ideal (but possibly degenerate and relativistic) gas. Probably this should be changed to using partial derivatives of the actual equation of state routines.

Rho-T Plot (Plot Type 5)¶

This plot should by self-explanatory. See above for an explanation of its header.

BURN Isotope Plot (Plot Type 6)¶

This plot should by self explanatory. The specific isotopes plotted are set by command setiso. See above for an explanation the meaning of its header and the zone sentinels shown. A sample list of isotopes included in the BURN network is given here.

Zonal Edit Plots (MONGO Plots)¶

Zonal variable plots made by the mongo and mon commands are titled by the first line of the standard graphics header described above. Otherwise they should be self-explanatory.

Time Plots¶

Deprecated Since Version <15.0.0.

Time plots made using the timeplot or tp commands are titled by the problem name, the current cycle number, and the reference time in seconds used in making the plot. To avoid truncation errors, any offsets involved in calculating the reference time (see subroutine plotit) are shown as additive quantities. Otherwise these plots should be self-explanatory.

Cycle Edits¶

Examples of many of the more complex edits that KEPLER can make to show the status of the calculation at a given cycle are listed below, together with an explanation of their formats.

These edits are produced automatically in the ASCII output fIle at intervals during the calculation depending on the settings of various edit control parameters. They can also be commanded to appear immediately on the terminal or in the ASCII output file by ASCII edit commands.

Cycle Edit Header¶

Wed Oct 31 03:42:03 2012     box a02       s1         weaver      53170    3.90529289311040E+17    ne20( 140) = 1.56E-17   iter =137

                                                      dtr( 602) = 2.56E+04
                                                      dtt( 134) = 6.29E+03                                               1200 zones
                                                      dtd( 134) = 1.58E+04
                                                      dtq( 736) = 5.87E+06
                                                     ne20( 140) = 1.87E+03
                                                      dtl( 669) = 3.67+103
                                                                                                                    mass    0.97640
    energy    source    kinetic    internal  potential    total     produced   escaped     check              luminosity  0.000E+00
   initial             0.000E+00  1.652E+47 -2.889E+47 -1.237E+47                                           neutrino lum  1.285E+34
   present  0.000E+00  1.986E+39  1.425E+49 -2.798E+49 -1.373E+49  4.174E+51  4.188E+51  1.810E+45              eff temp  0.000E+00
     rates  0.000E+00 -5.915E+33 -2.347E+43  4.499E+43  2.153E+43  2.140E+43 -1.247E+41 -1.401E+31                radius  1.241E+13

  reaction        p+p    he3+he3    he3+he4    c12(pg)    n14(pg)    o16(pg)     3a-c12    n14(ag)    c12(ag)    c12+c12    c12+o16
     moles  9.620E+31  4.098E+31  1.399E+31  1.216E+32  1.217E+32  6.014E+29  7.952E+29  6.214E+28  2.525E+25  3.032E+00  0.000E+00
   moles/s  7.472E+11  4.385E+13  9.478E+11  4.151E+15  2.795E+15  1.154E+13  2.228E+24  4.477E+23  2.513E+19  4.624E-05  0.000E+00
   maxzone        673        691        687        671        671        674        140        140        140        140          0

  reaction    o16+o16    o16(ag)   ne20(ag)   mg24(ag)   si28(ag)    s32(ag)   ar36(ag)   ca40(ag)   ti44(ag)   cr48(ag)   fe52(ag)
     moles  7.875E-31  1.239E+19  2.910E+19  3.392E+17  1.698E+13  1.266E+10  1.273E+06  4.111E-03  6.052E-04  2.987E-10  1.876E-56
   moles/s  1.067E-35  1.212E+13  4.071E+14  3.065E+12  1.488E+08  1.498E+05  1.588E+01  6.332E-08  8.335E-09  4.866E-15  6.127E-61
   maxzone        140        140        140        140        140        140        140        140        140        140        140

  reaction   fe52:2ng   fe52:a2p   fe54:2pg   he4-2n2p   p(e-nu)n   n(e+nu)p       weak
     moles  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
   moles/s  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
   maxzone          0          0          0          0          0          0          0

In the first line:

value	result
”`Wed Oct 03:42:03 2012`”	the current time and date.
”`box a02`”	the output box number, `craybox`.
”`s1`”	the name of the problem, `nameprob`.
”`weaver`”	the ID-word specified in the generator, `idword`.
”`53170`”	the current KEPLER cycle number, ncyc (q 4).
”`1.43960996754210E+07`”	the current problem time time (p 2) (`sec`)
”`ne20`”	the symbol of the ion or other quantity whose rate of change is deternlining the timestep, `idtcsym` $($ idtcon (q 13) $)$ .
”`140`”	the number of the zone where this rate of change is most time-step-limiting `jdtc` $($ idtcon (q 13) $)$
”`1.52E+03`”	the value of the current timestep, dt (q 1) (`sec`)
”`iter = 10`”	the number of iterations that KEPLER took to converge during the latest cycle, iter (q 7)

The next six lines give the limiting timesteps (DTC $(i)$ , in sec) due to

value	result
”`dtr`”	fractional changes in zonal radius
”`dtt`”	fractional changes in temperature
”`dtd`”	fractional changes in density
”`dtq`”	linear contraction
”`ne20`”	network ion abundance (“`ne20`” or other ion symbol)
”`dtl`”	fractional changes in luminosity

In KEPLER these are stired in idtcsym $(1\ldots6)$ . The zone limiting the timestep for each quantity, jdtc $(i)$ , is given in parentheses.

Along the right-hand margin is written the current number of zones in the problem (“1200’), followed by the current total mass in the problem, totm (q 17) $/$ scalem (p 273), in scalem (p 273) units (usually solar masses), the surface luminosity, xlum (q 46) (erg/sec), the surface neutrino luminosity, xlumn (q 25) (erg/sec), the effective surface temperature, teff (q 48) (K), and the photospheric radius, radius (q 47) (cm).

The energy edit given to the left of these quantities represents various kinds of energies and their rate of change summed over the whole star. The initial and present energies are in erg and the rates are in erg/sec. This edit summarizes in tabular form edit parameters enini (q 26) through encd (q 45).

There follows a table of all editable nuclear reaction rates, summed over the whole star. The reaction symbols, SYMR $(i)$ , in the row labeled reactions are defined (to the extent they aren’t obvious) in Section Zonal Nuclear Reaction Rates. The next row, labeled ‘moles’ gives the total net number of times in moles, totalr $(i)$ , that each reaction has taken place in the star since the beginning of the problem.

Note

Here the total number of backward reactions has been subtracted from the total number of forward reactions.

The row labeled “moles/s” gives the current net rate of each reaction in molels, rater $(i)$ , and the last row lists the zone where the rate of each reaction is currently the largest.

Primary Zonal Cycle Edit¶

      zone      mass    zmass    radius     density     temp      energy      sdot      snuc     sburn    abar    lumin    velocity
 conv    1    0.00000  0.00200  7.378E+08  2.365E+03  2.850E+08  4.211E+16  5.52E+05  5.53E+05  0.00E+00  14.8  2.241E+36 -2.760E-04
 conv    2    0.00200  0.00208  9.367E+08  2.355E+03  2.846E+08  4.203E+16  5.38E+05  5.39E+05  0.00E+00  14.8  4.517E+36 -3.503E-04
 conv    3    0.00408  0.00262  1.106E+09  2.346E+03  2.842E+08  4.195E+16  5.25E+05  5.26E+05  0.00E+00  14.8  7.314E+36 -4.135E-04
 conv    4    0.00670  0.00330  1.265E+09  2.337E+03  2.837E+08  4.187E+16  5.12E+05  5.13E+05  0.00E+00  14.8  1.074E+37 -4.727E-04
 conv    5    0.01000  0.00500  1.449E+09  2.325E+03  2.832E+08  4.177E+16  4.96E+05  4.97E+05  0.00E+00  14.8  1.579E+37 -5.416E-04
 conv    6    0.01500  0.00500  1.597E+09  2.313E+03  2.826E+08  4.166E+16  4.80E+05  4.81E+05  0.00E+00  14.8  2.067E+37 -5.966E-04
 conv    7    0.02000  0.00500  1.721E+09  2.302E+03  2.821E+08  4.156E+16  4.65E+05  4.67E+05  0.00E+00  14.8  2.542E+37 -6.431E-04
 conv    8    0.02500  0.00500  1.831E+09  2.292E+03  2.816E+08  4.147E+16  4.52E+05  4.54E+05  0.00E+00  14.8  3.003E+37 -6.838E-04
 conv    9    0.03000  0.00750  1.975E+09  2.280E+03  2.810E+08  4.137E+16  4.37E+05  4.38E+05  0.00E+00  14.8  3.672E+37 -7.373E-04
 conv   10    0.03750  0.00750  2.101E+09  2.266E+03  2.803E+08  4.125E+16  4.21E+05  4.22E+05  0.00E+00  14.8  4.317E+37 -7.842E-04
 conv   11    0.04500  0.00750  2.214E+09  2.254E+03  2.797E+08  4.114E+16  4.06E+05  4.07E+05  0.00E+00  14.8  4.940E+37 -8.262E-04
 conv   12    0.05250  0.00750  2.317E+09  2.242E+03  2.792E+08  4.103E+16  3.93E+05  3.94E+05  0.00E+00  14.8  5.542E+37 -8.645E-04
 conv   13    0.06000  0.00800  2.418E+09  2.230E+03  2.786E+08  4.093E+16  3.79E+05  3.81E+05  0.00E+00  14.8  6.164E+37 -9.020E-04
 conv   14    0.06800  0.00800  2.512E+09  2.219E+03  2.780E+08  4.082E+16  3.67E+05  3.68E+05  0.00E+00  14.8  6.765E+37 -9.368E-04
 conv   15    0.07600  0.00800  2.600E+09  2.208E+03  2.775E+08  4.072E+16  3.55E+05  3.56E+05  0.00E+00  14.8  7.347E+37 -9.693E-04
 conv   16    0.08400  0.00800  2.683E+09  2.197E+03  2.769E+08  4.063E+16  3.44E+05  3.45E+05  0.00E+00  14.8  7.912E+37 -9.999E-04
 conv   17    0.09200  0.00800  2.761E+09  2.186E+03  2.764E+08  4.053E+16  3.33E+05  3.34E+05  0.00E+00  14.8  8.460E+37 -1.029E-03
 conv   18    0.10000  0.00800  2.835E+09  2.176E+03  2.759E+08  4.044E+16  3.23E+05  3.24E+05  0.00E+00  14.8  8.991E+37 -1.056E-03
 conv   19    0.10800  0.00800  2.906E+09  2.166E+03  2.754E+08  4.035E+16  3.14E+05  3.15E+05  0.00E+00  14.8  9.508E+37 -1.082E-03
 conv   20    0.11600  0.00800  2.974E+09  2.157E+03  2.750E+08  4.026E+16  3.05E+05  3.06E+05  0.00E+00  14.8  1.001E+38 -1.107E-03

...

This edit lists the current values of each quantity named in the header for each zone currently in the problem.

The first column, labeled “zone” contains the zones number and a convection sentinel icon $($ j $)$ , that describes the status of convection between zones j and j $+$ 1. Possible values of the convection sentinel are describe with the description of icon.

The second column, labeled “mass” gives the mass interior to the inner boundary of each zone in units of scalem (p 273) (usually solar masses).

Third column, labeled “zmass” gives the mass of each zone, also in scalem (p 273) units. (These are converted from the internal KEPLER array ym $(j)$ , the mass exterior to the outer boundary of each zone in units of g) and from the zonal mass, xm $(j)$ in units of g.

The next ten columns correspond to the variables listed below:

value	result
”`radius`”	rn, Radius of outer zone boundary (`cm`).
”`density`”	dn, Density (`g/cc`).
”`temp`”	tn, Temperature (`K`).
”`energy`”	en, Specific energy density (`erg/g`).
”`sdot`”	sn, Total energy generation rate including neutrino losses (`erg/g/sec`).
”`snuc`”	snn, Nuclear energy generation rate excluding neutrino losses (`erg/g/sec`)
”`iosburn`”	sburn, Nuclear energy generation rate, excluding neutrino losses (`erg/g/sec`).
”`abar`”	abar, Mean atomic weight (`g/mol`).
”`lumin`”	xln, Calculated non-neutrino luminosity at the outer zone boundary (`erg/sec`).
”`velocity`”	un, Velocity of the outer zone boundary (`cm/sec`).

Secondary Zonal Cycle Edit¶

      zone        mass    entropy    mixtime      eta     pressure    stress     opacity     ye     yeburn  limnuc(cyczb)    dtburn
         1 n    0.00000  9.381E-01  8.070+115  5.754E+00  5.405E+25  0.000E+00  1.140E-02  0.45156  0.47897  mg25(112278)  2.759E-05
         2 n    0.00200  9.413E-01  5.027+115  5.720E+00  5.297E+25  1.655E+36  1.147E-02  0.45175  0.47897  mg25(112278)  3.901E-05
         3 n    0.00408  9.442E-01  3.628+115  5.689E+00  5.199E+25  3.989E+36  1.154E-02  0.45193  0.47897  mg25(112278)  5.425E-05
         4 n    0.00670  9.473E-01  2.794+115  5.657E+00  5.096E+25  8.132E+36  1.162E-02  0.45212  0.47897  mg25(112278)  8.481E-05
         5 n    0.01000  9.512E-01  2.149+115  5.617E+00  4.972E+25  1.930E+37  1.172E-02  0.45236  0.47897  mg25(112279)  1.624E-04
         6 n    0.01500  9.552E-01  1.783+115  5.576E+00  4.843E+25  2.809E+37  1.182E-02  0.45263  0.47897  mg25(112279)  4.555E-04
         7 n    0.02000  9.582E-01  1.621+115  5.546E+00  4.750E+25  2.215E+37  1.190E-02  0.45283  0.47897  mg25(112279)  1.149E-03
         8 n    0.02311  9.603E-01  1.511+115  5.525E+00  4.687E+25  2.181E+37  1.195E-02  0.45296  0.47897  mg25(112279)  2.832E-03
         9 n    0.02577  9.622E-01  1.416+115  5.507E+00  4.631E+25  2.582E+37  1.200E-02  0.45309  0.47897  mg25(112279)  7.275E-03
        10 n    0.02855  9.643E-01  1.314+115  5.486E+00  4.567E+25  3.787E+37  1.205E-02  0.45324  0.47897   p30(112279)  3.274E-02
        11 n    0.03213  9.665E-01  1.232+115  5.464E+00  4.501E+25  4.126E+37  1.211E-02  0.45339  0.47897   p30(112279)  1.324E-01
        12 n    0.03556  9.687E-01  1.158+115  5.443E+00  4.436E+25  4.941E+37  1.217E-02  0.45355  0.47897  ca43(112279)  2.332E-01
        13 n    0.03921  9.709E-01  1.099+115  5.422E+00  4.375E+25  5.074E+37  1.223E-02  0.45370  0.47897  ca43(112279)  3.058E-01
        14 n    0.04258  9.729E-01  1.043+115  5.402E+00  4.315E+25  6.236E+37  1.228E-02  0.45385  0.47897  ca43(112279)  1.299E+00
        15 n    0.04633  9.751E-01  9.940+114  5.382E+00  4.255E+25  6.745E+37  1.234E-02  0.45400  0.47897  ca43(112279)  4.912E+00
        16 n    0.05000  9.775E-01  9.375+114  5.359E+00  4.187E+25  1.020E+38  1.241E-02  0.45418  0.47897   p30(112280)  2.634E+00
        17 n    0.05500  9.796E-01  9.106+114  5.340E+00  4.131E+25  5.488E+37  1.246E-02  0.45433  0.47897  ar39(112280)  1.616E+01
        18 n    0.05747  9.808E-01  8.890+114  5.328E+00  4.096E+25  5.278E+37  1.250E-02  0.45442  0.47897  si32(112280)  2.377E+02
        19 n    0.05973  9.825E-01  8.566+114  5.313E+00  4.053E+25  9.492E+37  1.254E-02  0.45454  0.47897  si32(112280)  4.075E+02
        20 n    0.06355  9.844E-01  8.281+114  5.295E+00  4.001E+25  9.550E+37  1.259E-02  0.45468  0.47897  si32(112280)  3.416E+02
        21 n    0.06712  9.862E-01  8.032+114  5.279E+00  3.953E+25  9.740E+37  1.264E-02  0.45481  0.47897  si32(112280)  5.788E+02
        22 n    0.07053  9.880E-01  7.787+114  5.262E+00  3.906E+25  1.118E+38  1.269E-02  0.45495  0.47897  cl36(112280)  4.764E+02
        23 n    0.07421  9.898E-01  7.565+114  5.245E+00  3.858E+25  1.168E+38  1.274E-02  0.45508  0.47897   o16(112280)  7.308E+02
        24 n    0.07782  9.917E-01  7.353+114  5.229E+00  3.812E+25  1.280E+38  1.280E-02  0.45521  0.47897   o16(112280)  7.678E+02
        25 n    0.08155  9.935E-01  7.159+114  5.213E+00  3.766E+25  1.334E+38  1.285E-02  0.45535  0.47897   o16(112280)  1.136E+02

...

       229 i    0.97000  1.367E+00  2.145+114  2.676E+00  2.472E+24  7.342E+39  2.642E-02  0.46611  0.48953   nt0(115800)  6.613E+04
       230 i    0.97500  1.369E+00  2.139+114  2.666E+00  2.420E+24  7.435E+39  2.654E-02  0.46613  0.49030   nt0(116000)  7.738E+04
       231 i    0.98000  1.370E+00  2.133+114  2.657E+00  2.369E+24  7.529E+39  2.666E-02  0.46616  0.49108  si31(116300)  5.268E+04
 neut  232 i    0.98500  1.372E+00  2.128+114  2.648E+00  2.318E+24  7.626E+39  2.678E-02  0.46622  0.49186   nt0(117000)  2.205E+04
 semi  233 i    0.99000  1.374E+00  4.646E+12  2.642E+00  2.268E+24  7.732E+39  2.690E-02  0.46647  0.49283   nt0(118500)  1.268E+04
       234      0.99500  1.656E+00  2.235+114  3.520E+00  2.219E+24  1.065E+40  2.365E-02  0.50000  0.49313   s31(125000)  7.132E+02
       235      1.00000  1.688E+00  2.255+114  3.440E+00  2.170E+24  1.163E+40  2.383E-02  0.50000  0.49406   s31(125000)  7.427E+02
       236      1.00500  1.726E+00  2.284+114  3.334E+00  2.122E+24  1.373E+40  2.409E-02  0.50000  0.49477   s31(125000)  8.436E+02
       237      1.01000  1.773E+00  2.332+114  3.180E+00  2.075E+24  2.020E+40  2.448E-02  0.50000  0.49545   s35(125000)  6.500E+02
 osht  238      1.01500  1.788E+00  2.965E+12  3.129E+00  2.045E+24  6.237E+39  2.465E-02  0.50000  0.49606   nt0(125000)  1.249E+02
 conv  239      1.01638  2.188E+00  3.725E+02  1.345E+00  2.035E+24  2.829E+39  2.988E-02  0.50000  0.49694   he4(125000)  1.728E-09
 conv  240      1.01722  2.188E+00  3.188E+02  1.344E+00  2.029E+24  1.720E+39  2.989E-02  0.50000  0.49694   he4(125000)  1.859E-09

...

      1207     11.99999  2.182E+01  1.802+121 -3.452E+01  1.139E+01  0.000E+00  3.399E-04  0.80442  0.80301   nt0(125000)  9.437E+01
      1208     11.99999  2.189E+01  1.798+121 -3.460E+01  1.074E+01  0.000E+00  3.363E-04  0.80519  0.80380   nt0(125000)  9.437E+01
      1209     12.00000  2.192E+01  1.780+121 -3.473E+01  1.020E+01  0.000E+00  3.564E-04  0.80519  0.80380   nt0(125000)  9.437E+01

The columns labeled ‘zone’ and ‘mass’ list the same values as in the primary zonal edit, except that a network sentinel is also included in an unlabeled column between them. The network sentinel is

value	result
”`i`”	the zone is currently using the ISE network, netnum $(j) = 2$
”`n`”	the NSE network is being used netnum $(j) = 3$
’_’	(blank) if the APPROX network is being used netnum $(j) = 1$

The column labeled “limnuc(cyczb)” lists the isotope controlling the BURN timestep for each zone, along with the KEPLER cycle when this zone was last burned. This information is calculated from the KEPLER array limnuc, defined as:

limnuc $=$ nalim $+$ 100 $*$ nzlim $+$ 10000 $*$ ncyczb

where:

value	result
`nalim`	Atomic weight of the nucleus limiting the BURN zonal timestep.
`nzlim`	Atomic number of the nucleus limiting the BURN zonal timestep.
`ncyczb`	KEPLER cycle at which this zone was last coprocessed by BURN.

The remaining columns correspond to the variables listed below:

value	result
”`entropy`”	stot, Total zonal entropy (`kb/baryon`).
”`mixtime`”	taumix, Convective mixing time (`sec`), calculated from the convective mixing coefficient, difi $(j)$ and the zonal convective mixing length, mixlen by the equation `taumix` $=$ mixlen $^2$ $/$ difi $(j)$
”`eta`”	etan Electron degeneracy parameter ( $\mu_\mathrm{e}/k_\mathrm{B}T$ ).
”`pressure`”	pn Pressure (`erg/cc`).
”`stress`”	zn Viscous stress including linear and quadratic artificial viscosity (`erg`).
”`opacity`”	xkn Rosseland mean opacity (`cm**2/g`).
”`ye`”	`ye` Electron abundance in any current (non-BURN) network (`mol/g`).
”`yeburn`”	yeburn Electron abundance calculated for this zone by the BURN coprocessor (`mol/g`), derived from the BURN neutron excess, etab $(j)$ , by the equation: yeburn $=$ `0.5` $$ $($ `1.` $-$ etab* $(j)$ $)$ .
”`dtburn`”	dtimen..Current timestep for BURN coprocessing in this zone (`sec`).

Ion Abundance Cycle Edit¶

       zone      mass      nt1        ye         pn1        yq         he4        yf         eb0        o16        fe56       mg24
          1    0.00000  3.633E-09  4.513E-01  8.332E-10  1.755E-02  2.828E-06  1.000E-08  8.786E+00  0.000E+00  1.188E-01  1.561E-13
          2    0.00200  3.500E-09  4.514E-01  8.509E-10  1.755E-02  2.843E-06  1.000E-08  8.786E+00  0.000E+00  1.241E-01  1.562E-13
          3    0.00408  3.380E-09  4.516E-01  8.666E-10  1.755E-02  2.854E-06  1.000E-08  8.786E+00  0.000E+00  1.290E-01  1.560E-13
          4    0.00670  3.252E-09  4.517E-01  8.831E-10  1.755E-02  2.865E-06  1.000E-08  8.786E+00  0.000E+00  1.346E-01  1.557E-13
          5    0.01000  3.096E-09  4.519E-01  9.028E-10  1.755E-02  2.875E-06  1.000E-08  8.787E+00  0.000E+00  1.417E-01  1.550E-13

...

         96    0.72400  1.680E-12  4.650E-01  1.998E-09  1.785E-02  7.498E-07  1.000E-08  8.786E+00  0.000E+00  7.964E-01  2.524E-14
         97    0.73200  1.387E-12  4.652E-01  2.173E-09  1.784E-02  7.436E-07  1.000E-08  8.786E+00  0.000E+00  7.878E-01  2.749E-14
         98    0.74000  1.147E-12  4.654E-01  2.340E-09  1.784E-02  7.339E-07  1.000E-08  8.785E+00  0.000E+00  7.779E-01  2.959E-14
         99    0.74800  9.482E-13  4.657E-01  2.492E-09  1.785E-02  7.202E-07  1.000E-08  8.784E+00  0.000E+00  7.667E-01  3.143E-14
        100    0.75600  7.888E-13  4.658E-01  2.505E-09  1.785E-02  6.875E-07  1.000E-08  8.784E+00  0.000E+00  7.585E-01  3.094E-14
      total    0.76400  1.953E-11  0.000E+00  3.672E-11  0.000E+00  5.952E-08  0.000E+00  0.000E+00  0.000E+00  2.116E-02  1.996E-15

       zone      mass      si28       s32        ar36       ca40       ti44       cr48       'fe'       fe54       ni56       1-tot
          1    0.00000  1.151E-11  2.027E-08  1.642E-08  2.921E-09  1.607E-06  2.285E-02  8.579E-01  4.499E-04  3.257E-10  4.410E-13
          2    0.00200  1.222E-11  2.074E-08  1.694E-08  2.936E-09  1.560E-06  2.226E-02  8.532E-01  4.786E-04  3.634E-10  5.237E-13
          3    0.00408  1.289E-11  2.116E-08  1.742E-08  2.948E-09  1.517E-06  2.172E-02  8.487E-01  5.064E-04  4.013E-10  6.407E-13
          4    0.00670  1.365E-11  2.160E-08  1.794E-08  2.961E-09  1.471E-06  2.115E-02  8.437E-01  5.377E-04  4.461E-10  4.204E-13
          5    0.01000  1.462E-11  2.213E-08  1.858E-08  2.975E-09  1.416E-06  2.045E-02  8.373E-01  5.788E-04  5.080E-10  5.208E-13

...

         96    0.72400  1.208E-09  1.314E-08  3.900E-08  5.891E-09  3.884E-07  2.962E-04  1.392E-01  6.407E-02  1.403E-05 -2.676E-12
         97    0.73200  1.464E-09  1.360E-08  4.083E-08  6.748E-09  4.220E-07  3.179E-04  1.369E-01  7.496E-02  1.970E-05 -2.649E-12
         98    0.74000  1.746E-09  1.406E-08  4.244E-08  7.712E-09  4.552E-07  3.423E-04  1.346E-01  8.709E-02  2.708E-05 -2.613E-12
         99    0.74800  2.047E-09  1.450E-08  4.378E-08  8.781E-09  4.870E-07  3.687E-04  1.324E-01  1.004E-01  3.643E-05 -2.325E-12
        100    0.75600  2.199E-09  1.435E-08  4.337E-08  9.351E-09  4.974E-07  3.835E-04  1.306E-01  1.104E-01  4.396E-05 -2.364E-12
      total    0.76400  7.125E-12  6.880E-10  9.862E-10  1.032E-10  1.691E-08  1.899E-04  1.640E-02  4.503E-04  7.102E-08  9.618E-01

       zone      mass      nt1        ye         pn1        yq         he4        yf         eb0        o16        fe56       mg24
        101    0.76400  4.315E-13  4.661E-01  1.383E-09  1.825E-02  9.755E-08  1.820E-02  8.779E+00  2.247E-07  6.131E-01  5.888E-08
        102    0.77200  3.515E-13  4.664E-01  1.358E-09  1.827E-02  8.721E-08  1.820E-02  8.778E+00  2.533E-07  5.956E-01  8.640E-08
        103    0.78000  2.851E-13  4.668E-01  1.317E-09  1.830E-02  7.725E-08  1.819E-02  8.776E+00  2.860E-07  5.775E-01  1.277E-07

...

 thal   191    1.42800  2.341E-22  4.980E-01  2.994E-15  3.239E-02  3.177E-13  1.003E-08  8.464E+00  2.703E-02  2.569E-07  2.253E-08
 thal   192    1.43600  2.360E-22  4.980E-01  3.423E-15  3.214E-02  3.593E-13  1.002E-08  8.463E+00  3.099E-02  2.450E-07  2.028E-08
 thal   193    1.44400  2.399E-22  4.980E-01  3.856E-15  3.190E-02  4.028E-13  1.001E-08  8.462E+00  3.481E-02  2.354E-07  1.841E-08
        194    1.45200  2.428E-22  4.980E-01  4.429E-15  3.158E-02  4.604E-13  1.000E-08  8.461E+00  3.993E-02  2.239E-07  1.624E-08
      total    1.46000  8.746E-16  0.000E+00  7.180E-12  0.000E+00  3.162E-10  0.000E+00  0.000E+00  1.437E-04  6.968E-03  1.379E-08

       zone      mass      si28       s32        ar36       ca40       ti44       cr48       'fe'       fe54       ni56       1-tot
        101    0.76400  8.961E-04  5.523E-04  1.901E-04  1.000E-05  3.472E-05  4.057E-03  2.165E-01  1.646E-01  4.609E-05  2.298E-14
        102    0.77200  1.339E-03  7.289E-04  2.288E-04  1.273E-05  4.037E-05  4.540E-03  2.149E-01  1.825E-01  5.479E-05 -3.886E-14
        103    0.78000  2.011E-03  9.755E-04  2.782E-04  1.630E-05  4.661E-05  5.055E-03  2.134E-01  2.007E-01  6.386E-05  6.961E-14
        104    0.78800  3.010E-03  1.312E-03  3.398E-04  2.084E-05  5.326E-05  5.587E-03  2.118E-01  2.189E-01  7.300E-05 -7.772E-14

...

 thal   191    1.42800  5.585E-01  3.172E-01  7.610E-02  2.115E-02  5.271E-07  1.263E-10  1.391E-08  2.798E-07  2.598E-12 -6.586E-12
 thal   192    1.43600  5.387E-01  3.258E-01  7.940E-02  2.506E-02  6.214E-07  1.229E-10  1.270E-08  2.923E-07  2.954E-12 -6.491E-12
 thal   193    1.44400  5.202E-01  3.331E-01  8.285E-02  2.909E-02  7.271E-07  1.192E-10  1.177E-08  3.021E-07  3.277E-12 -5.648E-12
        194    1.45200  4.961E-01  3.415E-01  8.744E-02  3.498E-02  8.915E-07  1.144E-10  1.073E-08  3.138E-07  3.717E-12 -5.576E-12
      total    1.46000  1.064E-02  4.638E-03  9.791E-04  2.114E-04  9.289E-07  1.228E-04  2.324E-03  8.773E-03  2.511E-06  9.652E-01

       zone      mass      nt1        h1         pn1        he3        he4        c12        n14        o16        ne20       mg24
        195    1.46000  5.202E-32  0.000E+00  2.706E-17  0.000E+00  3.006E-14  1.420E-09  0.000E+00  4.530E-02  1.119E-07  2.092E-05
        196    1.46800  3.833E-32  0.000E+00  1.513E-17  0.000E+00  2.766E-14  1.531E-09  0.000E+00  5.127E-02  1.511E-07  2.326E-05
        197    1.47600  3.239E-32  0.000E+00  7.758E-18  0.000E+00  2.534E-14  1.738E-09  0.000E+00  5.798E-02  2.116E-07  2.580E-05

...

       1238   19.99994  0.000E+00  7.512E-01  0.000E+00  1.531E-05  2.488E-01  2.925E-14  4.609E-13  2.407E-15  9.980E-24  4.180E-27
       1239   19.99996  0.000E+00  7.512E-01  0.000E+00  1.531E-05  2.488E-01  2.925E-14  4.609E-13  2.407E-15  5.065E-25  1.837E-28
       1240   19.99997  0.000E+00  7.512E-01  0.000E+00  1.531E-05  2.488E-01  2.926E-14  4.609E-13  2.407E-15  7.537E-26  1.821E-30
       1241   19.99999  0.000E+00  7.512E-01  0.000E+00  1.531E-05  2.488E-01  2.926E-14  4.609E-13  2.406E-15  7.104E-26  9.317E-33
      total   20.00000  2.846E-30  4.675E-01  1.291E-17  9.509E-07  2.994E-01  6.537E-03  3.685E-03  1.050E-01  7.873E-03  5.612E-03

       zone      mass      si28       s32        ar36       ca40       ti44       cr48       fe52       fe54       ni56       1-tot
        195    1.46000  4.551E-01  3.978E-01  6.448E-02  3.733E-02  3.234E-06  1.268E-08  9.355E-12  7.307E-16  8.451E-21 -2.764E-08
        196    1.46800  4.503E-01  3.995E-01  6.540E-02  3.348E-02  2.895E-06  9.371E-09  5.919E-12  4.086E-16  5.194E-21 -2.764E-08
        197    1.47600  4.457E-01  4.003E-01  6.619E-02  2.977E-02  2.559E-06  6.767E-09  3.545E-12  2.095E-16  2.968E-21 -2.764E-08
        198    1.48400  4.412E-01  4.001E-01  6.678E-02  2.625E-02  2.236E-06  4.787E-09  1.981E-12  9.372E-17  1.552E-21 -2.764E-08

...

       1237   19.99992  2.695E-32  3.082E-42  6.098E-54  1.315E-68  1.618E-89  0.000E+00  0.000E+00  0.000E+00  0.000E+00  2.887E-15
       1238   19.99994  1.106E-32  1.265E-42  2.503E-54  5.397E-69  6.641E-90  0.000E+00  0.000E+00  0.000E+00  0.000E+00 -4.441E-16
       1239   19.99996  4.861E-34  5.559E-44  1.100E-55  2.372E-70  2.918E-91  0.000E+00  0.000E+00  0.000E+00  0.000E+00 -1.776E-15
       1240   19.99997  4.820E-36  5.512E-46  1.091E-57  2.352E-72  2.893E-93  0.000E+00  0.000E+00  0.000E+00  0.000E+00 -1.998E-15
       1241   19.99999  2.466E-38  2.820E-48  5.580E-60  1.203E-74  1.466E-95  0.000E+00  0.000E+00  0.000E+00  0.000E+00 -3.775E-15
      total   20.00000  1.800E-02  1.098E-02  1.908E-03  5.639E-04  2.477E-07  3.475E-10  5.692E-13  3.486E-16  1.034E-20  7.300E-02

The columns labeled ‘zone’ and ‘mass’ list the same values as in the primary zonal edit.

The following columns list the zonal mass fraction of the first ten ions in the current (non-BURN) nuclear network being employed in each zone, with each column labeled with the appropriate network ion symbol. The special cases where the ion quantity listed is not an ionic mass fraction are as follows:

value	result
”`ye`”	`ye`, electron abundance in any current (non-BURN) network (`mol/g`).
”`yq`”	`yq`, total abundance of ISE or NSE isotopes with $A \ge 24$ (`mol/g`).
”`yf`”	`yf`, total abundance of ISE or NSE isotopes with $A > 46$ plus $^{46}\mathrm{Ti}$ (`mol/g`).
”`eb0`”	`eb0`, mean binding energy per nucleon of all ISE or NSE isotopes (`MeV`).

After all the zones with a (non-BURN) network the same as the first zone are listed, a fmal row labeled “total” in the “zone” column summarizes the total mass fraction of each ion (as summed over the entire same-network block of zones) relative to the current total mass, totm (q 17) of the all the zones in the star. The column labeled “mass” in this row gives the total mass in the star interior to the outer boundary of this block of zones in scalem (p 273) units (usually solar masses).

When the edit described above for the first ten ions of the network employed by the first zone is complete, a essentially similar edit table is begun containing the same “zone” and “mass” columns, together with the next ten ions in that network, and so on until all the ions have been listed. A final column labeled “code:1-tot is appended that lists the amount by which the sum of the mass fractions of all ions in each zone differs from unity. Note that the total mass fraction given in the “1-tot” column for the “total” row includes all mass not accounted for by this same-network region, either due to mass-nonconservation in nuclear burning (generally due to zeroing small negative abundances) or to mass being present in other same-network regions.

When the edit for the first same-network block of zones is complete, edits in the same format are made for any other same-network region(s) that may be present.

ISE (or NSE) Zone-Cycle Edit¶

 zone  masscoord  cycle         time             dt        temp     density       ye        dye         yq        dyq       iternse
   1    0.000000  125000  2.135910012304E+06  3.557E+01  3.504E+09  2.015E+08  4.516E-01 -1.875E-06  1.744E-02  2.310E-15          1

    yp         yn        ysi        y56         f0         g0         h0         q0       qxc12      qxo16       qxne       betaq
 5.034E-11  2.015E-10  1.000E+00  1.000E+00 -1.611E-13 -7.244E-14  0.000E+00  0.000E+00  2.025E-16  7.216E-17  2.229E-20  1.000E+00

  wrate      rectot     rpdtot     redtot      snuw      snubps     eectot     epdtot     eedtot     ebind0     ebind     eexcite
 5.272E-08  1.091E-07  8.823E-12  5.637E-08  2.352E+11  2.081E+11  1.473E+11  0.000E+00  8.792E+10  8.783E+00  8.783E+00  3.395E-03

           k  isotope     xk         wk        rtot       rec        rpd        red        eec        etot       eex
           1     n1    2.015E-10  2.000E+00 -1.603E-16 -7.958E-07  0.000E+00  0.000E+00  2.332E-06  4.531E+02  0.000E+00
           2     h1    5.034E-11  2.000E+00  3.424E-12  6.802E-02  0.000E+00  0.000E+00  1.004E-01  4.876E+06  0.000E+00
          68     ti50  9.927E-03  1.032E+00 -1.585E-12  4.110E-10  1.803E-15  8.396E-09  5.019E-10  9.180E+05  4.653E-02
          69     v50   3.158E-06  2.666E+01  8.226E-11  1.295E-03  9.597E-06  2.851E-06  4.113E-03  2.510E+08  2.249E-01
          71     ti51  1.926E-05  4.156E+00 -1.869E-10  6.958E-10  1.263E-15  4.949E-04  8.065E-10  9.213E+07  5.254E-02
          72     v51   3.157E-03  1.032E+01  1.053E-09  1.713E-05  6.027E-11  1.071E-07  2.065E-05  1.236E+09  1.085E-01
          73     cr51  4.172E-06  8.776E+00  4.103E-10  5.015E-03  8.445E-07  1.910E-10  1.218E-02  9.617E+08  8.791E-02
          76     v52   1.489E-04  4.515E+01 -6.218E-09  3.211E-05  2.702E-11  2.203E-03  3.743E-05  3.215E+09  1.257E-01
          77     cr52  3.392E-02  1.049E+00  2.595E-09  3.978E-06  1.757E-11  1.564E-11  4.416E-06  2.779E+09  7.313E-02
          78     mn52  1.725E-09  1.696E+01  7.091E-13  2.081E-02  5.675E-04  2.679E-10  8.462E-02  2.717E+06  1.699E-01
          80     v53   2.741E-05  1.524E+01 -2.146E-09  1.193E-08  1.986E-14  4.150E-03  1.402E-08  1.058E+09  1.568E-01
          81     cr53  1.321E-02  4.721E+00  2.753E-09  1.141E-05  1.859E-11  3.633E-07  1.241E-05  3.029E+09  1.418E-01
          82     mn53  7.173E-05  9.875E+00  1.447E-09  1.069E-03  2.746E-07  3.738E-11  1.880E-03  2.454E+09  1.020E-01
          85     cr54  1.157E-01  1.322E+00 -2.898E-10  7.630E-10  1.496E-15  1.360E-07  8.441E-10  1.454E+08  2.204E-01
          86     mn54  2.761E-04  1.609E+01  5.452E-09  1.068E-03  8.603E-07  2.088E-06  2.580E-03  1.273E+10  1.577E-01
          89     cr55  1.973E-04  7.138E+00 -3.226E-09  6.568E-10  7.307E-15  8.991E-04  7.283E-10  1.590E+09  2.558E-01
          90     mn55  3.520E-02  1.221E+01  5.773E-09  9.436E-06  5.329E-11  4.152E-07  1.087E-05  6.840E+09  1.497E-01
          91     fe55  1.244E-04  4.962E+00  4.264E-09  1.885E-03  2.471E-07  3.557E-10  3.889E-03  8.487E+09  1.410E-01
          94     mn56  7.234E-04  2.675E+01 -2.826E-08  5.692E-05  9.158E-11  2.244E-03  7.350E-05  1.521E+10  2.092E-01
          95     fe56  1.361E-01  1.301E+00  1.660E-08  6.832E-06  2.392E-11  3.871E-11  7.670E-06  1.798E+10  2.080E-01
          96     co56  7.919E-09  1.500E+01  1.667E-12  1.161E-02  1.856E-04  8.616E-10  2.975E-02  4.071E+06  1.681E-01
          98     mn57  1.696E-04  6.932E+00 -5.200E-09  1.156E-09  5.367E-14  1.747E-03  1.323E-09  2.564E+09  1.263E-01
          99     fe57  3.970E-02  1.166E+01  2.072E-08  2.985E-05  4.930E-11  1.105E-07  3.360E-05  2.262E+10  1.594E-01
         100     co57  3.313E-05  8.298E+00  6.421E-10  1.104E-03  4.346E-07  2.573E-10  2.273E-03  1.275E+09  5.198E-02
         102     mn58  5.784E-07  1.233E+01 -6.736E-10  2.086E-08  1.214E-13  6.754E-02  2.302E-08  3.321E+08  3.960E-01
         103     fe58  4.243E-01  1.358E+00 -1.106E-10  5.697E-09  1.978E-14  2.081E-08  6.181E-09  1.187E+08  2.371E-01
         104     co58  1.364E-04  4.877E+01  2.326E-09  9.887E-04  8.571E-07  7.292E-07  1.886E-03  4.282E+09  1.163E-01
         106     fe59  2.306E-03  7.728E+00 -9.295E-10  3.761E-08  6.811E-14  2.382E-05  4.092E-08  4.606E+08  2.903E-01
         107     co59  8.405E-03  8.458E+00  4.580E-09  3.230E-05  1.073E-09  1.526E-07  4.565E-05  6.284E+09  7.187E-02
         108     ni59  3.647E-06  6.860E+00  2.742E-10  4.436E-03  3.702E-07  1.992E-12  1.006E-02  5.998E+08  2.251E-01
         110     fe60  7.783E-04  1.504E+00 -5.610E-11  4.861E-13  1.595E-17  4.325E-06  4.995E-13  2.766E+07  2.580E-01
         111     co60  6.277E-04  3.238E+01 -3.960E-09  6.907E-05  5.829E-09  4.476E-04  1.399E-04  3.721E+09  2.660E-01
         112     ni60  4.462E-03  1.069E+00  2.235E-09  3.005E-05  4.231E-10  4.004E-13  3.467E-05  2.488E+09  9.270E-02
         115     co61  1.351E-03  8.343E+00 -1.033E-09  4.711E-07  1.953E-12  4.714E-05  5.676E-07  5.268E+08  5.290E-02
         116     ni61  4.209E-03  1.746E+01  3.523E-08  5.106E-04  1.165E-09  1.276E-09  7.523E-04  5.008E+10  1.278E-01
         118     co62  3.740E-06  1.758E+01 -1.712E-09  5.367E-06  8.136E-12  2.839E-02  6.099E-06  8.447E+08  3.389E-01
         119     ni62  1.528E-01  1.109E+00  3.803E-10  1.544E-07  4.596E-13  8.660E-11  1.719E-07  4.087E+08  1.221E-01
         120     co63  7.278E-07  8.184E+00 -1.669E-10  1.807E-09  5.961E-16  1.445E-02  1.942E-09  2.469E+07  3.343E-02
         121     ni63  3.620E-03  1.174E+01  8.940E-11  1.761E-06  2.705E-12  2.049E-07  1.927E-06  1.126E+08  2.713E-01
         123     ni64  8.241E-03  1.061E+00 -4.958E-12  5.372E-10  3.634E-16  3.904E-08  5.535E-10  2.548E+06  7.875E-02

   snucq     snucq1      sdotq     tauqse       yf        yfqse      yffac       dyf        qj        yo160     dyo16q       eta
-3.332E+09 -2.115E+11 -4.466E+11  1.271E+02  1.742E-02  1.742E-02  1.074E+11  0.000E+00 -6.496E-17  0.000E+00  0.000E+00  5.754E+00

     sig       sigi       sige       sigp       sigr     sigion
 9.381E-01  3.332E-01  5.785E-01  2.552E-05  2.590E-02  5.266E-04

The variables listed in this edit and their units are defined on the following pages. They are drwan from editcom variables and internal KEPLER variables.

First Row:

value	result
”`zone`”	Index of the zone being edited, $j$
”`masscoord`”	Mass interior to the inner boundary of this zone in scalem (p 273) units (usually solar masses).
”`cycle`”	Current KEPLER cycle, ncyc (q 4)
“`time`”	Current problem time, time (p 2) (`sec`)
”`dt`”	Current timestep, q:dt (`sec`)
”`temp`”	Zonal temperature, tn $(j)$ (`K`)
”`density`”	Zonal density, dn $(j)$ (`g/cc`)
”`ye`”	Electron abundance (`mol/g`) [ionye].
”`dye`”	Total change in $Y_\mathrm{e}$ during the current timestep (`mol/g`).
”`yq`”	Total abundance of ISE or NSE isotopes with $A \ge 24$ (`mol/g`), `ionyq`
”`dyq`”	Total change in $Y_\mathrm{q}$ during the current timestep (`mol/g`).
”`iternse`”	Total number of iterations it took for ISE / NSE to converge in this zone, or `-1` if convergence was not achieved.

Second Row:

value	result
”`yp`”	Abundance of protons in the ISE / NSE network (`mol/g`) `yp`
”`yn`”	Abundance of neutrons in the ISF/NSE network (`mol/g`), `YN`
”`ysi`”	Abundance of $^{28}\mathrm{Si}$ in the ISE / NSE network (`mol/g`), `ysi`
”`y56`”	Abundance of $^{56}\mathrm{Fe}$ in the ISE / NSE network (`mol/g`), `y56`
”`f0`”	Convergence error for the ISE / NSE mass conservation equation, `f0`
”`g0`”	Convergence error for the ISE / NSE charge conservation equation (`mol/g`), `g0`
”`h0`”	Convergence error for the ISE / NSE heavy ( $A \ge 24$ ) isotope abundance conservation equation (`mol/g`), `h0`
”`q0`”	Convergence error for the ISE / NSE iron peak ( $A > 46$ plus $^{46}\mathrm{Ti}$ ) isotope abundance conservation equation (`mol/g`), `q0`
”`qxc12`”	Photodisintegration-flow-deternlined mass fraction of $^{12}\mathrm{C}$ in the ISE calculation.
”`qxo16`”	Photodisintegration-flow-deternlined mass fraction of $^{16}\mathrm{O}$ in the ISE calculation.
”`qxne`”	Photodisintegration-flow-deternlined mass fraction of $^{20}\mathrm{Ne}$ in the ISE calculation.
”`betaq`”	The ratio of the actual current abundance of $^{28}\mathrm{Si}$ to its “NSE” value based on the current abundance of neutron and protons.

Third Row:

value	result
”`wrate`”	Total zonal weak (i.e., neutronization) rate (`mol/g/sec`).
”`rectot`”	Total zonal electron capture rate (`mol/g/sec`).
”`rpdtot`”	Total zonal positron decay rate (`mol/g/sec`).
”`redtot`”	Total zonal electron decay rate (`mol/g/sec`).
”`snuw`”	Total zonal specific reutrino energy loss rate due to neutronization-related processes (`erg/g/sec`).
”`snubps`”	Specific neutrino energy loss rate due to plasma-related processes (`erg/g/sec`).
”`eectot`”	Total zonal electron capture plus positron decay neutrino loss rate (`erg/g/sec`).
”`epdtot`”	Formerly, the total zonal positron decay neutrino loss rate (`erg/g/sec`). This quantity is now combined into eectot, and epdtot is set to `0`.
”`eedtot`”	Total zonal electron decay neutrino loss rate (`erg/g/sec`).
”`ebind0`”	Average zonal nucleon binding energy at the beginning of the timestep (`MeV/nucleon`), including contribution from energy in nuclear excited states if iexciteh (p 171) $>$ `0`
”`ebind`”	Average zonal nucleon binding energy at the end of the timestep (`MeV/nucleon`), including contribution from energy in excited states if iexciteh (p 171) $>$ `0`
”`eexcite`”	Average zonal energy in nuclear excited states (`MeV/nucleon`).

Second Last Row:

value	result
”`snucq`”	ISE nuclear energy generation rate (`erg/g/sec`).
”`snuclq`”	ISE nuclear energy generation rate minus plasma neutrino loss rate (`erg/g/sec`).
”`sdotq`”	Net ISE energy generation rate including all neutrino losses (`erg/g/sec`).
”`yf`”	Total abundance of ISE or NSE isotopes with $A > 46$ plus $^{46}\mathrm{Ti}$ (`mol/g`) `ionyf`.
”`tauqse`”	Timescale for the $\mathrm{Fe}$ and $\mathrm{Si}$ quasiequilibrium blocks to equilibrate (`sec`).
”`yfqse`”	Quasiequilibrium value of $Y_\mathrm{q}$ based on the current abundance of neutron and protons and $^{28}\mathrm{Si}$ (`mol/g`).
”`yffac`”	The actual value of $Y_\mathrm{q}$ is less than `yfqse` by the factor `1` $/$ `yffac`.
”`dyf`”	Total change in $Y_\mathrm{f}$ during the current timestep (`mol/g`).
”`qj`”	Net photodisintegration flow downward from $^{24}\mathrm{Mg}$ (`mol/g/sec`).
”`y016O`”	Abundance of unburned $^{16}\mathrm{O}$ in the ISE / NSE network (`mol/g`), `yo160`.
”`dyo16q`”	Change in the $^{16}\mathrm{O}$ abundance during this timestep (`mol/g`)
”`eta`”	Electron (negatron) degeneracy parameter $\mu_\mathrm{e} / k_\mathrm{B}T$ , etan $(j)$ .

Last Row:

value	result
”`sig`”	total entropy sig
“`sigi`”	the ion entropy, sigi
“`sige`”	the electron entropy, sige
“`sigp`”	the positron entropy, sigp
“`sigr`”	radiation entropy, sigr
“`sigion`”	the entropy due to ionization, sigion

Detailed ISE / NSE Isotope Edit

The middle section of the ISE / NSE edit gives a list of all the isotopes in the edited zone whose mass fractions exceed xthres (p 183), together with quantities specific to each isotope. The edits for neutrons and protons, however, are always given, no matter what their mass fraction. Each isotope is identified by its isotope number “k” and by its isotope symbol, “isotope”.

The specific isotopic quantities listed are defined as follows:

value	result
”`xk`”	Isotopic mass raction
”`wk`”	Isotopic statistical weight
”`rtot`”	Total weak interaction rate due to this isotope (`mol/g/sec`)
”`rec`”	Electron ( $\mathrm{e}^-$ ) capture rate due to this isotope (`mol/g/sec`)
”`rpd`”	Positron emission rate due to this isotope (`mol/g/sec`)
”`red`”	Beta ( $\mathrm{e}^-$ ) decay rate due to this isotope (`mol/g/sec`)
”`eec`”	Neutrino energy loss due to $\mathrm{e}^-$ capture and $\mathrm{e}^+$ decay on this isotope (`erg/g/sec`)
”`etot`”	Total neutrino energy loss due to this isotope (`erg/g/sec`)
”`eex`”	Energy stored in the excited states of this isotope (`MeV/nucleus`)

BURN IsotoDe Abundance Cycle Edit¶

  burn ion edit:


       zone       mass     nt1        h1         he4        h2         h3         he3        li6        li7        li8        be7
          1    0.00000  6.677E-15  1.881E-13  2.084E-11  2.640E-25  4.451E-23  3.767E-24  9.887E-31  5.262E-29  2.201E-39  2.304E-30
          2    0.00200  6.623E-15  1.868E-13  2.071E-11  2.600E-25  4.401E-23  3.709E-24  9.674E-31  5.165E-29  2.141E-39  2.258E-30
          3    0.00408  6.567E-15  1.855E-13  2.058E-11  2.561E-25  4.353E-23  3.654E-24  9.479E-31  5.075E-29  2.087E-39  2.215E-30
          4    0.00670  6.504E-15  1.839E-13  2.043E-11  2.520E-25  4.300E-23  3.596E-24  9.279E-31  4.981E-29  2.035E-39  2.171E-30
          5    0.01000  4.544E-15  1.371E-13  1.579E-11  1.214E-25  2.662E-23  1.761E-24  3.131E-31  2.053E-29  4.656E-40  8.094E-31
          6    0.01500  4.523E-15  1.365E-13  1.574E-11  1.204E-25  2.647E-23  1.748E-24  3.098E-31  2.035E-29  4.596E-40  8.015E-31
          7    0.02000  4.500E-15  1.359E-13  1.567E-11  1.194E-25  2.632E-23  1.734E-24  3.066E-31  2.017E-29  4.540E-40  7.938E-31
          8    0.02311  4.483E-15  1.354E-13  1.562E-11  1.186E-25  2.621E-23  1.724E-24  3.043E-31  2.004E-29  4.502E-40  7.882E-31
          9    0.02577  4.465E-15  1.350E-13  1.558E-11  1.179E-25  2.610E-23  1.715E-24  3.022E-31  1.991E-29  4.466E-40  7.830E-31
         10    0.02855  4.445E-15  1.344E-13  1.552E-11  1.171E-25  2.598E-23  1.704E-24  2.998E-31  1.977E-29  4.427E-40  7.771E-31
         11    0.03213  4.423E-15  1.338E-13  1.546E-11  1.163E-25  2.584E-23  1.693E-24  2.973E-31  1.962E-29  4.387E-40  7.710E-31
         12    0.03556  4.401E-15  1.332E-13  1.539E-11  1.155E-25  2.571E-23  1.681E-24  2.950E-31  1.948E-29  4.349E-40  7.651E-31
         13    0.03921  4.379E-15  1.327E-13  1.533E-11  1.147E-25  2.558E-23  1.671E-24  2.928E-31  1.934E-29  4.315E-40  7.596E-31

...

 conv  1180   11.95768  1.359E-28  5.933E-01  3.733E-01  5.520E-13  1.958E-26  1.240E-06  8.707E-22  6.009E-10  2.951E-20  5.374E-23
 conv  1181   11.96475  1.359E-28  5.933E-01  3.733E-01  5.520E-13  1.956E-26  1.240E-06  8.707E-22  6.008E-10  2.951E-20  4.508E-23
 conv  1182   11.97112  1.359E-28  5.933E-01  3.733E-01  5.520E-13  1.953E-26  1.240E-06  8.708E-22  6.007E-10  2.951E-20  3.564E-23
 conv  1183   11.97899  1.359E-28  5.933E-01  3.733E-01  5.520E-13  1.950E-26  1.240E-06  8.708E-22  6.005E-10  2.951E-20  2.733E-23
 conv  1184   11.98553  1.360E-28  5.933E-01  3.733E-01  5.521E-13  1.949E-26  1.240E-06  8.708E-22  6.005E-10  2.951E-20  2.335E-23
 conv  1185   11.98920  1.360E-28  5.933E-01  3.733E-01  5.521E-13  1.947E-26  1.240E-06  8.708E-22  6.004E-10  2.951E-20  1.979E-23
 conv  1186   11.99225  1.360E-28  5.933E-01  3.733E-01  5.521E-13  1.946E-26  1.240E-06  8.709E-22  6.003E-10  2.951E-20  1.678E-23
 osht  1187   11.99455  1.360E-28  5.933E-01  3.733E-01  5.521E-13  1.944E-26  1.240E-06  8.709E-22  6.002E-10  2.951E-20  1.382E-23
       1188   11.99614  1.360E-28  5.933E-01  3.733E-01  5.521E-13  1.940E-26  1.240E-06  8.709E-22  6.001E-10  2.951E-20  9.228E-24
       1189   11.99700  1.361E-28  5.933E-01  3.733E-01  5.523E-13  1.923E-26  1.240E-06  8.711E-22  5.993E-10  2.951E-20  2.849E-24
       1190   11.99789  3.077E-28  5.941E-01  3.733E-01  7.180E-13  1.928E-36  1.247E-06  1.057E-21  3.256E-10  2.958E-20  9.254E-87
       1191   11.99852  4.465E-28  5.941E-01  3.733E-01  8.176E-13  5.269E-59  1.249E-06  1.155E-21  2.252E-10  2.958E-20  0.000E+00
       1192   11.99881  4.777E-28  5.941E-01  3.733E-01  8.231E-13  1.713E-80  1.249E-06  1.156E-21  2.253E-10  2.958E-20  0.000E+00
       1193   11.99913  4.920E-28  5.941E-01  3.733E-01  8.249E-13  4.113E-91  1.249E-06  1.156E-21  2.254E-10  2.958E-20  0.000E+00
       1194   11.99938  5.014E-28  5.941E-01  3.733E-01  8.253E-13  1.430E-96  1.249E-06  1.156E-21  2.255E-10  2.958E-20  0.000E+00
       1195   11.99954  5.118E-28  5.941E-01  3.733E-01  8.255E-13  0.000E+00  1.249E-06  1.157E-21  2.255E-10  2.958E-20  0.000E+00
       1196   11.99968  5.256E-28  5.941E-01  3.733E-01  8.258E-13  0.000E+00  1.249E-06  1.157E-21  2.255E-10  2.958E-20  0.000E+00
       1197   11.99977  5.422E-28  5.941E-01  3.733E-01  8.260E-13  0.000E+00  1.249E-06  1.157E-21  2.255E-10  2.958E-20  0.000E+00
       1198   11.99984  5.665E-28  5.941E-01  3.733E-01  8.264E-13  0.000E+00  1.249E-06  1.157E-21  2.256E-10  2.958E-20  0.000E+00
       1199   11.99987  5.967E-28  5.941E-01  3.733E-01  8.269E-13  0.000E+00  1.249E-06  1.157E-21  2.256E-10  2.958E-20  0.000E+00
       1200   11.99990  6.324E-28  5.941E-01  3.733E-01  8.276E-13  0.000E+00  1.249E-06  1.157E-21  2.256E-10  2.958E-20  0.000E+00
       1201   11.99992  6.800E-28  5.941E-01  3.733E-01  8.283E-13  0.000E+00  1.249E-06  1.157E-21  2.256E-10  2.958E-20  0.000E+00
       1202   11.99994  7.450E-28  5.941E-01  3.733E-01  8.295E-13  0.000E+00  1.249E-06  1.157E-21  2.257E-10  2.958E-20  0.000E+00
       1203   11.99995  8.678E-28  5.941E-01  3.733E-01  8.311E-13  0.000E+00  1.249E-06  1.157E-21  2.257E-10  2.958E-20  0.000E+00
       1204   11.99996  1.107E-27  5.941E-01  3.733E-01  8.342E-13  0.000E+00  1.249E-06  1.158E-21  2.258E-10  2.958E-20  0.000E+00
       1205   11.99998  3.031E-27  5.949E-01  3.728E-01  1.995E-12  0.000E+00  1.278E-06  2.764E-21  3.867E-10  3.024E-20  0.000E+00
       1206   11.99998  1.064E-26  5.979E-01  3.704E-01  7.002E-12  0.000E+00  1.394E-06  9.700E-21  4.629E-10  3.297E-20  0.000E+00
       1207   11.99999  4.032E-26  6.062E-01  3.642E-01  2.654E-11  0.000E+00  1.694E-06  3.677E-20  4.842E-10  4.024E-20  0.000E+00
       1208   11.99999  4.903E-26  6.077E-01  3.630E-01  3.227E-11  0.000E+00  1.751E-06  4.471E-20  4.781E-10  4.163E-20  0.000E+00
       1209   12.00000  4.903E-26  6.077E-01  3.630E-01  3.227E-11  0.000E+00  1.751E-06  4.471E-20  4.781E-10  4.163E-20  0.000E+00
      total   12.00000  4.675E-17  5.114E-01  3.219E-01  4.700E-13  4.576E-23  1.068E-06  1.325E-21  5.732E-10  2.542E-20  3.137E-11

...

       zone       mass     bi204      bi205      bi206      bi207      bi208      bi209      bi210       1-tot
          1    0.00000  1.122E-74  1.457E-77  4.331E-82  1.432E-85  1.151E-90  3.152E-95  0.000E+00 -2.252E-10
          2    0.00200  1.123E-74  1.459E-77  4.336E-82  1.433E-85  1.153E-90  3.156E-95  0.000E+00 -2.252E-10

...

       1203   11.99995  0.000E+00  0.000E+00  0.000E+00  0.000E+00  8.088E-79  8.996E-18  0.000E+00 -7.350E-14
       1204   11.99996  0.000E+00  0.000E+00  0.000E+00  0.000E+00  8.088E-79  8.996E-18  0.000E+00 -7.505E-14
       1205   11.99998  0.000E+00  0.000E+00  0.000E+00  0.000E+00  8.042E-79  8.946E-18  0.000E+00 -2.753E-14
       1206   11.99998  0.000E+00  0.000E+00  0.000E+00  0.000E+00  7.853E-79  8.738E-18  0.000E+00 -1.776E-14
       1207   11.99999  0.000E+00  0.000E+00  0.000E+00  0.000E+00  7.351E-79  8.183E-18  0.000E+00 -1.288E-14
       1208   11.99999  0.000E+00  0.000E+00  0.000E+00  0.000E+00  7.257E-79  8.077E-18  0.000E+00 -1.155E-14
       1209   12.00000  0.000E+00  0.000E+00  0.000E+00  0.000E+00  7.255E-79  8.077E-18  0.000E+00 -1.199E-14
      total   12.00000  1.923E-28  5.799E-28  5.150E-30  7.884E-27  2.105E-27  7.749E-18  2.071E-25 -3.167E-11

The columns labeled ‘zone’ and ‘mass’ list the same values as in the primary zonal edit.

The following columns list the zonal mass fraction of the first ten isotopes in the current BURN network being employed in each zone, with each column labeled with the appropriate BURN isotope symbol.

After all the zones in the BURN network are listed, a final row labeled “total” in the ‘zone’ colum summarizes the total mass fraction of each isotope (as summed over all BURN zones) relative to the current total mass totm (q 17) of the all the zones in the star. The column labeled “mass” in this row gives the total mass in the star interior to the outer boundary of this block of BURN zones in scalem (p 273) units (usually solar masses).

When the edit described above for the first ten BURN isotopes is complete, an essentially similar edit table is begun containing the same “zone” and “mass” columns, together with the next ten BURN isotopes, and so on until all the ions have been listed. A final column labeled “1-tot” is appended that lists the amount by which the sum of the mass fractions of all isotopes in each zone differs from unity.

Note

The total mass fraction given in the “1-tot” column for the “total” row includes all mass not accounted for by this block of BURN zones, either due to mass-nonconservation in nuclear burning (generally due to zeroing small negative abundances) or to mass being present in the problem that is not being co-processed by the BURN coprocessor.

Changeable Parameter Edit¶

parameter     name    value      parameter     name    value      parameter     name    value      parameter     name    value

       1     dtnew  3.680E+01          51  thickfac  1.000E-03         101   sneutmt  1.000E+00         151     fmax1  2.000E-03
       2      time  2.136E+06          52    maxbak         20         102    snucmt  1.000E+00         152     fmax2  6.000E-03
       3    extrap  1.000E+00          53     dtcut  1.000E-01         103     dsnum  1.000E-06         153   ddsfrac  1.000E-02
       4     q2fac  4.000E+00          54    tfcrbu  1.000E+01         104  tcorefac  1.000E+00         154   dtsfrac  1.000E-03
       5     maxit         40          55    tfctbu  1.000E+01         105   tqsemin  2.500E+09         155   iudflag          0
       6      dtcr  2.000E-02          56     xipot  1.361E+01         106      jqse          0         156    nsdump        100
       7      dtct  2.000E-02          57     dzer0  1.000E-01         107     artv1  1.000E+00         157  niterbar          1
       8      dtcd  2.000E-02          58    npflag          1         108     artv2  1.000E+00         158  iterbarm     999999
       9      dtcq  5.000E-02          59  fracneut  5.000E-02         109     artv3  1.000E+00         159  iautoout          0
   dtcdt  9.900E-01          60   radius0  0.000E+00         110    jartv1       1000         160  iflgabar          0
  fcrmax  1.000E-08          61     summ0  0.000E+00         111    jartv2       1000         161  fraccore  1.100E+00
  fctmax  1.000E-08          62     xlum0  0.000E+00         112    jartv3       1000         162    jmcalc       1209
   q1fac  1.000E-01          63     irnet          0         113  ipixtype         31         163  iupdflag          0
   nstop     150000          64  npixedit          1         114  nplotsym          8         164  niondump         40
   tstop  1.000E+30          65   tnucmin  0.000E+00         115   tscalem  1.000E+07         165  nisodump        400
   nedit    1000000          66   setparm  1.673E+01         116    igridt        370         166  nzondump         40
  dtedit  0.000E+00          67  iytsflag          1         117    igridm        498         167   iflagye          1
   ndump         10          68    tbound  0.000E+00         118    igridv        496         168   njeditq         10
   xmlen  1.000E+00          69    pbound  1.000E+01         119       jp0          1         169    dtqnum  1.000E-04
   fudgc  1.000E-02          70    fclmax  1.000E+99         120       jp1          0         170    ddqnum  1.000E-04
   difim  1.000E+00          71    fclext  1.000E-05         121    rpmin0  1.000E+07         171  iexciteh          1
   fcrbu  1.000E-02          72      dtcl  1.000E+99         122    rpmax0  1.000E+14         172   iterqmh       1000
   fctbu  1.000E-02          73     fclbu  1.000E+99         123    ymintd  1.000E+00         173   ypconvh  1.000E-08
  drmult  1.000E-01          74    tfclbu  1.000E+99         124    ymaxtd  1.000E+08         174   ynconvh  1.000E-08
   dtmax  1.750E+11          75   dtsmult  1.000E+99         125   sscalem  1.000E-03         175  ysiconvh  1.000E-08
  fcrext  1.000E+01          76  rnratmax  2.500E-01         126   pscalem  1.000E-01         176     cnseh  1.000E+00
  fctext  1.000E+01          77  rnratmin  1.000E-01         127  itvstart          0         177      fyph  5.000E-02
    ipup          1          78  tnratmax  2.000E-01         128   abunlim  1.000E-03         178      fynh  5.000E-02
   xk1mt  1.000E+00          79  tnratmin  8.000E-02         129    ivplot          1         179     fysih  1.500E-01
   xk2mt  1.000E+00          80  dnratmax  2.000E-01         130   vscalem  1.000E-05         180    neditq         20
   xk3mt  1.000E+00          81  dnratmin  8.000E-02         131   rscalem  1.000E-05         181  xiter1qe  5.000E-01
   xk4mt  1.000E+00          82     rnmin  2.000E+08         132    irtype          1         182   iter1qe         10
  rxkcmt  1.000E+00          83     tnmin  1.000E+04         133   convord  3.000E+00         183    xthres  1.000E-04
  pimult  1.000E+00          84     dnmin  1.000E-04         134  yplotmin  0.000E+00         184   tqselim  1.500E+09
  prmult  1.000E+00          85     jmmin         10         135  yplotmax  8.333E-03         185    o16lim  4.000E-02
  pemult  1.000E+00          86    izonef          1         136    etacut  1.000E+01         186   qn56lim  0.000E+00
  transm  1.000E+00          87   idzonef          1         137    nadapb          1         187  snuwmult  1.000E+00
  timex0 -1.000E+99          88     rnmax  1.000E+14         138   fracrz1  2.500E-01         188      jnse          0
 dyemult  2.000E+02          89    dypmin  1.000E-10         139   fracrz2  3.333E-01         189  siqselim  1.000E-03
 dyqmult  2.500E+00          90   dymgmin  1.000E-10         140   rzmult0  1.000E+00         190    jlcalc          0
   yemin  2.000E-02          91    zbound  1.000E-01         141   rzmult1  1.000E+00         191   vlimset -1.000E+00
iwinsize   10240750          92   etaconv  1.000E-08         142   rzmult2  1.000E+00         192    jpause         -1
  hstatm  1.000E+00          93   jshell0          1         143  abarrat0  2.000E+00         193   fracrz0  2.000E-01
 lenqmax    6000000          94   jshell1     842836         144  abarratm  1.300E+00         194   rzmultm  1.000E+00
 iwinloc          0          95  eexplode  0.000E+00         145  fracmlim  1.000E-02         195     fmaxm  1.000E-03
    dtcp  1.500E-01          96  texplode  1.000E+99         146  frcsound  1.000E+99         196  eionmult  1.000E+00
 yfloorx  3.000E-03          97    tauexp  1.000E-02         147   convlim  1.000E+00         197  nnewoutf       5000
    cenu  1.000E+00          98    p0core  0.000E+00         148  woversht  1.000E-02         198   neditq1          5
  t7peek  1.000E+50          99    d0core  3.000E+11         149   xmimult  1.000E+00         199  wilsonmt -1.000E+00
   xkmin  1.000E-10         100   d0power  1.000E+00         150     fmax0  1.000E-03         200    t11cut  1.000E+00

parameter     name    value      parameter     name    value      parameter     name    value      parameter     name    value

y56gessm  1.000E-04         251     ncent     999999         301   idtmaxl         40         351    binbet  1.000E+00
isi30brn          1         252      nedt     999999         302   idtlook         10         352     bina0  2.000E+00
 dqselim  1.000E+05         253      nzro          1         303  backfacq  5.000E-01         353    binmdt  1.000E-03
abunminx -1.000E-06         254   amaglim  0.000E+00         304  tempstop  1.000E+99         354    rocher  1.000E+99
  tfcybu  3.000E-01         255      ninv   ********         305   denstop  1.000E+99         355  xlanger1  0.000E+00
yfloorbx  3.000E-03         256     nbkup    5275740         306   vinstop  9.000E+07         356  xlanger2  0.000E+00
iqerrflg          1         257      ntty          6         307   o16stop -1.000E+00         357   iburnye          0
c12agmlt  1.200E+00         258     ileqs          1         308   timezms  1.000E+99         358   relmult  0.000E+00
xltaucon  0.000E+00         259   negbkup     157413         309  izonezms          1         359   geemult  1.000E+00
noqsecon          0         260   bkupdiv  1.000E+01         310  q1faczms  1.000E-01         360   grbparm  0.000E+00
 accrate  0.000E+00         261  bkupmass  1.000E-13         311   tempcig  1.000E+99         361    swmult  0.000E+00
xmacrete  0.000E+00         262    bkupmp  1.000E-19         312  yflrxcig  3.000E-03         362    tsharp  0.000E+00
 denconv  1.000E-07         263  nbkupmax         10         313  fmaxmcig  1.000E-03         363   xmlossn  0.000E+00
flamerad  0.000E+00         264  nsubcycm          1         314  fmax0cig  1.000E-03         364   nangmix          1
coulmult  1.000E+00         265    mazful          1         315   toffset  4.460E+14         365    angfmu  5.000E-02
  irzopt          3         266  al26mult  1.000E+00         316  abunminb  1.000E-04         366     angfc  3.330E-02
  nsurfz        -10         267    inburn          1         317  abunmaxb  1.000E+00         367    angfjc  1.000E+00
 fracsz0  1.000E+00         268     npage         53         318    numiso          0         368  angrcrit  2.500E+03
 fracsz1  1.000E+00         269    nsavez         33         319   timeref -1.000E+99         369    angric  2.500E-01
 xmlossm  0.000E+00         270    nsaveb         10         320  tosetref  0.000E+00         370  angfjdsi  1.000E+00
 xmloss0  0.000E+00         271     vloss  1.000E+99         321  timecmin  0.000E+00         371  angfjshi  1.000E+00
   totm0  2.387E+34         272  abunlimb  1.000E-05         322  timecmax  0.000E+00         372  angfjssi  9.000E-01
 fracdez  5.000E-01         273    scalem  1.989E+33         323     yemax  4.980E-01         373   angfjez  9.000E-01
xmratbak  2.000E-01         274   neditz1         10         324  abarsemi  4.000E+00         374  angfjgsf  9.000E-01
 xfracml  1.000E-02         275   neditz2         20         325  drmultlo  1.000E-01         375      fmin  1.000E-01
 mlcompf          0         276     medit          7         326  woverslo  1.000E-02         376   ncnvout          1
e1mltc12  1.700E+00         277  meditfin          7         327   maptime          2         377    kaptab          1
e2mltc12  1.700E+00         278    neditp         50         328   vminmap  1.000E+99         378     fkapz  1.000E+00
 nupdate        100         279  c12flame  0.000E+00         329   vmaxmap -1.000E+99         379   zfakexp  5.000E-01
  dtfrac  1.000E-02         280  o16flame  0.000E+00         330   vratmap  1.000E-99         380    angsmt  2.000E-01
bmassmin -1.000E+00         281   timenew  0.000E+00         331  tempcdep  1.000E+99         381   nangsmg          2
bmassmax  1.000E+99         282    sharp1  0.000E+00         332   o16odep -1.000E+99         382    angsml  1.000E-03
btempmin  1.000E-03         283    sharpr  1.000E+00         333  tempchar  1.000E+99         383    angsmm  1.000E-03
btempmax  1.000E+99         284    sharp2  0.000E+00         334   denchar  1.000E+99         384  ipapsize    6120792
 snucmin -1.000E+99         285   newflam        259         335  abarchar  1.000E+99         385   hstatxm  1.000E+99
 bdenmin  1.000E-99         286     taunu  0.000E+00         336  zonemmin  2.000E+30         386   hstatym  1.000E+32
 bdenmax  1.000E+99         287     enu53  3.000E+00         337   icalcne          0         387   xmlossw  0.000E+00
 tchange  2.000E-02         288     tmunu  8.000E+00         338   xneconv  1.000E-05         388      iold          0
 dchange  5.000E-02         289      tenu  4.000E+00         339   ionflag          0         389  rhotrans  1.000E+07
  netmax          2         290   noiland          1         340    xnemin  1.000E-05         390   nwndout          1
  neditb          5         291  charsizg  1.000E+00         341   xkapgam  5.400E-02         391   kapverb          0
 edmassl  1.000E-04         292  charsizc  8.500E-01         342     egamp  6.000E-01         392   xl0limf  4.000E+00
  nedita         10         293  charsizz  8.500E-01         343    tshock  1.000E+99         393   xl0limk  0.000E+00
 jmeditb          1         294  charsizh  1.000E+00         344   tnucleo  1.000E+99         394   llimout          1
neditall        100         295   widthtd  6.000E+01         345    tenvel  1.000E+99         395    nenout          0
  chimin  1.000E-10         296  ibackgnd         -1         346   nfirstq          0         396   ipromin          1
  delchi  1.500E-01         297   dscalem  1.989E+33         347  xlmxmult  1.000E+00         397   ipromax      99999
    fdtn  2.000E+00         298  ngridmin         10         348    binm10  2.387E+34         398   iprownd          1
  dtbkup  5.000E-01         299    ncycqq    1000000         349    binm20  0.000E+00         399   iproyld          0
   ncomp     999999         300  lentrack      19840         350    binalp  1.000E+00         400   minapro      -1000

parameter     name    value      parameter     name    value

 maxapro      -1000         451     ddmin -1.000E+00
 proymin  1.000E-10         452   iazonef          1
 proymax  1.000E+01         453    dynfac  1.000E-03
  xkdmin  1.000E-20         454   h1hburn -1.000E+00
  h1hdep -1.000E+00         455  c12heign  1.000E+99
he4hedep -1.000E+00         456  he4hebrn -1.000E+00
 nenuout         10         457  zonermin -1.000E+00
 optconv  0.000E+00         458  zonermax  1.000E+99
   rloss  1.000E+99         459  xmixnova  0.000E+00
   tloss -1.000E+00         460   accmass  0.000E+00
 tapprox -1.000E+99         461  accmassf  1.414E+00
 semilan  0.000E+00         462  vinstopm  1.000E+99
 profmin  1.000E-04         463   lowamul          0
 profmax  1.000E+03         464    pdmult  1.000E+00
 proamin  1.000E-12         465    edmult  1.000E+00
 proamax  1.000E+00         466   ihwamul       9999
  yelimb  4.970E-01         467   kapburn          0
   irecb          1         468    fackap  1.000E+00
bmasslow  1.979E+33         469  awwkloss  1.000E+99
 si28dep -1.000E+99         470   lossrot          0
 ivspecl          6         471    ymjkep -1.000E+99
  ivrate          1         472   maxzone          0
  magnet          2         473   cfakexp  0.000E+00
   nosht          2         474   minnucg         -1
   alpth  1.000E+00         475  mingaind         21
abarstep  1.000E+99         476  minlossd         21
zbarstep  1.000E+99         477  minnucgd         21
 xmustep  1.000E+99         478  minnucld         21
  netmin          1         479  tweakmin  5.000E-01
 awcorot -1.000E+00         480  centmult  0.000E+00
 ymcorot -1.000E+00         481    mixout          0
 nstrout          0         482    irprox          1
 mixcycl          2         483     n14pg          0
   lburn          0         484   r3amult  1.000E+00
  lbbkup          0         485    ibwarn          0
rlossmin -1.000E+00         486   ifallbk          0
   lcout          0         487   xnumu12  0.000E+00
 xmagfmu  1.000E+00         488    nzsave          1
  xmagft  1.000E+00         489     axion  0.000E+00
 xmagfnu  1.000E+00         490  zmhiconv  1.000E+99
xmagfdif  1.000E+00         491  rnhiconv  1.000E+99
dxncbkup  1.000E-13         492   icutbin          1
  iplotb          0         493  nconvers      10400
 minzone          1         494  nwndvers      20100
zonemmax  1.000E+99         495    h1hign  1.000E+99
 tenubar -1.000E+00         496      wimp  0.000E+00
  levcnv          1         497   ipdtmin          0
 mingain         -1
 minloss         -1
 minnucl         -1

Timeing information¶

>time
  timing information (seconds):

  activity    cpu        i/o       system     total
  iterate  1.322E+04  0.000E+00  1.782E+03  1.500E+04
  convect  8.282E+01  0.000E+00  4.000E-01  8.322E+01
  opacity  7.583E+03  0.000E+00  1.697E+03  9.280E+03
  nuclear  9.217E+04  0.000E+00  1.811E+03  9.399E+04
      eos  1.181E+04  0.000E+00  1.703E+03  1.351E+04
   rezone  8.170E+01  0.000E+00  1.300E-01  8.183E+01
     edit  1.242E+02  0.000E+00  4.379E+01  1.680E+02
    other  6.627E+03  0.000E+00  1.278E+02  6.755E+03
 convectb  7.755E+03  0.000E+00  1.040E+01  7.765E+03
   ratesb  1.793E+04  0.000E+00  2.066E+01  1.795E+04
  invertb  7.407E+05  0.000E+00  2.926E+02  7.410E+05
    miscb  2.901E+03  0.000E+00  1.541E+01  2.917E+03
 qq dumps  2.000E-02  0.000E+00  8.000E-02  1.000E-01
 rotation  2.046E+01  0.000E+00  2.900E-01  2.075E+01
    total  9.010E+05  0.000E+00  7.504E+03  9.085E+05

The information is identical to that returned by the time command.

Nuclear Reaction Edit¶

>tn
 reaction  total,moles  rate,moles/s  maxzone  mcoord,msun
 p+p        4.2728E+31    2.1411E+09       655    3.007
 he3+he3    1.2499E+30    2.3189E+01       655    3.007
 he3+he4    1.1924E+32    4.3323E+06       651    2.973
 c12(pg)    1.1178E+33    6.0616E+07       645    2.925
 n14(pg)    1.2425E+33    2.1904E+07       645    2.925
 o16(pg)    1.0152E+30    3.4449E+03       645    2.925
 3a-c12     3.1998E+32   -9.7031E+27       470    2.086
 n14(ag)    1.0932E+28    1.7131E+10       508    2.263
 c12(ag)    2.0685E+32   -3.4637E+29       356    1.551
 c12+c12    3.8854E+31    3.3053E+28       325    1.415
 c12+o16    3.8998E+29    2.5744E+29       325    1.415
 o16+o16    7.7836E+31    1.5220E+29       325    1.415
 o16(ag)   -8.0862E+32    4.6808E+28       326    1.420
 ne20(ag)   2.2400E+32    2.9440E+33       222   0.9531
 mg24(ag)   9.8828E+30    1.9631E+29       325    1.415
 si28(ag)  -3.3697E+30   -1.9381E+30       327    1.425
 s32(ag)   -3.6397E+30   -9.2583E+29       327    1.425
 ar36(ag)   2.5832E+28   -2.7917E+29       327    1.425
 ca40(ag)  -7.6219E+27   -5.2968E+27       327    1.425
 ti44(ag)   7.2435E+25    5.3662E+26       325    1.415
 cr48(ag)   4.3981E+25    4.4452E+26       325    1.415
 fe52(ag)   3.6627E+21    6.9082E+22       325    1.415
 fe52:2ng  -2.8732E+22   -6.0565E+23       326    1.420
 fe52:a2p   1.1556E+25    1.8346E+26       325    1.415
 fe54:2pg  -2.8673E+22   -1.0553E+24       329    1.434
 he4-2n2p   1.5740E+07    2.7317E+08       325    1.415
 p(e-nu)n   4.0351E+29    9.7311E+29       222   0.9531
 n(e+nu)p   3.1399E+25    2.8948E+25       222   0.9531
 weak       9.4775E+31    9.4940E+30         1   2.0000E-03

For each reaction (first column) give the total number of mols that have reacted since start of the problem (second column), the the current rate (third column), the zone where it is fastest (fourth column) and its (outer) mass coolrdionate (last column).

Core Information¶

>core
 core            j      zm         rn      ybind       stot      angjt
       center    1   0.000  0.000D+00 -5.372D+50  5.863D-01  0.000D+00
      ye core  294   1.275  1.036D+08 -3.839D+50  2.123D+00  0.000D+00
  approx netw  325   1.411  1.683D+08 -1.816D+50  3.707D+00  0.000D+00
      O shell  325   1.411  1.683D+08 -1.816D+50  3.707D+00  0.000D+00
    iron core  306   1.326  1.220D+08 -3.074D+50  2.589D+00  0.000D+00
      Si core  356   1.547  4.056D+08 -1.233D+50  4.909D+00  0.000D+00
 Ne/Mg/O core  434   1.912  1.646D+09 -4.118D+49  5.222D+00  0.000D+00
     C/O core  448   1.978  2.163D+09 -3.179D+49  6.937D+00  0.000D+00
      He core  651   2.965  4.038D+10 -1.632D+48  1.267D+01  0.000D+00
         star 1114  11.000  2.267D+13  0.000D+00  2.612D+01  0.000D+00

BURN Network Information¶

>burnamax
 iso    max massfrac      at mass     cycle
nt1      5.86249D-08        1.27944   38731    0    1
h1       7.59882D-01        0.00200     311    1    1
he4      9.99997D-01        0.00000    1218    2    4
h2       9.09974D-05       10.99814      59    1    2
h3       1.42741D-06        0.00970     104    1    3
he3      1.67836D-04       10.99968     766    2    3

...

pb207    0.00000D+00        0.00000      -1   82  207
pb208    0.00000D+00        0.00000      -1   82  208
bi209    0.00000D+00        0.00000      -1   83  209
  839 isotopes

BURN Network Matrix Diagnosis¶

>flowb 325 - - - + matrix
   marked ( h  3, c 12) but not used
   marked (he  3, c 12) but not used
   marked (li  9,be  9) but not used

 ...

   marked ( h  1, c 10) but not used
   marked ( h  1, n 11) but not used
   marked (he  4, n 11) but not used
   298116         matrix elements
   288631 zero    matrix elements
     9485 nonzero matrix elements
     9528 marked  matrix elements
       43 unused  matrix elements
        0 missing matrix elements

BURN Flow Analysis¶

>flowb 325 - - - +

 FLOWS for user-defined T, rho -- VERSION 10100


 cycle number =   38731
  complete flows for zone 325  last step = 6.503E-03   reference time = 0.00000E+00   update time = 3.4812559702403E+07

  t9 = 3.548995                rho = 3.632502E+06      eb =  3.148E+17,   3.028E+17   neg. y back-ups = 0
  nucleus changes fastest is     202 z=  21 n=  27     dth = 2.8047E-03               total inversions =       23
  mass non-conservation =   8.92069E-11                eta = 6.7657575E-04            total backups =        8
  plasma nu loss = 5.76E+13    beta nu loss = 4.27E+11 ecap+epos = 4.274E+11          e- decay = 4.675E+06
  weak rate = 8.647E-08        ecap = 2.768E-08        beta = 1.451E-12               rn= 1.718E+08



 MASS FRACTIONS

      p 4.698E-06      n 1.236E-11      a 2.210E-05   h  2 1.955E-17   h  3 2.945E-20  he  3 9.278E-16  li  6 1.598E-20
  li  7 1.216E-24  li  8 2.488E-34  li  9 5.344E-44  be  7 4.744E-20  be  9 1.862E-21  be 10 1.395E-25  be 11 2.692E-39
   b  8 2.324E-29   b 10 3.046E-18   b 11 7.998E-17   b 12 1.410E-27   b 13 6.637E-36   b 14 8.080E-50   c 10 2.160E-28
   c 11 6.215E-13   c 12 1.203E-04   c 13 3.069E-12   c 14 2.370E-14   c 15 1.262E-27   n 11 0.000E+00   n 12 5.456E-22

...

  zr 91 1.000E+00  zr 92 1.000E+00  zr 93 1.390E+00  zr 94 1.000E+00  zr 95 1.410E+00  zr 96 1.000E+00  nb 90 3.766E+00
  nb 91 1.300E+00  nb 92 2.369E+00  nb 93 1.000E+00  nb 94 5.510E+00  nb 95 1.503E+00  mo 92 1.000E+00  mo 93 1.000E+00
  mo 94 1.338E+00  mo 95 1.813E+00  mo 96 1.524E+00  mo 97 1.952E+00  mo 98 1.613E+00  mo 99 1.201E+01  mo100 2.459E+00

 ========= FLOWS =========


     3a     c12photo   c12+c12   c12+o16   o16+o16   na23+p    ne20+a    mg23+n    p31+p     si28+a    s31+n
  7.16E-14  7.63E-04  2.98E-03  1.23E-02  2.84E-03  6.51E-11  6.00E-11  1.43E-10  1.76E-18  6.89E-19  3.94E-19
  1.06E+00            1.19E+00  1.26E+00  1.36E+00  3.18E-01  6.36E-01  4.61E-02  6.00E-01  2.50E-01  1.00E-01

                                                     mg24+a    al27+p    si27+n              p(e,nu)n  n(e+,nub)p
                                                    3.83E-15  5.08E-15  5.03E-14             1.14E-09  1.90E-14
                                                    4.00E-01  5.00E-01  1.00E-01             1.31E+09  4.77E+04

  z/a             ng                pn                pg                ap                an                ag       weak     b

  h  2-9.50E-17 3.03E-20 0.00E+00 0.00E+00-8.25E-10 2.59E-13 0.00E+00 0.00E+00 0.00E+00 0.00E+00-5.64E-12 5.10E-16 0.0E+00 6.08E-17-8.31E-10
  h  3 0.00E+00 0.00E+00-1.17E-13 1.04E-11 2.66E-15 2.66E-15 0.00E+00 0.00E+00-3.07E-17 1.25E-17-9.56E-15 4.96E-16 1.7E-29 4.03E-19-1.24E-13
 he  3 5.51E-20 5.51E-20 0.00E+00 0.00E+00 0.00E+00 0.00E+00-1.70E-12 7.17E-13 0.00E+00 0.00E+00-8.23E-10 1.87E-12 0.0E+00 1.35E-14 1.87E-13
 li  6-2.92E-21 6.08E-22 0.00E+00 0.00E+00-7.42E-15 5.70E-17-4.08E-13 3.73E-15 0.00E+00 0.00E+00-7.09E-12 2.61E-17 0.0E+00 3.53E-19-1.57E-13

...

 mo 98 4.42E-43 4.42E-43 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.0E+00-2.51E-41-3.64E-41
 mo 99 6.15E-45 6.15E-45 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.0E+00-2.41E-43-4.36E-43
 mo100 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.0E+00-3.53E-43-6.15E-45
 total-9.96E-02 8.43E-01-8.02E-02 8.59E-01-1.33E+00 3.56E+02-1.99E+00 3.89E+02-9.38E-02 7.83E-01 1.53E+00 2.55E+00-7.9E-08-8.80E-08

  z/a             ng                pn                pg                ap                an                ag       weak     b


 ========= RATES =========


  z/a             ng                pn               beta               pg                ap                an                ag

  h  2 2.51E+08 9.68E+03 0.00E+00 0.00E+00 0.00E+00 0.00E+00 5.63E+09 2.67E+06 0.00E+00 0.00E+00 0.00E+00 0.00E+00 9.36E+06 2.12E+09
  h  3 0.00E+00 0.00E+00 2.25E+14 2.73E+15 1.78E-09 0.00E+00 5.79E+10 2.03E-13 0.00E+00 0.00E+00 2.28E+08 1.32E+15 9.06E+09 5.79E+10
 he  3 1.44E+07 4.08E-18 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 4.11E+08 1.96E+14 0.00E+00 0.00E+00 1.07E+09 1.22E+11
 li  6 1.85E+10 2.03E+04 0.00E+00 0.00E+00 0.00E+00 0.00E+00 4.60E+09 1.10E+06 2.46E+11 4.16E+14 0.00E+00 0.00E+00 1.72E+09 2.33E+07

...

 mo 98 6.09E+13 5.63E+08 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00
 mo 99 8.85E+13 5.23E+07 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00
 mo100 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00

BURN Statistics and Timing Information¶

::

>bstat: LAST ninvl = 288 nburnz = 35 TOTAL ninv = 69948603 nburnzt= 6417090 TIME ratesb= 4.151D+03 invertb= 1.952D+05 NETW imaxb = 546 nbmax = 839

Pulsar Data Estimate¶

> pulsar 2.
--------- PULSAR DATA ---------
 assuming pulsar parameters
 R         = 1.20000D+06 cm
 M_grav    = 1.40000D+00 M_sun
 M_baryonic= 2.00000D+00 M_sun
 I_PSR     = 3.60000D-01 M R^2

 J_PSR     = 6.85451D+41 erg s
 j_PSR     = 2.46133D+08 cm^2/s
 omega_PSR = 4.74794D-04 rad/s
 P_PSR     = 1.32335D+04 s
 -------------------------------

BURN and Themodynmaic Info Edit¶

>enucb 325 - - -
              temperature                  density                      eps                  deps/dT                deps/drho                    kappa                dkappa/dT              dkappa/drho
   3.5489954626487255E+09   3.6325018194859717E+06   0.0000000000000000E+00   0.0000000000000000E+00   0.0000000000000000E+00   4.5222822983423339E-02   2.1365107640083146E-11  -5.5622313963354002E-09