From 40543168aebff63ce17cd9996f26ebefb9df7155 Mon Sep 17 00:00:00 2001
From: nhorelik <nhorelik@mit.edu>
Date: Wed, 10 Apr 2013 16:38:35 -0400
Subject: [PATCH] Changed universe on new coord for lattice outside to be the
 same as previously

---
 src/geometry.F90 | 24 ++++++++++++------------
 1 file changed, 12 insertions(+), 12 deletions(-)

diff --git a/src/geometry.F90 b/src/geometry.F90
index 9dafd75f61..7f5e28ef3c 100644
--- a/src/geometry.F90
+++ b/src/geometry.F90
@@ -245,26 +245,26 @@ contains
             end if
             p % coord % next % uvw = p % coord % uvw
 
-            ! Move particle to next level
-            p % coord => p % coord % next
-
             ! set particle lattice indices
-            p % coord % lattice   = c % fill
-            p % coord % lattice_x = i_x
-            p % coord % lattice_y = i_y
-            p % coord % lattice_z = i_z
+            p % coord % next% lattice   = c % fill
+            p % coord % next% lattice_x = i_x
+            p % coord % next% lattice_y = i_y
+            p % coord % next% lattice_z = i_z
             if (.not. outside_lattice) then
-              p % coord % universe = lat % universes(i_x,i_y,i_z)
+              p % coord % next % universe = lat % universes(i_x,i_y,i_z)
             else
 
-              ! If we define a void universe, we could 
-              p % coord % universe = NONE
+              ! Set universe as the same for subsequent calls to find_cell
+              p % coord % next % universe = p % coord % universe
 
-              ! Set coord cell for distance_to_boundary
-              p % coord % cell = index_cell
+              ! Set coord cell for calls to distance_to_boundary
+              p % coord % next % cell = index_cell
 
             end if
 
+            ! Move particle to next level
+            p % coord => p % coord % next
+
           end if
 
           if (.not. outside_lattice) then