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.
Please let me (bovy -at- ias.edu) know if you make use of galpy in a publication.
Uses what later became the orbit integration routines and Dehnen and Shu disk distribution functions.
Employs galpy orbit integration in galpy.potential.MWPotential to characterize the orbits in the SEGUE G dwarf sample.
Uses galpy.potential force and density routines to characterize the difference between the vertical force and the surface density at large heights above the MW midplane.
Utilizes the Dehnen distribution function to inform a simple model of the velocity distribution of APOGEE stars in the Milky Way disk and to create mock data.
Makes use of potential models, the adiabatic and Staeckel actionAngle modules, and the quasiisothermal DF to model the dynamics of the SEGUE G dwarf sample in mono-abundance bins.
Employs galpy to perform orbit integrations in galpy.potential.MWPotential to characterize the orbits of stars in the APOGEE sample.
Please link back to http://github.com/jobovy/galpy . When using the galpy.actionAngle modules, please cite 2013ApJ...779..115B in addition to the papers describing the algorithm used. When orbit integrations are used, you could cite 2010ApJ...725.1676B (first galpy paper).