dynamics

The analytic dynamics module of planetplanet. These are a collection of Julia scripts and are separate from the rest of the code.

offset_diagram2.jl

Creates a diagram showing the analytic means to compute a timing offset. Requires installation of the ExoJulia package for carrying out Kepler solver before calling. Also requires PyPlot, which calls matplotlib, and CGS.jl and regress.jl. For example, from Julia prompt:

julia> include("/PathToExoJulia/exojulia.jl")

where PathToExoJulia indicates the path where it is installed.

This is followed by:

julia> include("offset_diagram2.jl")

timing_offset2.jl

Creates simulated timing offsets for planet-planet events, and then fits with analytic formulae to recover the eccentricities and make a plot for the paper. Same dependencies as offset_diagram2.jl.

monte_carlo_circular.jl

Computes the probability of density and coplanarity of TRAPPIST-1 based upon the measured transit periods and the observed durations plus uncertainties (as well as the measured planet-star radius ratios). Uses CGS.jl and JLD, which saves variables in HDF files.

plot_circular.jl

Creates plots of the posterior distribution of density and coplanarity parameter for TRAPPIST-1. Uses the Optim package for optimizing a fit to the distribution, as well as PyPlot and JLD. Requires first running monte_carlo_circular.jl and then reading in the jld file (the name of the jld file contains the number of simulations carried out; this needs to be modified in plot_circular.jld).