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()