docs: def _xys(date):
      """Get The X, Y and s coordinates
  
      Args:
          date (Date):
      Return:
          3-tuple of float: Values of X, Y and s, in radians
      """
  
      X, Y, s_xy2 = _xysxy2(date)
  
      # convert milli-arcsecond to arcsecond
      dX, dY = date.eop.dx / 1000., date.eop.dy / 1000.
  
      # Convert arcsecond to degrees then to radians
      X = np.radians((X + dX) / 3600.)
      Y = np.radians((Y + dY) / 3600.)
      s = np.radians(s_xy2 / 3600.) - (X * Y / 2)
  
      return X, Y, s