fixes for ReactornNull with eosparm

AMReX-Combustion · Jul 19, 2024 · bb1c21e · bb1c21e
1 parent d3e2b87
commit bb1c21e
Show file tree

Hide file tree

Showing 3 changed files with 10 additions and 2 deletions.
diff --git a/Source/Reactions/ReactorBase.H b/Source/Reactions/ReactorBase.H
@@ -89,7 +89,8 @@ public:
   {
     auto eos = pele::physics::PhysicsType::eos(eosparm);
     if (eos.identifier() == "Manifold") {
-      amrex::Abort("Manifold EOS only supported with ReactorRK64 for now");
+      amrex::Abort(
+        "Manifold EOS only supported with ReactorRK64 and ReactorNull for now");
     } else {
       m_eosparm = eosparm;
     }

diff --git a/Source/Reactions/ReactorNull.H b/Source/Reactions/ReactorNull.H
@@ -86,6 +86,13 @@ public:
       rhoE, frcEExt, FC_in, y_vect, vect_energy, FCunt, dt);
   }
 
+  void set_eos_parm(
+    const pele::physics::eos::EosParm<pele::physics::PhysicsType::eos_type>*
+      eosparm) override
+  {
+    m_eosparm = eosparm;
+  }
+
 private:
   utils::FlattenOps<Ordering> flatten_ops;
   int m_reactor_type{0};

diff --git a/Source/Reactions/ReactorNull.cpp b/Source/Reactions/ReactorNull.cpp
@@ -61,7 +61,7 @@ ReactorNull::react(
     }
     amrex::Real energy_loc = renergy_loc / rho_loc;
     amrex::Real T_loc = T_in(i, j, k, 0);
-    auto eos = pele::physics::PhysicsType::eos(m_eosparm);
+    auto eos = pele::physics::PhysicsType::eos(leosparm);
     if (captured_reactor_type == ReactorTypes::e_reactor_type) {
       eos.REY2T(rho_loc, energy_loc, Y_loc, T_loc);
     } else if (captured_reactor_type == ReactorTypes::h_reactor_type) {