anchor.diff:                  call write_qmc_var(io_unit, time_avg_psip_list_sq(2,i))
                   write (io_unit,'()')
               end do
  +            !write (io_unit, '(1X, "Selections",/)')
  +            !do i = 1, nstates_select
  +            !    call write_qmc_var(io_unit, selections(1,i))
  +            !    call write_qmc_var(io_unit, selections(2,i))
  +            !    call write_qmc_var(io_unit, selections(3,i))
  +            !    write (io_unit,'()')
  +            !end do
           end if
   
           call dSFMT_t_to_dSFMT_state_t(rng(0), qs%rng_state)
  
positive.code:         cluster%nexcitors = -1
          deexcit_count = 0
          do i = 0, max_size-min_size-1
              psize = psize + 1.0_p/2_int_64**(i+1)
              if (rand < psize) then
                  ! Found size!
                  cluster%nexcitors = i+min_size
                  cluster%pselect = cluster%pselect/2_int_64**(i+1)
                  exit
              end if
          end do
  
negative.code:                 proj_energy_cycle = cmplx(0.0, 0.0, p)
                  D0_population_cycle = cmplx(0.0, 0.0, p)
                  D0_population_ucc_cycle = 0.0_p
                  !$omp do schedule(dynamic,200) reduction(+:D0_population_cycle,proj_energy_cycle, D0_population_ucc_cycle, nattempts_spawn,ndeath)
                  do iattempt = 1, selection_data%nclusters
                      ! For OpenMP scalability, have this test inside a single loop rather
  