docs: def sample(self, wave):
          """Sample the spectrum.
  
          This uses :meth:`resample` to do the actual computation.
  
          Parameters
          ----------
          wave : number or array_like
              Wavelength set for sampling, given in user unit.
  
          Returns
          -------
          throughput : number or array_like
              Sampled throughput.
  
          """
          angwave = self.waveunits.ToAngstrom(wave)
          return self.__call__(angwave)