docs: def _obsolete_plot_iso_abund_marco(directory, name_h5_file, mass_range,
                                    cycle, logic_stable, i_decay,
                                    file_solar, solar_factor):
      """
      Interface to plot average over mass_range.
  
      Parameters
      ----------
      directory : string
          Location of h5 file to plot.  Needed for plot_tools.
      name_h5_file : string
          Name of h5 file.  Needed for plot_tools.
      mass_range : list
          A 1x2 array required to plot data in a certain mass range.  Needed for
          _read_iso_abund_marco.
      cycle : integer
          which cycle from the h5 file?.  Needed for _read_iso_abund_marco.
      logic_stable : boolean
          Do you want to plot only stable or not.
      i_decay : integer
          If i_decay is 1, then plot not decayed.  If i_decay is 2, then
          plot decayed.  Make sense only if stable is true.
      file_solar : string
          File where to take solar abundances.
      solar_factor : float
          value to correct initial abundances to solar, e.g. for Z=0.01
          and AG89 solar_factor = 2.
  
      See Also
      --------
      se._read_iso_abund_marco()
  
      """
  
  
      # provide library for Z versus element names, and Z for elements
      u.give_zip_element_z_and_names()
      # solar abundances are read here
      u.solar(file_solar,solar_factor)
      # from here I have average abundances in mass_range to plot
      average_iso_abund_marco(mass_range,cycle,logic_stable,i_decay)
  
      fig = pl.figure()            # Figure object
      ax = fig.add_subplot(1,1,1)     # Axes object: one row, one column, first plot (one plot!)
      # Tick marks
      xminorlocator = MultipleLocator(1)
      xmajorlocator = MultipleLocator(10)
      ax.xaxis.set_major_locator(xmajorlocator)
      ax.xaxis.set_minor_locator(xminorlocator)
      yminorlocator = MultipleLocator(0.1)
      ymajorlocator = MultipleLocator(1)
      ax.yaxis.set_major_locator(ymajorlocator)
      ax.yaxis.set_minor_locator(yminorlocator)
  
      ax.set_yscale('log')
  
      if not logic_stable:
          for i in range(len(u.spe)):
              pl.plot(amass_int[u.cl[spe[i]]],average_mass_frac[u.cl[spe[i]]],'ko')
  
          pl.xlabel('$Mass$ $number$', fontsize=20)
          pl.ylabel('$X_{i}$', fontsize=20)
  
          pl.ylim(1.0e-10,10.)
          pl.xlim(55,110)
  
      elif logic_stable:
      # plot stable
          #for i in range(len(stable)):
          #    pl.plot(amass_int[cl[stable[i].capitalize()]],average_mass_frac[cl[stable[i].capitalize()]]/u.solar_abundance[stable[i].lower()],'ko')
  
          if i_decay == 2:
              for j in range(len(stable)):
                      #print cl[stable[j].capitalize()],stable[j].capitalize(),amass_int[cl[stable[j].capitalize()]]
                  pl.plot(amass_int[u.cl[stable[j].capitalize()]],old_div(u.mass_fractions_array_decayed[u.back_ind[stable[j]]],u.solar_abundance[stable[j].lower()]),'Dk')
  
          for i in range(len(stable)):
              for j in range(len(stable)):
                  if stable[i][:2] == stable[j][:2]:
                      if stable[i] == stable[j-1]:
                          adum  =[amass_int[u.cl[stable[i].capitalize()]],amass_int[u.cl[stable[j].capitalize()]]]
                          mfdum =[old_div(float(average_mass_frac[u.cl[stable[i].capitalize()]]),float(u.solar_abundance[stable[i].lower()])),old_div(float(average_mass_frac[u.cl[stable[j].capitalize()]]),float(u.solar_abundance[stable[j].lower()]))]
                          mfddum=[old_div(float(u.average_mass_frac_decay[u.back_ind[stable[i]]]),float(u.solar_abundance[stable[i].lower()])),old_div(float(u.average_mass_frac_decay[u.back_ind[stable[j]]]),float(u.solar_abundance[stable[j].lower()]))]
                          #pl.plot(adum,mfdum,'k-')
                          # I had to add this try/except...why? I guess is someone related to H2, that I spotted that was wrong in stable_raw...
                          # should deal without this. Have to be solved when I have time Marco (June 7 2011)
                          if i_decay == 2:
                              try:
                                  pl.plot(adum,mfddum,'k-')
                              except UnboundLocalError:
                                  continue
  
      pl.xlabel('$Mass$ $number$', fontsize=20)
      pl.ylabel('$X_{i}/X_{sun}$', fontsize=20)
  
      pl.ylim(1.0e-3,5000.)
      pl.xlim(55,210)
  
  
  
      pl.grid()
      pl.show()