anchor.diff:        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(foutput_unit,'(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
  
positive.code: !!GH      close(14)
  !!GH   end if
  
     close(12);
     deallocate(spin,x,gIndx,aBas)
  enddo ! structure loop
  if (file_exists) close(14)
  
negative.code:    call map_enumStr_to_real_space(k,n,HNF,labeling,p,dvec,eps,sLV,aBas,spin,gIndx,x,L,diag,minkowskiReduce=.false.)
     write(strname,'("vasp.",i6.6)') strN
     write(strNstring,*) strN
     if (foutput_unit /= 6) then
        open(foutput_unit, file=strname)
     endif
     write(foutput_unit,'(a80)') trim(adjustl(title)) // " str #: " // adjustl(strNstring)
     !GH write(foutput_unit,'("scale factor")')
     write(foutput_unit,'("1.00")')
     do i = 1,3
        write(foutput_unit,'(3f12.8)') sLV(:,i)
     enddo
     call matrix_inverse(sLV,sLVinv)
     write(13,'("New inverse after reduction",/,3(3(f7.3,1x),/))') (sLVinv(i,:),i=1,3) 
  