anchor.diff:                  !$omp do schedule(dynamic,200) private(nspawned, connection, junk) reduction(+:D0_population_cycle,proj_energy)
                   do iattempt = 1, nattempts
   
  -                    call select_cluster(rng(it), nattempts, D0_pos, cumulative_abs_pops, tot_abs_pop, max_cluster_size, &
  -                                        cdet(it), cluster(it))
  +                    call select_cluster(rng(it), nattempts, D0_normalisation, D0_pos, cumulative_abs_pops, &
  +                                        tot_abs_pop, max_cluster_size, cdet(it), cluster(it))
   
                       if (cluster(it)%excitation_level <= truncation_level+2) then
   
  
positive.code:         use sort, only: insert_sort
  
          integer(lint), intent(in) :: nattempts
          integer, intent(in) :: D0_pos
          integer, intent(in) :: cumulative_excip_pop(:), tot_excip_pop, max_size
          type(dSFMT_t), intent(inout) :: rng
          type(det_info), intent(inout) :: cdet
  
negative.code: 
          integer, allocatable :: cumulative_abs_pops(:)
          integer :: D0_pos, max_cluster_size,tot_abs_pop
          integer :: D0_normalisation
          logical :: hit
          type(bloom_stats_t) :: bloom_stats
  
  