예제 #1
0
파일: Merc.py 프로젝트: RJWorth/AlphaCen 
	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))
예제 #2
0
파일: Merc.py 프로젝트: RJWorth/AlphaCen 
	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))
예제 #3
0
파일: Merc.py 프로젝트: RJWorth/AlphaCen 
	def v(self):
		return(Merc.mag(self.vel))
예제 #4
0
파일: Merc.py 프로젝트: RJWorth/AlphaCen 
	def dr(self,other):
		'''Distance between two objects.'''
		x1 = self.pos
		x2 = other.pos
		dr = Merc.mag(x1-x2)
		return(dr)
예제 #5
0
파일: Merc.py 프로젝트: RJWorth/AlphaCen 
	def r(self):
		return(Merc.mag( self.pos))