add an adaptive-NSE network with Fe-group nuclei

Docs/source/networks.rst

@@ -410,6 +410,19 @@ The full network appears as:
    :align: center
 
 
+``ase-iron``
+------------
+
+As with ``ase``, this network is constructed to have reverse rates for all forward rates, allowing
+it to be used with the :ref:`self_consistent_nse` solver.  It builds off of ``ase`` by including
+more iron-group nuclei (using the same nuclei as ``he-burn-28amnp``).
+
+The full network appears as:
+
+.. figure:: ../../networks/he-burn/ase-iron/ase-iron.png
+   :align: center
+
+
 ``cno_he_burn_33a``
 -------------------
 

networks/he-burn/ase-iron/Make.package

@@ -0,0 +1,15 @@
+CEXE_headers += network_properties.H
+
+ifeq ($(USE_REACT),TRUE)
+  CEXE_sources += actual_network_data.cpp
+  CEXE_headers += actual_network.H
+  CEXE_headers += tfactors.H
+  CEXE_headers += partition_functions.H
+  CEXE_sources += partition_functions_data.cpp
+  CEXE_headers += actual_rhs.H
+  CEXE_headers += reaclib_rates.H
+  CEXE_headers += table_rates.H
+  CEXE_sources += table_rates_data.cpp
+  USE_SCREENING = TRUE
+  USE_NEUTRINOS = TRUE
+endif

networks/he-burn/ase-iron/README.md

@@ -0,0 +1,19 @@
+# `ase`
+
+This network is specifically tailored for self-consistent NSE evolution.
+This network achieves machine precision accuracy when comparing equilibrium
+composition from direct integration and NSE composition.
+This can be demonstrated using the `nse_compatibility` script.
+
+This is similar to `he-burn-19a` except:
+
+* neutrinos are included.
+
+* N14 is not included.
+
+* Every forward rate has a corresponding reverse rate calculated from
+  detailed balance. This includes the reverse rates for C12+C12,
+  C12+O16, and O16+O16 as they're removed in other `he-burn` nets.
+
+* Q-value is recomputed in `DerivedRate` for better consistency
+  with NSE.

networks/he-burn/ase-iron/_parameters

@@ -0,0 +1,2 @@
+@namespace: network
+