Check for elastic scattering with isotropic mu.

This apparently happens for Sm150 in JEFF 3.2.
This commit is contained in:
Paul Romano 2017-04-25 15:55:29 -05:00
parent fa82543697
commit fb9f7e2d42

View file

@ -481,7 +481,11 @@ contains
! Sample scattering angle
select type (dist => rxn % products(1) % distribution(1) % obj)
type is (UncorrelatedAngleEnergy)
mu_cm = dist % angle % sample(E)
if (allocated(dist % angle % energy)) then
mu_cm = dist % angle % sample(E)
else
mu_cm = TWO*prn() - ONE
end if
end select
! Determine direction cosines in CM