anchor.diff:  
                   ! Find the spawning rate and add to the running
                   ! total.
  -                rspawn = rspawn + real(sum(abs(qmc_spawn%sdata(1,:))))/nattempts
  -
  -                call simple_annihilation(qmc_spawn, psip_list%walker_population)
  +                rspawn = rspawn + real(sum(abs(spawn%sdata(1,:))))/nattempts
  +
  +                call simple_annihilation(spawn, psip_list%walker_population)
   
               end do
   
  
positive.code:         !    psip_list: walker_t object for storing psips sampling the Hilbert
          !       space and related information.  Allocated and initialised on
          !       output.
  
          use parallel, only: nprocs, parent
          use checking, only: check_allocate
  
negative.code:         integer(i0), allocatable :: dets(:,:)
          real(p) :: H0i, Hii
          type(walker_t) :: psip_list
          type(spawn_t) :: spawn
  
          call init_simple_fciqmc(sys, qmc_in, reference, restart_in%read_restart, ndets, dets, ref_det, psip_list, spawn)
  
          if (parent) call rng_init_info(qmc_in%seed+iproc)
          call dSFMT_init(qmc_in%seed+iproc, 50000, rng)
  