Figure 10

To make Figure 10, follow these steps to run a new calculation for explosion energy 1.0 B but for zone 603. Be sure to create the currents XML file for this calculation, called currents_603_1.0e51_reference.xml, in the directory projects/nucnet-tools-code/data_pub/s15a28a28/. Once these are done, change into the projects/nucnet-tools-code/examples/misc/ directory and compile the codes there:

cd ~/projects/nucnet-tools-code/examples/misc/

make all_misc

Next, run the code to compute the difference in the currents between the zone 603 calculation and the zone 608 calculation:

./diff_current ../../data_pub/s15a28a28/currents_603_1.0e51_reference.xml ../../data_pub/s15a28a28/currents_608_1.0e51_reference.xml ../../data_pub/s15a28a28/currents_603_608_1.0e51_reference.xml

To make Figure 10, change into the projects/nucnet-tools-code/examples/graph/ directory and run the currents graph script on the currents difference XML file:

cd ../graph/

./current_graph.sh ../../data_pub/s15a28a28/currents_603_608_1.0e51_reference.xml "" "" "[z >= 6 and z <= 12 and a - z >= 6 and a - z <= 14]" linear currents_603_608_1_0e51_reference.pdf

Figure 10 is the output pdf file currents_603_608_1_0e51_reference.pdf. The code places the output file in the projects/nucnet-tools-code/examples/graph directory. You may choose to move it to the projects/nucnet-tools-code/data_pub/s15a28a28 directory, where it will reside with the other data files.

[next] [prev] [prev-tail] [front] [up]