anchor.diff:      ! Array to hold the spawned times of the walkers
   real(p), allocatable :: spawn_times(:) ! (spawned_walker_length)
   
  +
   contains
  -    subroutine do_tc_fciqmc()
  +    subroutine do_tc_fciqmc(spawner)
           
           integer :: nspawned, ireport, idet, iparticle
           integer, allocatable :: current_pos(:) ! (0:max(1,nprocs-1))
  
positive.code:                 ! loop over each walker on the determinant
                  do iparticle = 1, abs(walker_population(1,idet))
  
                      time = 0
                      do
                          time = time + timestep(R)
                          if ( time > t_barrier ) exit
  
negative.code:                 if(all(current_pos == spawning_head+1)) exit
                  
              end do
  
              do_annihilation
              if(vary_shift) update_shift
              print_out
  
          end do
  
      end subroutine do_tc_fciqmc
  
      
      subroutine ct_spawn_real(d, K_ii, parent_sgn, matel, nspawned, connection)
      
          ! randomly select a (valid) excitation from the current determinant
          ! stored in "d" for the Hubbard model in realspace
          !
          ! In: 
          !    d: info on current determinant that we will spawn from.
          !    R_ii: the diagonal Hamiltonian matrix element for the determinant d
          !    parent_sgn: sgn on the parent determinant (i.e. +ve or -ve integer)
          !
          ! Out:
  