THE EGS5 CODE SYSTEM

! locally (in fact EDEP = particles kinetic energy).
! For IARG=3, we are discarding the particle since it is in
! region 1 or 3, so score its energy.
!
!***********************************************************************
subroutine ausgab(iarg)
implicit none
include ’include/egs5_h.f’
! Main EGS "header" file
include ’include/egs5_epcont.f’ ! epcont contains edep
include ’include/egs5_stack.f’ ! stack contains x, y, z, u, v,
! w, ir and np
common/score/escore(3)
real*8 escore
integer iarg
integer irl
! Arguments
! Local variables
if (iarg.le.4) then
irl=ir(np)
! Pick up current region number
escore(irl)=escore(irl)+edep
end if
return
end
!--------------------------last line of ausgab.f------------------------
AUSGAB is still very simple since all we need to do is to keep track of the energy deposited
in the three regions. The variable EDEP (available through COMMON/EPCONT/) contains the energy
deposited during a particular step for a variety of different IARG-situations, as described in the
comments above and further elaborated upon in Appendix B. In this example, but not always,
we can sum EDEP for any value of IARG up to 4. The following is the output provided by tutor2
(named tutor2.out in distribution file).
PEGS5-call comes next
inseed= 1 (seed for generating unique sequences of Ranlux)
ranlux luxury level set by rluxgo : 1 p= 48
ranlux initialized by rluxgo from seed 1
Start tutor2
Call hatch to get cross-section data
151