galpy.potential.Potential.lindbladR

Potential.lindbladR(OmegaP, m=2, **kwargs)

NAME:

lindbladR

PURPOSE:

calculate the radius of a Lindblad resonance

INPUT:

OmegaP - pattern speed

m= order of the resonance (as in m(O-Op)=kappa (negative m for outer)
use m=’corotation’ for corotation +scipy.optimize.brentq xtol,rtol,maxiter kwargs

OUTPUT:

radius of Linblad resonance, None if there is no resonance

HISTORY:

2011-10-09 - Written - Bovy (IAS)