aboutsummaryrefslogtreecommitdiff
path: root/src
diff options
context:
space:
mode:
authorshawley <shawley@89daf98e-ef62-4674-b946-b8ff9de2216c>2002-03-13 16:09:58 +0000
committershawley <shawley@89daf98e-ef62-4674-b946-b8ff9de2216c>2002-03-13 16:09:58 +0000
commit3a9f3270333c886f62b9dc8e77ab9aaa82fb7eda (patch)
tree4302fce3bc4bc3c88f4a4d67142521d2481506df /src
parente8e3d0345ccb0b4a3ad64c36490851e7c886235c (diff)
added possiblity to set ahf_centroid using parameter centroid_horizon
to decide which horizon to use when setting the centroid git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@288 89daf98e-ef62-4674-b946-b8ff9de2216c
Diffstat (limited to 'src')
-rw-r--r--src/AHFinder_gau.F17
1 files changed, 10 insertions, 7 deletions
diff --git a/src/AHFinder_gau.F b/src/AHFinder_gau.F
index 4baad67..8ef92a9 100644
--- a/src/AHFinder_gau.F
+++ b/src/AHFinder_gau.F
@@ -590,13 +590,16 @@
avgx = avgx/dble(npoints)
avgy = avgy/dble(npoints)
- ahf_centroid_x = avgx
- ahf_centroid_y = avgy
-
- if (verbose) then
- write(1,*) 'AHFinder_gau: ahf_centroid_x = ', ahf_centroid_x
- write(1,*) 'AHFinder_gau: ahf_centroid_y = ', ahf_centroid_y
- end if
+ if ((.not.find3).or.
+ . (find3.and.(mfind-1.eq.centroid_horizon))) then
+ ahf_centroid_x = avgx
+ ahf_centroid_y = avgy
+
+ if (verbose) then
+ write(1,*) 'AHFinder_gau: ahf_centroid_x = ', ahf_centroid_x
+ write(1,*) 'AHFinder_gau: ahf_centroid_y = ', ahf_centroid_y
+ end if
+ endif
! Other processors.