Skip to content

update CNO_extras #2051

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
zingale wants to merge 8 commits into
base: development
Choose a base branch
from
Open

update CNO_extras #2051

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
6210a9b
update CNO_extras
zingale Aug 28, 2026
9102160
update
zingale Aug 28, 2026
d8c2d97
update
zingale Aug 28, 2026
0f555ff
update
zingale Sep 1, 2026
8553b47
more updates
zingale Sep 2, 2026
fcf37b3
update
zingale Sep 2, 2026
b9ff61b
update
zingale Sep 4, 2026
609a5ad
remove old fig
zingale Sep 4, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 2 additions & 0 deletions Docs/source/networks.rst
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -166,6 +166,8 @@ burning. This network is managed by pynucastro.
.. figure:: ../../networks/CNO_extras/cno_extras_hide_alpha.png
:align: center

Overall there are 22 nuclei and 74 rates explicitly connecting the nuclei.

.. note::

We add ${}^{56}\mathrm{Fe}$ as an inert nucleus to allow this to be
Expand Down
26 changes: 8 additions & 18 deletions networks/CNO_extras/CNO_extras.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3,15 +3,8 @@

def create_network():

net = pyna.network_helper(["h1", "he4",
"c12", "c13",
"n13", "n14", "n15",
"o14", "o15", "o16", "o17", "o18",
"f17", "f18", "f19",
"ne18", "ne19", "ne20",
"mg22", "mg24"],
tabular_ordering=["ffn", "langanke", "oda"],
inert_nuclei=["fe56"], network_type="amrex")
net = pyna.common_networks.cno(network_type="amrex")
net.add_inert_nucleus("fe56")

return net

Expand All @@ -22,20 +15,17 @@ def doit():

net.write_network()

comp = pyna.Composition(net.get_nuclei())
comp.set_solar_like()
comp = pyna.Composition(net.get_nuclei(), init="solar")

rho = 1.e6
rho = 1.e4
T = 1.e8

net.plot(rho, T, comp, outfile="cno_extras.png",
Z_range=[1, 13], N_range=[1, 13])

net.plot(outfile="cno_extras_hide_alpha.png",
Z_range=[1, 13], N_range=[1, 13],
rotated=True,
hide_xalpha=True)
node_size=500, node_font_size="9",
Z_range=[0, 13], N_range=[-0.5, 12.5],
always_show_p=True)

net.summary()

if __name__ == "__main__":
doit()
2 changes: 2 additions & 0 deletions networks/CNO_extras/Make.package
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -10,9 +10,11 @@ ifeq ($(USE_REACT),TRUE)
CEXE_headers += rate_type.H
CEXE_headers += reaclib_rates.H
CEXE_headers += approximate_rates.H
CEXE_headers += branched_rates.H
CEXE_headers += modified_rates.H
CEXE_headers += table_rates.H
CEXE_headers += temperature_table_rates.H
CEXE_headers += tabular_weak_rates.H
CEXE_headers += derived_rates.H
USE_SCREENING = TRUE
USE_NEUTRINOS = TRUE
Expand Down
Loading
Loading