anchor.diff:  !logical                   :: conc_check
   integer, pointer          :: degen_list(:)
   integer, pointer          :: rotProdLab(:)
  -
  +integer, pointer          :: labPerms(:,:)
   
   SNF = 0
   allocate(HNFin(3,3,1))
  
positive.code:       write(*,'("Concentrations match for structures ",2(i6,1x))') strN, strN2
        write(13,'("Concentrations match for structures ",2(i6,1x))') strN, strN2
        if (any(L1/=L2)) then ! We must find the automorphism that connects the two different groups and adjust ilabeling2 accordingly
           print*,"hardwired change of labeling2 -- this is just temporary while testing/developing"
           write(*,'("labeling before: ",20(i1))') ilabeling2
           ilabeling2 = ilabeling2((/1,7,6,5,4,3,2,8,14,13,12,11,10,9/))
           write(*,'("labeling after:  ",20(i1))') ilabeling2
           if (associated(g1)) deallocate(g1)
           if (associated(g2)) deallocate(g2)
           !1) build the member list
  
           call make_member_list(diag,g1) ! This routine allocates g1
           do i = 1,3
              write(*,'("g1: ",20(i2,1x))') g1(i,:)
           enddo
           print*
  
           !2) find the group automorphism
           call matrix_inverse(real(L2(:,:,1),dp),L2inv,err)
           if (err) stop "L1inv is not defined"
           T = matmul(L1,L2inv(:,:,1))
           ! do i = 1,3
           !    write(*,'("L1:  ",3(i2,1x))') L1(i,:,1)
           ! enddo; print*
           ! do i = 1,3
           !    write(*,'("L2:  ",3(i2,1x))') L2(i,:,1)
           ! enddo
           do i = 1,3
              write(*,'(" T:  ",3(f7.3,1x))') T(i,:)
           enddo
  
           allocate(g2(3,size(g1,2)))
           g2 = -1
           g2 = nint(matmul(T,g1))
           !if (any(g2==-1)) stop "Bug in group member mapping"
           do i = 1,3
              write(*,'("g2: ",20(i2,1x))') g2(i,:)
           enddo
           print*
           do i = 1,size(qLabel,1)
              write(*,'("qL: ",20(i1))') qLabel(i,:)
           enddo
  
negative.code:    match = 0
     ! Read in each structure from the second file and see if it matches the current structures from
     ! file 1.
     f2 = 14; ! Define a unit number for the second file.
     open(f2,file=f2name,status="old")
     lc = 0 ! Count the number of lines
     do ! read the second file until the structure list begins
  