diff --git a/Source/PeleLMeX.H b/Source/PeleLMeX.H index 8f07e385..975d352e 100644 --- a/Source/PeleLMeX.H +++ b/Source/PeleLMeX.H @@ -316,6 +316,7 @@ public: * parameters. */ void checkRunParams(); + void checkSetupParams(); //----------------------------------------------------------------------------- // REGRID / LOAD BALANCE diff --git a/Source/PeleLMeX_Setup.cpp b/Source/PeleLMeX_Setup.cpp index 292030ce..d77a4dc5 100644 --- a/Source/PeleLMeX_Setup.cpp +++ b/Source/PeleLMeX_Setup.cpp @@ -209,6 +209,9 @@ PeleLM::Setup() // Initialize active control initActiveControl(); + + // Check setup parameters + checkSetupParams(); } void @@ -748,6 +751,18 @@ PeleLM::readParameters() } #endif + // ----------------------------------------- + // External Sources + // ----------------------------------------- + m_user_defined_ext_sources = false; + m_ext_sources_SDC = false; // TODO: add capability to update ext_srcs in SDC + pp.query("user_defined_ext_sources", m_user_defined_ext_sources); +} + +void +PeleLM::checkSetupParams() +{ + BL_PROFILE("PeleLMeX::checkSetupParams()"); // Ensure unsupported physics is not used with manifiold models if (pele::physics::PhysicsType::eos_type::identifier() == "Manifold") { if (m_closed_chamber != 0) { @@ -772,15 +787,18 @@ PeleLM::readParameters() #endif #ifdef PELE_USE_EFIELD amrex::Abort("Efield models are not yet supported for Manifold EOS"); +#endif +#ifdef USE_MANIFOLD_EOS + if ( + std::abs( + (0.1 * eos_parms.host_parm().Pnom_cgs - prob_parm->P_mean) / + prob_parm->P_mean) > 1e-6) { + amrex::Abort("For Manifold EOS, pressure in manifold model " + "(manifold.nominal_pressure_cgs) and pressure in PeleLMeX " + "(prob.Pmean) must match"); + } #endif } - - // ----------------------------------------- - // External Sources - // ----------------------------------------- - m_user_defined_ext_sources = false; - m_ext_sources_SDC = false; // TODO: add capability to update ext_srcs in SDC - pp.query("user_defined_ext_sources", m_user_defined_ext_sources); } void