These output files roughly model [rp-process nucleosynthesis](https://en.wikipedia.org/wiki/Rp-process). *out.xml* was generated with the single-zone network code with the accompanying master.h, Reaclib V2.2, and the command: **../single_zone_network --t9_0 2 --rho_0 1.e9 --tau 30 --network_xml ../nucnet-tools-code/data_pub/my_net.xml --nuc_xpath "[z <= 60]" --iterative_solver gmres --iterative_t9 2 --output_xml out.xml --init_mass_frac "{h1; 0.7}" "{he4; 0.2}" "{c12; 0.1}" --xml_steps 50 --t_end 1000.** *out2.xml* was generated with the single-zone network code with the accompanying master.h, Reaclib V2.2, and the command: **./single_zone_network --t9_0 2 --rho_0 1.e9 --tau 30 --network_xml ../nucnet-tools-code/data_pub/my_net.xml --nuc_xpath "[z <= 60]" --iterative_solver gmres --iterative_t9 2 --output_xml out2.xml --init_mass_frac "{h1; 0.7}" "{he4; 0.2}" "{c12; 0.1}" --xml_steps 50 --t_end 1000. --rate_mod "{[.]; [count(reactant[. = 'he4']) = 3]; 0.1}** *out3.xml* was generated with the single-zone network code with the accompanying master.h and net_rp.xml, and the command: **./single_zone_network --t9_0 2 --rho_0 1.e9 --tau 30 --network_xml net_rp.xml --nuc_xpath "[z <= 60]" --iterative_solver gmres --iterative_t9 2 --output_xml out3.xml --init_mass_frac "{h1; 0.7}" "{he4; 0.2}" "{c12; 0.1}" --xml_steps 50 --t_end 1000.**