Nuclear Data Libraries

General Overview
REACLIB, etc.

Vital.F90
This file contains very important reaction rates such as 22Ne(a,n), 12C(a,g), etc. I changed this file for my 22Ne(a,g)&(a,n) sensitivity purpose. Note! when you change the rate here, you MUST make sure that the setting in ppn_physics.input: the reaction rate of interest is set to "T" (not "F"). (Check after that, networksetup.txt shows the rate is "T" there too)

How to change nuclear reaction rate in NuGrid
In ppn_physics.input, there is a parameter named INDEX_REACLIB, which gives which Reaclib libraries are used (0 = Basel, 1 = JINA, 2 = PURE (no idea what it means). For changing reaction rates of specific reactions, set "RATE_INDEX" and "RATE_FACTOR". You need to find "RATE_INDEX" from the networksetup.txt (which is created when you run PPN first time in the directory). Find the first column of the reaction you want to change. For example, 59Cu(p,g)60Zn in XRB example is 10694. For details of networksetup.txt, see PPN Physics wiki https://github.com/NuGrid/NuPPN/wiki/physics of Nugrid. I believe you can set number of reactions up to 10 by default. If you want to set more, you need to change it in Physics_Knob.F90 (see below).

NuPPN/physics/source directory
- physics_knobs.F90

* "rate_index" * "rate_factor" will change a reaction rate(s) simply by a factor set here. * Note the two variables above are array, whose number of elements are set by "num_rate_factors (10 by default)". * Also INDEX_REACLIB in pph_physics.input is linked to indenx_reaclib here.

- reaclib.F90

* This shows what files are used for rate calculations.

case(0) reacfile = '../NPDATA/REACLIB/reaclib.nosmo'

case(1) reacfile = '../NPDATA/REACLIB/20081109ReaclibV0.5'

case(2) reacfile = '../NPDATA/REACLIB/20120510ReaclibV1.1'

case(3) reacfile = '../NPDATA/REACLIB/results01111258'   <- This is Reaclib ver.2 (see Physics wiki https://github.com/NuGrid/NuPPN/wiki/physics)