anchor.diff:              cluster%pselect = cluster%pselect * select_proportion(choice)
   
               first = .true.
  +            allowed = .true.
               accumulate_excitors : do i = 1, size(select_info, dim=2)
                   if (select_info(choice, i) /= 0) then
   
  
positive.code:         integer :: i, pos, prev_pos
          real(p) :: pops(1:nexcitors), tot_level_pop
          logical :: hit
          logical, intent(out) :: allowed
          complex(p) :: excitor_pop
  
          ! Calculate total population of excips on states with this excit level.
  
negative.code:         end if
  
          if (abs(excitor_pop) <= initiator_pop) cdet%initiator_flag = 3
          call remove_ex_level_bit_string(sys%basis%string_len, cdet%f)
          cluster%excitation_level = get_excitation_level(f0, cdet%f)
          cluster%amplitude = excitor_pop
  
          ! Sign change due to difference between determinant
          ! representation and excitors and excitation level.
  