From 92b58c2b5663c408a345bac337d0d030dac8a5b5 Mon Sep 17 00:00:00 2001 From: Sterling Harper Date: Mon, 18 Apr 2016 19:40:23 -0400 Subject: [PATCH] Avoid divide-by-zero errors with score==ZERO check --- src/tally.F90 | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/src/tally.F90 b/src/tally.F90 index 360bb35bc..43e10986e 100644 --- a/src/tally.F90 +++ b/src/tally.F90 @@ -3039,6 +3039,8 @@ contains logical :: scoring_diff_nuclide real(8) :: dsigT, dsigA, dsigF, cum_dsig + if (score == ZERO) return + associate(deriv => tally_derivs(t % deriv)) select case (deriv % variable) @@ -3347,7 +3349,7 @@ contains end if end associate score = score * (deriv % flux_deriv & - + dsigF * mat % atom_density(l) / material_xs % fission & + + dsigF * mat % atom_density(l) / material_xs % nu_fission& * micro_xs(p % event_nuclide) % nu_fission & / micro_xs(p % event_nuclide) % fission) end associate