only do new gridding if max_grid_size not set (erf-model#1854)

ewquon · Oct 7, 2024 · 444a187 · 444a187
1 parent 0f89138
commit 444a187
Showing 1 changed file with 6 additions and 2 deletions.
diff --git a/Source/ERF_make_new_level.cpp b/Source/ERF_make_new_level.cpp
@@ -25,10 +25,14 @@ void ERF::MakeNewLevelFromScratch (int lev, Real time, const BoxArray& ba_in,
 {
     BoxArray ba;
     DistributionMapping dm;
-    if (lev == 0 && ba_in.size() != ParallelDescriptor::NProcs())
+    Box domain(Geom(0).Domain());
+    if (lev == 0 &&
+        (max_grid_size[0][0] >= domain.length(0)) &&
+        (max_grid_size[0][1] >= domain.length(1)) &&
+        ba_in.size() != ParallelDescriptor::NProcs())
     {
         // We only decompose in z if max_grid_size_z indicates we should
-        bool decompose_in_z = (max_grid_size[0][2] < Geom(0).Domain().length(2));
+        bool decompose_in_z = (max_grid_size[0][2] < domain.length(2));
 
         ba = ERFPostProcessBaseGrids(Geom(0).Domain(),decompose_in_z);
         dm = DistributionMapping(ba);