Commit 59de56d0 authored by Max Capelle's avatar Max Capelle

add theta for asteroid belt

parent 9795d80a
......@@ -221,6 +221,13 @@ class Resonance_astroids:
self.T_ref = 12*365.25*24*3600/Y
print("Belt reference period not given, set to period of Jupiter {}".format(self.T_ref))
# Initial angle
try:
self.theta = par['theta']
except:
self.theta = 0
print("Resonance belt initial angle not given, set to {}".format(self.theta))
# Determine periods of asteroids
self.T = np.linspace(1,self.N,self.N)*self.T_ref
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment