Fix for when a new level is created

cgilet · Apr 15, 2024 · b6aa3be · b6aa3be
1 parent b1c00f7
commit b6aa3be
Showing 1 changed file with 3 additions and 4 deletions.
diff --git a/src/incflo_regrid.cpp b/src/incflo_regrid.cpp
@@ -43,18 +43,17 @@ void incflo::MakeNewLevelFromCoarse (int lev,
     m_leveldata[lev] = std::move(new_leveldata);
     m_factory[lev] = std::move(new_fact);
 
-    //make_mixedBC_mask(lev, ba, dm);
-
     m_diffusion_tensor_op.reset();
     m_diffusion_scalar_op.reset();
 
+    // Note: finest_level has not yet been updated and so we use lev
 #ifdef AMREX_USE_EB
-    macproj = std::make_unique<Hydro::MacProjector>(Geom(0,finest_level),
+    macproj = std::make_unique<Hydro::MacProjector>(Geom(0,lev),
                       MLMG::Location::FaceCentroid,  // Location of mac_vec
                       MLMG::Location::FaceCentroid,  // Location of beta
                       MLMG::Location::CellCenter  ); // Location of solution variable phi
 #else
-    macproj = std::make_unique<Hydro::MacProjector>(Geom(0,finest_level));
+    macproj = std::make_unique<Hydro::MacProjector>(Geom(0,lev));
 #endif
 }