My implementation of orbit calculation with AGAMA.
This code uses Astropy, pyGaia, and AGAMA (Vasiliev 2019) to integrate orbits in an axisymmetric Galactic model potential and, then, extract several kinematic/dynamical quantities, including angular momenta, actions, and orbital energies.
The standard model potental I adopt is McMillan 2017, which is described in the file "McMillan17.ini" in this repository.
If you have any questions, get in contact with me: guilherme.limberg [at] usp.br.