Skip to content

Commit

Permalink
fixes for ReactornNull with eosparm
Browse files Browse the repository at this point in the history
  • Loading branch information
baperry2 committed Jul 19, 2024
1 parent d3e2b87 commit 2a6ea71
Show file tree
Hide file tree
Showing 2 changed files with 8 additions and 1 deletion.
7 changes: 7 additions & 0 deletions Source/Reactions/ReactorNull.H
Original file line number Diff line number Diff line change
Expand Up @@ -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};
Expand Down
2 changes: 1 addition & 1 deletion Source/Reactions/ReactorNull.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -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) {
Expand Down

0 comments on commit 2a6ea71

Please sign in to comment.