Edit Quantities by Group

KEPLER Saved Zonal Interface Variables

Note

These variables are defined as staring from Zone 0 to allow specifying the lower boundary of the innermost zone.

ym

Exterior mass coordinate.

name

type

unit

centering

storage

ym

float

g

boundary

dump

rn

Radius of outer zone boundary.

name

type

unit

centering

storage

rn

float

cm

boundary

dump

rd

Calculated change in radius (during a cycle) and radius at last timestep (at the end of a cycle).

name

type

unit

centering

storage

rd

float

cm

boundary

dump

un

Velocity of the outer zone boundary.

name

type

unit

centering

storage

un

float

cm/sec

boundary

dump

xln

Calculated non-neutrino luminosity at the outer zone boundary.

name

type

unit

centering

storage

xln

float

erg/sec

boundary

dump

qln

Converged non-neutrino luminosity at the outer zone boundary.

name

type

unit

centering

storage

qln

float

erg/sec

boundary

dump

qld

Calculated change in converged luminosity (during a cycle) and converged luminosity at last timestep (at the end of a cycle).

name

type

unit

centering

storage

qld

float

erg/sec

boundary

dump

difi

Diffusion coefficient for material mixing due to convection.

name

type

unit

centering

storage

difi

float

cm**2/sec

boundary

dump

angdg

Total rotationally induced diffusion coefficients.

name

type

unit

centering

storage

angdg

float

cm**2/sec

boundary

dump

Added In Version 15.2.

Saved Array of Rotational mixing coefficients

angd

Rotational diffusion coefficients separated for processes.

name

type

unit

centering

storage

angd

float

cm**2/sec

boundary

dump

2D-array, with second dimension 1nangmd with the following layout

value

result

1

dynamical shear instability

2

Solberg-H{o}iland instability

3

secular shear instability

4

Eddington-Sweet circulation

5

Goldreich-Schubert-Fricke instability

These can be accessed from loadbuf using angd1 $-$ angd5

Added In Version 15.2.

bfdiff

Diffision coefficient due to magnetic fields, $D_B$ .

name

type

unit

centering

storage

bfdiff

float

cm**2/sec

boundary

dump

Added In Version 16.11.

Other entry:
Magnetic Quantities

bfvisc

Viscosity due to magnetic fields, $D_B$ .

name

type

unit

centering

storage

bfvisc

float

cm**2/sec

boundary

dump

Added In Version 16.11.

Other entry:
Magnetic Quantities

Saved Arrays Evaluated at Zone Centers

Note

These are alos defined intenally to start a 0, but the 0 value usually is not used.

netnum

Number of the nuclear network currently used in the zone.

name

type

unit

centering

storage

netnum

integer

-

center

loadbuf

value

result

1

APPROX

2

ISE

3

NSE

xm

Zonal mass.

name

type

unit

centering

storage

xm

float

g

center

dump

dn

Density.

name

type

unit

centering

storage

dn

float

g/cc

center

dump

tn

Temperature.

name

type

unit

centering

storage

tn

float

K

center

dump

td

Calculated change in temperature (during a cycle) and temperature at last timestep (at the end of a cycle).

name

type

unit

centering

storage

td

float

K

center

dump

en

Specific energy density.

name

type

unit

centering

storage

en

float

erg/g

center

loadbuf

pn

Pressure.

name

type

unit

centering

storage

pn

float

erg/cc

center

dump

zn

Viscous stress, including linear and quadratic artificial viscosity.

name

type

unit

centering

storage

zn

float

erg

center

dump

etan

Electron degeneracy parameter ( $\mu_\mathrm{e} / k_\mathrm{B} T$ ).

name

type

unit

centering

storage

etan

float

-

center

dump

sn

Total energy generation rate – including neutrino losses.

name

type

unit

centering

storage

sn

float

erg/g/sec

center

dump

snn

Nuclear energy generation rate excluding neutrino losses.

name

type

unit

centering

storage

snn

float

erg/g/sec

center

dump

abar

Mean atomic weight.

name

type

unit

centering

storage

abar

float

g/mol

center

dump

zbar

Mean atomic charge.

name

type

unit

centering

storage

zbar

float

?

center

dump

xkn

Rosseland mean opacity.

name

type

unit

centering

storage

xkn

float

cm**2/g

center

dump

xnei

Number density of ionized electrons.

name

type

unit

centering

storage

xnei

float

1/cc

center

dump

stot

Total zonal entropy ( $k_\mathrm{B}$ $/$ baryon).

name

type

unit

centering

storage

stot

float

-

center

dump

angjx

Specific angular momentum x-component.

name

type

unit

centering

storage

angjx

float

cm**2/sec

center

dump

Added In Version 17.10.

angjy

Specific angular momentum y-component.

name

type

unit

centering

storage

angjy

float

cm**2/sec

center

dump

Added In Version 17.10.

angjz

Specific angular momentum z-component.

name

type

unit

centering

storage

angjz

float

cm**2/sec

center

dump

Added In Version 17.10.

Saved (Zone-Centered) Nuclear Burn Derivative Arrays (APPROX and ISE)

dsold

Density at which the energy generation rate was last calculated.

name

type

unit

centering

storage

dsold

float

g/cc

center

dump

tsold

Temperature at which the energy generation rate was last calculated.

name

type

unit

centering

storage

tsold

float

K

center

dump

snold

Last calculated energy generation rate.

name

type

unit

centering

storage

snold

float

erg/g/sec

center

zdum

snbd

Last calculated partial derivative of the energy generation rate with respect to density.

name

type

unit

centering

storage

snbd

float

cm**5/g/sec**3

center

zdum

snbt

Last calculated partial derivative of the energy generation rate with respect to temperature.

name

type

unit

centering

storage

snbt

float

erg/g/sec/K

center

zdum

abarold

Last calculated value of $\bar{A}$ .

name

type

unit

centering

storage

abarold

float

g/mol

center

zdum

abarnbd

Last calculated partial derivative of $\bar{A}$ with respect to density (cc/g).

name

type

unit

centering

storage

abarnbd

float

cc/mol

center

zdum

abarnbt

Last calculated partial derivative of $\bar{A}$ with respect to temperature.

name

type

unit

centering

storage

abarnbt

float

g/mol/K

center

zdum

ypbtime

Last calculated total time derivative of the proton abundance.

name

type

unit

centering

storage

ypbtime

float

mol/g/sec

center

zdum

ynbtime

Last calculated total time derivative of the neutron abundance.

name

type

unit

centering

storage

ynbtime

float

mol/g/sec

center

zdum

Saved Array of Ion Abundances (Generally Elemental)

ppn
Ion Abundance Array.

name

type

unit

centering

storage

ppn

float

mol/g

center

dump
ppn(i,j) = ppn((j-1)*imax+i)
This stored the APPROX. ISE, and NSE abundaces but is not directly accessible on the command line

Saved Zonal BURN Co-Processor Arrays

netnumb

Number of the BURN isotopic nuclear network used in this zone.

name

type

unit

centering

storage

netnumb

integer

-

center

dump

limnuc

Integer containing the Z and A of the isotope controlling the BURN timestep along with the KEPLER cycle that this zone was last burned coded in the form:.

name

type

unit

centering

storage

limnuc

integer

-

center

loadbuf
LIMNUC  = A + 100*Z + 10000* NCYC
Note

I suspect this will not work properly now ot in the future - as we have done runs with > 10M time steps. NEEDS UPDATE

timen

Time to which the BURN coprocessing in this zone has been updated.

name

type

unit

centering

storage

timen

float

sec

center

dump

dtimen

Current timestep for BURN coprocessing in this zone.

name

type

unit

centering

storage

dtimen

float

sec

center

dump

dnold

Density at which BURN coprocessing was last done in this zone.

name

type

unit

centering

storage

dnold

float

g/cc

center

dump

tnold

Temperature at which BURN coprocessing was last done in this zone.

name

type

unit

centering

storage

tnold

float

K

center

dump

ymb

Mass interior to the outer boundary of this zone, as last considered by the BURN coprocessor.

name

type

unit

centering

storage

ymb

float

g

center

dump

sburn

BURN nuclear energy generation rate, excluding neutrino losses.

name

type

unit

centering

storage

sburn

float

erg/g/sec

center

dump

etab

BURN neutron excess.

name

type

unit

centering

storage

etab

float

mol/g

center

dump
The corresponding BURN electron abundance is given by:
yeburn = 0.5 * (1. - etab)

sadv

Advection specofoc energy generation rate.

name

type

unit

centering

storage

sadv

float

erg/g/sec

center

dump

Added In Version 16.82.

sv

Shear viscous heating rate.

name

type

unit

centering

storage

sv

float

erg/g/sec

center

dump

Added In Version 16.98.

pbuf

Abundance buffer (internal working array - see subroutine CONVECTB).

name

type

unit

centering

storage

pbuf

float

mol/g

center

dump

Zonal ‘Elemental’ Mass Fractions And Related Variables

The edit or dump variables corresponding to the elemental mass fractions of the ‘ions’ involved in the APPROX, NSE, and ISE networks take the form, ‘ionfe54‘, where fe54 exemplifies any valid network ‘ion’ symbol. Valid ion symbols are shown below in bold-face letters in the table that follows. The symbol with which each ion’s curve is labeled in abundance plots is shown in square brackets, followed by (in the case of ISE or NSE) the isotopes whose mass fractions are explicitly summed to get each ‘elemental’ mass fraction. A list of the precise isotopes included is given further down.

APPROX ion list

ION#	Network Symbol	Plot Symbol
1	nt1	n
2	h1	H
3	pn1	p
4	he3	$^3\mathrm{He}$
5	he4	He
6	c12	C
7	n14	N
8	o16	O
9	ne20	Ne
10	mg24	Mg
11	si28	Si
12	s32	S
13	ar36	Ar
14	ca40	Ca
15	ti44	Ti
16	cr48	Cr
17	fe52	$^{52}\mathrm{Fe}$
18	fe54	Fe
19	ni56	Ni

ISE and NSE Networks ion list

ION#	Network symbol	Plot Symbol	description
1	nt1	n	neutron only
2	ye	-	electron abundance
3	pn1	p	proton only
4	yq	-	$A \ge$ `24` abundances
5	he4	He	`2` $\le A \le$ `5`
6	yf	-	$A >$ `46` plus $^{46}\mathrm{Ti}$ abundance
7	eb0	-	binding energy / nucleon
8	o16	O	‘unburned’ $^{16}\mathrm{O}$ only
9	fe56	$^{56}\mathrm{Fe}$	$^{56}\mathrm{Fe}$ only
10	mg24	Mg	`23` $<= A <=$ `28`, excluding $^{28}\mathrm{Si}$
11	si28	Si	$^{28}\mathrm{Si}$ only
12	s32	S	`29` $\le A \le$ `35`
13	ar36	Ar	`36` $\le A \le$ `39`
14	ca40	Ca	`40` $\le A \le$ `43`
15	ti44	Ti	`44` $\le A \le$ `47`
16	cr48	Cr	`48` $\le A \le$ `51`
17	‘fe’	‘Fe’	$A \ge$ `2` * $Z +$ `4` Iron Peak, excluding $^{56}\mathrm{Fe}$
18	fe54	Fe	$A \approx$ `2` * $Z +$ `2` Iron Peak
19	ni56	Ni	$A \le$ `2` * $Z +$ `1` Iron Peak

The distinction between h1 and pn1 in the APPROX network is that the former represents unburned hydrogen, while the latter represents protons resulting from photodisintegration. This separation makes inverting the resulting reaction matrix more efficient. Note also that, while not explicitly included in the APPROX network, flows proceeding through the isotopes $^{27}\mathrm{\!Al}$ , $^{31}\mathrm{P}$ , $^{35}\mathrm{\!Cl}$ , $^{39}\mathrm{K}$ , $^{43}\mathrm{Sc}$ , $^{47}\mathrm{V}$ , $^{51}\mathrm{Mn}$ , $^{55}\mathrm{\!Co}$ , $^{53}\mathrm{Fe}$ , and $^{2}\mathrm{H}$ are included implicitly by assuming that their abundances are in steady state with their neighboring nuclei (see [WZW78]).

For zones using the ISE or NSE networks, the exact list of isotopes whose mass fractions are summed to get the edited values for each ‘ion’ symbol are as follows:

Symbol	Ion List
ionnt1	$\mathrm{n}$
ionpn1	$^{1}\mathrm{H}$
ionhe4	$^{2}\mathrm{H}$ + $^{3}\mathrm{H}$ + $^{3}\mathrm{He}$ + $^{4}\mathrm{He}$ + $^{5}\mathrm{He}$ + $^{5}\mathrm{Li}$
ionmg24	$^{23}\mathrm{Na}$ + $^{24}\mathrm{Mg}$ + $^{25}\mathrm{Mg}$ + $^{26}\mathrm{Mg}$ + $^{26}\mathrm{\!Al}$ + $^{27}\mathrm{\!Al}$ + $^{28}\mathrm{\!Al}$
ion016	$^{16}\mathrm{O}$
ionsi28	$^{28}\mathrm{Si}$
ions32	$^{29}\mathrm{Si}$ + $^{30}\mathrm{Si}$ + $^{30}\mathrm{P}$ + $^{31}\mathrm{P}$ + $^{31}\mathrm{S}$ + $^{32}\mathrm{P}$ + $^{32}\mathrm{S}$ + $^{33}\mathrm{P}$ + $^{33}\mathrm{S}$ + $^{34}\mathrm{S}$ + $^{35}\mathrm{\!Cl}$
ionar36	$^{36}\mathrm{\!Cl}$ + $^{36}\mathrm{\!Ar}$ + $^{37}\mathrm{\!Cl}$ + $^{37}\mathrm{\!Ar}$ + $^{38}\mathrm{\!Ar}$ + $^{39}\mathrm{\!Ar}$ + $^{39}\mathrm{K}$
ionca40	$^{40}\mathrm{\!Ar}$ + $^{40}\mathrm{K}$ + $^{40}\mathrm{\!Ca}$ + $^{41}\mathrm{K}$ + $^{41}\mathrm{\!Ca}$ + $^{42}\mathrm{K}$ + $^{42}\mathrm{\!Ca}$ + $^{43}\mathrm{K}$ + $^{43}\mathrm{\!Ca}$ + $^{43}\mathrm{Sc}$
ionti44	$^{44}\mathrm{\!Ca}$ + $^{44}\mathrm{Sc}$ + $^{44}\mathrm{Ti}$ + $^{45}\mathrm{Sc}$ + $^{45}\mathrm{Ti}$ + $^{46}\mathrm{\!Ca}$ + $^{46}\mathrm{Sc}$ + $^{46}\mathrm{Ti}$ + $^{47}\mathrm{Sc}$ + $^{47}\mathrm{Ti}$ + $^{47}\mathrm{V}$
ioncr48	$^{48}\mathrm{\!Ca}$ + $^{48}\mathrm{Sc}$ + $^{48}\mathrm{Ti}$ + $^{48}\mathrm{V}$ + $^{48}\mathrm{\!Cr}$ + $^{49}\mathrm{Sc}$ + $^{49}\mathrm{Ti}$ + $^{49}\mathrm{V}$ + $^{49}\mathrm{\!Cr}$ + $^{50}\mathrm{Ti}$ + $^{50}\mathrm{V}$ + $^{50}\mathrm{\!Cr}$ + $^{51}\mathrm{Ti}$ + $^{51}\mathrm{V}$ + $^{51}\mathrm{\!Cr}$ + $^{51}\mathrm{Mn}$
ion’fe’	$^{52}\mathrm{Ti}$ + $^{52}\mathrm{V}$ + $^{52}\mathrm{\!Cr}$ + $^{53}\mathrm{V}$ + $^{53}\mathrm{\!Cr}$ + $^{54}\mathrm{V}$ + $^{54}\mathrm{\!Cr}$ + $^{54}\mathrm{Mn}$ + $^{55}\mathrm{\!Cr}$ + $^{55}\mathrm{Mn}$ + $^{56}\mathrm{\!Cr}$ + $^{56}\mathrm{Mn}$ + $^{57}\mathrm{Mn}$ + $^{57}\mathrm{Fe}$ + $^{58}\mathrm{Mn}$ + $^{58}\mathrm{Fe}$ + $^{58}\mathrm{\!Co}$ + $^{59}\mathrm{Fe}$ + $^{59}\mathrm{\!Co}$ + $^{60}\mathrm{Fe}$ + $^{60}\mathrm{\!Co}$ + $^{60}\mathrm{Ni}$ + $^{61}\mathrm{Fe}$ + $^{61}\mathrm{\!Co}$ + $^{61}\mathrm{Ni}$ + $^{62}\mathrm{Fe}$ + $^{62}\mathrm{\!Co}$ + $^{62}\mathrm{Ni}$ + $^{63}\mathrm{\!Co}$ + $^{63}\mathrm{Ni}$ + $^{64}\mathrm{\!Co}$ + $^{64}\mathrm{Ni}$ + $^{65}\mathrm{Ni}$ + $^{66}\mathrm{Ni}$
ionfe54	$^{52}\mathrm{Mn}$ + $^{53}\mathrm{Mn}$ + $^{54}\mathrm{Fe}$ + $^{55}\mathrm{\!Co}$ + $^{56}\mathrm{\!Co}$ + $^{57}\mathrm{\!Co}$ + $^{58}\mathrm{Ni}$ + $^{59}\mathrm{Ni}$
ionni56	$^{52}\mathrm{Fe}$ + $^{53}\mathrm{Fe}$ + $^{54}\mathrm{\!Co}$ + $^{55}\mathrm{\!Co}$ + $^{56}\mathrm{Ni}$ + $^{57}\mathrm{Ni}$ + $^{59}\mathrm{\!Cu}$ + $^{60}\mathrm{Zn}$
ionfe56	$^{56}\mathrm{Fe}$

Note

The special cases where the edit variables formed from the ‘ion’ symbols listed do not correspond to ‘elemental’ mass fractions are as follows:

Symbol	Description
ionye	electron abundance in any current (non-BURN) network (`mol/g`)
ionyq	total abundance of ISE or NSE isotopes with $A\ge$ `24` (`mol/g`)
ionyf	total abundance of ISE or NSE isotopes with $A>$ `46` plus $^{46}\mathrm{Ti}$ (`mol/g`)
ioneb0	mean binding energy per nucleon of all ISE or NSE isotopes (`MeV`)

Note

if lburn (p 434) is enabled, ye is set to the BURN value.

Ise Edit Quantities

Zonal edit variables for ISE network isotopic mass-fractions, weak rates, and related quantities can be referenced in the form: ‘insefe56‘ where ‘nse‘ is a fixed flag, ‘fe56‘ represents the symbol of an isotope in the ISE network (stored in array nucsymd in subroutine qnse - see below), and ‘r‘ represents a variable flag depending on which isotopic quantity or rate is desired. ‘i‘ may take the values to get:

value	result
’i‘	isotopic mass fraction,
’w‘	isotopic statistical weight,
’r‘	total weak interaction rate due to this isotope (`mol/g/sec`),
’c‘	electron capture rate due to this isotope (`mol/g/sec`)
’p‘	positron emission rate due to this isotope (`mol/g/sec`)
’b‘	beta decay rate due to this isotope (`mol/g/sec`)
’e‘	total neutrino energy loss due to this isotope (`erg/g/sec`)
’k‘	neutrino energy loss due to e- capture and e+ decay on this isotope (`erg/g/sec`)
’x‘	energy stored in the excited states of this isotope (`MeV/nucleus`)

The symbols of the 125 isotopes comprising the ISE (and NSE) network are:

n1   h1   h2   h3   he3  he4  he5  li5  c12  o16   (  1- 10)
ne20 na23 mg24 mg25 mg26 al26 al27 al28 si28 si29  ( 11- 20)
si30 p30  p31  s31  p32  s32  p33  s33  s34  cl35  ( 21- 30)
cl36 ar36 cl37 ar37 ar38 ar39 k39  ar40 k40  ca40  ( 31- 40)
k41  ca41 k42  ca42 k43  ca43 sc43 ca44 sc44 ti44  ( 41- 50)
sc45 ti45 ca46 sc46 ti46 sc47 ti47 v47  ca48 sc48  ( 51- 60)
ti48 v48  cr48 sc49 ti49 v49  cr49 ti50 v50  cr50  ( 61- 70)
ti51 v51  cr51 mn51 ti52 v52  cr52 mn52 fe52 v53   ( 71- 80)
cr53 mn53 fe53 v54  cr54 mn54 fe54 co54 cr55 mn55  ( 81- 90)
fe55 co55 cr56 mn56 fe56 co56 ni56 mn57 fe57 co57  ( 91-100)
ni57 mn58 fe58 co58 ni58 fe59 co59 ni59 cu59 fe60  (101-110)
co60 ni60 zn60 fe61 co61 ni61 fe62 co62 ni62 co63  (111-120)
ni63 co64 ni64 ni65 ni66                           (120-125)

Zonal Nuclear Reaction Rates

The edit or dump variables corresponding to editable nuclear reaction rates are listed below. The values edited are the zonal net rates (forward-back) in mol/g/sec, and take into account flows through intermediate nuclei that are only implicitly included in the APPROX network (see [WZW78] and sdot). The indices correspond to location in the rr array.

number

name

reaction

networks

1

p+p

p+p

APPROX

2

he3+he3

3He+3He

APPROX

3

he3+he4

3He+4He

APPROX

4

c12(pg)

12C(p,g)13N

APPROX

5

n14(pg)

14N(p,g)15O

APPROX

6

o16(pg)

16O(p,g)17F

APPROX

7

3a-c12

4He(2a,g)12C

APPROX

8

n14(ag)

14N(a,g)18F

APPROX

9

c12(ag)

12C(a,g)16O

APPROX

10

c12+c12

12C+12C

APPROX

11

c12+o16

12C+16O

APPROX

12

o16+o16

16O+16O

APPROX, ISE

13

o16(ag)

16O(a,g)20Ne

APPROX

14

ne20(ag)

20Ne(a,g)24Mg

APPROX, ISE

15

mg24(ag)

24Mg(a,g)28Si

APPROX

16

si28(ag)

28Si(a,g)32S

APPROX

17

s32(ag)

32S(a,g)36Ar

APPROX

18

ar36(ag)

36Ar(a,g)40Ca

APPROX

19

ca40(ag)

40Ca(a,g)44Ti

APPROX

20

ti44(ag)

44rri(a,g)48Cr

APPROX

21

cr48(ag)

48Cr(a,g)52Fe

APPROX

22

fe52(ag)

52Fe(a,g)56Ni

APPROX

23

fe52:2ng

52Fe(2n,g)54Fe

APPROX

24

fe52:a2p

52Fe(a,2p)54Fe

APPROX

25

fe54:2pg

54Fe(2p,g)56Ni

APPROX

26

he4-2n2p

4He(g,2p2n)

APPROX

27

p(e-nu)n

p(e-,v)n

APPROX, ISE

28

n(e+nu)p

n(e+,v)p

APPROX, ISE

29

weak

neutronization rate

ISE

In detail:

p+p

(APPROX only) $^{1}\mathrm{H}$ + $^{1}\mathrm{H}$ $\rightarrow$ $^{2}\mathrm{H}$ + ${\mathrm{e}}^{+}$ + ${\nu}_{\mathrm{e}}$ followed by $^{2}\mathrm{H}$ + $^{1}\mathrm{H}$ $\rightarrow$ $^{3}\mathrm{He}$ .

name

type

unit

centering

storage

p+p

float

-

center

loadbuf

he3+he3

(APPROX only) $^{3}\mathrm{He}$ + $^{3}\mathrm{He}$ $\rightarrow$ $^{4}\mathrm{He}$ + $^{1}\mathrm{H}$ + $^{1}\mathrm{H}$ .

name

type

unit

centering

storage

he3+he3

float

-

center

loadbuf

he3+he4

(APPROX only) $^{3}\mathrm{He}$ + $^{4}\mathrm{He}$ $\rightarrow$ $^{7}\mathrm{Be}$ followed by PPII and PPIII.

name

type

unit

centering

storage

he3+he4

float

-

center

loadbuf

c12(pg)

(APPROX only) $^{12}\mathrm{\!C}$ + $^{1}\mathrm{H}$ $\rightarrow$ $^{13}\mathrm{N}$ followed by $^{13}\mathrm{N}$ $\rightarrow$ $^{13}\mathrm{\!C}$ + ${\mathrm{e}}^{+}$ + ${\nu}_{\mathrm{e}}$ and $^{13}\mathrm{\!C}$ + $^{1}\mathrm{H}$ $\rightarrow$ $^{14}\mathrm{N}$ .

name

type

unit

centering

storage

c12(pg)

float

-

center

loadbuf

n14(pg)

14N(p,g)15O (APPROX only).

name

type

unit

centering

storage

n14(pg)

float

-

center

loadbuf

o16(pg)

160(p,g)17F (APPROX only).

name

type

unit

centering

storage

o16(pg)

float

-

center

loadbuf

3a-c12

4He(2a,g)12C (APPROX only).

name

type

unit

centering

storage

3a-c12

float

-

center

loadbuf

n14(ag)

14N(a,g)18F (APPROX only).

name

type

unit

centering

storage

n14(ag)

float

-

center

loadbuf

c12(ag)

12C(a,g)16O (APPROX only).

name

type

unit

centering

storage

c12(ag)

float

-

center

loadbuf

c12+c12

12C+12C (APPROX only).

name

type

unit

centering

storage

c12+c12

float

-

center

loadbuf

c12+o16

12C+16O (APPROX only).

name

type

unit

centering

storage

c12+o16

float

-

center

loadbuf

o16+o16

16O+16O (APPROX and ISE).

name

type

unit

centering

storage

o16+o16

float

-

center

loadbuf

o16(ag)

16O(a,g)20Ne (APPROX only).

name

type

unit

centering

storage

o16(ag)

float

-

center

loadbuf

ne20(ag)

20Ne(a,g)24Mg (APPROX and ISE).

name

type

unit

centering

storage

ne20(ag)

float

-

center

loadbuf

mg24(ag)

24Mg(a,g)28Si (APPROX only).

name

type

unit

centering

storage

mg24(ag)

float

-

center

loadbuf

si28(ag)

28Si(a,g)32S (APPROX only).

name

type

unit

centering

storage

si28(ag)

float

-

center

loadbuf

s32(ag)

32S(a,g)36Ar (APPROX only).

name

type

unit

centering

storage

s32(ag)

float

-

center

loadbuf

ar36(ag)

36Ar(a,g)40Ca (APPROX only).

name

type

unit

centering

storage

ar36(ag)

float

-

center

loadbuf

ca40(ag)

4OCa(a,g)44Ti (APPROX only).

name

type

unit

centering

storage

ca40(ag)

float

-

center

loadbuf

ti44(ag)

44rri(a,g)48Cr (APPROX only).

name

type

unit

centering

storage

ti44(ag)

float

-

center

loadbuf

cr48(ag)

48Cr(a,g)52Fe (APPROX only).

name

type

unit

centering

storage

cr48(ag)

float

-

center

loadbuf

fe52(ag)

52Fe(a,g)56Ni (APPROX only).

name

type

unit

centering

storage

fe52(ag)

float

-

center

loadbuf

fe52:2ng

52Fe(2n,g)54Fe (APPROX only).

name

type

unit

centering

storage

fe52:2ng

float

-

center

loadbuf

fe52:a2p

52Fe(a,2p)54Fe (APPROX only).

name

type

unit

centering

storage

fe52:a2p

float

-

center

loadbuf

fe54:2pg

54Fe(2p,g)56Ni (APPROX only).

name

type

unit

centering

storage

fe54:2pg

float

-

center

loadbuf

he4-2n2p

4He(g,2p2n) (APPROX only).

name

type

unit

centering

storage

he4-2n2p

float

-

center

loadbuf

p(e-nu)n

P(e-,v)n (APPROX and ISE).

name

type

unit

centering

storage

p(e-nu)n

float

-

center

loadbuf

n(e+nu)p

N(e+,v)p (APPROX and ISE).

name

type

unit

centering

storage

n(e+nu)p

float

-

center

loadbuf

weak

Neutronization rate (ISE only).

name

type

unit

centering

storage

weak

float

-

center

loadbuf

Here p denotes a proton; n, a neutron; a, an alpha particle; g, a ganuna-ray; nu, a neutrino; and ecap, the rate for electron capture (ISE zones only).

Zonal BURN Isotope Mass Fractions

The edit or dump variables corresponding to the zonal isotopic mass fractions of the isotopes involved in the BURN coprocessor network take the form ‘isoal26‘, where al26 exemplifies any valid BURN isotope symbol. Valid isotope symbols depend on the BURN network specified in the BURN generator.

A commonly used set of 149 isotopes is given by:

ntl  h1   h2   h3   he3  he4  li7  be7  b10  bl1   (  1- 10)
cll  c12  c13  c14  n13  n14  n15  o15  o16  o17   ( 11- 20)
o18  f17  f18  f19  ne19 ne20 ne21 ne22 ne23 na22  ( 21- 30)
na23 na24 mg23 mg24 mg25 mg26 mg27 al26 al27 aI28  ( 31- 40)
si27 si28 si29 si30 si31 p30  p31  p32  p33  s31   ( 41- 50)
s32  s33  s34  s35  s36  s37  cl35 cl36 cI37 cl38  ( 51- 60)
ar36 ar37 ar38 ar39 ar40 ar41 k39  k40  k41  k42   ( 61- 70)
ca40 ca41 ca42 ca43 ca44 ca45 ca46 ca47 sc43 sc44  ( 71- 80)
sc45 sc46 sc47 ti44 ti45 ti46 ti47 ti48 ti49 ti50  ( 81- 90)
ti51 v47  v48  v49  v50  v51  v52  cr48 cr49 cr50  ( 91-100)
cr51 cr52 cr53 cr54 crS5 mn51 mn52 mn53 mn54 mn55  (101-110)
mn56 fe52 fe53 fe54 fe55 fe56 feS7 fe58 feS9 fe60  (111-120)
fe61 co55 co56 co57 co58 co59 co60 co61 niS6 ni57  (121-130)
ni58 ni59 ni60 ni61 ni62 ni63 ni64 ni65 cu59 cu60  (131-140)
cu61 cu62 cu63 cu64 zn60 zn61 zn62 zn63 zn64       (141-149)

Editcom Zonal Variables

Another set of editable zonal quantities are those listed in the editcom common block with an ‘ed‘ appended to the beginning of their names to identify them as edit quantities. Note that the ASCII names of these variables are specified in the the nameedit array and an ASCII string specifying their units is contained in the unitedit array (see block data KEPDAT). This structure makes it easy to add additional zonal edit variables. The current set of editcom zonal variables is as follows:

Derivative `editcom` Variables

gamma

Adiabatic exponent defined by $\gamma=\rho\,(dP/d\rho)_\mathrm{ad}/P=\Gamma_1$ .

name

type

unit

centering

storage

gamma

float

-

center

loadbuf

gammae

Adiabatic exponent for electrons defined by $\gamma_\mathrm{e}=\rho\,(\mathrm{d}P_\mathrm{e}/\mathrm{d}\rho)_\mathrm{ad}/P_\mathrm{e}=\Gamma_\mathrm{1e}$ .

name

type

unit

centering

storage

gammae

float

-

center

loadbuf

gammaf

Adiabatic exponent defined by $\gamma_\mathrm{f}=P\,(\mathrm{d}T/\mathrm{d}P)_\mathrm{ad}/T=(\Gamma_2-1)/\Gamma_2$ .

name

type

unit

centering

storage

gammaf

float

-

center

loadbuf

eit

Partial derivative of the ion specific energy with respect to temperature.

name

type

unit

centering

storage

eit

float

erg/g/K

center

loadbuf

eet

Partial derivative of the electron specific energy with respect to temperature including pairs.

name

type

unit

centering

storage

eet

float

erg/g/K

center

loadbuf

ept

Partial derivative of the positron specific energy with respect to temperature.

name

type

unit

centering

storage

ept

float

erg/g/K

center

loadbuf

ert

Partial derivative of the radiative specific energy with respect to temperature.

name

type

unit

centering

storage

ert

float

erg/g/K

center

loadbuf

eid

Partial derivative of the ion specific energy with respect to density.

name

type

unit

centering

storage

eid

float

erg/cc/g**2

center

loadbuf

eed

Partial derivative of the electron specific energy with respect to density including pairs.

name

type

unit

centering

storage

eed

float

erg/cc/g**2

center

loadbuf

epd

Partial derivative of the positron specific energy with respect to density.

name

type

unit

centering

storage

epd

float

erg/cc/g**2

center

loadbuf

erd

Partial derivative of the radiative specific energy with respect to density.

name

type

unit

centering

storage

erd

float

erg/cc/g**2

center

loadbuf

Note

Original manual says: (Not implemented yet).

pit

Partial derivative of the ion pressure with respect to temperature.

name

type

unit

centering

storage

pit

float

erg/cc/K

center

loadbuf

pet

Partial derivative of the electron pressure with respect to temperature including pairs.

name

type

unit

centering

storage

pet

float

erg/cc/K

center

loadbuf

ppl

Partial derivative of the positron pressure with respect to temperature.

name

type

unit

centering

storage

ppl

float

erg/cc/K

center

loadbuf

prt

Partial derivative of the radiation pressure with respect to temperature.

name

type

unit

centering

storage

prt

float

erg/cc/K

center

loadbuf

pid

Partial derivative of the ion pressure with respect to density.

name

type

unit

centering

storage

pid

float

erg/g

center

loadbuf

ped

Partial derivative of the electron pressure with respect to density including pairs.

name

type

unit

centering

storage

ped

float

erg/g

center

loadbuf

ppd

Partial derivative of the positron pressure with respect to density.

name

type

unit

centering

storage

ppd

float

erg/g

center

loadbuf

prd

Partial derivative of the radiation pressure with respect to density.

name

type

unit

centering

storage

prd

float

erg/g

center

loadbuf

xnebt

Partial derivative of the electron number density with respect to temperature including pairs.

name

type

unit

centering

storage

xnebt

float

1/cc/K

center

loadbuf

xnebd

Partial derivative of the electron number density with respect to density including pairs.

name

type

unit

centering

storage

xnebd

float

1/g

center

loadbuf

etabt

Partial derivative of the negatron degeneracy parameter with respect to temperature.

name

type

unit

centering

storage

etabt

float

1/K

center

loadbuf

etabd

Partial derivative of the negatron degeneracy parameter with respect to density.

name

type

unit

centering

storage

etabd

float

cc/g

center

loadbuf

ent

Partial derivative of the total specific energy with respect to temperature.

name

type

unit

centering

storage

ent

float

erg/g/K

center

loadbuf

end

Partial derivative of the total specific ,energy with respect to density.

name

type

unit

centering

storage

end

float

erg/cc/g**2

center

loadbuf

pnt

Partial derivative of the total pressure with respect to temperature.

name

type

unit

centering

storage

pnt

float

erg/cc/K

center

loadbuf

pnd

Partial derivative of the total pressure with respect to density.

name

type

unit

centering

storage

pnd

float

erg/g

center

loadbuf

xknt

Partial derivative of the total opacity with respect to temperature.

name

type

unit

centering

storage

xknt

float

cm**2/g/K

center

loadbuf

xknd

Partial derivative of the total opacity with respect to density.

name

type

unit

centering

storage

xknd

float

cm**5/g**2

center

loadbuf

snt

Partial derivative of the energy generation rate with respect to temperature.

name

type

unit

centering

storage

snt

float

erg/g/K

center

loadbuf

snd

Partial derivative of the energy generation rate with respect to density.

name

type

unit

centering

storage

snd

float

erg*cc/g**2

center

loadbuf

Note

Reserved variable spaces in editcom common block:

derivl  Reserved for future derivative edit.
deriv2  Reserved for future derivative edit.
deriv3  Reserved for future derivative edit.
deriv4  Reserved for future derivative edit.
derivS  Reserved for future derivative edit.

EOS `editcom` Variables:

sig

Total entropy.

name

type

unit

centering

storage

sig

float

kb/baryon

center

loadbuf

sigi

Ion entropy.

name

type

unit

centering

storage

sigi

float

kb/baryon

center

loadbuf

sige

Electron entropy (negatrons only).

name

type

unit

centering

storage

sige

float

kb/baryon

center

loadbuf

sigp

Positron entropy.

name

type

unit

centering

storage

sigp

float

kb/baryon

center

loadbuf

sigr

Radiation entropy.

name

type

unit

centering

storage

sigr

float

kb/baryon

center

loadbuf

sigion

Ionizatiqn entropy.

name

type

unit

centering

storage

sigion

float

kb/baryon

center

loadbuf

ei

Specific ion energy.

name

type

unit

centering

storage

ei

float

erg/g

center

loadbuf

ee

Specific electron energy including pairs.

name

type

unit

centering

storage

ee

float

erg/g

center

loadbuf

ep

Specific positron energy.

name

type

unit

centering

storage

ep

float

erg/g

center

loadbuf

er

Specific radiation energy.

name

type

unit

centering

storage

er

float

erg/g

center

loadbuf

pi

Ion pressure.

name

type

unit

centering

storage

pi

float

erg/cc

center

loadbuf

pe

Electron pressure including pairs.

name

type

unit

centering

storage

pe

float

erg/cc

center

loadbuf

pp

Positron pressure.

name

type

unit

centering

storage

pp

float

erg/cc

center

loadbuf

pr

Radiation pressure.

name

type

unit

centering

storage

pr

float

erg/cc

center

loadbuf

xni

Ion number density.

name

type

unit

centering

storage

xni

float

1/cc

center

loadbuf

xne

Electron number density including pairs.

name

type

unit

centering

storage

xne

float

1/cc

center

loadbuf

xnemc

Electron number density without contribution from pairs (#/cc).

name

type

unit

centering

storage

xnemc

float

-

center

loadbuf

xnp

Positron (or pair) number density.

name

type

unit

centering

storage

xnp

float

1/cc

center

loadbuf

eta

Electron (negatron) degeneracy parameter ( $\mu_\mathrm{e}/k_\mathrm{B}T$ ).

name

type

unit

centering

storage

eta

float

-

center

loadbuf

etap

Posltron degeneracy parameter ( $\mu_\mathrm{positron}/k_\mathrm{B}T$ ).

name

type

unit

centering

storage

etap

float

-

center

loadbuf

xkibenl

IBEN1 Rosselandmean opacity.

name

type

unit

centering

storage

xkibenl

float

cm**2/g

center

loadbuf

xkiben2

IBEN2 Rosseland mean opacity.

name

type

unit

centering

storage

xkiben2

float

cm**2/g

center

loadbuf

xkchrsty

Christy Rosseland mean opacity.

name

type

unit

centering

storage

xkchrsty

float

cm**2/g

center

loadbuf

xkcompt

Rosseland mean opacity from Compton scattering.

name

type

unit

centering

storage

xkcompt

float

cm**2/g

center

loadbuf

Note

Original manual says: (Not implemented yet).

xkr

Total radiative Rosseland mean opacity.

name

type

unit

centering

storage

xkr

float

cm**2/g

center

loadbuf

xkcond

Effective Rosseland mean opacity due to electron heat conduction.

name

type

unit

centering

storage

xkcond

float

cm**2/g

center

loadbuf

xkcomp

Compton scattering opacity.

name

type

unit

centering

storage

xkcomp

float

cm**2/g

center

loadbuf

xnem

Matter electron number density.

name

type

unit

centering

storage

xnem

float

1/cc

center

loadbuf

xnembt

Temperature derivative of Matter electron number density.

name

type

unit

centering

storage

xnembt

float

1/cc/K

center

loadbuf

xnembd

Density derivative of Matter electron number density.

name

type

unit

centering

storage

xnembd

float

1/g

center

loadbuf

pcoul

Pressure contribution due to Coulomb corrections.

name

type

unit

centering

storage

pcoul

float

erg/cc

center

loadbuf

ecoul

Specific internal energy contribution due to Coulomb corrections.

name

type

unit

centering

storage

ecoul

float

erg/g

center

loadbuf

Note

Reserved variable spaces in editcom common block:

eosl    Reserved for future EOS quantity.
eos2    Reserved for future EOS quantity.

ISE EDITCOM Variables:

dye

Total change in $Y_\mathrm{e}$ during the current timestep.

name

type

unit

centering

storage

dye

float

mol/g

center

loadbuf

dyq

Total change in $Y_\mathrm{q}$ during the current timestep.

name

type

unit

centering

storage

dyq

float

mol/g

center

loadbuf

dyf

Total change in $Y_\mathrm{f}$ during the current timestep.

name

type

unit

centering

storage

dyf

float

mol/g

center

loadbuf

iternse

Total number of iterations it took for ISE/NSE to converge in this zone, or to -1 if convergence was not achieved.

name

type

unit

centering

storage

iternse

integer

-

center

loadbuf

qxc12

Photodisintegration-flow-deterrriined mass fraction of $^{12}\mathrm{\!C}$ in the ISE calculation.

name

type

unit

centering

storage

qxc12

float

-

center

loadbuf

qxo16

Photodisintegration-flow-deterrnined mass fraction of $^{16}\mathrm{O}$ in the ISE calculation.

name

type

unit

centering

storage

qxo16

float

-

center

loadbuf

qxne

Photodisintegration-flow-deterrnined mass fraction of $^{20}\mathrm{Ne}$ in the ISE calculation.

name

type

unit

centering

storage

qxne

float

-

center

loadbuf

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.

name

type

unit

centering

storage

betaq

float

-

center

loadbuf

wrate

Total zonal weak (i.e., neutronization) rate.

name

type

unit

centering

storage

wrate

float

mol/g/sec

center

loadbuf

rectot

Total zonal electron capture rate.

name

type

unit

centering

storage

rectot

float

mol/g/sec

center

loadbuf

rpdtot

Total zonal positron decay rate.

name

type

unit

centering

storage

rpdtot

float

mol/g/sec

center

loadbuf

redtot

Total zonal electron decay rate.

name

type

unit

centering

storage

redtot

float

mol/g/sec

center

loadbuf

snuw

Total zonal specific neutrino energy loss rate due to neutronization-related processes.

name

type

unit

centering

storage

snuw

float

erg/g/sec

center

loadbuf

snubps

Specific neutrino energy loss rate due to plasma-related processes.

name

type

unit

centering

storage

snubps

float

erg/g/sec

center

loadbuf

eectot

Total zonal electron capture plus positron decay neutrino loss rate.

name

type

unit

centering

storage

eectot

float

erg/g/sec

center

loadbuf

epdtot

Formerly, the total zonal positron decay neutrino loss rate.

name

type

unit

centering

storage

epdtot

float

erg/g/sec

center

loadbuf

This quantity is now combined into eectot, andepdtot is set to 0..

eedtot

Total zonal electron decay neutrino loss rate.

name

type

unit

centering

storage

eedtot

float

erg/g/sec

center

loadbuf

ebind0

Average zonal nucleon binding energy at the beginning of the timestep including contribution from energy in nuclear excited states if iexciteh (p 171) $>$ 0.

name

type

unit

centering

storage

ebind0

float

MeV/nucleon

center

loadbuf

ebind

Average zonal nucleon binding energy at the end of the timestep including contribution from energy in excited states if iexciteh (p 171) $>$ 0.

name

type

unit

centering

storage

ebind

float

MeV/nucleon

center

loadbuf

eexcite

Average zonal energy in nuclear excited states.

name

type

unit

centering

storage

eexcite

float

MeV/nucleon

center

loadbuf

snucq

ISE nuclear energy generation rate.

name

type

unit

centering

storage

snucq

float

erg/g/sec

center

loadbuf

snuclq

ISE nuclear energy generation rate minus plasma neutrino loss rate.

name

type

unit

centering

storage

snuclq

float

erg/g/sec

center

loadbuf

sdotq

Net ISE energy generation rate including all neutrino losses.

name

type

unit

centering

storage

sdotq

float

erg/g/sec

center

loadbuf

tauqse

Timescale for the $\mathrm{Fe}$ and $\mathrm{Si}$ quasiequilibrium blocks to equilibrate.

name

type

unit

centering

storage

tauqse

float

sec

center

loadbuf

yfqse

Quasiequilibrium value of $Y_\mathrm{q}$ based on the current abundance of neutron and protons and $^{28}\mathrm{Si}$ .

name

type

unit

centering

storage

yfqse

float

mol/g

center

loadbuf

yffac

The actual value of $Y_\mathrm{q}$ is less than yfqse by the factor 1 $/$ yffac.

name

type

unit

centering

storage

yffac

float

-

center

loadbuf

qj

Net photodisintegration flow downward from “i:24Mg.

name

type

unit

centering

storage

qj

float

mol/g/sec

center

loadbuf

dyo16q

Change in the $^{16}\mathrm{O}$ abundance during this timestep.

name

type

unit

centering

storage

dyo16q

float

mol/g

center

loadbuf

`loadbuf` Edit Quantities

The following miscellaneous zonal edit variables are defined in subroutine loadbuf (and more may easily be added):

Miscellaneous

rhor

Integral of $\rho\,\mathrm{d}\,r$ for mass exterior to the outer boundary of this zone.

name

type

unit

centering

storage

rhor

float

g/cm**2

boundary

loadbuf

Mixing and Transport

mixlen

Mixing length used in calculating convection.

name

type

unit

centering

storage

mixlen

float

cm

boundary

loadbuf

taumix

Convective mixing time.

name

type

unit

centering

storage

taumix

float

sec

boundary

loadbuf

convect

Convection flag.

name

type

unit

centering

storage

convect

float

-

boundary

loadbuf

value

meaning

0.

radiative

1.

neutral

2.

overshooting

3.

semiconvection

4.

converction

5.

thermohaline

Rotation

angj

Specific angular momentum magnitude.

name

type

unit

centering

storage

angj

float

cm**2/sec

center

loadbuf

Added In Version 15.2.

Changed In Version 17.10: now derived quantity, components are stored as primary data instead.

angja

Specific angular momentum magnitude.

name

type

unit

centering

storage

angja

float

cm**2/sec

center

loadbuf

currently same as angj

Added In Version 17.10.

angi

Specific moment of inertia.

name

type

unit

centering

storage

angi

float

cm**2

center

cycle

Added In Version 15.2.

angw

Angular velocity magnitude.

name

type

unit

centering

storage

angw

float

rad/sec

center

cycle

Added In Version 15.2.

angwa

Angular velocity magnitude.

name

type

unit

centering

storage

angwa

float

rad/sec

center

cycle

currently same as angw

Added In Version 17.10.

angwx

Angular velocity x-component.

name

type

unit

centering

storage

angwx

float

rad/sec

center

cycle

Added In Version 17.10.

angwy

Angular velocity y-component.

name

type

unit

centering

storage

angwy

float

rad/sec

center

cycle

Added In Version 17.10.

angwz

Angular velocity z-component.

name

type

unit

centering

storage

angwz

float

rad/sec

center

cycle

Added In Version 17.10.

angwcst

Rotation constant magnitude $\bar{\omega}\,\bar{\rho}^{-2/3}$ .

name

type

unit

centering

storage

angwcst

float

?

center

loadbuf

where

$\bar{\rho}=\frac{3\,m}{4\pi r^3(m)}$

$\bar{\omega}=J(m)/I(m)$

$J(m) = \int_0^m j(m^\prime)\,\mathrm{d}m^\prime$

$I(m) = \int_0^m i(m^\prime)\,\mathrm{d}m^\prime$

$i$ : specific moment of inertia (cm**2)

$j$ : specific angular momentum (cm**2/sec)

Changed In Version 17.10: renamed from angwx to angwcst

angwcsta

Rotation constant magnitude.

name

type

unit

centering

storage

angwcsta

float

?

center

loadbuf

Same as angwcst

Added In Version 17.10.

angwcstx

Rotation constant x-component.

name

type

unit

centering

storage

angwcstx

float

?

center

loadbuf

See angwcst.

Added In Version 17.10.

angwcsty

Rotation constant y-component.

name

type

unit

centering

storage

angwcsty

float

?

center

loadbuf

See angwcst.

Added In Version 17.10.

angwcstz

Rotation constant z-component.

name

type

unit

centering

storage

angwcstz

float

?

center

loadbuf

See angwcst.

Added In Version 17.10.

angr

Zone gyration-average radius.

name

type

unit

centering

storage

angr

float

cm

center

loadbuf

$\sqrt{1.5\,i}$

Added In Version 17.10.

angv

Equatorial rotational velocity magnitude.

name

type

unit

centering

storage

angv

float

cm/sec

center

loadbuf

Added In Version 15.2.

angva

Equatorial rotational velocity magnitude.

name

type

unit

centering

storage

angva

float

cm/sec

center

loadbuf

currently same as angv

Added In Version 17.10.

angvx

Equatorial rotational velocity x-component.

name

type

unit

centering

storage

angvx

float

cm/sec

center

loadbuf

Added In Version 17.10.

angvy

Equatorial rotational velocity y-component.

name

type

unit

centering

storage

angvy

float

cm/sec

center

loadbuf

Added In Version 17.10.

angvz

Equatorial rotational velocity z-component.

name

type

unit

centering

storage

angvz

float

cm/sec

center

loadbuf

Added In Version 17.10.

angri

Rotation/shear flow Richardson Number.

name

type

unit

centering

storage

angri

float

-

center

cycle

Added In Version 15.2.

ange

Specific rotational energy.

name

type

unit

centering

storage

ange

float

erg/g

center

loadbuf

angf

Specific rotational energy.

name

type

unit

centering

storage

angf

float

erg/g

center

loadbuf

angp

Rotation frequency.

name

type

unit

centering

storage

angp

float

Hz

center

loadbuf

angvk

Keplerian rotational velocity.

name

type

unit

centering

storage

angvk

float

cm/sec

center

loadbuf

angwc

Critical angular velocity.

name

type

unit

centering

storage

angwc

float

rad/sec

center

loadbuf

angwwc

Angular velocity / critical velocity.

name

type

unit

centering

storage

angwwc

float

-

center

loadbuf

Same as as angvvc

angvvc

Rotational velocity / critical velocity.

name

type

unit

centering

storage

angvvc

float

-

center

loadbuf

Same as as angwwc

angwk

Keplerian angular velocity.

name

type

unit

centering

storage

angwk

float

rad/sec

center

loadbuf

angjk

Keplerian specific angular momentum.

name

type

unit

centering

storage

angjk

float

cm**2/sec

center

loadbuf

angek

Specific Keplerian energy.

name

type

unit

centering

storage

angek

float

erg/g

center

loadbuf

angwwk

Ratio of angular velocity to Keplerian velocity.

name

type

unit

centering

storage

angwwk

float

-

center

loadbuf

Same as angvvk.

angvvk

Ratio of rotational velocity to Keplerian velocity.

name

type

unit

centering

storage

angvvk

float

-

center

loadbuf

Same as angwwk.

angeek

Ratio of equatorial rotational energy to Keplerian energy.

name

type

unit

centering

storage

angeek

float

-

center

loadbuf

Same as angcg

angcg

Ratio of centripetal force to gravity.

name

type

unit

centering

storage

angcg

float

-

center

loadbuf

Same as angeek

angjjk

Ratio of shell average specific angular momentum to Keplerian specific angular momentum.

name

type

unit

centering

storage

angjjk

float

-

center

loadbuf

angjt

Magnitude integrated angular momentum.

name

type

unit

centering

storage

angjt

float

erg*sec

boundary

loadbuf

angjta

Magnitude integrated angular momentum.

name

type

unit

centering

storage

angjta

float

erg*sec

boundary

loadbuf

Same as angjt

Added In Version 17.10.

angjtx

Integrated angular momentum x-component.

name

type

unit

centering

storage

angjtx

float

erg*sec

boundary

loadbuf

Added In Version 17.10.

angjty

Integrated angular momentum y-component.

name

type

unit

centering

storage

angjty

float

erg*sec

boundary

loadbuf

Added In Version 17.10.

angjtz

Integrated angular momentum z-component.

name

type

unit

centering

storage

angjtz

float

erg*sec

boundary

loadbuf

Added In Version 17.10.

angit

Integrated moment of inertia momentum.

name

type

unit

centering

storage

angit

float

g*cm**2

boundary

loadbuf

angd0

Diffusion coefficient for convective processes.

name

type

unit

centering

storage

angd0

float

cm**2/sec

boundary

loadbuf

same as difi

angd1

Diffusion coefficient for dynamical shear instability.

name

type

unit

centering

storage

angd1

float

cm**2/sec

boundary

loadbuf

from angd(:,:n:1)

angd2

Diffusion coefficient for Solberg-Hoiland instability.

name

type

unit

centering

storage

angd2

float

cm**2/sec

boundary

loadbuf

from angd(:,:n:2)

angd3

Diffusion coefficient for secular shear instability.

name

type

unit

centering

storage

angd3

float

cm**2/sec

boundary

loadbuf

from angd(:,:n:3)

angd4

Diffusion coefficient for Eddington-Sweet circulation.

name

type

unit

centering

storage

angd4

float

cm**2/sec

boundary

loadbuf

from angd(:,:n:4)

angd5

Diffusion coefficient for Goldreich-Schubert-Fricke instability (cm:math:^2/s).

name

type

unit

centering

storage

angd5

float

cm**2/sec

boundary

loadbuf

from angd(:,:n:5)

Mass Coordinate

zm

Interior mass coordinate.

name

type

unit

centering

storage

zm

float

g

center

cycle

zm_sun

Interior mass coordinate.

name

type

unit

centering

storage

zm_sun

float

Msun

center

cycle

..versionadded:: 17.10.0

xm_sun

Zone mass.

name

type

unit

centering

storage

xm_sun

float

Msun

center

cycle

..versionadded:: 17.10.0

ym_sun

Outer mass coordinate.

name

type

unit

centering

storage

ym_sun

float

Msun

center

cycle

..versionadded:: 17.10.0

Miscellaneous

gamed

Eddington Gamma.

name

type

unit

centering

storage

gamed

float

-

center

loadbuf

mu

Mean molecular weight.

name

type

unit

centering

storage

mu

float

g/mol

center

loadbuf

mue

Mean molecular weight per electron.

name

type

unit

centering

storage

mue

float

g/mol

center

loadbuf

dr

Zone thickness.

name

type

unit

centering

storage

dr

float

cm

center

loadbuf

Added In Version 17.11.7.

tau

Optical depth.

name

type

unit

centering

storage

tau

float

-

center

loadbuf

dtau

Optical thickness of zone.

name

type

unit

centering

storage

dtau

float

-

center

loadbuf

Added In Version 17.0.13.

etot

Specific binding energy of zone.

name

type

unit

centering

storage

etot

float

erg/g

center

loadbuf

Added In Version 18.15.6.

xbind

Binding energy of zone.

name

type

unit

centering

storage

xbind

float

erg

center

loadbuf

ybind

Exterior binding energy.

name

type

unit

centering

storage

ybind

float

erg

center

loadbuf

zbind

Interior binding energy.

name

type

unit

centering

storage

zbind

float

erg

center

loadbuf

uesc

Local escape velocity w/r layers below.

name

type

unit

centering

storage

uesc

float

cm/sec

center

loadbuf

coldepth

Column depth.

name

type

unit

centering

storage

coldepth

float

g/cm**2

center

loadbuf

visc

Kinematic viscosity.

name

type

unit

centering

storage

visc

float

cm**2/sec

center

loadbuf

dynvisc

Dynamic viscosity.

name

type

unit

centering

storage

dynvisc

float

g/cm/sec

center

loadbuf

gr

Enhancement factor of gravity due to GR.

name

type

unit

centering

storage

gr

float

-

center

loadbuf

y

Column depth.

name

type

unit

centering

storage

y

float

g/cm**2

boundary

loadbuf

logy

$log$ ( column depth / g/cm**2 ).

name

type

unit

centering

storage

logy

float

-

boundary

loadbuf

rhor3

$\rho\,r^2$ .

name

type

unit

centering

storage

rhor3

float

g

center

loadbuf

hp

Pressure scale height.

name

type

unit

centering

storage

hp

float

cm

boundary

loadbuf

vconv

Convective velocity.

name

type

unit

centering

storage

vconv

float

cm/sec

boundary

loadbuf

rossby

Convective Rossby Number.

name

type

unit

centering

storage

rossby

float

-

boundary

loadbuf

difieff

Effective convective diffusion coefficient.

name

type

unit

centering

storage

difieff

float

cm**2/sec

boundary

loadbuf

angdgeff

Effective rotational diffusion coefficient.

name

type

unit

centering

storage

angdgeff

float

cm**2/sec

boundary

loadbuf

bfdiffef

Effective magnetic diffusion coefficient.

name

type

unit

centering

storage

bfdiffef

float

cm**2/sec

boundary

loadbuf

gam

Adiabatic exponent, $d\,\ln\,P/d\,\ln\,\rho$ .

name

type

unit

centering

storage

gam

float

-

boundary

loadbuf

ciso

Isothermal sound speed.

name

type

unit

centering

storage

ciso

float

cm/sec

center

cycle

cad

Adiabatic sound speed.

name

type

unit

centering

storage

cad

float

cm/sec

center

cycle

n2

Square of brunt-vaisala frequency.

name

type

unit

centering

storage

n2

integer

1/sec**2

boundary

loadbuf

gbuoy

Buoyancy effective acceleration.

name

type

unit

centering

storage

gbuoy

float

cm/sec**2

boundary

loadbuf

Thermodynmaic Derivatives

stotd

Old specific entropy.

name

type

unit

centering

storage

stotd

float

erg/g/K

center

loadbuf

gamma1

Adiabatic exponent $\Gamma_{\!1}$ .

name

type

unit

centering

storage

gamma1

float

-

center

cycle

gamma2

Adiabatic exponent $\Gamma_{\!2}$ .

name

type

unit

centering

storage

gamma2

float

-

center

cycle

gamma3

Adiabatic exponent $\Gamma_{\!3}$ .

name

type

unit

centering

storage

gamma3

float

-

center

loadbuf

Computed from gamma1 and gamma2 using $1 + \Gamma_{\!1}\,(\Gamma_{\!2}-1)/\Gamma_{\!2}$

enbtn

( $\partial$ en / $\partial$ tn ):math:_{mathrm{dn}}.

name

type

unit

centering

storage

enbtn

float

erg/g/K

center

cycle

enbdn

( $\partial$ en / $\partial$ dn ):math:_{mathrm{tn}}.

name

type

unit

centering

storage

enbdn

float

erg*cc/g**2

center

cycle

pnbtn

( $\partial$ pn / $\partial$ tn ):math:_{mathrm{dn}}.

name

type

unit

centering

storage

pnbtn

float

dyn/K

center

cycle

pnbdn

( $\partial$ pn / $\partial$ dn ):math:_{mathrm{tn}}.

name

type

unit

centering

storage

pnbdn

float

cm**2/sec**2

center

cycle

h

Specific enthalpy.

name

type

unit

centering

storage

h

float

erg/g

center

loadbuf

htot

Enthalpy.

name

type

unit

centering

storage

htot

float

erg

center

loadbuf

alpha

Isotheremal compressibility.

name

type

unit

centering

storage

alpha

float

-

center

loadbuf

delta

Thermal expansion coefficient.

name

type

unit

centering

storage

delta

float

-

center

loadbuf

cp

Specific heat capacity at constant pressure.

name

type

unit

centering

storage

cp

float

erg/g/K

center

loadbuf

Energy Generation

sneut

Neutrino loss rate where a positive value implies net loss.

name

type

unit

centering

storage

sneut

float

erg/g/sec

center

loadbuf

sgain

Maximum of the net nuclear energy generation rate and zero.

name

type

unit

centering

storage

sgain

float

erg/g/sec

center

loadbuf

sloss

Total energy ‘loss’ rate, incl. neutrino loss and net nuclear loss (if any).

name

type

unit

centering

storage

sloss

float

erg/g/sec

center

loadbuf

eg

Specific gravitational energy generation rate.

name

type

unit

centering

storage

eg

float

erg/g/sec

center

loadbuf

snt

$\partial$ sn / $\partial$ $T$ ) $_\mathrm{dn}$ .

name

type

unit

centering

storage

snt

float

erg/g/sec/K

center

loadbuf

snd

$\partial$ sn / $\partial$ $\rho$ ) $_\mathrm{tn}$ .

name

type

unit

centering

storage

snd

float

cm**5/g/sec**3

center

loadbuf

sneutbt

$\partial$ sneut / $\partial$ $T$ ) $_\mathrm{dn}$ .

name

type

unit

centering

storage

sneutbt

float

erg/g/sec/K

center

loadbuf

sneutbd

$\partial$ sneut / $\partial$ $\rho$ ) $_\mathrm{tn}$ .

name

type

unit

centering

storage

sneutbd

float

cm**5/g/sec**3

center

loadbuf

snuc

Nuclear energy generation rate.

name

type

unit

centering

storage

snuc

float

erg/g/sec

center

loadbuf

snucbt

$\partial$ snuc / $\partial$ $T$ ) $_\mathrm{dn}$ .

name

type

unit

centering

storage

snucbt

float

erg/g/sec/K

center

loadbuf

snucbd

$\partial$ snuc / $\partial$ $\rho$ ) $_\mathrm{tn}$ .

name

type

unit

centering

storage

snucbd

float

cm**5/g/sec**3

center

loadbuf

snuw

Energy loss/gain due to weak processes.

name

type

unit

centering

storage

snuw

float

erg/g/sec

center

loadbuf

snuwbt

$\partial$ snuw / $\partial$ $T$ ) $_\mathrm{dn}$ .

name

type

unit

centering

storage

snuwbt

float

erg/g/sec/K

center

loadbuf

snuwbd

$\partial$ snuw / $\partial$ $\rho$ ) $_\mathrm{tn}$ .

name

type

unit

centering

storage

snuwbd

float

cm**5/g/sec**3

center

loadbuf

snubps

.total “pair” neutrino energy loss rate.

name

type

unit

centering

storage

snubps

float

erg/g/sec

center

loadbuf

snubpsbt

( $\partial$ snubps / $\partial$ $T$ ) $_\mathrm{dn}$ .

name

type

unit

centering

storage

snubpsbt

float

erg/g/sec/K

center

loadbuf

snubpsbd

( $\partial$ snubps / $\partial$ $\rho$ ) $_\mathrm{tn}$ .

name

type

unit

centering

storage

snubpsbd

float

cm**5/g/sec**3

center

loadbuf

snlt

( $\partial$ ln sn / $\partial$ ln $T$ ) $_\mathrm{dn}$ .

name

type

unit

centering

storage

snlt

float

-

center

loadbuf

snld

( $\partial$ ln sn / $\partial$ ln $\rho$ ) $_\mathrm{tn}$ .

name

type

unit

centering

storage

snld

float

-

center

loadbuf

sneutlt

( $\partial$ ln sneut / $\partial$ ln $T$ ) $_\mathrm{dn}$ .

name

type

unit

centering

storage

sneutlt

float

-

center

loadbuf

sneutld

( $\partial$ ln sneut / $\partial$ ln $\rho$ ) $_\mathrm{tn}$ .

name

type

unit

centering

storage

sneutld

float

-

center

loadbuf

snuclt

( $\partial$ ln snuc / $\partial$ ln $T$ ) $_\mathrm{dn}$ .

name

type

unit

centering

storage

snuclt

float

-

center

loadbuf

snucld

( $\partial$ ln snuc / $\partial$ ln $\rho$ ) $_\mathrm{tn}$ .

name

type

unit

centering

storage

snucld

float

-

center

loadbuf

snuwlt

( $\partial$ ln snuw / $\partial$ ln $T$ ) $_\mathrm{dn}$ .

name

type

unit

centering

storage

snuwlt

float

-

center

loadbuf

snuwld

( $\partial$ ln snuw / $\partial$ ln $\rho$ ) $_\mathrm{tn}$ .

name

type

unit

centering

storage

snuwld

float

-

center

loadbuf

snubpslt

( $\partial$ ln snubps / $\partial$ ln $T$ ) $_\mathrm{dn}$ .

name

type

unit

centering

storage

snubpslt

float

-

center

loadbuf

snubpsld

( $\partial$ ln snubps / $\partial$ ln $\rho$ ) $_\mathrm{tn}$ .

name

type

unit

centering

storage

snubpsld

float

-

center

loadbuf

Magnetic Quantities

magvcc

arsinh $($ mag. VC criterion $\times$ 1.17 $)$ according to Henk Spruit.

name

type

unit

centering

storage

magvcc

float

-

center

loadbuf

Note

Ask Henk Spruit for details.

eb

Specific magnetic energy density.

name

type

unit

centering

storage

eb

float

erg/g

center

loadbuf

ebv

Magnetic energy density per volume.

name

type

unit

centering

storage

ebv

float

erg/cc

center

loadbuf

bfdiff

Diffision coefficient due to magnetic fields, $D_B$ .

name

type

unit

centering

storage

bfdiff

float

cm**2/sec

boundary

dump

Added In Version 16.11.

Main entry:
Saved Array of Rotational mixing coefficients

bfvisc

Viscosity due to magnetic fields, $D_B$ .

name

type

unit

centering

storage

bfvisc

float

cm**2/sec

boundary

dump

Added In Version 16.11.

Main entry:
Saved Array of Rotational mixing coefficients

bfbr

Radial component of magnetic field, $B_\mathrm{r}$ .

name

type

unit

centering

storage

bfbr

float

G

boundary

loadbuf

Added In Version 16.11.

bfbt

Toroidal component of magnetic field, $B_\mathrm{t}$ .

name

type

unit

centering

storage

bfbt

float

G

boundary

loadbuf

Added In Version 16.11.

The following magnetic quantities are only available if magnet (p 423) $=$ `-1` or `-2`

magaid

Old $I_\mathrm{tot}$ .

name

type

unit

centering

storage

magaid

float

g*cm**2

center

loadbuf

magai

New $I_\mathrm{tot}$ .

name

type

unit

centering

storage

magai

float

g*cm**2

center

loadbuf

magdidt

$\mathrm{d}I/\mathrm{d}t$ .

name

type

unit

centering

storage

magdidt

float

g*cm**2/sec

center

loadbuf

magtau

$\mathrm{d} t /\mathrm{d \ln I}$ .

name

type

unit

centering

storage

magtau

float

sec

center

loadbuf

magb2

$B^2$ .

name

type

unit

centering

storage

magb2

float

g/cm/sec**2

center

loadbuf

mags

Stress $S$ .

name

type

unit

centering

storage

mags

float

g/cm/sec**2

center

loadbuf

magqw

$\mathrm{d}\Omega/d \ln r$ .

name

type

unit

centering

storage

magqw

float

rad/sec

center

loadbuf

magnu

Effective viscosity $\nu$ .

name

type

unit

centering

storage

magnu

float

cm**2/sec

center

loadbuf

The following magnetic quantities are only available if magnet (p 423) $\ge$ `0`

See [Spr02] for definitions of the magnetic quantities.

magn

$N^2$ .

name

type

unit

centering

storage

magn

float

1/sec**2

boundary

loadbuf

magnmu

$N_\mu^2$ .

name

type

unit

centering

storage

magnmu

float

1/sec**2

boundary

loadbuf

magnt

$N_T^2$ .

name

type

unit

centering

storage

magnt

float

1/sec**2

boundary

loadbuf

magcp

Specific heat capacity $c_P$ .

name

type

unit

centering

storage

magcp

float

erg/g/K

boundary

loadbuf

magcap

Effective radiative opacity.

name

type

unit

centering

storage

magcap

float

cm**2/g

boundary

loadbuf

mageta

Viscosity $\eta$ ?.

name

type

unit

centering

storage

mageta

float

?

boundary

loadbuf

magdel

” $\delta$ ”.

name

type

unit

centering

storage

magdel

float

-

boundary

loadbuf

magalp

” $\alpha$ ”.

name

type

unit

centering

storage

magalp

float

-

boundary

loadbuf

maglam

” $\Lambda$ ” - logarithm of Coulmb length cutoff?.

name

type

unit

centering

storage

maglam

float

-

boundary

loadbuf

magq

Shear rate.

name

type

unit

centering

storage

magq

float

-

boundary

loadbuf

magwa0

$\omega_\mathrm{A}$ (0).

name

type

unit

centering

storage

magwa0

float

rad/sec

boundary

loadbuf

magq0

Critical shear rate (0).

name

type

unit

centering

storage

magq0

float

-

boundary

loadbuf

mags0

Magentic stress (0).

name

type

unit

centering

storage

mags0

float

?

boundary

loadbuf

magbphi0

$B_\phi$ (0).

name

type

unit

centering

storage

magbphi0

float

cgs

boundary

loadbuf

magbr0

$B_r$ (0).

name

type

unit

centering

storage

magbr0

float

cgs

boundary

loadbuf

mageta0

Magnetic “viscousity” (0).

name

type

unit

centering

storage

mageta0

float

?

boundary

loadbuf

magwa1

$\omega_\mathrm{A}$ (1).

name

type

unit

centering

storage

magwa1

float

rad/sec

boundary

loadbuf

magq1

Critical shear rate (1).

name

type

unit

centering

storage

magq1

float

-

boundary

loadbuf

mags1a

Magentic stress (1a).

name

type

unit

centering

storage

mags1a

float

?

boundary

loadbuf

mags1b

Magentic stress (1b).

name

type

unit

centering

storage

mags1b

float

?

boundary

loadbuf

mags1

Magentic stress (1).

name

type

unit

centering

storage

mags1

float

?

boundary

loadbuf

magbphi1

$B_\phi$ (1).

name

type

unit

centering

storage

magbphi1

float

cgs

boundary

loadbuf

magbr1

$B_r$ (1).

name

type

unit

centering

storage

magbr1

float

cgs

boundary

loadbuf

mageta1a

Magnetic “viscousity” (1a).

name

type

unit

centering

storage

mageta1a

float

?

boundary

loadbuf

mageta1b

Magnetic “viscousity” (1b).

name

type

unit

centering

storage

mageta1b

float

?

boundary

loadbuf

mageta1

Magnetic “viscousity” (1).

name

type

unit

centering

storage

mageta1

float

?

boundary

loadbuf

magsm

Minimum stress.

name

type

unit

centering

storage

magsm

float

?

boundary

loadbuf

magqm

Maximum of minimum shear values.

name

type

unit

centering

storage

magqm

float

-

boundary

loadbuf

magsf

$(1-q_\mathrm{min}/q)$ .

name

type

unit

centering

storage

magsf

float

-

boundary

loadbuf

magnu

Magnetic viscousity.

name

type

unit

centering

storage

magnu

float

?

boundary

loadbuf

magdiff

Magnetically-induced componennt of diffusion.

name

type

unit

centering

storage

magdiff

float

cm**2/sec

boundary

loadbuf

magvisc

Magnetically-induced componennt of viscosity.

name

type

unit

centering

storage

magvisc

float

cm**2/sec

boundary

loadbuf

magtmix

Magnetically-induced mixing time-scale.

name

type

unit

centering

storage

magtmix

float

sec

boundary

loadbuf

magtrot

Magnetically-induced spin-down time-scale.

name

type

unit

centering

storage

magtrot

float

sec

boundary

loadbuf

mags

Effective magnetic stress.

name

type

unit

centering

storage

mags

float

?

boundary

loadbuf

magbphi

$B_\phi$ (0).

name

type

unit

centering

storage

magbphi

float

cgs

boundary

loadbuf

magbr

$B_r$ .

name

type

unit

centering

storage

magbr

float

cgs

boundary

loadbuf

Edit Quantities by Group

KEPLER Saved Zonal Interface Variables

Saved Array of Rotational mixing coefficients

Saved Arrays Evaluated at Zone Centers

Saved (Zone-Centered) Nuclear Burn Derivative Arrays (APPROX and ISE)

Saved Array of Ion Abundances (Generally Elemental)

Saved Zonal BURN Co-Processor Arrays

Zonal ‘Elemental’ Mass Fractions And Related Variables

APPROX ion list

ISE and NSE Networks ion list

Ise Edit Quantities

Zonal Nuclear Reaction Rates

Zonal BURN Isotope Mass Fractions

Editcom Zonal Variables

Derivative editcom Variables

EOS editcom Variables:

ISE EDITCOM Variables:

loadbuf Edit Quantities

BURN Related Edit Quantities

Miscellaneous

Mixing and Transport

Rotation

Mass Coordinate

Miscellaneous

Thermodynmaic Derivatives

Energy Generation

Magnetic Quantities

The following magnetic quantities are only available if magnet (p 423) -1 or -2

The following magnetic quantities are only available if magnet (p 423) 0

Derivative `editcom` Variables

EOS `editcom` Variables:

`loadbuf` Edit Quantities

The following magnetic quantities are only available if magnet (p 423) $=$ `-1` or `-2`

The following magnetic quantities are only available if magnet (p 423) $\ge$ `0`