Add the shift applied to particles to the snapshot headers

Implement #852 (closed).

• Store the initial shift applied to all the particles in the space object.
• Write the shift to the output headers.
• Document this.