rdf script to compute the radial distribution functions from the trajectories