Dear Justin Thanks for your quick answer. I want to obtain *radial density profile* relative to COM and not radial distribution function (RDF) relative to COM. Are you sure g_rdf -com is appropriate for this aim? Thanks