anchor.diff:         end if
   c      the lattice
          if(file_lattice.ne." ") then
  +         if(overb) write(*,*) "reading lattive from
  +     +       file"//trim(file_lattice)
            call lattice_read_file(file_lattice,lat)
          else
            do k=1,3
  
positive.code:          write(*,*) "Output file",infile(nfil)
         end if
  c       
         rmax = 4.0
         nb = 100
         atom1_tag = " "
         atom1_id = 3910
         atom2_tag = "O"
         latv = 10.0014453
         file_in = "test.xyz"
         file_in = "m.trj"
         file_out = "test.out"
         file_lattice = "lat.dat"
  c       nfrm = 2
  c      ---------------------------      
  c      start checks/balances
  c      ---------------------------      
  c      files
         if(nfil.lt.2) then
           message = "please provide input/output files"
           goto 911
         end if
  c      atom tags
  
negative.code:        fn_out = 11
         nfrm = 0
         aformat = " "
         overb = .false.
         file_lattice = " "
         latv = -1.0
         lat = 0.0
  