diff options
| author | pollney <pollney@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2001-12-19 11:26:36 +0000 |
|---|---|---|
| committer | pollney <pollney@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2001-12-19 11:26:36 +0000 |
| commit | 9db55c0af91cfa53b74046528e8bd4f7d3552c3d (patch) | |
| tree | 9a59a5a3a3c22fc26432a7ea10025037be9e1824 | |
| parent | fb336778cf035af842437aa7f381e7c1f67047ff (diff) | |
A quick fix to get around a pgf90 compiler bug that seems to have
been triggered by the last set of commits. I commented out the
ahf_radvel_corr grid function and the reference to it in the code,
and this gets around the compiler problem. This is necessary for
us to do our runs on platinum. Note that you can still use Scott's
new drift correction parameters, you just can't see analyse their
effects by looking at the ahf_radvel_corr grid function until we
work out a more permanent fix.
Thanks to Gab and ThomasR for helping to diagnose this.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@271 89daf98e-ef62-4674-b946-b8ff9de2216c
| -rw-r--r-- | interface.ccl | 2 | ||||
| -rw-r--r-- | src/AHFinder_gau.F | 5 |
2 files changed, 3 insertions, 4 deletions
diff --git a/interface.ccl b/interface.ccl index 89afd9f..08fcdd6 100644 --- a/interface.ccl +++ b/interface.ccl @@ -99,7 +99,7 @@ real drift_correct_vars type=SCALAR ahf_centroid_y, ahf_omega_corr, ahf_omega_cum - ahf_radvel_corr, +# ahf_radvel_corr, } "Variables for drift correction" diff --git a/src/AHFinder_gau.F b/src/AHFinder_gau.F index 063dd22..021d612 100644 --- a/src/AHFinder_gau.F +++ b/src/AHFinder_gau.F @@ -602,7 +602,7 @@ write(0,*) 'AHFinder_gau: xi_drift = ',xi_drift - omega_corr = - drift_correct_factor * xi_drift + omega_corr = - 1.1 * xi_drift . /(dble(ahf_findevery)*cctk_delta_time) omega_cum = omega_cum + omega_corr @@ -629,7 +629,6 @@ ! Do drift correction in the radial direction if (rad_drift_correct.eq.1) then - r_new = sqrt(avgx**2 + avgy**2) if (ahf_ncall.gt.1) then @@ -701,7 +700,7 @@ radvel_corr = radvel_corr_all - ahf_radvel_corr = radvel_corr +! ahf_radvel_corr = radvel_corr if (ahf_ncall.gt.1) then |