Neutron sources - Floating-point exception

Dear FLUKA experts,

I am reaching out for assistance with a problem I am encountering while using a custom source.f file to simulate AmBe, Cf-252 and moderated Cf sources. The file compiles successfully without any errors or warnings, but when I attempt to use it in my simulation, the run crashes with an error that appears to be linked to a division by zero.

I suspect that the issue might be related to the energy values provided in the input. Could it be that extremely low energy values are being interpreted as zero, causing this problem?

I have reviewed the code to ensure no explicit divisions by zero, but I might be overlooking something.

Has anyone suggestions on how to debug and resolve this issue?

Thank you in advance for your help!
Best regards,
Antonella

phantom test .flair (46.5 KB)
MT96_2024_cf_ambe_cfmod_WIP.f (13.4 KB)
AmBe_0901_ph_water.inp (5.4 KB)

Dear Antonella,

the problem is that you load the cumulative distribution functions into the arrays CF, AMBE and CFMOD, but you try to sample from the uninitialized arrays CUMUCF, CUMUAM and CUMCFMOD.

Cheers,
David

Dear David,
thank you very much for your help!

Cheers,
Antonella