diff --git a/src/physics.F90 b/src/physics.F90 index 940239fc62..71a5cc6cc3 100644 --- a/src/physics.F90 +++ b/src/physics.F90 @@ -105,6 +105,10 @@ contains end if if (.not. p % alive) return + ! Sample a scattering reaction and determine the secondary energy of the + ! exiting neutron + call scatter(p, i_nuclide) + ! Play russian roulette if survival biasing is turned on if (survival_biasing) then @@ -112,11 +116,6 @@ contains if (.not. p % alive) return end if - ! Sample a scattering reaction and determine the secondary energy of the - ! exiting neutron - - call scatter(p, i_nuclide) - end subroutine sample_reaction !=============================================================================== @@ -297,6 +296,7 @@ contains p % last_wgt = p % wgt else p % wgt = ZERO + p % last_wgt = ZERO p % alive = .false. end if end if