application comis quit
      subroutine newnt
      integer nh,part
      real x,y,z,u,v,w,p,en,edep
      character*100 chform
      call hropen(42,'NTUPLE','out.rz','N',1024,istat)
      open(unit=41,file='pions.dat',status='old')
      call hbnt(1,'pions',' ')
      chform = 'part:I, x, y, z, u, v, w, p, en, edep'
      call hbname(1,'pions',part,chform)
10    read(41,fmt=30,end=20) nh,part,x,y,z,u,v,w,p,en,edep
30    format(2I3,9E13.6)
      call hfnt(1)
      goto 10
20    close(41)
      call hrout(1,icycle,' ')
      call hrend('NTUPLE')
      close(42)
      end
      call newnt
      end
quit