anchor.diff:     foundHNF = .false.
      do iHNF = 1, nHNF
         if(.not. all(HNF==HNFin(:,:,iHNF))) cycle
  +      foundHNF = .true.
      enddo
      if(.not. foundHNF) then
         write(13,'("HNF doesn''t match for str #:",i8," ---Skipping")') StrN
  
positive.code: nAt = size(aBas,2)
  
  write(18,'("Diagonal entries of HNF: ",3(i3,1x))') diag
  g = 0
  do iAt = 1, nAt ! Map each real space vector (atom position) into the group
     g(:,iAt)=matmul(matmul(L,Ainv),aBas(:,iAt))
  
negative.code: do ip = 1, nP
     write(17,'("Perm #",i3,":",1x,200(i2,1x))') ip,LattRotList(1)%perm(ip,:)
  enddo
  call find_labeling_from_atom_basis(L(:,:,1),pLV,aBas,aTyp,SNF,labeling,eps)
  write(17,'("Input atom labels: ",200(i2,1x))') aTyp
  write(17,'("Labeling:          ",200(i2,1x))') labeling
  
  