Post regrid clean

Source/PeleLMBC.cpp

@@ -96,6 +96,12 @@ InterpBase* PeleLM::getInterpolator(int a_method) {
                                            : &eb_mf_cell_cons_interp;
 #else
     mapper = &mf_cell_cons_interp;
+#endif
+  } else if (a_method == 2) {
+#ifdef AMREX_USE_EB
+    Abort("Regrid interpolation method = 2 not available with EB !");
+#else
+    mapper = &mf_linear_slope_minmax_interp;
 #endif
   } else {
     Abort("Unknown interpolation method");

Source/PeleLMRegrid.cpp

@@ -81,6 +81,9 @@ void PeleLM::MakeNewLevelFromCoarse( int lev,
          setRhoToSumRhoY(lev, AmrNewTime);
       }
 
+      // Recompute temperature from interpolated rhoYs and rhoH
+      setTemperature(lev, AmrNewTime);
+
       // Initialize thermodynamic pressure
       setThermoPress(lev, AmrNewTime);
    }
@@ -176,6 +179,15 @@ void PeleLM::RemakeLevel( int lev,
    m_resetCoveredMask = 1;
 
    if (!m_incompressible) {
+      // Enforce density / species density consistency
+      // only usefull when using cell cons interp
+      if (m_regrid_interp_method == 1) {
+         setRhoToSumRhoY(lev, AmrNewTime);
+      }
+
+      // Recompute temperature from interpolated rhoYs and rhoH
+      setTemperature(lev, AmrNewTime);
+
       // Initialize thermodynamic pressure
       setThermoPress(lev, AmrNewTime);
    }