Address some of Paul's comments

This commit is contained in:
guillaume 2019-03-17 17:11:12 -04:00
parent 66f7b559e8
commit 707df81242

View file

@ -9,14 +9,14 @@ GROUP_STRUCTURES = {}
"""Dictionary of commonly used energy group structures:
- "CASMO-X" (where X is 2, 4, 8, 16, 25, 40 or 70) from the CASMO_ lattice
physics code
- "XMAS-172" designed for LWR analysis (Sartori 1990, Santamarina et al. 2004)
- "SHEM-361" designed for LWR analysis to eliminate self-shielding calculations
- "XMAS-172_" designed for LWR analysis (Sartori_ 1990, Santamarina et al. 2004)
- "SHEM-361_" designed for LWR analysis to eliminate self-shielding calculations
of thermal resonances (Hfaiedh et al. 2005, Santamarina 2007, Hebert et al. 2008)
- activation_ energy group structures "VITAMIN-J-175", "TRIPOLI-315",
"CCFE-709_" and "UKAEA-1102_"
.. _CASMO: https://www.studsvik.com/SharepointFiles/CASMO-5%20Development%20and%20Applications.pdf
.. _XMAS: https://www-nds.iaea.org/wimsd/energy.htm
.. _XMAS-172: https://www-nds.iaea.org/wimsd/energy.htm
.. _SHEM-361: https://www.polymtl.ca/merlin/libraries.htm
.. _activation: https://fispact.ukaea.uk/wiki/Keyword:GETXS
.. _CCFE-709: https://fispact.ukaea.uk/wiki/CCFE-709_group_structure
@ -24,14 +24,14 @@ of thermal resonances (Hfaiedh et al. 2005, Santamarina 2007, Hebert et al. 2008
"""
GROUP_STRUCTURES['CASMO-2'] = np.array([
0., 6.25e-1, 2.e7])
0., 6.25e-1, 2.e7])
GROUP_STRUCTURES['CASMO-4'] = np.array([
0., 6.25e-1, 5.53e3, 8.21e5, 2.e7])
0., 6.25e-1, 5.53e3, 8.21e5, 2.e7])
GROUP_STRUCTURES['CASMO-8'] = np.array([
0., 5.8e-2, 1.4e-1, 2.8e-1, 6.25e-1, 4., 5.53e3, 8.21e5, 2.e7])
0., 5.8e-2, 1.4e-1, 2.8e-1, 6.25e-1, 4., 5.53e3, 8.21e5, 2.e7])
GROUP_STRUCTURES['CASMO-16'] = np.array([
0., 3.e-2, 5.8e-2, 1.4e-1, 2.8e-1, 3.5e-1, 6.25e-1, 8.5e-1,
9.72e-1, 1.02, 1.097, 1.15, 1.3, 4., 5.53e3, 8.21e5, 2.e7])
0., 3.e-2, 5.8e-2, 1.4e-1, 2.8e-1, 3.5e-1, 6.25e-1, 8.5e-1,
9.72e-1, 1.02, 1.097, 1.15, 1.3, 4., 5.53e3, 8.21e5, 2.e7])
GROUP_STRUCTURES['CASMO-25'] = np.array([
0., 3.e-2, 5.8e-2, 1.4e-1, 2.8e-1, 3.5e-1, 6.25e-1, 9.72e-1, 1.02, 1.097,
1.15, 1.855, 4., 9.877, 1.5968e1, 1.4873e2, 5.53e3, 9.118e3, 1.11e5, 5.e5,
@ -166,79 +166,79 @@ GROUP_STRUCTURES['TRIPOLI-315,'] = np.array([
1.105e7, 1.162e7, 1.221e7, 1.284e7, 1.350e7, 1.384e7, 1.419e7,
1.455e7, 1.492e7, 1.568e7, 1.649e7, 1.691e7, 1.733e7, 1.964e7])
GROUP_STRUCTURES['SHEM-361'] = np.array([
2.49990e-03, 4.55602e-03, 7.14526e-03, 1.04505e-02, 1.48300e-02,
2.00104e-02, 2.49394e-02, 2.92989e-02, 3.43998e-02, 4.02999e-02,
4.73019e-02, 5.54982e-02, 6.51999e-02, 7.64969e-02, 8.97968e-02,
1.04298e-01, 1.19995e-01, 1.37999e-01, 1.61895e-01, 1.90005e-01,
2.09610e-01, 2.31192e-01, 2.54997e-01, 2.79989e-01, 3.05012e-01,
3.25008e-01, 3.52994e-01, 3.90001e-01, 4.31579e-01, 4.75017e-01,
5.20011e-01, 5.54990e-01, 5.94993e-01, 6.24999e-01, 7.19999e-01,
8.00371e-01, 8.80024e-01, 9.19978e-01, 9.44022e-01, 9.63960e-01,
9.81959e-01, 9.96501e-01, 1.00904e+00, 1.02101e+00, 1.03499e+00,
1.07799e+00, 1.09198e+00, 1.10395e+00, 1.11605e+00, 1.12997e+00,
1.14797e+00, 1.16999e+00, 1.21397e+00, 1.25094e+00, 1.29304e+00,
1.33095e+00, 1.38098e+00, 1.41001e+00, 1.44397e+00, 1.51998e+00,
1.58803e+00, 1.66895e+00, 1.77997e+00, 1.90008e+00, 1.98992e+00,
2.07010e+00, 2.15695e+00, 2.21709e+00, 2.27299e+00, 2.33006e+00,
2.46994e+00, 2.55000e+00, 2.59009e+00, 2.62005e+00, 2.64004e+00,
2.70012e+00, 2.71990e+00, 2.74092e+00, 2.77512e+00, 2.88405e+00,
3.14211e+00, 3.54307e+00, 3.71209e+00, 3.88217e+00, 4.00000e+00,
4.21983e+00, 4.30981e+00, 4.41980e+00, 4.76785e+00, 4.93323e+00,
5.10997e+00, 5.21008e+00, 5.32011e+00, 5.38003e+00, 5.41025e+00,
5.48817e+00, 5.53004e+00, 5.61979e+00, 5.72015e+00, 5.80021e+00,
5.96014e+00, 6.05991e+00, 6.16011e+00, 6.28016e+00, 6.35978e+00,
6.43206e+00, 6.48178e+00, 6.51492e+00, 6.53907e+00, 6.55609e+00,
6.57184e+00, 6.58829e+00, 6.60611e+00, 6.63126e+00, 6.71668e+00,
6.74225e+00, 6.75981e+00, 6.77605e+00, 6.79165e+00, 6.81070e+00,
6.83526e+00, 6.87021e+00, 6.91778e+00, 6.99429e+00, 7.13987e+00,
7.38015e+00, 7.60035e+00, 7.73994e+00, 7.83965e+00, 7.97008e+00,
8.13027e+00, 8.30032e+00, 8.52407e+00, 8.67369e+00, 8.80038e+00,
8.97995e+00, 9.14031e+00, 9.50002e+00, 1.05793e+01, 1.08038e+01,
1.10529e+01, 1.12694e+01, 1.15894e+01, 1.17094e+01, 1.18153e+01,
1.19795e+01, 1.21302e+01, 1.23086e+01, 1.24721e+01, 1.26000e+01,
1.33297e+01, 1.35460e+01, 1.40496e+01, 1.42505e+01, 1.44702e+01,
1.45952e+01, 1.47301e+01, 1.48662e+01, 1.57792e+01, 1.60498e+01,
1.65501e+01, 1.68305e+01, 1.74457e+01, 1.75648e+01, 1.77590e+01,
1.79591e+01, 1.90848e+01, 1.91997e+01, 1.93927e+01, 1.95974e+01,
2.00734e+01, 2.02751e+01, 2.04175e+01, 2.05199e+01, 2.06021e+01,
2.06847e+01, 2.07676e+01, 2.09763e+01, 2.10604e+01, 2.11448e+01,
2.12296e+01, 2.13360e+01, 2.14859e+01, 2.17018e+01, 2.20011e+01,
2.21557e+01, 2.23788e+01, 2.25356e+01, 2.46578e+01, 2.78852e+01,
3.16930e+01, 3.30855e+01, 3.45392e+01, 3.56980e+01, 3.60568e+01,
3.64191e+01, 3.68588e+01, 3.73038e+01, 3.77919e+01, 3.87874e+01,
3.97295e+01, 4.12270e+01, 4.21441e+01, 4.31246e+01, 4.41721e+01,
4.52904e+01, 4.62053e+01, 4.75173e+01, 4.92591e+01, 5.17847e+01,
5.29895e+01, 5.40600e+01, 5.70595e+01, 5.99250e+01, 6.23083e+01,
6.36306e+01, 6.45923e+01, 6.50460e+01, 6.55029e+01, 6.58312e+01,
6.61612e+01, 6.64929e+01, 6.68261e+01, 6.90682e+01, 7.18869e+01,
7.35595e+01, 7.63322e+01, 7.93679e+01, 8.39393e+01, 8.87741e+01,
9.33256e+01, 9.73287e+01, 1.00594e+02, 1.01098e+02, 1.01605e+02,
1.02115e+02, 1.03038e+02, 1.05646e+02, 1.10288e+02, 1.12854e+02,
1.15480e+02, 1.16524e+02, 1.17577e+02, 1.20554e+02, 1.26229e+02,
1.32701e+02, 1.39504e+02, 1.46657e+02, 1.54176e+02, 1.63056e+02,
1.67519e+02, 1.75229e+02, 1.83295e+02, 1.84952e+02, 1.86251e+02,
1.87559e+02, 1.88877e+02, 1.90204e+02, 1.93078e+02, 1.95996e+02,
2.00958e+02, 2.12108e+02, 2.24325e+02, 2.35590e+02, 2.41796e+02,
2.56748e+02, 2.68297e+02, 2.76468e+02, 2.84888e+02, 2.88327e+02,
2.95922e+02, 3.19928e+02, 3.35323e+02, 3.53575e+02, 3.71703e+02,
3.90760e+02, 4.19094e+02, 4.53999e+02, 5.01746e+02, 5.39204e+02,
5.77146e+02, 5.92941e+02, 6.00099e+02, 6.12834e+02, 6.46837e+02,
6.77287e+02, 7.48517e+02, 8.32218e+02, 9.09681e+02, 9.82494e+02,
1.06432e+03, 1.13467e+03, 1.34358e+03, 1.58620e+03, 1.81183e+03,
2.08410e+03, 2.39729e+03, 2.70024e+03, 2.99618e+03, 3.48107e+03,
4.09735e+03, 5.00451e+03, 6.11252e+03, 7.46585e+03, 9.11881e+03,
1.11377e+04, 1.36037e+04, 1.48997e+04, 1.62005e+04, 1.85847e+04,
2.26994e+04, 2.49991e+04, 2.61001e+04, 2.73944e+04, 2.92810e+04,
3.34596e+04, 3.69786e+04, 4.08677e+04, 4.99159e+04, 5.51656e+04,
6.73794e+04, 8.22974e+04, 9.46645e+04, 1.15624e+05, 1.22773e+05,
1.40000e+05, 1.64999e+05, 1.95008e+05, 2.30014e+05, 2.67826e+05,
3.20646e+05, 3.83884e+05, 4.12501e+05, 4.56021e+05, 4.94002e+05,
5.78443e+05, 7.06511e+05, 8.60006e+05, 9.51119e+05, 1.05115e+06,
1.16205e+06, 1.28696e+06, 1.33694e+06, 1.40577e+06, 1.63654e+06,
1.90139e+06, 2.23130e+06, 2.72531e+06, 3.32871e+06, 4.06569e+06,
4.96585e+06, 6.06530e+06, 6.70319e+06, 7.40817e+06, 8.18730e+06,
9.04836e+06, 9.99999e+06, 1.16183e+07, 1.38403e+07, 1.49182e+07,
1.96403e+07])
0.00000e+00, 2.49990e-03, 4.55602e-03, 7.14526e-03, 1.04505e-02,
1.48300e-02, 2.00104e-02, 2.49394e-02, 2.92989e-02, 3.43998e-02,
4.02999e-02, 4.73019e-02, 5.54982e-02, 6.51999e-02, 7.64969e-02,
8.97968e-02, 1.04298e-01, 1.19995e-01, 1.37999e-01, 1.61895e-01,
1.90005e-01, 2.09610e-01, 2.31192e-01, 2.54997e-01, 2.79989e-01,
3.05012e-01, 3.25008e-01, 3.52994e-01, 3.90001e-01, 4.31579e-01,
4.75017e-01, 5.20011e-01, 5.54990e-01, 5.94993e-01, 6.24999e-01,
7.19999e-01, 8.00371e-01, 8.80024e-01, 9.19978e-01, 9.44022e-01,
9.63960e-01, 9.81959e-01, 9.96501e-01, 1.00904e+00, 1.02101e+00,
1.03499e+00, 1.07799e+00, 1.09198e+00, 1.10395e+00, 1.11605e+00,
1.12997e+00, 1.14797e+00, 1.16999e+00, 1.21397e+00, 1.25094e+00,
1.29304e+00, 1.33095e+00, 1.38098e+00, 1.41001e+00, 1.44397e+00,
1.51998e+00, 1.58803e+00, 1.66895e+00, 1.77997e+00, 1.90008e+00,
1.98992e+00, 2.07010e+00, 2.15695e+00, 2.21709e+00, 2.27299e+00,
2.33006e+00, 2.46994e+00, 2.55000e+00, 2.59009e+00, 2.62005e+00,
2.64004e+00, 2.70012e+00, 2.71990e+00, 2.74092e+00, 2.77512e+00,
2.88405e+00, 3.14211e+00, 3.54307e+00, 3.71209e+00, 3.88217e+00,
4.00000e+00, 4.21983e+00, 4.30981e+00, 4.41980e+00, 4.76785e+00,
4.93323e+00, 5.10997e+00, 5.21008e+00, 5.32011e+00, 5.38003e+00,
5.41025e+00, 5.48817e+00, 5.53004e+00, 5.61979e+00, 5.72015e+00,
5.80021e+00, 5.96014e+00, 6.05991e+00, 6.16011e+00, 6.28016e+00,
6.35978e+00, 6.43206e+00, 6.48178e+00, 6.51492e+00, 6.53907e+00,
6.55609e+00, 6.57184e+00, 6.58829e+00, 6.60611e+00, 6.63126e+00,
6.71668e+00, 6.74225e+00, 6.75981e+00, 6.77605e+00, 6.79165e+00,
6.81070e+00, 6.83526e+00, 6.87021e+00, 6.91778e+00, 6.99429e+00,
7.13987e+00, 7.38015e+00, 7.60035e+00, 7.73994e+00, 7.83965e+00,
7.97008e+00, 8.13027e+00, 8.30032e+00, 8.52407e+00, 8.67369e+00,
8.80038e+00, 8.97995e+00, 9.14031e+00, 9.50002e+00, 1.05793e+01,
1.08038e+01, 1.10529e+01, 1.12694e+01, 1.15894e+01, 1.17094e+01,
1.18153e+01, 1.19795e+01, 1.21302e+01, 1.23086e+01, 1.24721e+01,
1.26000e+01, 1.33297e+01, 1.35460e+01, 1.40496e+01, 1.42505e+01,
1.44702e+01, 1.45952e+01, 1.47301e+01, 1.48662e+01, 1.57792e+01,
1.60498e+01, 1.65501e+01, 1.68305e+01, 1.74457e+01, 1.75648e+01,
1.77590e+01, 1.79591e+01, 1.90848e+01, 1.91997e+01, 1.93927e+01,
1.95974e+01, 2.00734e+01, 2.02751e+01, 2.04175e+01, 2.05199e+01,
2.06021e+01, 2.06847e+01, 2.07676e+01, 2.09763e+01, 2.10604e+01,
2.11448e+01, 2.12296e+01, 2.13360e+01, 2.14859e+01, 2.17018e+01,
2.20011e+01, 2.21557e+01, 2.23788e+01, 2.25356e+01, 2.46578e+01,
2.78852e+01, 3.16930e+01, 3.30855e+01, 3.45392e+01, 3.56980e+01,
3.60568e+01, 3.64191e+01, 3.68588e+01, 3.73038e+01, 3.77919e+01,
3.87874e+01, 3.97295e+01, 4.12270e+01, 4.21441e+01, 4.31246e+01,
4.41721e+01, 4.52904e+01, 4.62053e+01, 4.75173e+01, 4.92591e+01,
5.17847e+01, 5.29895e+01, 5.40600e+01, 5.70595e+01, 5.99250e+01,
6.23083e+01, 6.36306e+01, 6.45923e+01, 6.50460e+01, 6.55029e+01,
6.58312e+01, 6.61612e+01, 6.64929e+01, 6.68261e+01, 6.90682e+01,
7.18869e+01, 7.35595e+01, 7.63322e+01, 7.93679e+01, 8.39393e+01,
8.87741e+01, 9.33256e+01, 9.73287e+01, 1.00594e+02, 1.01098e+02,
1.01605e+02, 1.02115e+02, 1.03038e+02, 1.05646e+02, 1.10288e+02,
1.12854e+02, 1.15480e+02, 1.16524e+02, 1.17577e+02, 1.20554e+02,
1.26229e+02, 1.32701e+02, 1.39504e+02, 1.46657e+02, 1.54176e+02,
1.63056e+02, 1.67519e+02, 1.75229e+02, 1.83295e+02, 1.84952e+02,
1.86251e+02, 1.87559e+02, 1.88877e+02, 1.90204e+02, 1.93078e+02,
1.95996e+02, 2.00958e+02, 2.12108e+02, 2.24325e+02, 2.35590e+02,
2.41796e+02, 2.56748e+02, 2.68297e+02, 2.76468e+02, 2.84888e+02,
2.88327e+02, 2.95922e+02, 3.19928e+02, 3.35323e+02, 3.53575e+02,
3.71703e+02, 3.90760e+02, 4.19094e+02, 4.53999e+02, 5.01746e+02,
5.39204e+02, 5.77146e+02, 5.92941e+02, 6.00099e+02, 6.12834e+02,
6.46837e+02, 6.77287e+02, 7.48517e+02, 8.32218e+02, 9.09681e+02,
9.82494e+02, 1.06432e+03, 1.13467e+03, 1.34358e+03, 1.58620e+03,
1.81183e+03, 2.08410e+03, 2.39729e+03, 2.70024e+03, 2.99618e+03,
3.48107e+03, 4.09735e+03, 5.00451e+03, 6.11252e+03, 7.46585e+03,
9.11881e+03, 1.11377e+04, 1.36037e+04, 1.48997e+04, 1.62005e+04,
1.85847e+04, 2.26994e+04, 2.49991e+04, 2.61001e+04, 2.73944e+04,
2.92810e+04, 3.34596e+04, 3.69786e+04, 4.08677e+04, 4.99159e+04,
5.51656e+04, 6.73794e+04, 8.22974e+04, 9.46645e+04, 1.15624e+05,
1.22773e+05, 1.40000e+05, 1.64999e+05, 1.95008e+05, 2.30014e+05,
2.67826e+05, 3.20646e+05, 3.83884e+05, 4.12501e+05, 4.56021e+05,
4.94002e+05, 5.78443e+05, 7.06511e+05, 8.60006e+05, 9.51119e+05,
1.05115e+06, 1.16205e+06, 1.28696e+06, 1.33694e+06, 1.40577e+06,
1.63654e+06, 1.90139e+06, 2.23130e+06, 2.72531e+06, 3.32871e+06,
4.06569e+06, 4.96585e+06, 6.06530e+06, 6.70319e+06, 7.40817e+06,
8.18730e+06, 9.04836e+06, 9.99999e+06, 1.16183e+07, 1.38403e+07,
1.49182e+07, 1.96403e+07])
GROUP_STRUCTURES['CCFE-709'] = np.array([
1.e-5, 1.0471e-5, 1.0965e-5, 1.1482e-5, 1.2023e-5,
1.2589e-5, 1.3183e-5, 1.3804e-5, 1.4454e-5, 1.5136e-5,