Explore the results

While a simulation is in progress or when it has ended, the historic snapshot file can be interpreted to generate a plain text file and a plot file with the history of the simulation:

python scripts/explore_history.py target/case3.json   target/case3_history.bin
python scripts/explore_history.py target/case4.json   target/case4_history.bin
python scripts/explore_history.py target/case7.json   target/case7_history.bin
python scripts/explore_history.py target/example.json target/example_history.bin

TRAPPIST-1 explore history example

The plain text file follows a tab separated format and contains the following columns:

  • current_time: Current year.
  • planet: Identifier 1, 2, 3... for bodies orbiting the central body.
  • semi-major_axis_AU: Semi-major axis (AU).
  • corrotation_radius_AU: Distance at which the body should be to be in pseudo-synchronization given its current rotation period.
  • planet_obliquity_deg: Obliquity (degrees).
  • eccentricity: Eccentricity (degrees).
  • inclination_deg: Inclination (degrees).
  • energy_lost_due_to_tides_W_per_m2: Instantaneous tidal heat flux (W/m^2).
  • mean_energy_lost_due_to_tides_W_per_m2: Mean tidal heat flux (W/m^2).
  • planet_rotation_period_hours: Rotation period (hours).
  • planet_pseudo_synchronization_period: Period at which the body should rotate to be in pseudo-synchronization given its current orbital period.
  • energy_lost_due_to_tides_W: Energy lost due to tides (W).
  • mean_energy_lost_due_to_tides_W: Mean energy lost due to tides (W).
  • star_obliquity_deg: Central body obliquity with respect to the current body (degrees).
  • planet_precession_angle_deg: Precession angle (degrees).
  • conservation_of_angular_momentum: Conservation of angular momentum (ΔL/L0).
  • star_rotation_period_days: Central body rotation period (days).
  • conservation_of_energy: Conservation of (kinetic+potential) energy (ΔE/E0).


To explore what possible resonances might be present in the system:

python scripts/explore_timed_resonances.py target/case3_history.bin
python scripts/explore_timed_resonances.py target/case4_history.bin
python scripts/explore_timed_resonances.py target/case7_history.bin
python scripts/explore_timed_resonances.py target/example_history.bin

TRAPPIST-1 explore timed resonances example

Finally, to study a given resonance (e.g., 3:2) between planet one and two:

python scripts/explore_single_resonance.py target/case3_history.bin 1 2 3 2
python scripts/explore_single_resonance.py target/case4_history.bin 1 2 3 2
python scripts/explore_single_resonance.py target/case7_history.bin 1 2 3 2
python scripts/explore_single_resonance.py target/example_history.bin 1 2 3 2

TRAPPIST-1 explore resonance 3:2 between planet 6 and 7 example