diff options
| -rw-r--r-- | src/puncture_tracker.c | 25 |
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]); + } } } } |