Also integrate the mesh acceleration when computing the stars' and BHs'...
Also integrate the mesh acceleration when computing the stars' and BHs' half-dt kick to write the velocities in snapshots.
@jch you noted that this was inconsistent with the hydro and gravity code in the current version. Are you happy with this fix?