From c6123dee5fb164e78bf522d4bcf473e078bcc619 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Thu, 9 Nov 2023 07:26:09 -0500 Subject: [PATCH] more updates --- nse_tabular/README.md | 13 ++++++++++++- nse_tabular/nse_table_data.cpp | 2 +- 2 files changed, 13 insertions(+), 2 deletions(-) diff --git a/nse_tabular/README.md b/nse_tabular/README.md index fa9d2d609b..ffbc002f39 100644 --- a/nse_tabular/README.md +++ b/nse_tabular/README.md @@ -9,7 +9,7 @@ elsewhere. This requires compiling with ``` -USE_TABULAR_NSE=TRUE +USE_NSE_TABLE=TRUE ``` This will change the equation of state to work in terms of (Y_e, abar, @@ -39,6 +39,8 @@ The table provides: * dAbar/dt : evolution of Abar in 1/s +* d/dt : evolution of in MeV/s + * e_nu : weak rate neutrino loss energy in erg/g/s * X(A), X(B), ... : the reduced composition mass fractions. They are @@ -92,3 +94,12 @@ The data is ordered such that rho varies the slowest (from low to high), T varies the next slowest (from low to high), and Ye varies the fastest (from high to low). +## Outputting for a different network + +At the moment, the script is configured for ``aprox19``. To change it +to output to a different network, a new function needs to be added in +the same form as the ``get_aprox19_comp`` function. The job of that +run is to reduce the composition down to that of the network on the +grid. The main thing that would need to be done is to change the list +of nuclei and update the list ``X[]`` to output them in the proper +order. diff --git a/nse_tabular/nse_table_data.cpp b/nse_tabular/nse_table_data.cpp index 4f357af787..fe25777209 100644 --- a/nse_tabular/nse_table_data.cpp +++ b/nse_tabular/nse_table_data.cpp @@ -10,6 +10,6 @@ namespace nse_table AMREX_GPU_MANAGED amrex::Array1D dbeadttab; AMREX_GPU_MANAGED amrex::Array1D enutab; - AMREX_GPU_MANAGED amrex::Array2D massfractab; + AMREX_GPU_MANAGED amrex::Array2D massfractab; } #endif