anchor.diff:         grs = 0.0d0
   c      get reciprocal lattice vectors 
          call smd_lat_invrt(lat,rlat)
  -       write(*,*) "Inverse Lattice vectors"
  -       write(*,*) (rlat(i,1),i=1,3)
  -       write(*,*) (rlat(i,2),i=1,3)
  -       write(*,*) (rlat(i,3),i=1,3)
  +c       write(*,*) "Inverse Lattice vectors"
  +c       write(*,*) (rlat(i,1),i=1,3)
  +c       write(*,*) (rlat(i,2),i=1,3)
  +c       write(*,*) (rlat(i,3),i=1,3)
   c      get lattice volume
          call smd_latt_vol(lat,vol)
   c
  
positive.code: 11     continue
         message = "looking for number of atoms "//pname
         read(fn_in,20,end=135) bigbuf
         write(*,*) bigbuf(1:len_trim(bigbuf)),is_integer(bigbuf)
  c      increment line number
         l=l+1
         if(.not.is_integer(bigbuf)) goto 11
  
negative.code:           call get_command_argument(i,buffer,l,istatus)
            if(istatus.ne.0) goto 18
            if(is_number(buffer)) then
  c            message = "Only atom tags are supported now"
  c            goto 911
              read(buffer,*) atom2_id
            else
              atom2_tag = buffer
  