pyoculus.solvers.resonance_zone

class pyoculus.solvers.resonance_zone.ResonanceZone(manifold1: Manifold, manifold2: Manifold)

Class representing the resonnance zone formed by 2 stable and 2 unstable manifolds, intersected by 2 heteroclinic point and 2 fixed points.

Parameters:

• manifold1 (-) – The first manifold.

• manifold2 (-) – The second manifold.

- area

Compute the area of the resonance zone.

- plot

Plot the stable and unstable manifolds.

- contour

Compute the closed contour of the resonance zone.

area(whole_chain=False)

Compute the area of the resonance zone using Meiss’s action principle.

whole_chain: if True, compute the area of the whole chain of islands, otherwise only one island.

plot_manifolds(which='both', stepsize_limit=None, with_clinics=False, **kwargs)

Plot the 2 stable and the 2 unstable manifolds.

which (str): which manifold to plot. Can be ‘stable’, ‘unstable’ or ‘both’. stepsize_limit = threshold to remove big steps in the manifolds. with_clinics (bool): if True, plot the clinic points.

plot_reso_zone()

Compute the closed contour of the resonnance zone, formed by the stable and unstable manifolds of the two fixed points. Returns the contour points as a Nx2 array.

flux_approx_SL()

Approximate the flux through the resonnance zone using area computed with shoelace method, multiplied by the mean magnetic field inside the contour