This python class provides resonance information using semianalytic methods developed in Gallardo (2020, CMDA 132:9)
Only restricted case has been developed so far. Only works in linux systems so far. User must have installed numpy, matplotlib and gfortran compiler as the core code was written in FORTRAN <3
User must initialize a resonance as follows:
reso1 = semianalyticResonance(star_mass, a_pla, e_pla, pla_mass, k, kp, e, inc, argper, lonNod)
The class semianalyticResonance() currently has six methods
- which_resonance() prints the resonance
- which_units() prints the units used in the code
- width() wich gives stable width
- sigmas() which gives equilibrium points
- period() which gives libration period for each equilibrium point
- hamiltonian(a_width, levels) which plots hamiltonian level curves of the theory
¡Check the example!
Further detailed documentation and examples are in development.
#SuperAtlas (only for windows)
In adittion to the mentioned class we provide a simple python version of Superatlas, an easy way to plot resonance shapes. As an example check funcion_varios_atlas.py which generates the next plot: