From 46acb854a93a1696cb41bf5e64fb76ff717908e6 Mon Sep 17 00:00:00 2001
From: Weiqun Zhang <weiqunzhang@lbl.gov>
Date: Mon, 27 May 2024 07:51:04 -0700
Subject: [PATCH] Fix bug introduced in #3959

---
 Src/Boundary/AMReX_InterpBndryData.H | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/Src/Boundary/AMReX_InterpBndryData.H b/Src/Boundary/AMReX_InterpBndryData.H
index 0a561fae71b..ed611738c83 100644
--- a/Src/Boundary/AMReX_InterpBndryData.H
+++ b/Src/Boundary/AMReX_InterpBndryData.H
@@ -171,7 +171,7 @@ InterpBndryDataT<MF>::setBndryValues (BndryRegisterT<MF> const& crse, int c_star
                                       int f_start, int bnd_start, int num_comp, const IntVect& ratio,
                                       int max_order, int max_width)
 {
-    AMREX_ASSERT(this->grids == fine.boxArray());
+    AMREX_ASSERT(fine.empty() || this->grids == fine.boxArray());
 
     const Box& fine_domain = this->geom.Domain();
 
@@ -182,7 +182,8 @@ InterpBndryDataT<MF>::setBndryValues (BndryRegisterT<MF> const& crse, int c_star
 #ifdef AMREX_USE_OMP
 #pragma omp parallel if (Gpu::notInLaunchRegion())
 #endif
-        for (MFIter mfi(fine,info); mfi.isValid(); ++mfi)
+        for (MFIter mfi(this->grids, this->bndry[0].DistributionMap(), info);
+             mfi.isValid(); ++mfi)
         {
             const Box& fine_bx  = mfi.validbox();
             for (OrientationIter fi; fi; ++fi) {