diff options
| author | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2001-06-07 08:40:52 +0000 |
|---|---|---|
| committer | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2001-06-07 08:40:52 +0000 |
| commit | a559c62181ca54d696da394063e71be075942521 (patch) | |
| tree | e0e2c02771bd7da48a131323b805bb61b7ac5875 /src/AHFinder.F | |
| parent | 18c1ce2107381948c1aaad4a3a049172d9d8b40b (diff) | |
Adding parameters "ahf_mask_0" (etc) to control what horizons will be excised
when looking for more than one. This allowas the possibility of looking for
all three horizons and excising only one of them (the common one, say).
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@203 89daf98e-ef62-4674-b946-b8ff9de2216c
Diffstat (limited to 'src/AHFinder.F')
| -rw-r--r-- | src/AHFinder.F | 14 |
1 files changed, 13 insertions, 1 deletions
diff --git a/src/AHFinder.F b/src/AHFinder.F index 6c3d501..8e45be1 100644 --- a/src/AHFinder.F +++ b/src/AHFinder.F @@ -1698,7 +1698,19 @@ if ((horizon.and.(CCTK_EQUALS(ahf_mask,'strong'))).or. . (status.and.(mtype.eq.1).and.(CCTK_EQUALS(ahf_mask,'weak')))) then - call AHFinder_mask(CCTK_FARGUMENTS,rhor) + rhor = 0.0D0 + + if (find3) then + if ((mfind.eq.1).and.(ahf_mask_0.eq.1)) then + call AHFinder_mask(CCTK_FARGUMENTS,rhor) + else if ((mfind.eq.2).and.(ahf_mask_1.eq.1)) then + call AHFinder_mask(CCTK_FARGUMENTS,rhor) + else if ((mfind.eq.3).and.(ahf_mask_1.eq.1)) then + call AHFinder_mask(CCTK_FARGUMENTS,rhor) + end if + else + call AHFinder_mask(CCTK_FARGUMENTS,rhor) + end if ! Here I alter the parameters for SimpleExcision. Notice that at ! the moment, I only allow the excised cube to grow, and never |