anchor.diff:  use numerical_utilities
   use enumeration_utilities ! This maps structures from real space into the group
   use enumeration_types, only: maxLabLength
  +use io_utils, only: co_ca
   implicit none
   character(80) fname, title, strname, strNstring
   character(maxLabLength) :: labeling
  -integer ioerr, iline, ic, i, ilab, pgOps, nD, hnfN, iAt, jAt, iSpec, nSpec, idx
  +integer ioerr, iline, ic, i, ilab, pgOps, nD, hnfN, iAt, jAt, iSpec, nSpec, idx, foutput_unit
   integer k, strN, istrN, strNi, strNf, sizeN, n, diag(3), a,b,c,d,e,f, nType
   integer HNF(3,3), L(3,3), hnf_degen, lab_degen, tot_degen
   integer, pointer :: gIndx(:)
  
positive.code:       do iAt = 1, n*(nD+nSpec)
           if (labeling(gIndx(iAt):gIndx(iAt))==char(ilab+48)) then ! We have a match to the current
              ! label so print it
             write(12,'(3(f12.8,1x))') aBas(:,iAt)
             write(13,'("At. #: ",i2," position:",3(f7.3,1x),"<",a1,">")') iAt, aBas(:,iAt), labeling(gIndx(iAt):gIndx(iAt))
           endif
        enddo
  
negative.code: !!GH      close(14)
  !!GH   end if
  
     close(foutput_unit);
     deallocate(spin,x,gIndx,aBas)
  enddo ! structure loop
  if (file_exists) close(14)
  