From 59532bb00f9d4a73b22ee2a87c8019d8e8f236cd Mon Sep 17 00:00:00 2001 From: Paul Romano Date: Sun, 27 Nov 2011 22:55:29 -0500 Subject: [PATCH] Small change to nu-fission analog tallies. --- src/tally.f90 | 15 +++++++-------- 1 file changed, 7 insertions(+), 8 deletions(-) diff --git a/src/tally.f90 b/src/tally.f90 index 0342619405..b4fe31cc9f 100644 --- a/src/tally.f90 +++ b/src/tally.f90 @@ -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