application comis quit
      subroutine newnt
      real en,u,v,w
      character*100 chform
      COMMON/QUEST/IQUEST(100)                         ! Declare IQUEST communication vector
      IQUEST(10) = 256000                               ! I require 65000 records
      call hropen(42,'NTUPLE','out.rz','NQE',4096,istat)
      open(unit=41,file='in.dat',status='old')
      call hbnt(1,'photons',' ')
      chform = 'en, u, v, w'
      call hbname(1,'photons',en,chform)
10    read(41,fmt=30,end=20) en,u,v,w
30    format(E15.6,E14.6,E14.6,E14.6)
      call hfnt(1)
      goto 10
20    close(41)
      call hrout(1,icycle,' ')
      call hrend('NTUPLE')
      close(42)
      end
      call newnt
      end
quit