aboutsummaryrefslogtreecommitdiff
path: root/src
diff options
context:
space:
mode:
authordiener <diener@89daf98e-ef62-4674-b946-b8ff9de2216c>2001-12-17 15:12:07 +0000
committerdiener <diener@89daf98e-ef62-4674-b946-b8ff9de2216c>2001-12-17 15:12:07 +0000
commit4ae20305102bd9ad0cb706279a7360e127cc2142 (patch)
tree2c23e3245507d627a5310506fd1e5cfa13d3c46f /src
parentd29c6c5b8c401c0ec15fb3310be5fd9dbe0836c8 (diff)
Moved the assignment of ahf_centroid_x and ahf_centroid_y to make sure
that they are output even when drift correction is not on. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@265 89daf98e-ef62-4674-b946-b8ff9de2216c
Diffstat (limited to 'src')
-rw-r--r--src/AHFinder_gau.F6
1 files changed, 3 insertions, 3 deletions
diff --git a/src/AHFinder_gau.F b/src/AHFinder_gau.F
index bdd3119..803c7da 100644
--- a/src/AHFinder_gau.F
+++ b/src/AHFinder_gau.F
@@ -581,6 +581,9 @@
avgx = avgx/dble(npoints)
avgy = avgy/dble(npoints)
+ ahf_centroid_x = avgx
+ ahf_centroid_y = avgy
+
! Calculate correction in omega to stop "drift" of apparent horizon.
if (drift_correct_on.eq.1) then
@@ -588,9 +591,6 @@
if ((.not.find3).or.
. (find3.and.(mfind-1.eq.drift_correct_horizon))) then
- ahf_centroid_x = avgx
- ahf_centroid_y = avgy
-
xi_new = atan2(avgx,avgy)
write(0,*) 'AHFinder_gau: ahf_centroid_x = ',ahf_centroid_x