anchor.diff:         character*30 pname 
   c
          pname = "xyz_read_coords_byname"
  +       c = 0.0d0
   c      ------------------------------------------
   c      get number of atoms (skipping empty lines)
   c      ------------------------------------------
  
positive.code:             goto 911
            end if
            go to 16
          else if(buffer.eq."-nbins") then
            i = i+1
            call get_command_argument(i,buffer,l,istatus)
  
negative.code:        character*30 pname 
  c
         pname = "xyz_read_coords_byname"
         c1 = 0.0d0
         c2 = 0.0d0
         write(*,*) "in "//pname
         j1 = 0
         j2 = 0
  