crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 22:26:08 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Make sure isotropic angular distribution gets created from ACE files

Browse source
This commit is contained in:
Paul Romano 2020-01-14 07:45:09 -06:00
parent d4ad3cfb75
commit d0d1b1db12
1 changed files with 8 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

12
openmc/data/reaction.py
View file

@ -1073,11 +1073,15 @@ class Reaction(EqualityMixin):
if i_reaction < ace.nxs[5] + 1:
# Check if angular distribution data exist
loc = int(ace.xss[ace.jxs[8] + i_reaction])
if loc <= 0:
# Angular distribution is either given as part of a product
# angle-energy distribution or is not given at all (in which
# case isotropic scattering is assumed)
if loc < 0:
# Angular distribution is given as part of a product
# angle-energy distribution
angle_dist = None
elif loc == 0:
# Angular distribution is isotropic
energy = [0.0, grid[-1]]
mu = Uniform(-1., 1.)
angle_dist = AngleDistribution(energy, [mu, mu])
else:
angle_dist = AngleDistribution.from_ace(ace, ace.jxs[9], loc)