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:    ! Use the permutations effected by the rotations that fix the superlattice to generate labelings
     ! that are equivalent. The list of equivalent labelings will be used when we look for a match in the
     ! struct_enum file
     allocate(ilabeling(n*nD1),ilabeling2(n*nD1))!,permutedLab(n*nD1)) 26 July, permutedLab not used
     read(labeling,'(500i1)') ilabeling
     !write(*,'("ilabeling:",/,30(i1),/)') ilabeling
     !print *,"size",size(ilabeling)
     !print *,"strN",strN
     call find_equivalent_labelings(ilabeling,LattRotList1,pLabel,rotProdLab)
     ! do i = 1,size(plabel,1)
     !    write(*,'("equiv_labelings <> ",100(i1))') pLabel(i,:)
     ! enddo
  
negative.code:       write(*,'("Structure number:",i8," is a match to # ",i9," in file 2.")') strN, match
     endif
     deallocate(ilabeling2)
     deallocate(ilabeling1)
  enddo
  close(13)
  
  