Option to calculate extra values (e.g. potential) for output in snapshots

Not urgent obviously, but in principle it'd be nice to run with the faster non-potential gravity for the simulation but to still calculate the potential for each particle when about to write each snapshot.