Action-angle formalism for axisymmetric potentials using the adiabatic approximation, grid-based interpolation
INPUT:
pot= potential or list of potentials
zmax= zmax for building Ez grid
Rmax = Rmax for building grids
gamma= (default=1.) replace Lz by Lz+gamma Jz in effective potential
nEz=, nEr=, nLz, nR= grid size
numcores= number of cpus to use to parallellize
c= if True, use C to calculate actions
+scipy.integrate.quad keywords
OUTPUT: HISTORY:
2012-07-27 - Written - Bovy (IAS@MPIA)