fix comp

AMReX-Astro · Nov 4, 2023 · bd33b5c · bd33b5c
1 parent e1ff5fc
commit bd33b5c
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 6 deletions.
diff --git a/integration/nse_update_strang.H b/integration/nse_update_strang.H
@@ -41,6 +41,7 @@ void nse_burn(BurnT& state, const Real dt) {
   Real dq_out;
   Real dyedt;
   Real dabardt;
+  Real dbeadt;
   Real enu;
   Real X[NumSpec];
 
@@ -58,7 +59,7 @@ void nse_burn(BurnT& state, const Real dt) {
   // in terms of the auxiliary composition, Ye, abar, and B/A
 
   nse_interp(T_in, state.rho, state.aux[iye],
-             abar_out, dq_out, dyedt, dabardt, enu, X);
+             abar_out, dq_out, dyedt, dabardt, dbeadt, enu, X);
 
   // update Ye
 
@@ -67,7 +68,7 @@ void nse_burn(BurnT& state, const Real dt) {
   // now get the composition from the table using the updated Ye
 
   nse_interp(T_in, state.rho, state.aux[iye],
-             abar_out, dq_out, dyedt, dabardt, enu, X);
+             abar_out, dq_out, dyedt, dabardt, dbeadt, enu, X);
 
 
   // this is MeV / nucleon -- here aux has not yet been updated, so we

diff --git a/unit_test/nse_table_cell/nse_cell.H b/unit_test/nse_table_cell/nse_cell.H
@@ -19,14 +19,19 @@ void nse_cell_c()
     Real abar{};
     Real dq{};
     Real dyedt{};
+    Real dabardt{};
+    Real dbeadt{};
+    Real e_nu{};
     Real X[NumSpec] = {};
 
-    nse_interp(density, temperature, ye,
-               abar, dq, dyedt, X);
+    nse_interp(temperature, density, ye,
+               abar, dq, dyedt, dabardt, dbeadt, e_nu, X);
 
-    std::cout << "abar = " << abar << std::endl;
-    std::cout << "B/A = " << dq << std::endl;
+    std::cout << "Abar = " << abar << std::endl;
+    std::cout << "<B/A> = " << dq << std::endl;
     std::cout << "dYe/dt = " << dyedt << std::endl;
+    std::cout << "dAbar/dt = " << dabardt << std::endl;
+    std::cout << "d<B/A>/dt = " << dbeadt << std::endl;
     for (int n = 0; n < NumSpec; ++n) {
         std::cout << "X(" << short_spec_names_cxx[n] << ") = " << X[n] << std::endl;
     }