To make the data for Figure 33, you need the results from the E = 1.0 B multi-zone nucleosynthesis calculation. Once those are ready, compile the codes in the nuclear_decay project:
First, decay the results from the E = 1.0 B calculation for some time, for example, 1010 years, and then compute the species overabundances:
In the calls to compute_species_overabundances, the input pre-supernova file provides the full list of available species. Also, the calculations that generated s15a28a28_presn.xml used initial Anders and Grevesse Solar System abundances. For this reason, we use the Solar mass fractions in Anders.xml, not Lodders.xml. Plot the appropriate columns of these text files against the interior mass coordinate from column 2 of mass_pre.txt (divided by 1.9891 × 1033 to convert grams to solar masses) that you previously created. Note that the overabundances text files contain all three zone labels (“0” for the other two labels), so be careful in reading these files.