From e7b98e72d65a1c51fca3354371cbf6229cb3b11e Mon Sep 17 00:00:00 2001 From: Jo Bovy Date: Tue, 2 Jun 2015 21:15:21 -0400 Subject: [PATCH] directly specify degree unit to avoid confusion with kpc unit --- galpy/orbit_src/OrbitTop.py | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/galpy/orbit_src/OrbitTop.py b/galpy/orbit_src/OrbitTop.py index 304364c6d..c2b39a1c8 100644 --- a/galpy/orbit_src/OrbitTop.py +++ b/galpy/orbit_src/OrbitTop.py @@ -771,9 +771,10 @@ def SkyCoord(self,*args,**kwargs): """ radec= self._radec(*args,**kwargs) tdist= self.dist(*args,**kwargs) - return coordinates.SkyCoord(radec[:,0],radec[:,1], + return coordinates.SkyCoord(radec[:,0]*units.degree, + radec[:,1]*units.degree, distance=tdist*units.kpc, - frame='icrs',unit='deg') + frame='icrs') def _radec(self,*args,**kwargs): """Calculate ra and dec"""