Add improved method for phasespace output