anchor.diff:             gIndx(iAt), iAt, labeling(gIndx(iAt):gIndx(iAt)), aBas(:,iAt)
      enddo
      
  -   write(strname,'("vasp.",i4.4)') strN
  +   write(strname,'("vasp.",i3.3)') strN
      write(strNstring,*) strN
  -   open(12,file=strname)
  +   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)
  
positive.code: enddo
  write(13,'("Skipped to the ",i10,"nth line in ",a80)') strNi,fname 
  
  ! 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
     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
  
negative.code: 
  enddo ! structure loop
  close(11)
  if (append) close(12)
  close(13)
  
  END PROGRAM makeStr
  