njoy / NJOY2016

Nuclear data processing with legacy NJOY
https://www.njoy21.io/NJOY2016
Other
96 stars 85 forks source link

Issue processing MF34 with multiple subsections #102

Closed whaeck closed 5 years ago

whaeck commented 5 years ago

The following email was received from @luca-fiorito-11.

NJOY sets by default the number of subsections ni=1 for MF34. This is not specified in the ENDF-6 format manual, which, on the contrary, allows using more than 1 subsection also for MF34.

Line 1942 of error.f90 in NJOY2016.42 if (mfcov.eq.34) then nc=0 ni=1 endif

For Fe56 in JEFF-33, the covariance between P1/P1 is given in 3 subsections (ni=3)…therefore NJOY fails to run. By commenting out the lines reported above, NJOY completes the calculation.

Still, I think unit tests should be run to see if the change does not disrupt the NJOY regular behaviour.

Cheers, Luca

whaeck commented 5 years ago

For some reason, NJOY seems to reset the values of ni and nc from what it read before.

whaeck commented 5 years ago

Merged in #104