anchor.diff:      integer C(size(conc)), X(size(conc))
       integer iK, n
       integer(li) p
  -    
  +
       n = size(conc)
   
       call get_Cs(conc,C)
  
positive.code:                 write(*,'("Multiplier: ",20(i3,1x))') multiplier
                  stop
               endif
               
               if (lab(idx)=='') then
                  degeneracy_list(nUniq) = degeneracy_list(nUniq) + 1
                  lab(idx) = 'D'  ! Mark as a duplicate
  
negative.code:        enddo
         if (j < 1) exit ! We're done counting (hit all possible numbers), exit
         digCnt(j) = digCnt(j) + 1
         a(j) = label(digCnt(j),j)
         c(a(j)) = c(a(j)) + 1     ! Add 1 to the number of digits of the j+1-th kind
         c(label(digCnt(j)-1,j)) = c(label(digCnt(j)-1,j)) - 1 ! subtract
         ! 1 from the number of digits of the j-th kind
  