anchor.diff:  character(80) dummy, specAtoms 
   integer, pointer :: spin(:) ! Occupation variable of the positions
   real(dp), pointer :: x(:) ! Concentration of each component
  -logical :: file_exists
  +logical :: file_exists, co_ca_error
   
   open(13,file="readcheck_makestr.out")
   write(13,'(5(/),"This file echoes the input for the makestr.x program and")')
  
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: write(13,'("be useful for debugging and fixing bad input")')
  write(13,'("***************************************************")')
  
  if(iargc()/=2 .and. (iargc()/=3 .and. iargc()/=4)) stop "makestr.x requires two (optional three or four)  arguments: &
                 & the filename to read from and the number of the structure (to start and to stop at)"
  call getarg(1,dummy)
  read(dummy,'(a80)') fname
  call getarg(2,dummy)
  read(dummy,*) strNi
  if (iargc()>=3) then
    call getarg(3,dummy)
    read(dummy,*) strNf
  else
    strNf=strNi
  endif
  if (iargc()==4) then
    call getarg(4, dummy)
    read(dummy, *) foutput_unit
  else
    foutput_unit = 12
  endif
  
  ! Modifications made by GLWH 03/05/2014.
  