crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 14:15:42 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Small change to nu-fission analog tallies.

Browse source
This commit is contained in:
Paul Romano 2011-11-27 22:55:29 -05:00
parent bf1bbb9845
commit 59532bb00f
1 changed files with 7 additions and 8 deletions
Download patch file Download diff file Expand all files Collapse all files

15
src/tally.f90
View file

@ -515,6 +515,13 @@ contains
bin_energyout = bins(T_ENERGYOUT)
score_index0 = score_index
! Since the creation of fission sites is weighted such that it
! is expected to create n_particles sites, we need to multiply
! the score by keff to get the true nu-fission rate. Otherwise,
! the sum of all nu-fission rates would be ~1.0.
score = keff
! loop over number of particles banked
do k = 1, p % n_bank
! determine outgoing energy from fission bank
@ -526,14 +533,6 @@ contains
! determine scoring index
score_index = sum((bins - 1) * t % stride) + 1
! Since the creation of fission sites is weighted such that
! it is expected to create n_particles sites, we need to
! multiply the score by keff to get the true nu-fission
! rate. Otherwise, the sum of all nu-fission rates would be
! ~1.0.
score = keff
! Add score to tally
call add_to_score(t % scores(score_index, j), score)
end do