From ae9ef677cb595968c6d15e840db425c573d98c5d Mon Sep 17 00:00:00 2001
From: Paul Romano <romano7@gmail.com>
Date: Sat, 15 Oct 2011 22:53:45 -0400
Subject: [PATCH] Fixed bug in fission events when running an analog
 calculation (no survival biasing). This caused some fission events to not
 kill particles.

---
 src/physics.f90 | 4 +---
 1 file changed, 1 insertion(+), 3 deletions(-)

diff --git a/src/physics.f90 b/src/physics.f90
index 5bfd00fccb..27e4c60063 100644
--- a/src/physics.f90
+++ b/src/physics.f90
@@ -472,6 +472,7 @@ contains
           scattered = .true.
        case (N_FISSION, N_F, N_NF, N_2NF, N_3NF)
           call create_fission_sites(p, index_nuclide, rxn, .true.)
+          p % alive = .false.
        case (N_GAMMA : N_DA)
           call n_absorption(p)
        case default
@@ -767,9 +768,6 @@ contains
     ! increment number of bank sites
     n_bank = min(n_bank + nu, 3*n_particles)
 
-    ! kill original neutron if no survival biasing
-    if (actual_event) p % alive = .false.
-
   end subroutine create_fission_sites
 
 !===============================================================================