anchor.diff:  character(1) bulksurf
   character(40) dummy 
   
  +logical append
  +
   ! Testing
   integer, pointer :: spin(:) ! Occupation variable of the positions
   real(dp), pointer :: x(:) ! Concentration of each component
  
positive.code:          endif
        enddo
     enddo
     close(12);
     
     call map_enumStr_to_real_space(k,nAt,HNF,labeling,p,dvec,eps,sLVorig,aBas,spin,x,L,diag)
     
  
negative.code:            gIndx(iAt), iAt, labeling(gIndx(iAt):gIndx(iAt)), aBas(:,iAt)
     enddo
     
     write(strname,'("vasp.",i3.3)') strN
     write(strNstring,*) strN
     if (.not. append) open(12,file=strname)
     write(12,'(a80)') trim(adjustl(title)) // " str #: " // adjustl(strNstring)
     write(12,'("scale factor")')
     ! Make "nice" superlattice vectors (maximally orthogonal, Minkowski reduced)
  