Page 1 of 2

Problems with OECD/NEA JEFF-32 ACE files

Posted: Fri Oct 27, 2017 6:00 pm
by 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)

Re: Problems with OECD/NEA JEFF-32 ACE files

Posted: Fri Oct 27, 2017 9:45 pm
by 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?

Re: Problems with OECD/NEA JEFF-32 ACE files

Posted: Fri Oct 27, 2017 10:36 pm
by 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.

Re: Problems with OECD/NEA JEFF-32 ACE files

Posted: Wed Nov 01, 2017 5:35 pm
by 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.

Re: Problems with OECD/NEA JEFF-32 ACE files

Posted: Wed Nov 01, 2017 8:40 pm
by Jaakko Leppänen
There's another similar call on line 347. Could you try the same changes there?

Re: Problems with OECD/NEA JEFF-32 ACE files

Posted: Wed Nov 01, 2017 9:29 pm
by emilf
Yes, I saw it and replaced it also. No success.
I got another die call. I will upload the error description

Re: Problems with OECD/NEA JEFF-32 ACE files

Posted: Wed Nov 01, 2017 9:37 pm
by 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.

Re: Problems with OECD/NEA JEFF-32 ACE files

Posted: Wed Nov 01, 2017 9:49 pm
by Jaakko Leppänen
Can you give the direct link to the data library that you used?

Re: Problems with OECD/NEA JEFF-32 ACE files

Posted: Wed Nov 01, 2017 9:57 pm
by emilf
Looks exactly like yours (except https :))
http://www.oecd-nea.org/dbforms/data/ev ... s/jeff_32/

Re: Problems with OECD/NEA JEFF-32 ACE files

Posted: Wed Nov 01, 2017 11:16 pm
by Jaakko Leppänen
Which file?