The three following problems have been found with the stats outputs:

• The time and step values are incorrect. Easy to show, run simulation with runtime = 1.0, then read the stats and you will see that the final time value is nowhere near 1.0. Similarly, the step number does not match when using data.attrs['output.stats.interval_time']. Current work around is using data.attrs['output.stats.interval_time'] and making new time array based off index of array, so seems to be outputting at right times but steps and time are wrong values.

• T00 in stats is not in units of m0 n0, seems to be related to the sum of particle Lorentz factor, but off by factor of 2 when testing?

• quantities labels does not work for some values, eg E2 and Npart give errors as "Unrecognized stats name" (while E^2 works...)