docs: def abup_se_plot(mod,species):
  
          """
          plot species from one ABUPP file and the se file.
  
          You must use this function in the directory where the ABP files
          are and an ABUP file for model mod must exist.
  
          Parameters
          ----------
          mod : integer
              Model to plot, you need to have an ABUPP file for that
              model.
          species : string
              The species to plot.
  
          Notes
          -----
          The species is set to 'C-12'.
  
          """
  
  # Marco, you have already implemented finding headers and columns in
  # ABUP files. You may want to transplant that into here?
          species='C-12'
  
          filename = 'ABUPP%07d0000.DAT' % mod
          print(filename)
          mass,c12=np.loadtxt(filename,skiprows=4,usecols=[1,18],unpack=True)
          c12_se=self.se.get(mod,'iso_massf','C-12')
          mass_se=self.se.get(mod,'mass')
  
          pyl.plot(mass,c12)
          pyl.plot(mass_se,c12_se,'o',label='cycle '+str(mod))
          pyl.legend()