docs: def _obsolete_plot_el_abund_marco(directory,name_h5_file,mass_range,cycle,logic_stable,i_decay,file_solar,solar_factor,symbol='ko'):
      """
      Interface to plot elements abundances averaged 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)
      # element abundances are calculated here
      mass_fractions_array_decayed = average_mass_frac_decay
      mass_fractions_array_not_decayed = average_mass_frac
      u.element_abund_marco(i_decay,stable,jjdum,mass_fractions_array_not_decayed,mass_fractions_array_decayed)
  
  
      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(u.z_bismuth):
              pl.plot(z_for_elem[i],elem_prod_fac[i],symbol,markersize=10.)
  
          pl.xlabel('$Atomic$ $number$', fontsize=20)
          pl.ylabel('$X_{i}/X_{sun}$', fontsize=20)
  
          pl.ylim(1.0e-2,1000.)
          pl.xlim(0,95)
  
      elif logic_stable:
          for i in range(u.z_bismuth):
              if index_stable[i] == 1:
                  continue
                  #pl.plot(z_for_elem[i],elem_prod_fac[i],'ko')
          if i_decay == 2:
              for i in range(u.z_bismuth):
                  if index_stable[i] == 1:
                      pl.plot(z_for_elem[i],elem_prod_fac_decayed[i],symbol,markersize=10.)
  
          pl.xlabel('$Atomic$ $number$', fontsize=20)
          pl.ylabel('$X_{i}/X_{sun}$', fontsize=20)
  
          pl.ylim(1.0e-2,1000.)
          pl.xlim(0,95)
  
  
  
      pl.grid()
      pl.show()