anchor.diff:  enddo
   write(13,'("Skipped to the ",i10,"nth line in ",a80)') strNi,fname 
   
  +if (append) open(12,file="vasp.POSCAR")
   ! Read in the info for the given structure
   do istrN=strNi,strNf
      read(11,*) strN, hnfN, sizeN, nAt, pgOps, diag, a,b,c,d,e,f, L, labeling
  +   if (append) print *, "current structure: ", strN
      write(13,'("strN, hnfN, sizeN, nAt, pgOps, diag, a,b,c,d,e,f, L, labeling")')  
      write(13,'(i11,1x,i7,1x,i8,1x,i2,1x,i2,1x,3(i3,1x),1x,6(i3,1x),1x,9(i3,1x),1x,a40)') &
           strN, hnfN, sizeN, nAt, pgOps, diag, a,b,c,d,e,f, L, labeling
  
positive.code: character(1) bulksurf
  character(40) dummy 
  
  ! Testing
  integer, pointer :: spin(:) ! Occupation variable of the positions
  real(dp), pointer :: x(:) ! Concentration of each component
  
negative.code:          endif
        enddo
     enddo
     if (.not. append) close(12);
     
     call map_enumStr_to_real_space(k,nAt,HNF,labeling,p,dvec,eps,sLVorig,aBas,spin,x,L,diag)
     
  