diff --git a/Exec/RegTests/CMakeLists.txt b/Exec/RegTests/CMakeLists.txt
index 88474eff..85d73eee 100644
--- a/Exec/RegTests/CMakeLists.txt
+++ b/Exec/RegTests/CMakeLists.txt
@@ -21,6 +21,7 @@ else()
   endif()
   if(PELE_DIM EQUAL 2)
     add_subdirectory(TripleFlame)
+    add_subdirectory(IsothermalSoretChannel)
   endif()
   if(PELE_DIM EQUAL 3)
     add_subdirectory(TurbInflow)
diff --git a/Exec/RegTests/IsothermalSoretChannel/input.inp b/Exec/RegTests/IsothermalSoretChannel/isothermal-soret.inp
similarity index 97%
rename from Exec/RegTests/IsothermalSoretChannel/input.inp
rename to Exec/RegTests/IsothermalSoretChannel/isothermal-soret.inp
index d7360016..16a75794 100644
--- a/Exec/RegTests/IsothermalSoretChannel/input.inp
+++ b/Exec/RegTests/IsothermalSoretChannel/isothermal-soret.inp
@@ -25,7 +25,7 @@ prob.P_mean = 101325.0
 prob.Vin = 10.0
 prob.Tlow = 400.0
 prob.Thigh = 800.0
-prob.phi = 0.5
+prob.phi = 1.0
 #---------------------- Transport --------------------------------
 transport.use_soret = 1
 
@@ -46,7 +46,7 @@ peleLM.do_species_balance = 1             # Compute species balance
 amr.max_step = 2000                       # Maximum number of time steps
 amr.stop_time = 4.00                      # final simulation physical time [s]
 amr.max_wall_time = 1                     # Maximum simulation run time [h]
-amr.cfl = 0.5                             # CFL number for hyperbolic system 
+amr.cfl = 0.5                             # CFL number for hyperbolic system
 amr.dt_shrink = 0.01                     # Scale back initial timestep
 amr.dt_change_max = 1.1                   # Maximum dt increase btw successive steps
 
@@ -86,6 +86,6 @@ amr.gradT.adjacent_difference_greater = 200
 amr.gradT.field_name = temp
 
 #---------------------- Debug/HPC CONTROL-------------------------
-amrex.fpe_trap_invalid = 1
-amrex.fpe_trap_zero = 1
-amrex.fpe_trap_overflow = 1
+#amrex.fpe_trap_invalid = 1
+#amrex.fpe_trap_zero = 1
+#amrex.fpe_trap_overflow = 1
diff --git a/Exec/RegTests/IsothermalSoretChannel/pelelmex_prob.H b/Exec/RegTests/IsothermalSoretChannel/pelelmex_prob.H
index 33c12a84..f824e29d 100644
--- a/Exec/RegTests/IsothermalSoretChannel/pelelmex_prob.H
+++ b/Exec/RegTests/IsothermalSoretChannel/pelelmex_prob.H
@@ -59,11 +59,8 @@ pelelmex_initdata(
   const amrex::Real* prob_hi = geomdata.ProbHi();
   const amrex::Real* dx = geomdata.CellSize();
 
-  const amrex::Real x = prob_lo[0] + (i + 0.5) * dx[0];
   const amrex::Real y = prob_lo[1] + (j + 0.5) * dx[1];
   const amrex::Real y_low = prob_lo[1];
-
-  const amrex::Real Lx = prob_hi[0] - prob_lo[0];
   const amrex::Real Ly = prob_hi[1] - prob_lo[1];
 
   auto eos = pele::physics::PhysicsType::eos();
@@ -98,19 +95,16 @@ AMREX_GPU_HOST_DEVICE
 AMREX_FORCE_INLINE
 void
 bcnormal(
-  const amrex::Real x[AMREX_SPACEDIM],
+  const amrex::Real* /*x[AMREX_SPACEDIM]*/,
   const int /*m_nAux*/,
   amrex::Real s_ext[NVAR],
   const int /*idir*/,
   const int sgn,
-  const amrex::Real time,
+  const amrex::Real /*time*/,
   amrex::GeometryData const& /*geomdata*/,
   ProbParm const& prob_parm,
   pele::physics::PMF::PmfData::DataContainer const* /*pmf_data*/)
 {
-  auto eos = pele::physics::PhysicsType::eos();
-  amrex::Real massfrac[NUM_SPECIES] = {0.0};
-
   if (sgn == 1) {
     s_ext[TEMP] = prob_parm.Tlow;
   } else if (sgn == -1) {
diff --git a/Exec/RegTests/IsothermalSoretChannel/pelelmex_prob.cpp b/Exec/RegTests/IsothermalSoretChannel/pelelmex_prob.cpp
index 53c005d5..d454a284 100644
--- a/Exec/RegTests/IsothermalSoretChannel/pelelmex_prob.cpp
+++ b/Exec/RegTests/IsothermalSoretChannel/pelelmex_prob.cpp
@@ -10,6 +10,7 @@ PeleLM::readProbParm()
   pp.query("Vin", PeleLM::prob_parm->Vin);
   pp.query("Thigh", PeleLM::prob_parm->Thigh);
   pp.query("Tlow", PeleLM::prob_parm->Tlow);
+  pp.query("phi", PeleLM::prob_parm->phi);
 
   PeleLM::prob_parm->bathID = N2_ID;
   PeleLM::prob_parm->oxidID = O2_ID;
diff --git a/Tests/CMakeLists.txt b/Tests/CMakeLists.txt
index 21fa5e81..7e42d225 100644
--- a/Tests/CMakeLists.txt
+++ b/Tests/CMakeLists.txt
@@ -236,6 +236,7 @@ if(NOT PELE_ENABLE_EB)
   endif()
   if(PELE_DIM EQUAL 2)
     add_test_r(tripleflame-${PELE_DIM}d TripleFlame)
+    add_test_rv(isothermal-soret IsothermalSoretChannel)
   endif()
   if(PELE_DIM EQUAL 3)
     add_test_r(hit-${PELE_DIM}d HITDecay)