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

moved conversion of decay_rate units from tally.F90 to openmc/data/product.py

Browse source
This commit is contained in:
Sam Shaner 2016-09-15 10:43:45 -04:00
parent cc8f36edd6
commit 3020800ba4
2 changed files with 7 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

4
openmc/data/product.py
View file

@ -98,7 +98,9 @@ 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)
self._decay_rate = decay_rate
# Convert the decay rate from units of inverse shakes to inverse seconds
self._decay_rate = decay_rate * 1.e-8
@distribution.setter
def distribution(self, distribution):

10
src/tally.F90
View file

@ -708,7 +708,7 @@ contains
score = p % absorb_wgt * yield * &
micro_xs(p % event_nuclide) % fission &
/ micro_xs(p % event_nuclide) % absorption &
* rxn % products(1 + d) % decay_rate * 1.e8_8
* rxn % products(1 + d) % decay_rate
end associate
! Tally to bin
@ -733,7 +733,7 @@ contains
! and not the MAX_DELAYED_GROUPS constant for this loop.
do d = 1, size(rxn % products) - 2
score = score + rxn % products(1 + d) % decay_rate * 1.e8_8 * &
score = score + rxn % products(1 + d) % decay_rate * &
p % absorb_wgt * micro_xs(p % event_nuclide) % fission *&
nuclides(p % event_nuclide) % nu(E, EMISSION_DELAYED, d)&
/ micro_xs(p % event_nuclide) % absorption
@ -769,11 +769,9 @@ contains
associate (rxn => nuclides(p % event_nuclide) % &
reactions(nuclides(p % event_nuclide) % index_fission(1)))
! determine score based on bank site weight and keff. Note that
! the units of the decay rate have been converted from inverse
! shakes to inverse seconds (1 shake = 1.e-8 seconds)
! determine score based on bank site weight and keff.
score = score + keff * fission_bank(n_bank - p % n_bank + k) &
% wgt * rxn % products(1 + g) % decay_rate * 1.e8_8
% wgt * rxn % products(1 + g) % decay_rate
end associate
! if the delayed group filter is present, tally to corresponding