Tibeuleu/KozaiLidov
1
0
Fork 0
Code Activity

display sma and ecc for 2 bodies

Browse Source
This commit is contained in:
Thibault Barnouin
2021-11-19 16:28:41 +01:00
parent 400e140c35
commit 1c20570200
2 changed files with 7 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
main.py
View File

@@ -27,8 +27,8 @@ def main():
v = np.array([v1, v2, v3])
#integration parameters
duration, step = 100*yr, np.array([1./(365.25*2.), 1./(365.25*2.), 1./365.25])*yr #integration time and step in years
duration, step = 100*yr, np.array([1./(365.25*2.), 1./(365.25*1.), 5./(365.25*1.)])*yr #integration time and step in years
step = np.sort(step)[::-1]
integrator = "leapfrog"
n_bodies = 2
display = False
@@ -52,7 +52,7 @@ def main():
parameters = [duration, step, dyn_syst, integrator]
display_parameters(E, L, sma, ecc, parameters=parameters, savename=savename)
print(sma/au)
return 0
if __name__ == '__main__':