the attached patch allows the puncture velocities to be output (simply moves the result of an interpolation into grid functions).
The patch looks fine, but why do you define pt_vel_t if the only thing it will ever be set to is cctk_time?
I only define pt_vel_t to do the very same thing that is done for pt_loc (maybe it also makes output nicer?).
I applied the patch (as is, the extra double value should not hurt anyone and it makes pt_vel the same as pt_loc).