code: def multitaper_spectrum(self, clm, k, convention='power', unit='per_l',
                              **kwargs):
          
          return self._multitaper_spectrum(clm, k, convention=convention,
                                           unit=unit, **kwargs)
docs: Return the multitaper spectrum estimate and standard error.
  
          Usage
          -----
          mtse, sd = x.multitaper_spectrum(clm, k, [convention, unit, lmax,
                                                    taper_wt, clat, clon,
                                                    coord_degrees])
  
          Returns
          -------
          mtse : ndarray, shape (lmax-lwin+1)
              The localized multitaper spectrum estimate, where lmax is the
              spherical-harmonic bandwidth of clm, and lwin is the
              spherical-harmonic bandwidth of the localization windows.
          sd : ndarray, shape (lmax-lwin+1)
              The standard error of the localized multitaper spectrum
              estimate.
  
          Parameters
          ----------
          clm : SHCoeffs class instance
              SHCoeffs class instance containing the spherical harmonic
              coefficients of the global field to analyze.
          k : int
              The number of tapers to be utilized in performing the multitaper
              spectral analysis.
          convention : str, optional, default = 'power'
              The type of output spectra: 'power' for power spectra, and
              'energy' for energy spectra.
          unit : str, optional, default = 'per_l'
              The units of the output spectra. If 'per_l', the spectra contain
              the total contribution for each spherical harmonic degree l. If
              'per_lm', the spectra contain the average contribution for each
              coefficient at spherical harmonic degree l.
          lmax : int, optional, default = clm.lmax
              The maximum spherical-harmonic degree of clm to use.
          taper_wt : ndarray, optional, default = None
              1-D numpy array of the weights used in calculating the multitaper
              spectral estimates and standard error.
          clat, clon : float, optional, default = 90., 0.
              Latitude and longitude of the center of the spherical-cap
              localization windows.
          coord_degrees : bool, optional, default = True
              True if clat and clon are in degrees.