def e(self): '''Calculate eccentricity of an object's orbit''' val = 1.+2.*self.eps()*Merc.mag(self.h())**2./self.gm()**2. if ((val < 0) & (-val < tol)): val = 0 e = sqrt(val) return(np.array(e))
def i(self): '''Calculate inclination of an object's orbit. Takes mks, returns degrees''' i = arccos(self.h()[2]/Merc.mag(self.h()))*180./pi return(np.array(i))
def v(self): return(Merc.mag(self.vel))
def dr(self,other): '''Distance between two objects.''' x1 = self.pos x2 = other.pos dr = Merc.mag(x1-x2) return(dr)
def r(self): return(Merc.mag( self.pos))