real coord(3,2000)
      real prec,time,potE,uconst,t_bath,qfrag(100)
      integer ihpb(100),jhpb(100)
      character*80 arg,ctime_start /"0"/,ctime_end /"1000000000"/,
     &  cfreq /"1"/
      integer*4 time_start,time_end,freq
c      real time_start,time_end,freq
      
      call getarg(1,arg)
      if (iargc().eq.0) stop
      if (iargc().gt.1) then
        call getarg(2,ctime_start)
        read (ctime_start,*) time_start
      else
        time_start=0
      endif
      if (iargc().gt.2) then
        call getarg(3,ctime_end)
        read (ctime_end,*) time_end
      else
        time_end=1000000000
      endif
      if (iargc().gt.3) then
        call getarg(4,cfreq)
        read (cfreq,*) freq
      else
        freq=1
      endif

      call xdrfopen(ixdrf,arg, "r", iret)
       
      licz=0
      do while(.true.) 
       call xdrffloat(ixdrf, time, iret)
c       print *,"time",time," iret",iret
       if(iret.eq.0) exit
       kk = kk + 1
       call xdrffloat(ixdrf, potE, iret)
       call xdrffloat(ixdrf, uconst, iret)
c       print *,"potE",potE," uconst",uconst
#ifdef NEWUNRES
       call xdrffloat(ixdrf, uconst_back, iret)
#endif
c       print *,"uconst_back",uconst_back
       call xdrffloat(ixdrf, t_bath, iret)
c       print *,"t_bath",t_bath
       call xdrfint(ixdrf, nss, iret) 
       do j=1,nss
        call xdrfint(ixdrf, ihpb(j), iret)
        call xdrfint(ixdrf, jhpb(j), iret)
       enddo
c       print *,"nss",nss
       call xdrfint(ixdrf, nfrag, iret)
       do i=1,nfrag
        call xdrffloat(ixdrf, qfrag(i), iret)
       enddo
       prec=10000.0

       isize=0
       call xdrf3dfcoord(ixdrf, coord, isize, prec, iret)

       licz=licz+1
c       print *,licz,time_start,time_end,freq,mod(licz,freq)
       if (licz.gt.time_end) exit
       if (licz.ge.time_start .and. mod(licz,freq).eq.0) then
       write (*,'(e15.8,2e15.5,f12.5,$)') time,potE,uconst,t_bath
       write (*,'(i4,$)') nss,(ihpb(j),jhpb(j),j=1,nss)
       write (*,'(i4,20f7.4)') nfrag,(qfrag(i),i=1,nfrag)
       write (*,'(8f10.5)') ((coord(k,j),k=1,3),j=1,isize)
       endif
      enddo
     
      end