crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-27 13:45:36 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

Ignore obsolete energy transformation matrix in MF=4.

Browse source
This commit is contained in:
Paul Romano 2016-10-05 15:10:41 -05:00
parent b3e351474c
commit 9079ae7b74
1 changed files with 11 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

11
openmc/data/angle_distribution.py
View file

@ -1,6 +1,7 @@
from collections import Iterable
from io import StringIO
from numbers import Real
from warnings import warn
import numpy as np
@ -224,13 +225,23 @@ class AngleDistribution(EqualityMixin):
# Read HEAD record
items = get_head_record(file_obj)
lvt = items[2]
ltt = items[3]
# Read CONT record
items = get_cont_record(file_obj)
li = items[2]
nk = items[4]
center_of_mass = (items[3] == 2)
# Check for obsolete energy transformation matrix. If present, just skip
# it and keep reading
if lvt > 0:
warn('Obsolete energy transformation matrix in MF=4 angular '
'distribution.')
for _ in range((nk + 5)//6):
file_obj.readline()
if ltt == 0 and li == 1:
# Purely isotropic
energy = np.array([0., ev.info['energy_max']])