Page 1 of 1

Fatal error in function ProcessPhotonProd

Posted: Tue Jul 09, 2019 12:10 am
by gtdrew
Dear all,

I am attempting to get coupled neutron-photon transport simulations working in Serpent, and I am encountering the following error: "Fatal error in function ProcessPhotonProd, No fission channels found, Simulation aborted". In addition, I get the warning: "Warning message from function ProcessPhotonProd: CDF(NP - 1) = 0.000000E+00". I have the MCNP photon data libraries in my acelib path, and I have the auxiliary photon data in my pdatadir path. I have attempted runs with both the 2.1.30 version (with source code changes from viewtopic.php?f=25&t=2839&p=8402#p8402) and 2.1.31 version, both yielding the same error. I feel like the issue is stemming from an issue with my pdatadir files, but if anyone could provide some wisdom regarding this issue and possible solutions, it would be much appreciated.

Re: Fatal error in function ProcessPhotonProd

Posted: Tue Jul 09, 2019 10:39 am
by Jaakko Leppänen
The subroutine actually handles the photon production reactions in neutron data. Which neutron data libraries are you using?

Re: Fatal error in function ProcessPhotonProd

Posted: Wed Aug 28, 2019 7:50 pm
by MohammedQah
I got this error when I use KERMA data:

Fatal error in function ReadPendfData:

Mismatch in lowest energy of ACE and KERMA datas for nuclide 90232.03c (mt=3318)
ACE: 9.99999999999999939497E-12 KERMA: 3.99441299999999960973E-03


and also many warning messages of negative xs points or negative non-fission KERMAs in data.
any suggestion on that error.

Re: Fatal error in function ProcessPhotonProd

Posted: Thu Aug 29, 2019 2:00 pm
by Riku Tuominen
This error has nothing to do with processphotonprod, so next time please use the correct section of the forum and make a new thread if necessary. In order to fix the error remove the following lines of code from readpendfdata.c:

lines 419-422:
else
Die(FUNCTION_NAME, "Mismatch in lowest energy of ACE and KERMA datas "
"for nuclide %s (mt=%ld) \nACE: %.20E KERMA: %.20E",
GetText(nuc + NUCLIDE_PTR_NAME), mt, E[0], E0[0]);

lines 426-429:
else
Die(FUNCTION_NAME, "Mismatch in highest energy of ACE and KERMA datas "
"for nuclide %s (mt=%ld) \nACE: %.20E KERMA: %.20E",
GetText(nuc + NUCLIDE_PTR_NAME), mt, E[NES - 1], E0[np - 1]);

The issue of negative kerma coefficients is briefly discussed in section 2.2 of the energy deposition paper.