To make the data for Figure 2, first create the explosion code output. Next, in the projects/simple_snII/ directory, run the command
(Typical execution time for an x86_64 type processor: less than a minute.) This creates a text file step_526_1.0e51.txt in the directory ../nucnet-tools-code/data_pub/s15a28a28/ giving the temperature T9 = T∕109 K and density (in g/cc) as a function of interior mass coordinate at the time (indicated on the second line of the file) the shock has just crossed zone 526 (as labeled in the original pre-supernova model XML file). The first line of the file gives the total number of zones. Similarly, run the command
Graph column 1 (starting at line 3) of step_526_1.0e51.txt and step_661_1.0e51.txt versus column 2 of mass_pre.txt (divided by 1.9891 × 1033 to convert grams to solar masses), which you previously created.
To plot the t = 0 part of the figure, make the data by typing
Graph column 3 (divided by 109 to convert from T in Kelvins to T9 = T∕109 K) versus column 2 (divided by 1.9891 × 1033) of step_0_1.0e51.txt