1 files changed, 10 insertions, 0 deletions
diff --git a/src/EHFinder_SetMask.F90 b/src/EHFinder_SetMask.F90
index f521ba2..4be4ac7 100644
--- a/src/EHFinder_SetMask.F90
+++ b/src/EHFinder_SetMask.F90
@@ -426,6 +426,16 @@ subroutine EHFinder_SetMask1(CCTK_ARGUMENTS)
       end where
     end if
 
+!   Now check if any remaining active points where actually excised in the
+!   numerical run generating the metric data.
+
+    if ( use_mask .gt. 0 ) then
+      where ( emask(ixl:ixr,jyl:jyr,kzl:kzr) .lt. quarter )
+        eh_mask(ixl:ixr,jyl:jyr,kzl:kzr) = -1
+        f(ixl:ixr,jyl:jyr,kzl:kzr) = ex_value
+      end where
+    end if
+
 !   Make sure to mark all points outside of the rectangular box as excised
 !   points and set the value of f to -ex_value.
     if ( use_outer_excision .gt. 0 ) then