Before output_00040 in LCDM, don’t use the field “boxlen”in hdf5 files which is not correct. The true boxlen is 2625 com Mpc/h.

In pfof hdf5 particle files, don’t use the potential which is very approximative (NGP). Instead the potential interpolated from the gravity cells with Magrathea is much more accurate (TSC).