mirror/rebound-collisions
1
0
Fork 0
mirror of https://github.com/Findus23/rebound-collisions.git synced 2024-09-19 15:53:48 +02:00
Code Releases Activity

fix inclination

Browse source
This commit is contained in:
Lukas Winkler 2021-03-23 13:21:20 +01:00
parent fb01d09404
commit c28b112b7a
Signed by: lukas
GPG key ID: 54DE4D798D244853
2 changed files with 5 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
timeplot.py
View file

@ -1,6 +1,6 @@
import argparse
from collections import namedtuple
from math import log10
from math import log10, degrees
import matplotlib
import matplotlib.animation as animation
@ -65,9 +65,9 @@ def update_plot(num: int, args: MyProgramArgs, sa: SimulationArchive, ed: ExtraD
if args.y_axis == "e":
bla = np.array([a, [p.e for p in sim.particles[1:]]])
elif args.y_axis == "i":
bla = np.array([a, [p.inc for p in sim.particles[1:]]])
bla = np.array([a, [degrees(p.inc) for p in sim.particles[1:]]])
elif args.y_axis == "Omega":
bla = np.array([a, [p.Omega for p in sim.particles[1:]]])
bla = np.array([a, [degrees(p.Omega) for p in sim.particles[1:]]])
else:
raise ValueError("invalid y-axis")
dots.set_offsets(bla.T)
@ -104,7 +104,7 @@ def main(args: MyProgramArgs):
plt.ylim(-0.1, 1) # e
plt.ylabel("e")
elif args.y_axis == "i":
plt.ylim(0, 4) # i
plt.ylim(0, 10) # i
plt.ylabel("i")
elif args.y_axis == "Omega":
plt.ylim(0, 360) # i

2
water_sim.py
View file

@ -124,7 +124,7 @@ def main(fn: Path, testrun=False):
else:
part = Particle(
m=columns[0], a=columns[1], e=columns[2],
inc=radians(columns[3]), omega=columns[4],
inc=columns[3], omega=columns[4],
Omega=columns[5], M=columns[6],
simulation=sim,
hash=hash,