def SkyCoord(self,*args,**kwargs):
        """
        NAME:
           SkyCoord
        PURPOSE:
           return the position as an astropy SkyCoord
        INPUT:
           t - (optional) time at which to get the position
           obs=[X,Y,Z] - (optional) position of observer (in kpc) 
                         (default=Object-wide default)
                         OR Orbit object that corresponds to the orbit
                         of the observer
                         Y is ignored and always assumed to be zero
           ro= distance in kpc corresponding to R=1. (default=Object-wide default)
        OUTPUT:
           SkyCoord(t)
        HISTORY:
           2015-06-02 - Written - Bovy (IAS)
        """
        kwargs.pop('quantity',None) # rm useless keyword to no conflict later
        _check_roSet(self,kwargs,'SkyCoord')
        radec= self._radec(*args,**kwargs)
        tdist= self.dist(quantity=False,*args,**kwargs)
        if not _APY3: # pragma: no cover
            return coordinates.SkyCoord(radec[:,0]*units.degree,
                                        radec[:,1]*units.degree,
                                        distance=tdist*units.kpc,
                                        frame='icrs')
        pmrapmdec= self._pmrapmdec(*args,**kwargs)
        vlos= self._lbdvrpmllpmbb(*args,**kwargs)[:,3]
        # Also return the Galactocentric frame used
        v_sun= coordinates.CartesianDifferential(\
            nu.array([-self._solarmotion[0],
                       self._solarmotion[1]+self._vo,
                       self._solarmotion[2]])*units.km/units.s)
        return coordinates.SkyCoord(radec[:,0]*units.degree,
                                    radec[:,1]*units.degree,
                                    distance=tdist*units.kpc,
                                    pm_ra_cosdec=pmrapmdec[:,0]\
                                        *units.mas/units.yr,
                                    pm_dec=pmrapmdec[:,1]*units.mas/units.yr,
                                    radial_velocity=vlos*units.km/units.s,
                                    frame='icrs',
                                    galcen_distance=\
                                        nu.sqrt(self._ro**2.+self._zo**2.)\
                                        *units.kpc,
                                    z_sun=self._zo*units.kpc,
                                    galcen_v_sun=v_sun)
NAME:
           SkyCoord
        PURPOSE:
           return the position as an astropy SkyCoord
        INPUT:
           t - (optional) time at which to get the position
           obs=[X,Y,Z] - (optional) position of observer (in kpc) 
                         (default=Object-wide default)
                         OR Orbit object that corresponds to the orbit
                         of the observer
                         Y is ignored and always assumed to be zero
           ro= distance in kpc corresponding to R=1. (default=Object-wide default)
        OUTPUT:
           SkyCoord(t)
        HISTORY:
           2015-06-02 - Written - Bovy (IAS)