anchor.diff:  
           if (nexcitors /= 0) then
               do i = 1, nexcitors
  -                pops(i) = get_rand_close_open(rng)*(tot_level_pop-1) + &
  +                pops(i) = get_rand_close_open(rng)*(tot_level_pop-depsilon) + depsilon + &
                                   ex_lvl_dist%cumulative_pop_ex_lvl(ex_level - 1)
               end do
               call insert_sort(pops)
  
positive.code:             cluster%pselect = cluster%pselect * select_proportion(choice)
  
              first = .true.
              accumulate_excitors : do i = 1, size(select_info, dim=2)
                  if (select_info(choice, i) /= 0) then
  
  
negative.code:         integer :: i, pos, prev_pos
          real(p) :: pops(1:nexcitors), tot_level_pop
          logical :: hit
          logical, intent(inout) :: allowed
          complex(p) :: excitor_pop
  
          ! Calculate total population of excips on states with this excit level.
  