anchor.diff:           cycle
         endif
         do i = 0, k-1
  -         if (count(ilabeling2==i)/=count(ilabeling==i)) then
  +         if (count(ilabeling2==i)/=count(ilabeling1==i)) then
               write(13,'("Concentrations don''t match for color ",i1)') i
               write(*,'("Concentrations don''t match for color ",i1)') i
               cycle compareloop
  
positive.code: integer :: iuq, nuq, iP, nP, lc, iStr2, HNFtest(3,3), hdgenfact, labdgen, totdgen
  integer :: strN2, hnfN2, n2, pgOps2, diag2(3), iL, nL, idx, f1, f2, iLP, nLP, jp
  !integer :: sizeN2
  integer, allocatable :: ilabeling(:), ilabeling2(:), atomType(:), permutedLab(:)
  real(dp) :: eps
  logical full, HNFmatch, foundLab, err, equiatomic
  character(maxLabLength)   :: labeling ! List, 0..k-1, of the atomic type at each site
  
negative.code:       write(*,'("Structure number:",i8," is a match to # ",i9," in file 2.")') strN, match
     endif
     deallocate(ilabeling2)
     deallocate(ilabeling1)
  enddo
  close(13)
  
  