CarpetReduce: Handle multi-patch boundaries correctly

Use MultiPatch_GetBoundarySpecification to determine the boundary
specifications if it is available; otherwise, use GetBoundarySpecification.

1 files changed, 9 insertions, 2 deletions

diff --git a/Carpet/CarpetReduce/src/mask_coords.c b/Carpet/CarpetReduce/src/mask_coords.c
index 2b7d706bd..6ea6d52ea 100644
--- a/Carpet/CarpetReduce/src/mask_coords.c
+++ b/Carpet/CarpetReduce/src/mask_coords.c
@@ -24,8 +24,15 @@ CoordBase_SetupMask (CCTK_ARGUMENTS)
   
   
   
-  ierr = GetBoundarySpecification
-    (6, nboundaryzones, is_internal, is_staggered, shiftout);
+  if (CCTK_IsFunctionAliased ("MultiPatch_GetBoundarySpecification")) {
+    int const m = MultiPatch_GetMap (cctkGH);
+    assert (m >= 0);
+    ierr = MultiPatch_GetBoundarySpecification
+      (m, 6, nboundaryzones, is_internal, is_staggered, shiftout);
+  } else {
+    ierr = GetBoundarySpecification
+      (6, nboundaryzones, is_internal, is_staggered, shiftout);
+  }
   if (ierr != 0) {
     CCTK_WARN (0, "Could not get boundary specification");
   }