docs: 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