diff options
author | diener <diener@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2001-12-17 15:12:07 +0000 |
---|---|---|
committer | diener <diener@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2001-12-17 15:12:07 +0000 |
commit | 4ae20305102bd9ad0cb706279a7360e127cc2142 (patch) | |
tree | 2c23e3245507d627a5310506fd1e5cfa13d3c46f /src | |
parent | d29c6c5b8c401c0ec15fb3310be5fd9dbe0836c8 (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.F | 6 |
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 |