code: def query(self, coords, order=1):
          
          return super(SFDQuery, self).query(coords, order=order)
docs: Returns E(B-V) at the specified location(s) on the sky. See Table 6 of
          Schlafly & Finkbeiner (2011) for instructions on how to convert this
          quantity to extinction in various passbands.
  
          Args:
              coords (`astropy.coordinates.SkyCoord`): The coordinates to query.
              order (Optional[int]): Interpolation order to use. Defaults to `1`,
                  for linear interpolation.
  
          Returns:
              A float array containing the SFD E(B-V) at every input coordinate.
              The shape of the output will be the same as the shape of the
              coordinates stored by `coords`.