SUBROUTINE ughbook
      implicit none

      integer istat, i
c------------------------------------------------------------------------------
      include 'ntpl.inc'
c------------------------------------------------------------------------------
      include 'edep1.inc'
c
c      open(53,file='edep.bin',form='UNFORMATTED')

c      call uzero(elowz, 1,17920)
c      call uzero(elum , 1,160000)
c      call uzero(eqd0 , 1,232320)
c      call uzero(eqfd1, 1,181632)

      etotal0 = 0.
      do i = 1, 300
         edepbuf(i) = 0.
      enddo
c
      call hropen(1,'nlchit','nlchit.rz','N',1024,istat)
      call hbnt( 1,'nlchit', ' ')
      call hbname( 1,'nlchit', nh,
     -  'nh[0,1000],irl(nh),ql(nh),el(nh),xl(nh),yl(nh),zl(nh),'//
     -  'ul(nh),vl(nh),wl(nh),ph(nh),fh(nh),evnt(nh)')
c
c--- neutron production
c
c      call hbook1(11,' Egam (MeV)', 50,-2.,3.,0.)
c      call hbook1(12,' Egam (MeV)', 50,0.,1000.,0.)

c      call hbook1(12,' step (cm)', 40,0., 10.,0.)
c      call hbook1(13,' Eneutron (MeV)', 40,0., 20.,0.)
c      call hbook1(14,' Egam (MeV) with neutrons', 40,0.,50.,0.)

      return
      end