eosparm in BC stuff

AMReX-Combustion · Sep 17, 2024 · a50c552 · a50c552
1 parent 0b22bbf
commit a50c552
Show file tree

Hide file tree

Showing 3 changed files with 6 additions and 2 deletions.
diff --git a/Exec/RegTests/FlameSheet/pelelmex_prob.H b/Exec/RegTests/FlameSheet/pelelmex_prob.H
@@ -38,7 +38,7 @@ pelelmex_initdata(
 
   constexpr amrex::Real Pi = 3.14159265358979323846264338327950288;
 
-  auto eos = pele::physics::PhysicsType::eos();
+  auto eos = pele::physics::PhysicsType::eos(prob_parm.eosparm);
   amrex::GpuArray<amrex::Real, NUM_SPECIES + 4> pmf_vals = {0.0};
   amrex::Real molefrac[NUM_SPECIES] = {0.0};
   amrex::Real massfrac[NUM_SPECIES] = {0.0};
@@ -140,7 +140,7 @@ bcnormal(
   amrex::Real molefrac[NUM_SPECIES] = {0.0};
   amrex::Real massfrac[NUM_SPECIES] = {0.0};
 
-  auto eos = pele::physics::PhysicsType::eos();
+  auto eos = pele::physics::PhysicsType::eos(prob_parm.eosparm);
   if (sgn == 1) {
     pele::physics::PMF::pmf(pmf_data, prob_lo[idir], prob_lo[idir], pmf_vals);
 

diff --git a/Exec/RegTests/FlameSheet/pelelmex_prob.cpp b/Exec/RegTests/FlameSheet/pelelmex_prob.cpp
@@ -12,5 +12,6 @@ PeleLM::readProbParm()
   pp.query("pertmag", PeleLM::prob_parm->pertmag);
   pp.query("pertlength", PeleLM::prob_parm->pertlength);
 
+  PeleLM::prob_parm->eosparm = PeleLM::eos_parms.device_parm();
   PeleLM::pmf_data.initialize();
 }
diff --git a/Exec/RegTests/FlameSheet/pelelmex_prob_parm.H b/Exec/RegTests/FlameSheet/pelelmex_prob_parm.H
@@ -3,6 +3,7 @@
 
 #include <AMReX_REAL.H>
 #include <AMReX_GpuMemory.H>
+#include <PelePhysics.H>
 
 using namespace amrex::literals;
 
@@ -13,5 +14,7 @@ struct ProbParm
   amrex::Real pertmag = 0.0004_rt;
   amrex::Real pertlength = 0.008_rt;
   int meanFlowDir = 1;
+  pele::physics::eos::EosParm<pele::physics::PhysicsType::eos_type> const*
+    eosparm;
 };
 #endif