From a0ac0d089bc851ab9fd05a8dc3862813befb33d9 Mon Sep 17 00:00:00 2001
From: Bhargav Sriram Siddani <siddanibhargav@gmail.com>
Date: Fri, 1 Nov 2024 15:31:57 -0700
Subject: [PATCH] Changed the location of query for epsilons of mu(I)

---
 src/rheology/incflo_read_rheology_parameters.cpp | 2 ++
 src/setup/init.cpp                               | 2 --
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/rheology/incflo_read_rheology_parameters.cpp b/src/rheology/incflo_read_rheology_parameters.cpp
index cbe8cd99..566b64e7 100644
--- a/src/rheology/incflo_read_rheology_parameters.cpp
+++ b/src/rheology/incflo_read_rheology_parameters.cpp
@@ -201,6 +201,8 @@ void incflo::ReadRheologyParameters()
         else if(fluid_model_s_snd == "mpmd")
         {
             m_fluid_model_second = FluidModel::DataDrivenMPMD;
+            pp_scnd.query("mu_p_eps_second",m_mu_p_eps_second);
+            pp_scnd.query("mu_sr_eps_second",m_mu_sr_eps_second);
             amrex::Print() << "Data-driven model through AMReX-MPMD."<<std::endl;
         }
 #endif
diff --git a/src/setup/init.cpp b/src/setup/init.cpp
index 27614752..2fd329ac 100644
--- a/src/setup/init.cpp
+++ b/src/setup/init.cpp
@@ -155,8 +155,6 @@ void incflo::ReadParameters ()
         pp.query("two_fluid_rho_harmonic",m_two_fluid_rho_harmonic);
         pp.query("two_fluid_eta_harmonic",m_two_fluid_eta_harmonic);
         pp.query("two_fluid_cc_rho_conc",m_two_fluid_cc_rho_conc);
-        pp.query("mu_p_eps_second",m_mu_p_eps_second);
-        pp.query("mu_sr_eps_second",m_mu_sr_eps_second);
 
         // Density (if constant)
         pp.query("ro_0", m_ro_0);