add ALaDyn_particles_hists.py - this subrotuine plots and save the phase space histogram