THE EGS5 CODE SYSTEM

! usersc contains emaxe
common/geom/zbound
real*8 zbound
! geom passes info to our howfar routine
real*8 ein,xin,yin,zin,
* uin,vin,win,wtin
integer iqin,irin
integer i,j
character*24 medarr(1)
! Arguments
! Local variables
! ----------
! Open files
! ----------
open(UNIT= 6,FILE=’egs5job.out’,STATUS=’unknown’)
! ====================
call counters_out(0)
! ====================
!-----------------------------------------------------------------------
! Step 2: pegs5-call
!-----------------------------------------------------------------------
! ==============
call block_set
! Initialize some general variables
! ==============
! ---------------------------------
! define media before calling PEGS5
! ---------------------------------
nmed=1
medarr(1)=’TA ’
do j=1,nmed
do i=1,24
media(i,j)=medarr(j)(i:i)
end do
end do
! nmed and dunit default to 1, i.e. one medium and we work in cm
chard(1) = 0.1d0 ! optional, but recommended to invoke
! automatic step-size control
! ---------------------------------------------
! Run KEK version of PEGS5 before calling HATCH
! (method was developed by Y. Namito - 010306)
! ---------------------------------------------
142