docs: def get_central_wave(wavl, resp):
      """Calculate the central wavelength or the central wavenumber,
      depending on what is input
      """
      return np.trapz(resp * wavl, wavl) / np.trapz(resp, wavl)