From 444a1877404d9355633edc28978f9f4ba6ee159a Mon Sep 17 00:00:00 2001 From: Ann Almgren Date: Mon, 7 Oct 2024 09:33:19 -0700 Subject: [PATCH] only do new gridding if max_grid_size not set (#1854) --- Source/ERF_make_new_level.cpp | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) diff --git a/Source/ERF_make_new_level.cpp b/Source/ERF_make_new_level.cpp index 76174516a..2de79276f 100644 --- 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);