THE EGS5 CODE SYSTEM

! ====================
call counters_out(0)
! ====================
!-----------------------------------------------------------------------
! Step 2: pegs5-call
!-----------------------------------------------------------------------
! ==============
call block_set
! Initialize some general variables
! ==============
! nmed and dunit default to 1, i.e. one medium and we work in cm
if(loop.eq.3) then
chard(1) = 0.20d0 ! optional, but recommended to invoke
chard(2) = 0.20d0 ! automatic step-size control
else
chard(1) = 0.00d0 ! optional, but recommended to invoke
chard(2) = 0.00d0 ! automatic step-size control
endif
write(6,100) loop, chard(1)
100 FORMAT(72(’*’),/,
*’Initializing EGS5, loop = ’,I1,’: charD = ’,f5.2,/,
*72(’*’),/)
if(loop.eq.1) then
! ---------------------------------
! define media before calling PEGS5
! ---------------------------------
nmed=2
medarr(1)=’SI with long steps ’
medarr(2)=’SI with short steps ’
do j=1,nmed
do i=1,24
media(i,j)=medarr(j)(i:i)
end do
end do
! ---------------------------------------------
! Run KEK version of PEGS5 before calling HATCH
! (method was developed by Y. Namito - 010306)
! ---------------------------------------------
write(6,110)
110 FORMAT(’ PEGS5-call comes next’/)
! ==========
call pegs5
! ==========
158