
      
      subroutine wroutput (iRun,tgen,texec,nout1)
c----------------------------------------------------------------------
c     Print out the final results.
c----------------------------------------------------------------------
      include 'params.inc'
      include 'geomz.inc'
      include 'globals.inc'
      include 'shared.inc'

      integer*4 iRun

      INTEGER*4 nescgp(negrps)
      real*8 enesc(negrps)
      real*8 wcut, wmin, wmax
      real*8 wlost, wesc, wrr
      real*8 wabs, wcen, epgain, etot
      integer*4 npart, nphtot, nploss
      integer*4 nlost, nesc, nrr, nabs, ncen
      integer*4 nscat, nsplt, ntrac
      common /mjcstuff/  nescgp, enesc, wcut, wmin, wmax,
     &     wlost, wesc, wrr, wabs, wcen, epgain, etot, 
     &     npart, nphtot, nploss, nlost, nesc, nrr, nabs,
     &     ncen, nscat, nsplt, ntrac

      INTEGER nout1
      
c------------------------------print out final results


      if (print_flag .gt. 1) then
      write(nout1,209) iRun
 209  format('Run number: ',i4,//)
      endif

      npart = nphtot

      if (print_flag .gt. 1) then
      write (nout1,210)
 210  format(/20x,'plankian source:',
     .     / 2x,'# emitted',
     .     3x,'# killed',
     .     3x,'energy emitted (kev)'
     .     ,3x,'energy gained via rr')
      
      write (nout1,211) npart, nploss, etot, epgain
 211  format(1x,i8,4x,i8,5x,e12.4,12x,e12.4)
      endif
      
      trials = dble(npart)
      do 498 i=1,12
         wesc = wesc + enesc (i)
         nesc = nesc + nescgp(i)
 498  continue
      ffesc = wesc/trials
      g_ffesc(iRun) = ffesc
      std   = sqrt(ffesc*(1.0 - ffesc)/trials)

      if (print_flag .gt. 1) then      
      write (nout1,251) nabs, nesc, ntrac, nlost, ncen, nscat, nrr, nsplt
 251  format(/'nabs   = ',i10,
     |     /'nesc   = ',i10,
     |     /'ntrac  = ',i10,
     |     /'nlost  = ',i10,
     |     /'ncen   = ',i10,
     |     /'nscat  = ',i10,
     |     /'nrr    = ',i10,
     |     /'nsplit = ',i10)
      
      write (nout1,252) ffesc, std
 252  format(/' escape prob = ',f10.6,
     .     /' std =         ',f10.6)
      
      write (nout1,258)
 258  format(/1x,'energy group',5x,'# bundles escaping',5x,'energy (kev)
     |     ')
      write (nout1,259) (ig,nescgp(ig),enesc(ig)*bwgt, ig=1,12)
 259  format(4x,i4,15x,i6,13x,e10.4)
      
      write (nout1,271)
 271  format(/20x,'energy balance (kev):',/3x,'emitted',4x,'lost via rr'
     |     ,4x,'absorbed',5x,'escaped',5x,'censused',7x,'lost')
      endif
     
      err = epgain + wrr*bwgt


      if (print_flag .gt. 1) then   
      write (nout1,272) etot, err, wabs*bwgt, wesc*bwgt, wcen*bwgt,
     |     wlost*bwgt
 272  format(5(e11.4,2x),e11.4/)
      endif

      
      ttrack = 1.0d+6 * texec / dble(ntrac)

      if (print_flag .gt. 1) then  
      write(nout1,253) tgen, texec, ntrac, ttrack
 253  format(/5x,'generation time (sec)  = ',e15.7
     +     /5x,'execution time  (sec)  = ',e15.7
     +     //5x,'number of tracks       = ',i10,
     +     /5x,'time/track (microsec)  = ',e15.7)

      write(nout1,9902)
 9902 format('-------------------------------------------------')
      endif

      
      return
      end