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 noneinclude ’include/egs5_h.f’! Main EGS "header" fileinclude ’include/egs5_epcont.f’ ! epcont contains edepinclude ’include/egs5_stack.f’ ! stack contains x, y, z, u, v,! w, ir and npcommon/score/escore(3)real*8 escoreinteger iarginteger irl! Arguments! Local variablesif (iarg.le.4) thenirl=ir(np)! Pick up current region numberescore(irl)=escore(irl)+edepend ifreturnend!--------------------------last line of ausgab.f------------------------AUSGAB is still very simple since all we need to do is to keep track of the energy depositedin the three regions. The variable EDEP (available through COMMON/EPCONT/) contains the energydeposited during a particular step for a variety of different IARG-situations, as described in thecomments 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 nextinseed= 1 (seed for generating unique sequences of Ranlux)ranlux luxury level set by rluxgo : 1 p= 48ranlux initialized by rluxgo from seed 1Start tutor2Call hatch to get cross-section data151