diff options
| author | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2001-09-24 15:36:34 +0000 |
|---|---|---|
| committer | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2001-09-24 15:36:34 +0000 |
| commit | 23ed1898df8ea0e029cb59ddb1f19c79d9e12bfe (patch) | |
| tree | 88d31c9e78bfde6892c93045741e7c61c9e98846 /src/AHFinder_gau.F | |
| parent | 505a935cb8fa38bcf971da2e663a0fed766a2a33 (diff) | |
Changing the drift correction so that it sees the correct power of
psi controlled by the parameter rotation_psipower.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@247 89daf98e-ef62-4674-b946-b8ff9de2216c
Diffstat (limited to 'src/AHFinder_gau.F')
| -rw-r--r-- | src/AHFinder_gau.F | 29 |
1 files changed, 16 insertions, 13 deletions
diff --git a/src/AHFinder_gau.F b/src/AHFinder_gau.F index 6534a56..3b1bc67 100644 --- a/src/AHFinder_gau.F +++ b/src/AHFinder_gau.F @@ -650,24 +650,27 @@ ! Actually apply drift correction to shift. if (drift_correct_on.eq.1) then - if (ahf_ncall.gt.1) then - if ((.not.find3).or. - . (find3.and.(mfind-1.eq.drift_correct_horizon))) then + if ((.not.find3).or. + . (find3.and.(mfind-1.eq.drift_correct_horizon))) then - call CCTK_ReduceLocalScalar(ierror,cctkGH,-1,sum_handle, - . omega_corr,omega_corr_all,CCTK_VARIABLE_REAL) - if (ierror.ne.0) then - call CCTK_WARN(1,"Reduction failed!") - end if + call CCTK_ReduceLocalScalar(ierror,cctkGH,-1,sum_handle, + . omega_corr,omega_corr_all,CCTK_VARIABLE_REAL) + if (ierror.ne.0) then + call CCTK_WARN(1,"Reduction failed!") + end if - omega_corr = omega_corr_all + omega_corr = omega_corr_all - ahf_omega_corr = omega_corr - ahf_omega_cum = omega_cum + ahf_omega_corr = omega_corr + ahf_omega_cum = omega_cum + + if (ahf_ncall.gt.1) then + + betax = betax - y*omega_corr/psi**rotation_psipower + betay = betay + x*omega_corr/psi**rotation_psipower - betax = betax - y*omega_corr/psi**2 - betay = betay + x*omega_corr/psi**2 end if + end if end if |