To make graphs of explosive nucleosynthesis yields versus time for a given zone, you first need to create explosion code outputs for energies 0.4 B, 0.6 B, 0.8 B, and 1.2 B, as done previously for energy 1.0 B. Once the explosion code output files are ready, retrieve the parameter data for the nucleosynthesis calculations (if you have not already done so):
Now you can run the single-zone code. To run the code for explosion energy 1.0 B, for example, for zone 608, type
(The execution time depends on the particular zone, but, for an x86_64 type processor: ~10-20 minutes.)
To make integrated currents graphs, you will need to indicate the name of the XML file that will contain the currents. To do this, edit the input zone XML file ../nucnet-tools-code/data_pub/simple_snII__zone.xml. Find the lines
and edit them so that they read
Rerun the above calculation, which will take longer because it now is computing and keeping track of the currents. This calculation will create both my_output_608_1.0e51_reference.xml and currents_608_1.0e51_reference.xml in the directory ../nucnet-tools-code/data_pub/s15a28a28/. Change the zone number and/or energy to run other single-zone calculations.
For the plots that include steady-state mass fractions, such as 15N and 18F, versus time, it might be necessary to change the “steps” property in the input zone XML file ../nucnet-tools-code/data_pub/simple_snII__zone.xml from 5 to 1 before the calculation, as shown below:
Saving the mass fractions to the output XML file every 5 steps may miss the sudden rise or fall in the steady-state 15N and 18F mass fractions due to shock passage. Choosing to save every timestep captures the sudden change more accurately at the cost of a larger output file.