Use MT=101 for "absorption" in ESZ block for IncidentNeutron.from_ace

This commit is contained in:
Paul Romano 2018-11-12 07:52:07 -06:00
parent a6a95b2d33
commit d1e427ed9a

View file

@ -641,10 +641,10 @@ class IncidentNeutron(EqualityMixin):
data.reactions[1] = total
if np.count_nonzero(absorption_xs) > 0:
absorption = Reaction(27)
absorption = Reaction(101)
absorption.xs[strT] = Tabulated1D(energy, absorption_xs)
absorption.redundant = True
data.reactions[27] = absorption
data.reactions[101] = absorption
# Read each reaction
n_reaction = ace.nxs[4] + 1