Commit c50964f0 authored by Olaf's avatar Olaf

minor changes

parent 51f3ed44
......@@ -7,7 +7,7 @@ M = 5.972e24 # Earth' mass [kg]
Au = 149597870 # Astronomical unit [km]
Y = 365.25*24*3600*11.862 # Year [s]
G = 6.67408e-11*M*(Au*10**(3))**(-3)*Y**2 # Gravitational constant [m**3 * kg**(-1) * s**(-2)]
massSun = 333000 # AU
massSun = 333000 # M
class Body:
......@@ -181,7 +181,7 @@ class Resonance_astroids:
self.mass = par['mass']
except:
self.mass = 0.012
print("Total mass of belt not given, set to mass of Moon: {} AU".format(self.mass))
print("Total mass of belt not given, set to mass of Moon: {} Earth masses".format(self.mass))
# Individual masses
self.masses = (self.mass/self.N*np.ones(self.N)).tolist()
......@@ -191,7 +191,7 @@ class Resonance_astroids:
self.masscenter = pas['masscenter']
except:
self.masscenter = massSun
print("Center mass of belt not given, set to mass of Sun: {} AU".format(self.masscenter))
print("Center mass of belt not given, set to mass of Sun: {} Earth masses".format(self.masscenter))
# Belt radius
try:
......@@ -281,7 +281,7 @@ class Asteroids:
self.mass = par['mass']
except:
self.mass = 0.012
print("Total mass of belt not given, set to mass of Moon: {} AU".format(self.mass))
print("Total mass of belt not given, set to mass of Moon: {} Earth masses".format(self.mass))
# Individual masses
self.masses = (self.mass/self.N*np.ones(self.N)).tolist()
......@@ -291,7 +291,7 @@ class Asteroids:
self.masscenter = pas['masscenter']
except:
self.masscenter = massSun
print("Center mass of belt not given, set to mass of Sun: {} AU".format(self.masscenter))
print("Center mass of belt not given, set to mass of Sun: {} Earth masses".format(self.masscenter))
# Inner belt radius
try:
......
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