Tibeuleu/KozaiLidov
1
0
Fork 0
Code Activity

revert with old dimensions

Browse Source
This commit is contained in:
Thibault Barnouin
2021-11-18 20:10:14 +01:00
parent ff1c46df71
commit e076127f6d
6 changed files with 10 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
main.py
View File

@@ -9,8 +9,8 @@ from lib.units import *
def main():
#initialisation
m = np.array([1., 1., 1e-5])*Ms/Ms # Masses in Solar mass
a = np.array([1., 1., 5.])*au/au # Semi-major axis in astronomical units
m = np.array([1., 1., 1e-5])*Ms # Masses in Solar mass
a = np.array([1., 1., 5.])*au # Semi-major axis in astronomical units
e = np.array([0., 0., 1./4.]) # Eccentricity
psi = np.array([0., 0., 0.])*np.pi/180. # Inclination of the orbital plane in degrees
@@ -25,10 +25,10 @@ def main():
v = np.array([v1, v2, v3])
#integration parameters
duration, step = 100*yr/yr, np.array([1./(365.25*2.), 1./365.25])*yr/yr #integration time and step in years
duration, step = 100*yr, np.array([1./(365.25*2.), 1./365.25])*yr #integration time and step in years
integrator = "leapfrog"
n_bodies = 2
display = True
display = False
savename = "{0:d}bodies_{1:s}".format(n_bodies, integrator)
#simulation start