galport.OrbitTools.classify_orbits

OrbitTools.classify_orbits(t_out=0.0, theta_p=None, time_resolution=None, family='ILR', time_around_res=False, amplitude_res=False)[source]
Parameters:

  • t_out ((M, ) float or numpy array) – array of times, in which we define the orbital type, by default 0.

  • theta_p ((N, ) numpy array, optional) – array of the perturbation (e.g. bar) rotation angle Default: None

  • time_resolution (float, optional) – time accuracy of series. Recommend don’t take too small Default: 5.

  • family (str, optional) – Default: ‘ILR’

  • time_around_res (bool, optional) – if True function estimate the resonance entry and exit times for resonant orbits, by default False

  • amplitude_res (bool, optional) – if True function estimate the maximum libration amplitude of the resonant angle, by default False

Returns:

  • types ((M, ) numpy array) – array of types (integer)

  • amplitude ((M, ) numpy array, optional) – array of angles amplitude for passage or resonant orbit.

  • times ((M, 2) numpy array, optional) – if time_around=True array of times for resonance and passage orbits, when they entered/left into resonance or began/end to pass through it.