anchor.diff:                  if(qs%psip_list%nstates-nD0_proc == 0) then
                       max_cluster_size = 0
                   else
  -                    max_cluster_size = uccmc_in%pow_trunc
  +                    max_cluster_size = min(sys%nel,qs%ref%max_ex_level+2)
                   end if
   
                   call init_mc_cycle(qs%psip_list, qs%spawn_store%spawn, qs%estimators(1)%nattempts, ndeath, &
  
positive.code:              ![todo] deal with restarting
          end if
  
          call init_particle_t(qmc_in%walker_length, 1, sys%basis%tensor_label_len, qmc_in%real_amplitudes, &
                               qmc_in%real_amplitude_force_32, time_avg_psip_list, io_unit=io_unit)
               time_avg_psip_list%pops(:,1) = qs%psip_list%pops(:,1)
               time_avg_psip_list%states(:,1) = qs%psip_list%states(:,1)
               time_avg_psip_list%nparticles = qs%psip_list%pops(:,1)/qs%psip_list%pop_real_factor
               time_avg_psip_list%nstates = 1
  
  
          allocate(time_avg_psip_list_sq(2,size(qs%psip_list%states(1,:))))
          time_avg_psip_list_sq(1,1) = qs%psip_list%states(1,1)
          time_avg_psip_list_sq(2,1) = (real(qs%psip_list%pops(1,1))/qs%psip_list%pop_real_factor)**2
          nstates_sq = 1
  
negative.code: 
                  nattempts_spawn=0
  
                  !call cumulative_population(qs%psip_list%pops, qs%psip_list%states(sys%basis%tot_string_len,:), &
                  !                           qs%psip_list%nstates, D0_proc, D0_pos, qs%psip_list%pop_real_factor, &
                  !                           .false., sys%read_in%comp, cumulative_abs_real_pops, &
                  !                           tot_abs_real_pop, D0_pop = D0_normalisation, p_ref=p_ref)
                  call update_bloom_threshold_prop(bloom_stats, nparticles_old(1))
  
                  ! Three options for evolution:
  