anchor.diff:  
          i=INDEX(file_in,".",.true.) 
          aformat=file_in(i+1:)
  -       write(*,*) "format is ",aformat
   
          if(aformat.eq."trj") then
            call trj_read_header(fn_in,nvm,nva,nu,nprec)
  
positive.code:        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
  
negative.code: c     -------------------------------------
    100 continue
        read(fn,1000,end=11,err=911) card
   1000 format(a)
        if(card(1:5).ne.'frame') goto 100
  c
  c     skip one lines here
  c     --------------------
        read(fn,1000,err=911,end=911) card
  c     read lattice
        message = "reading lattice"
        do i=1,3
          read(fn,1000,err=911,end=911) card
          i0=1
          do k=1,3
            message = message(1:len_trim(message)+1)//token
  