PunctureTracker: Correct bug introduced in last commit

Was setting the spherical surface variables from the CarpetRegrid2 variables
rather than the puncture tracker variables.


git-svn-id: http://svn.aei.mpg.de/numrel/AEIThorns/PunctureTracker/trunk@16 a2659f00-0f4f-0410-9214-a4596bbb8a4f

1 files changed, 13 insertions, 12 deletions

diff --git a/src/puncture_tracker.c b/src/puncture_tracker.c
index 786ed93..8c3341b 100644
--- a/src/puncture_tracker.c
+++ b/src/puncture_tracker.c
@@ -282,25 +282,26 @@ PunctureTracker_SetPositions (CCTK_ARGUMENTS)
 
   for (int n = 0; n < max_num_tracked; ++ n) {
     if (track[n]) {
-      if (verbose) {
-        CCTK_VInfo (CCTK_THORNSTRING,
-                    "Setting spherical surface centroid from puncture #%d to (%g,%g,%g)",
-                    n,
-                    (double)pt_loc_x[n],
-                    (double)pt_loc_y[n],
-                    (double)pt_loc_z[n]);
-      }
-    
       // store puncture location in spherical surface
       if (which_surface_to_store_info[n] != -1) {
 	int sn = which_surface_to_store_info[n];
 	
-	sf_centroid_x[sn] = position_x[n];
-	sf_centroid_y[sn] = position_y[n];
-	sf_centroid_z[sn] = position_z[n];
+	sf_centroid_x[sn] = pt_loc_x[n];
+	sf_centroid_y[sn] = pt_loc_y[n];
+	sf_centroid_z[sn] = pt_loc_z[n];
 	
 	sf_active[sn] = 1;
 	sf_valid[sn] = 1;
+
+        if (verbose) {
+          CCTK_VInfo (CCTK_THORNSTRING,
+                      "Setting spherical surface %d centroid from puncture #%d to (%g,%g,%g)",
+                      sn,
+                      n,
+                      (double)pt_loc_x[n],
+                      (double)pt_loc_y[n],
+                      (double)pt_loc_z[n]);
+        }
       }
     }
   }