Example #1
0
	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))
Example #2
0
	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))
Example #3
0
	def v(self):
		return(Merc.mag(self.vel))
Example #4
0
	def dr(self,other):
		'''Distance between two objects.'''
		x1 = self.pos
		x2 = other.pos
		dr = Merc.mag(x1-x2)
		return(dr)
Example #5
0
	def r(self):
		return(Merc.mag( self.pos))