Once the single-zone nucleosynthesis results are available, make the data for Figure 8. To do so, change into your data directory and retrieve the mass fractions of interest as a function of time dump:
This creates a text file n_h1_n14_n15_o18_f18_608_1.0e51_reference.txt, which contains in column 1 the time dump number and in columns 2-7 the final mass fractions of the requested species for the single-zone nucleosynthesis calculation of zone 608 in the E = 1.0 B explosion of s15a28. It is also necessary to get the time for each time dump:
To make the first panel of Figure 8, plot column 2 (neutrons) and column 3 (1H) of n_h1_n14_n15_o18_f18_608_1.0e51_reference.txt versus column 2 of time_608_1.0e51_reference.txt. To make the second panel of Figure 8, plot column 4 (14N), column 5 (15N), column 6 (18O), and column 7 (18F) of n_h1_n14_n15_o18_f18_608_1.0e51_reference.txt versus column 2 of time_608_1.0e51_reference.txt.