diff options
| author | diener <diener@2a26948c-0e4f-0410-aee8-f1d3e353619c> | 2003-05-13 17:55:20 +0000 |
|---|---|---|
| committer | diener <diener@2a26948c-0e4f-0410-aee8-f1d3e353619c> | 2003-05-13 17:55:20 +0000 |
| commit | 7ee9407fc3320e2ca0c428265f3bcd0eaa88beac (patch) | |
| tree | 26e0c4da5e850c08fcf4a501c1ce30416a271bfa | |
| parent | 11237939d932a29eb65ac50b7414f5e0dabed203 (diff) | |
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
| -rw-r--r-- | schedule.ccl | 2 | ||||
| -rw-r--r-- | src/EHFinder_SetMask.F90 | 10 |
2 files changed, 11 insertions, 1 deletions
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 |