Figure 2

      To make the data for Figure 2, first create the explosion code output. Next, in the projects/simple_snII/ directory, run the command

./t9_rho_vs_mr ../nucnet-tools-code/data_pub/s15a28a28/expl_1.0e51.xml "[@label1 = '526']" > ../nucnet-tools-code/data_pub/s15a28a28/step_526_1.0e51.txt

(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

./t9_rho_vs_mr ../nucnet-tools-code/data_pub/s15a28a28/expl_1.0e51.xml "[@label1 = '661']" > ../nucnet-tools-code/data_pub/s15a28a28/step_661_1.0e51.txt

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

../nucnet-tools-code/examples/analysis/print_properties ../nucnet-tools-code/data_pub/s15a28a28/s15a28a28_presn.xml "mass below" "cell temperature" > ../nucnet-tools-code/data_pub/s15a28a28/step_0_1.0e51.txt

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