From 891ad670b7d95e13bd6ea6672fd75b68e30a5cd7 Mon Sep 17 00:00:00 2001 From: Paul Romano Date: Thu, 15 May 2014 23:24:49 -0400 Subject: [PATCH] Fix bug when using all with a void material present. --- src/tally.F90 | 11 +++++++---- 1 file changed, 7 insertions(+), 4 deletions(-) diff --git a/src/tally.F90 b/src/tally.F90 index b731c1ae5..31ef3afd3 100644 --- a/src/tally.F90 +++ b/src/tally.F90 @@ -498,7 +498,9 @@ contains filter_index = sum((matching_bins(1:t%n_filters) - 1) * t % stride) + 1 if (t % all_nuclides) then - call score_all_nuclides(p, i_tally, flux, filter_index) + if (p % material /= MATERIAL_VOID) then + call score_all_nuclides(p, i_tally, flux, filter_index) + end if else NUCLIDE_BIN_LOOP: do k = 1, t % n_nuclide_bins @@ -1167,9 +1169,10 @@ contains filter_index = sum((matching_bins(1:t%n_filters) - 1) * t % stride) + 1 if (t % all_nuclides) then - ! Score reaction rates for each nuclide in material - call score_all_nuclides(p, i_tally, flux, filter_index) - + if (p % material /= MATERIAL_VOID) then + ! Score reaction rates for each nuclide in material + call score_all_nuclides(p, i_tally, flux, filter_index) + end if else NUCLIDE_BIN_LOOP: do b = 1, t % n_nuclide_bins ! Get index of nuclide in nuclides array