add read_phasespace.py script

[ci skip]