From 9079ae7b7478919d2cba685d3cdf491f6bc13738 Mon Sep 17 00:00:00 2001 From: Paul Romano Date: Wed, 5 Oct 2016 15:10:41 -0500 Subject: [PATCH] Ignore obsolete energy transformation matrix in MF=4. --- openmc/data/angle_distribution.py | 11 +++++++++++ 1 file changed, 11 insertions(+) diff --git a/openmc/data/angle_distribution.py b/openmc/data/angle_distribution.py index 0f0105c2a..669951059 100644 --- a/openmc/data/angle_distribution.py +++ b/openmc/data/angle_distribution.py @@ -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']])