def element_abund_marco(i_decay, stable_isotope_list,<EOL>stable_isotope_identifier,<EOL>mass_fractions_array_not_decayed,<EOL>mass_fractions_array_decayed):
<EOL>global elem_abund<EOL>elem_abund = np.zeros(z_bismuth)<EOL>global elem_abund_decayed<EOL>elem_abund_decayed = np.zeros(z_bismuth)<EOL>global elem_prod_fac<EOL>elem_prod_fac = np.zeros(z_bismuth)<EOL>global elem_prod_fac_decayed<EOL>elem_prod_fac_decayed = np.zeros(z_bismuth)<EOL>for i in range(z_bismuth):<EOL><INDENT>dummy = <NUM_LIT:0.><EOL>for j in range(len(spe)):<EOL><INDENT>if znum_int[j] == i+<NUM_LIT:1> and stable_isotope_identifier[j] > <NUM_LIT:0.5>:<EOL><INDENT>dummy = dummy + float(mass_fractions_array_not_decayed[j])<EOL><DEDENT><DEDENT>elem_abund[i] = dummy<EOL><DEDENT>for i in range(z_bismuth):<EOL><INDENT>if index_stable[i] == <NUM_LIT:1>:<EOL><INDENT>elem_prod_fac[i] = float(old_div(elem_abund[i],solar_elem_abund[i]))<EOL><DEDENT>elif index_stable[i] == <NUM_LIT:0>:<EOL><INDENT>elem_prod_fac[i] = <NUM_LIT:0.><EOL><DEDENT><DEDENT>if i_decay == <NUM_LIT:2>:<EOL><INDENT>for i in range(z_bismuth):<EOL><INDENT>dummy = <NUM_LIT:0.><EOL>for j in range(len(mass_fractions_array_decayed)):<EOL><INDENT>if znum_int[cl[stable_isotope_list[j].capitalize()]] == i+<NUM_LIT:1>:<EOL><INDENT>dummy = dummy + float(mass_fractions_array_decayed[j])<EOL><DEDENT><DEDENT>elem_abund_decayed[i] = dummy<EOL><DEDENT>for i in range(z_bismuth):<EOL><INDENT>if index_stable[i] == <NUM_LIT:1>:<EOL><INDENT>elem_prod_fac_decayed[i] = float(old_div(elem_abund_decayed[i],solar_elem_abund[i]))<EOL><DEDENT>elif index_stable[i] == <NUM_LIT:0>:<EOL><INDENT>elem_prod_fac_decayed[i] = <NUM_LIT:0.><EOL><DEDENT><DEDENT><DEDENT>
Given an array of isotopic abundances not decayed and a similar
array of isotopic abundances not decayed, here elements abundances,
and production factors for elements are calculated