docs: def get_abundance_iso_decay(self,cycle):
          """
          returns the decayed stable isotopes.
  
          Parameters
          ----------
          cycle : integer
              The cycle.
  
          """
  
          import nuutils as u
  
          masses_for_this_cycle = self.se.get(cycle,'mass')
          self._read_iso_abund_marco([min(masses_for_this_cycle),max(masses_for_this_cycle)],cycle)
  
          u.stable_specie()
          self.decay(self.mass_frac)
  
          self.index_for_all_species = u.cl
          self.index_for_stable_species = u.back_ind
  
          self.decayed_stable_isotopes_per_cycle = decayed_multi_d
  
          # from here read solar abundances
          solar_factor = 2.
          u.solar('iniab1.0E-02.ppn_GN93',solar_factor)
  
          self.stable_isotope_identifier=u.jjdum
          self.stable_isotope_list=u.stable
  
          self.isotopic_production_factors=[]
          for i in range(len(masses_for_this_cycle)):
              pf_dum=[]
              jj=0
              for j in range(len(self.stable_isotope_identifier)):
                  if self.stable_isotope_identifier[j] == 1:
                      pf_dum.append(float(old_div(self.mass_frac[i][self.index_for_all_species[self.stable_isotope_list
  [jj].capitalize()]],u.solar_abundance[self.stable_isotope_list[jj].lower()])))
                      jj=jj+1
                  #elif self.stable_isotope_identifier[j] == 0:
                  #       pf_dum.append(float(0.))
              self.isotopic_production_factors.append(pf_dum)
  
          self.isotopic_production_factors_decayed=[]
          for i in range(len(masses_for_this_cycle)):
              pf_dum_d=[]
              jj=0
              for j in range(len(self.stable_isotope_identifier)):
                  if self.stable_isotope_identifier[j] == 1:
                      pf_dum_d.append(float(old_div(self.decayed_stable_isotopes_per_cycle[i][self.index_for_stable_species[self.stable_isotope_list
  [jj].upper()]],u.solar_abundance[self.stable_isotope_list[jj].lower()])))
                      jj=jj+1
              self.isotopic_production_factors_decayed.append(pf_dum_d)