NUCGRAPH.md¶

NUCGRAPH numg floor uncert nopt(i), i=1,numg

This keyword creates data output files as .grn and .pln that contain the data and gnuplot plotting instructions for the total and dominant nuclide contributions for a set of requested output graphs. The number of graphs required (numg) is input, and for each graph an option number (nopt(i)) is read. Allowable values for the options are:

1. specific activity
2. specific total heat
3. gamma dose rate
4. specific gamma heat
5. specific beta heat
6. specific ingestion
7. specific inhalation

The topx and topxx inputs from the SORTDOMINANT keyword may be used to limit or expand the number of dominant contributions considered. The floor argument specifies the cut-off for a dominant contribution to be plotted. Only those curves that contribute more than the floor percentage of the total value are included in the graph data.

The uncert argument allows the user to plot the UNCERTAINTY pathways-based uncertainty data (for the total response only) using uncert=1 or the uncertainty may be ignored using uncert=0.

This keyword requires the user to specify a gnuplot style file as ng_style within the files file to generate the full plot script, which may be plotted using

gnuplot <fileroot>.pln

Example usage of the keyword:

 < -- Control phase -- >
 ...
 FISPACT
 * Title of the simulation
 < -- Initial phase -- > 
 ...
 NUCGRAPH 7 0.01 1  1 2 3 4 5 6 7
 ...
 < -- Inventory phase -- > 
 ...

The default .pln created by this keyword should be treated as a guiding template only - the user should modify it to produce visually acceptable output.

Example plots of .grn files (using modified versions of .pln):

Becquerel activity from pure iron after a 2 year irradiation under fusion power plant conditions

Decay heat from pure tungsten after a 2 year irradiation under fusion power plant conditions

Contact gamma dose rate from pure tin after a 2 year irradiation under fusion power plant conditions