use *global* xyz positions to interpolate from Cactus grid

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@557 f88db872-0e4f-0410-b76b-b9085cfa78c5

1 files changed, 10 insertions, 7 deletions

diff --git a/src/gr/horizon_function.cc b/src/gr/horizon_function.cc
index 68d2a94..c21bfad 100644
--- a/src/gr/horizon_function.cc
+++ b/src/gr/horizon_function.cc
@@ -136,13 +136,16 @@ CCTK_VInfo(CCTK_THORNSTRING, "      xyz positions and derivative coefficients");
 			= p.ghosted_gridfn(ghosted_gfns::gfn__h, irho,isigma);
 		const fp rho = p.rho_of_irho(irho);
 		const fp sigma = p.sigma_of_isigma(isigma);
-		fp x, y, z;
-		p.xyz_of_r_rho_sigma(r,rho,sigma, x,y,z);
-
-		// xyz positions of grid points
-		p.gridfn(nominal_gfns::gfn__xx, irho,isigma) = x;
-		p.gridfn(nominal_gfns::gfn__yy, irho,isigma) = y;
-		p.gridfn(nominal_gfns::gfn__zz, irho,isigma) = z;
+		fp local_x, local_y, local_z;
+		p.xyz_of_r_rho_sigma(r,rho,sigma, local_x,local_y,local_z);
+		const fp global_x = ps.global_x_of_local_x(local_x);
+		const fp global_y = ps.global_y_of_local_y(local_y);
+		const fp global_z = ps.global_z_of_local_z(local_z);
+
+		// global xyz positions of grid points
+		p.gridfn(nominal_gfns::gfn__xx, irho,isigma) = global_x;
+		p.gridfn(nominal_gfns::gfn__yy, irho,isigma) = global_y;
+		p.gridfn(nominal_gfns::gfn__zz, irho,isigma) = global_z;
 
 		// 1st derivative coefficient gridfns X_ud
 		p.gridfn(nominal_gfns::gfn__X_ud_11, irho,isigma)