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 numberk” 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
      10     dtcdt  9.900E-01          60   radius0  0.000E+00         110    jartv1       1000         160  iflgabar          0
      11    fcrmax  1.000E-08          61     summ0  0.000E+00         111    jartv2       1000         161  fraccore  1.100E+00
      12    fctmax  1.000E-08          62     xlum0  0.000E+00         112    jartv3       1000         162    jmcalc       1209
      13     q1fac  1.000E-01          63     irnet          0         113  ipixtype         31         163  iupdflag          0
      14     nstop     150000          64  npixedit          1         114  nplotsym          8         164  niondump         40
      15     tstop  1.000E+30          65   tnucmin  0.000E+00         115   tscalem  1.000E+07         165  nisodump        400
      16     nedit    1000000          66   setparm  1.673E+01         116    igridt        370         166  nzondump         40
      17    dtedit  0.000E+00          67  iytsflag          1         117    igridm        498         167   iflagye          1
      18     ndump         10          68    tbound  0.000E+00         118    igridv        496         168   njeditq         10
      19     xmlen  1.000E+00          69    pbound  1.000E+01         119       jp0          1         169    dtqnum  1.000E-04
      20     fudgc  1.000E-02          70    fclmax  1.000E+99         120       jp1          0         170    ddqnum  1.000E-04
      21     difim  1.000E+00          71    fclext  1.000E-05         121    rpmin0  1.000E+07         171  iexciteh          1
      22     fcrbu  1.000E-02          72      dtcl  1.000E+99         122    rpmax0  1.000E+14         172   iterqmh       1000
      23     fctbu  1.000E-02          73     fclbu  1.000E+99         123    ymintd  1.000E+00         173   ypconvh  1.000E-08
      24    drmult  1.000E-01          74    tfclbu  1.000E+99         124    ymaxtd  1.000E+08         174   ynconvh  1.000E-08
      25     dtmax  1.750E+11          75   dtsmult  1.000E+99         125   sscalem  1.000E-03         175  ysiconvh  1.000E-08
      26    fcrext  1.000E+01          76  rnratmax  2.500E-01         126   pscalem  1.000E-01         176     cnseh  1.000E+00
      27    fctext  1.000E+01          77  rnratmin  1.000E-01         127  itvstart          0         177      fyph  5.000E-02
      28      ipup          1          78  tnratmax  2.000E-01         128   abunlim  1.000E-03         178      fynh  5.000E-02
      29     xk1mt  1.000E+00          79  tnratmin  8.000E-02         129    ivplot          1         179     fysih  1.500E-01
      30     xk2mt  1.000E+00          80  dnratmax  2.000E-01         130   vscalem  1.000E-05         180    neditq         20
      31     xk3mt  1.000E+00          81  dnratmin  8.000E-02         131   rscalem  1.000E-05         181  xiter1qe  5.000E-01
      32     xk4mt  1.000E+00          82     rnmin  2.000E+08         132    irtype          1         182   iter1qe         10
      33    rxkcmt  1.000E+00          83     tnmin  1.000E+04         133   convord  3.000E+00         183    xthres  1.000E-04
      34    pimult  1.000E+00          84     dnmin  1.000E-04         134  yplotmin  0.000E+00         184   tqselim  1.500E+09
      35    prmult  1.000E+00          85     jmmin         10         135  yplotmax  8.333E-03         185    o16lim  4.000E-02
      36    pemult  1.000E+00          86    izonef          1         136    etacut  1.000E+01         186   qn56lim  0.000E+00
      37    transm  1.000E+00          87   idzonef          1         137    nadapb          1         187  snuwmult  1.000E+00
      38    timex0 -1.000E+99          88     rnmax  1.000E+14         138   fracrz1  2.500E-01         188      jnse          0
      39   dyemult  2.000E+02          89    dypmin  1.000E-10         139   fracrz2  3.333E-01         189  siqselim  1.000E-03
      40   dyqmult  2.500E+00          90   dymgmin  1.000E-10         140   rzmult0  1.000E+00         190    jlcalc          0
      41     yemin  2.000E-02          91    zbound  1.000E-01         141   rzmult1  1.000E+00         191   vlimset -1.000E+00
      42  iwinsize   10240750          92   etaconv  1.000E-08         142   rzmult2  1.000E+00         192    jpause         -1
      43    hstatm  1.000E+00          93   jshell0          1         143  abarrat0  2.000E+00         193   fracrz0  2.000E-01
      44   lenqmax    6000000          94   jshell1     842836         144  abarratm  1.300E+00         194   rzmultm  1.000E+00
      45   iwinloc          0          95  eexplode  0.000E+00         145  fracmlim  1.000E-02         195     fmaxm  1.000E-03
      46      dtcp  1.500E-01          96  texplode  1.000E+99         146  frcsound  1.000E+99         196  eionmult  1.000E+00
      47   yfloorx  3.000E-03          97    tauexp  1.000E-02         147   convlim  1.000E+00         197  nnewoutf       5000
      48      cenu  1.000E+00          98    p0core  0.000E+00         148  woversht  1.000E-02         198   neditq1          5
      49    t7peek  1.000E+50          99    d0core  3.000E+11         149   xmimult  1.000E+00         199  wilsonmt -1.000E+00
      50     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

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

parameter     name    value      parameter     name    value

     401   maxapro      -1000         451     ddmin -1.000E+00
     402   proymin  1.000E-10         452   iazonef          1
     403   proymax  1.000E+01         453    dynfac  1.000E-03
     404    xkdmin  1.000E-20         454   h1hburn -1.000E+00
     405    h1hdep -1.000E+00         455  c12heign  1.000E+99
     406  he4hedep -1.000E+00         456  he4hebrn -1.000E+00
     407   nenuout         10         457  zonermin -1.000E+00
     408   optconv  0.000E+00         458  zonermax  1.000E+99
     409     rloss  1.000E+99         459  xmixnova  0.000E+00
     410     tloss -1.000E+00         460   accmass  0.000E+00
     411   tapprox -1.000E+99         461  accmassf  1.414E+00
     412   semilan  0.000E+00         462  vinstopm  1.000E+99
     413   profmin  1.000E-04         463   lowamul          0
     414   profmax  1.000E+03         464    pdmult  1.000E+00
     415   proamin  1.000E-12         465    edmult  1.000E+00
     416   proamax  1.000E+00         466   ihwamul       9999
     417    yelimb  4.970E-01         467   kapburn          0
     418     irecb          1         468    fackap  1.000E+00
     419  bmasslow  1.979E+33         469  awwkloss  1.000E+99
     420   si28dep -1.000E+99         470   lossrot          0
     421   ivspecl          6         471    ymjkep -1.000E+99
     422    ivrate          1         472   maxzone          0
     423    magnet          2         473   cfakexp  0.000E+00
     424     nosht          2         474   minnucg         -1
     425     alpth  1.000E+00         475  mingaind         21
     426  abarstep  1.000E+99         476  minlossd         21
     427  zbarstep  1.000E+99         477  minnucgd         21
     428   xmustep  1.000E+99         478  minnucld         21
     429    netmin          1         479  tweakmin  5.000E-01
     430   awcorot -1.000E+00         480  centmult  0.000E+00
     431   ymcorot -1.000E+00         481    mixout          0
     432   nstrout          0         482    irprox          1
     433   mixcycl          2         483     n14pg          0
     434     lburn          0         484   r3amult  1.000E+00
     435    lbbkup          0         485    ibwarn          0
     436  rlossmin -1.000E+00         486   ifallbk          0
     437     lcout          0         487   xnumu12  0.000E+00
     438   xmagfmu  1.000E+00         488    nzsave          1
     439    xmagft  1.000E+00         489     axion  0.000E+00
     440   xmagfnu  1.000E+00         490  zmhiconv  1.000E+99
     441  xmagfdif  1.000E+00         491  rnhiconv  1.000E+99
     442  dxncbkup  1.000E-13         492   icutbin          1
     443    iplotb          0         493  nconvers      10400
     444   minzone          1         494  nwndvers      20100
     445  zonemmax  1.000E+99         495    h1hign  1.000E+99
     446   tenubar -1.000E+00         496      wimp  0.000E+00
     447    levcnv          1         497   ipdtmin          0
     448   mingain         -1
     449   minloss         -1
     450   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