From b894ad6ff35f26b739666285c413bec82aef3846 Mon Sep 17 00:00:00 2001 From: Sam Shaner Date: Thu, 15 Sep 2016 16:15:34 -0400 Subject: [PATCH] moved decay rate units conversion to reaction.py --- openmc/data/product.py | 4 +--- openmc/data/reaction.py | 4 +++- 2 files changed, 4 insertions(+), 4 deletions(-) diff --git a/openmc/data/product.py b/openmc/data/product.py index 50f1482213..753888f6df 100644 --- a/openmc/data/product.py +++ b/openmc/data/product.py @@ -98,9 +98,7 @@ class Product(EqualityMixin): def decay_rate(self, decay_rate): cv.check_type('product decay rate', decay_rate, Real) cv.check_greater_than('product decay rate', decay_rate, 0.0, True) - - # Convert the decay rate from units of inverse shakes to inverse seconds - self._decay_rate = decay_rate * 1.e8 + self._decay_rate = decay_rate @distribution.setter def distribution(self, distribution): diff --git a/openmc/data/reaction.py b/openmc/data/reaction.py index 45d668b92f..eb33e1395f 100644 --- a/openmc/data/reaction.py +++ b/openmc/data/reaction.py @@ -107,7 +107,9 @@ def _get_fission_products(ace): for group in range(n_group): delayed_neutron = Product('neutron') delayed_neutron.emission_mode = 'delayed' - delayed_neutron.decay_rate = ace.xss[idx] + + # Convert units of inverse shakes to inverse seconds + delayed_neutron.decay_rate = ace.xss[idx] * 1.e8 group_probability = Tabulated1D.from_ace(ace, idx + 1) if np.all(group_probability.y == group_probability.y[0]):