From bd33b5c50c097da562fefaec4d9225b3f361c9a6 Mon Sep 17 00:00:00 2001
From: Michael Zingale <michael.zingale@stonybrook.edu>
Date: Sat, 4 Nov 2023 14:16:26 -0400
Subject: [PATCH] fix comp

---
 integration/nse_update_strang.H     |  5 +++--
 unit_test/nse_table_cell/nse_cell.H | 13 +++++++++----
 2 files changed, 12 insertions(+), 6 deletions(-)

diff --git a/integration/nse_update_strang.H b/integration/nse_update_strang.H
index d6336896b3..4bbbc5d8e4 100644
--- 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
index 6d93cc21ae..88b21c1b02 100644
--- 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;
     }