Problems with OECD/NEA JEFF-32 ACE files

[emilf]

I’m trying to use JEFF-32 ace file from OECD/NEA. Serpent v29 gives the same error (so far for all Ni and some Fe nuclides ).
???

Nuclide 28060.03c -- nickel 60 at 300K (Ni-60)

***** Fri Oct 27 16:23:08 2017:

- MPI task = 1
- OpenMP thread = 0
- RNG parent seed = 1509114145

Fatal error in function ProcessMuDistributions:

Error in angular distribution (mt 2)

Simulation aborted.

Nuclide 28061.03c -- nickel 61 at 300K (Ni-61)

[Jaakko Leppänen]

I've made a note in the subroutine that in some machines NJOY-2012 processes some angular distributions in JEFF-3.2 in such way that the highest value of the PDF of the scattering cosine differs from 1.0. I've been doing some extensive validations with JEFF-3.2 libraries processed with NJOY-2012, so this is probably a numerical problem in NJOY. In any case, if the data is not completely screwed up, it would be good if Serpent could use it.

Try replacing line 380 in processmudistributions.c:

Code: Select all

   Die(FUNCTION_NAME, "Error in angular distribution (mt %ld)", mt);

with:

Code: Select all

   Warn(FUNCTION_NAME, "Error in angular distribution (mt %ld %E %E)", mt, RDB[l1 - NP], RDB[l1 - 1] - 1.0);

If this is related to the same problem, you should get a warning instead of an error. Could post the results?

[Jaakko Leppänen]

Which data files are you using? The 300K files from:

https://www.oecd-nea.org/dbforms/data/e ... s/jeff_32/

are read and processed without errors.

[emilf]

Yes, the library is from the OECD website. However, replacing line 380 retains the error in my case.

Nuclide 26058.03c -- iron 58 at 300K (Fe-58)

***** Wed Nov 1 15:31:31 2017:

- MPI task = 0
- OpenMP thread = 0
- RNG parent seed = 1509546671

Fatal error in function ProcessMuDistributions:

Error in angular distribution (mt 2)

Simulation aborted.

[Jaakko Leppänen]

There's another similar call on line 347. Could you try the same changes there?

[emilf]

Yes, I saw it and replaced it also. No success.
I got another die call. I will upload the error description

[emilf]

This is the output. So after two warnings from new "Warn" call comes fatal error:

Nuclide 26058.03c -- iron 58 at 300K (Fe-58)

***** Wed Nov 1 19:33:02 2017 (seed = 1509561154, MPI task = 1)
Warning message from function ProcessMuDistributions:

Error in angular distribution (mt 2 -1.000000E+00 -1.000000E+00)


***** Wed Nov 1 19:33:02 2017 (seed = 1509561154, MPI task = 1)
Warning message from function ProcessMuDistributions:

Error in angular distribution (mt 2 0.000000E+00 -1.000000E+00)

***** Wed Nov 1 19:33:02 2017:

- MPI task = 1
- OpenMP thread = 0
- RNG parent seed = 1509561154

Fatal error in function ProcessMuDistributions:

mixed mu interpolations (mt 2)

Simulation aborted.

[Jaakko Leppänen]

Can you give the direct link to the data library that you used?

[emilf]

Looks exactly like yours (except https )
http://www.oecd-nea.org/dbforms/data/ev ... s/jeff_32/

[Jaakko Leppänen]

Which file?