The current distances from the Sun, proper motions, and line-of-sight velocities, and the related uncertainties, of all 159 globular clusters are taken respectively from Baumgardt & Vasiliev (2021) and Vasiliev & Baumgardt (2021). These values are then converted into Galactocentric positions and velocities by making use of Astropy, and used as initial conditions to execute step (i) here.
Step (ii) here requires generating an N-body system, representing the globular cluster, whose initial total mass and half-mass radius are assigned on the basis of their current values, as given by Baumgardt & Hilker (2018) and reported in this Table. These values are also those adopted to generate the phase-space coordinates of each cluster, modeled as a Plummer distribution.