anchor.diff:             message = "no free file units"
              goto 911
          else
  -           message = "opening "//trim(solutefile)
  -           open(fns,file=trim(solutefile),
  +           message = "opening "//trim(infile)
  +           open(fns,file=trim(infile),
        $                form='formatted',status='old',err=911)
          end if 
          if(.not.util_get_io_unit(fno)) then
  
positive.code:        rewind(fns)
         call xyz_read(nat,c,atag,fns)
         write(*,*) "nat",nat
         do i=1,nid-2
           v  = c(:,id(nid-1)) 
           v1 = c(:,id(nid))-v
           x = c(:,id(i))
           write(*,*) "rotating by",angle
           angle = angle*pi/180.0
  
negative.code:          y = y1
           z = z1
         end subroutine
  
        subroutine my_get_command_argument(i,buffer,l,istatus)
        implicit none
        integer i
        integer l
        character*(*) buffer
        integer istatus
  c
        istatus = 0
        l = 0
        call getarg(i,buffer)
        if(buffer.eq." ") istatus = 1
  c      call get_command_argument(i,buffer,l,istatus)
        end subroutine
  
        function is_number(string)
        implicit none
        character*(*) string
        logical is_number
  c
        character(len=11) :: anumber = "0123456789."
        integer i,l
        is_number = .false.
        if(string.eq." ") return
        l = len_trim(string)
        do i=1,l
         if(string(i:i).ne." ") exit
        end do
        is_number = verify(string(i:l), anumber).eq.0
        end function
  
  
  