def set_get_abu_distr_decay_old(self,cycles=20*[-1],mass_range=20*[[0,0]],ylim=20*[[0,0]],isotopes=['all'],linestyle=[],save_dir=''): 
		'''
			Plots HRDs


			end_model - array, control how far in models a run is plottet, if -1 till end
			symbs_1  - set symbols of runs
		'''
		if len(linestyle)==0:
			linestyle=200*['-']
		import nugridse as mp
		import utils as u
		print self.runs_H5_restart
		massfrac_all=[]
		iso_all=[]
                for i in range(len(self.runs_H5_restart)):
        		sefiles=mp.se(self.runs_H5_restart[i])
			cycle=cycles[i]
                        if cycle==-1:
                                cycle=int(sefiles.se.cycles[-1])			
			if mass_range[i][0] ==0 and mass_range[i][1]==0:
				mass_range[i][1]=sefiles.get(cycle,'mass')[-1]				
			sefiles.read_iso_abund_marco(mass_range[i],cycle)
		    	u.stable_specie()
		    	sefiles.decay(sefiles.mass_frac)
			idx_species=[]
			massfrac=[]
			iso=[]
			if not isotopes[0]=='all':
				for k in range(len(isotopes)):
					other_name_scheme=isotopes[k].split("-")[0].upper()+(5-len(isotopes[k])+1)*" "+isotopes[k].split("-")[1]
					#other_name_scheme=other_name_scheme.capitalize()

					idx_specie=u.back_ind[other_name_scheme]
					idx_species.append(idx_specie)
					massfrac.append(average_massfrac_decay[idx_specie])
				iso=isotopes
			else:
				massfrac=mp.average_mass_frac_decay
				other_name_scheme=u.back_ind
				iso=[]
				import re
				for kk in range(len(other_name_scheme)):
					list1=re.split('(\d+)',other_name_scheme[kk])
					newname=list1[0].capitalize()+'-'+list1[1]
					iso.append(newname)
			massfrac_all.append(massfrac)
			iso_all.append(iso)

		return iso_all,massfrac_all