GitHub - garyForeman/LMCorbit: orbits calculation for two of the Milky Way's satellite galaxies

garyForeman / LMCorbit Public

Notifications You must be signed in to change notification settings
Fork 0
Star 0

orbits calculation for two of the Milky Way's satellite galaxies

Notifications

Name		Name	Last commit message	Last commit date
Latest commit History 5 Commits
README		README
accelerate.py		accelerate.py
constants.py		constants.py
driver.py		driver.py
leapfrogKDK.py		leapfrogKDK.py
plot.py		plot.py

Repository files navigation

Author: Gary Foreman
email: [email protected]
Reference: Bekki & Chiba (2005)
http://arxiv.org/abs/astro-ph/0412318

Implementation of the three body orbits calculation for the Milky Way, 
Large Magellanic Cloud (LMC), and Small Magellanic Cloud. I perform integration
over time using the KDK leapfrog algorithm. Initialization parameters are 
outlined in the paper by Bekki & Chiba referenced above, and can be found in the
constants.py file. accelerate.py contains functions used to compute the 
accelerations on the LMC and SMC based on both gravitational potential and
dynamical friction. leapfrogKDK.py contains the actual implementation of the 
leapfrog algorithm. driver.py runs the orbits calculation and outputs the
orbits.dat file. plot.py generates a plot from the orbits.dat file (see Figure 2
of Bekki & Chiba for comparison).

Required python libraries:
matplotlib.pyplot
numpy 
scipy.optimize