galpy is a python package for galactic dynamics. It supports orbit integration in a variety of potentials, evaluating and sampling various distribution functions, and the calculation of action-angle coordinates for all static potentials.
https://ui.adsabs.harvard.edu/abs/2015ApJS..216...29B