Added trace information on collision.

src/geometry.F90

@@ -137,7 +137,7 @@ contains
           p % coord % cell = index_cell
 
           ! Show cell information on trace
-          if (trace) then
+          if (verbosity >= 10 .or. trace) then
              message = "    Entering cell " // trim(to_str(c % id))
              call write_message()
           end if

src/particle_header.F90

@@ -52,7 +52,8 @@ module particle_header
      real(8)    :: last_E        ! last energy
 
      ! What event last took place
-     integer    :: event
+     integer    :: event         ! scatter, absorption, fission
+     integer    :: event_nuclide ! index in nuclides array
 
      ! Post-collision physical data
      integer    :: n_bank        ! number of fission sites banked

src/physics.F90

@@ -182,9 +182,18 @@ contains
 
     ! Display information about collision
     if (verbosity >= 10 .or. trace) then
-!!$       message = "    " // trim(reaction_name(MT)) // ". Energy = " // &
-!!$            trim(to_str(p % E * 1e6_8)) // " eV."
-!!$       call write_message()
+       select case (p % event)
+       case (EVENT_SCATTER)
+          message = "    Scattered off of"
+       case (EVENT_FISSION)
+          message = "    Fissioned in"
+       case (EVENT_ABSORB)
+          message = "    Absorbed in"
+       end select
+       message = trim(message) // " " // trim(nuclides(&
+            p % event_nuclide) % name) // ". Energy = " // &
+            trim(to_str(p % E * 1e6_8)) // " eV." 
+       call write_message()
     end if
 
     ! check for very low energy
@@ -262,6 +271,9 @@ contains
     IE  =  micro_xs(index_nuclide) % index_grid
     f   =  micro_xs(index_nuclide) % interp_factor
 
+    ! Save which nuclide particle had collision with
+    p % event_nuclide = index_nuclide
+
     ! ==========================================================================
     ! DISAPPEARANCE REACTIONS (ANALOG) OR IMPLICIT CAPTURE (SURVIVAL BIASING)