From 7ee9407fc3320e2ca0c428265f3bcd0eaa88beac Mon Sep 17 00:00:00 2001 From: diener Date: Tue, 13 May 2003 17:55:20 +0000 Subject: Added support for using emask. Needs to be tested more. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/EHFinder/trunk@102 2a26948c-0e4f-0410-aee8-f1d3e353619c --- schedule.ccl | 2 +- src/EHFinder_SetMask.F90 | 10 ++++++++++ 2 files changed, 11 insertions(+), 1 deletion(-) diff --git a/schedule.ccl b/schedule.ccl index 4ff5fe0..38fec6b 100644 --- a/schedule.ccl +++ b/schedule.ccl @@ -67,7 +67,7 @@ if (CCTK_Equals(metric_type,"static conformal") && CCTK_Equals(eh_metric_type,"n if (use_mask) { - schedule EHFinder_Read_Mask at CCTK_INITIAL + schedule EHFinder_Read_Mask at CCTK_INITIAL after MaskOne { LANG: Fortran } "Read in excision mask from file" 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 -- cgit v1.2.3