Once the single-zone nucleosynthesis results are available, make the data for Figure 9. To do so, if you have not already done so, compile the code that computes the steady-state:
Now you can run the code:
This creates a text file f18_steady_state_608_1.0e51_reference.txt, which contains in column 1 the time, in columns 2-7 various total flows for 18F, in column 8 the steady-state abundance of 18F, and in column 9 the network abundance of 18F for the single-zone nucleosynthesis calculation of zone 608 in the E = 1.0 B explosion of s15a28. Plot the product of column 8 and 18 and the product of column 9 and 18 (since mass fraction is ~ abundance times mass number) against column 1. Repeat for 15N by using n15 in place of f18 in the above command for the compute_species_flows code and 15 in place of 18 in converting to mass fraction.