From 6053b3e6cc01a4e2341f336aab6a85a8830a2a2d Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Fri, 26 Apr 2024 15:29:34 -0400 Subject: [PATCH] fix a check on whether the outer boundary is open in gravity we were testing physbc_lo[dim] != amrex::PhysBCType::interior when it should be physbc_hi[dim] --- Source/driver/Castro.cpp | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/Source/driver/Castro.cpp b/Source/driver/Castro.cpp index 702385e75c..4ceac43446 100644 --- a/Source/driver/Castro.cpp +++ b/Source/driver/Castro.cpp @@ -3551,12 +3551,14 @@ Castro::apply_tagging_restrictions(TagBoxArray& tags, [[maybe_unused]] Real time int boundary_buf = n_error_buf[dim] + blocking_factor[dim] / ref_ratio[dim]; - if ((physbc_lo[dim] != amrex::PhysBCType::symmetry && physbc_lo[dim] != amrex::PhysBCType::interior) && + if ((physbc_lo[dim] != amrex::PhysBCType::symmetry && + physbc_lo[dim] != amrex::PhysBCType::interior) && (idx[dim] <= domlo[dim] + boundary_buf)) { outer_boundary_test[dim] = true; } - if ((physbc_hi[dim] != amrex::PhysBCType::symmetry && physbc_lo[dim] != amrex::PhysBCType::interior) && + if ((physbc_hi[dim] != amrex::PhysBCType::symmetry && + physbc_hi[dim] != amrex::PhysBCType::interior) && (idx[dim] >= domhi[dim] - boundary_buf)) { outer_boundary_test[dim] = true; }