ucnaicgv (usercode̲NaI̲CG-volume) をシンプル体系にした版のファイル構成とその中身
ucnaicgv.f

ucnaicgv.inp!**********************************************************************!***************************** KEK, High Energy Accelerator Research!** u c n a i c g v _ base **Organization!*****************************
EGS5.0 USER CODE - 09 Oct 2015/1430!**********************************************************************!* This is a general User Code based on the cg geometry scheme.
!**********************************************************************! Original "ucnaicgv_vsp.f" by H. Hirayama and Y. Naito! the simplified version by H. Iwas!**********************************************************************implicit non
include 'include/egs5_h.f'
include 'include/egs5_bounds.f
include 'include/egs5_brempr.f
include 'include/egs5_edge.f
include 'include/egs5_media.f
include 'include/egs5_misc.f
include 'include/egs5_thresh.f
include 'include/egs5_uphiot.f
include 'include/egs5_useful.finclude 'include/egs5_usersc.f
include 'include/egs5_userxt.f
include 'include/randomm.f
include 'auxcommons/aux_h.f'
include 'auxcommons/edata.f
include 'auxcommons/etaly1.f
include 'auxcommons/instuf.f
include 'auxcommons/lines.f
include 'auxcommons/nfac.finclude 'auxcommons/watch.f
include 'auxcommons/geom_common.f' ! geom-common fil
common/totals/depe,deltae,spec(3,50),maxpic
real*8 depe,deltae,spe
real*8 totke,rnnow,etot,esumt,wti
integer maxpic
integer i,icases,idin,ie,ifti,ifto,ii,iiz,imed,ireg,isam* izn,nlist,j,k,n,ner,ntype,idu
character*24 medarr(MXMED!#####(6) User setting parameters-3 (EGS options
do i=1,nregif(med(i).ne.0) theiphter(i) = 1! Switches for PE-angle samplin
iedgfl(i) = 1! K & L-edge fluorescenc
iauger(i) = 0! K & L-Augeiraylr(i) = 0! Rayleigh scatterin
lpolar(i) = 0! Linearly-polarized photon scatterin
incohr(i) = 0! S/Z rejectio
iprofr(i) = 0! Doppler broadenin
impacr(i) = 0! Electron impact ionizatio
end iend d!

!

(1) Open fileopen(6, FILE='egs5job.out',STATUS='unknown')
open(4, FILE='egs5job.inp',STATUS='old')open(39,FILE='egs5job.pic',STATUS='unknown')

call source(wtin)! source generating her
irin = 0! set 0 for Autose
call cgcontrol(2,idum,idum) ! Autoset iri
etot = ekein + iabs(iqin)*RM
call shower (iqin,etot,xin,yin,zin,uin,vin,win,irin,wtin) ! main cal
ncount = ncount + 1
if(iwatch.gt.0) call swatch(-1,iwatch
end do!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MAIN LOOP EN!============================================================ MAIN PROGRAM EN!#####(3) User setting parameters-1
nmed=medarr(1)='AL(link PEGS5-material to EGS5 med!#####(10) Source descriptio
subroutine source(wtin
implicit non
include 'auxcommons/instuf.f
include 'include/randomm.f
real*8 wti!#####(4) User setting parametersncases= 1000! number of calculation
maxpict = 200! number of incident radiations in C
chard(1) = 7.62d0! set character dimension
chard(2) = 0.1d
chard(3) = 0.5dchard(4) = 5.0diqin = -1ekein = 1.253
xin= 0.0yin= 0.0zin= -5.0uin= 0.0vin= 0.0win= 1.0
wtin = 1.0! set medicDD********)

PGssst y,tt0 eg   sste'nDzyxent)

gg)))

0an'n

5nrt 

'' ''''' m e'''''

2'''''')))')   

c 2)

n)

endn)

ssn000Dt)) )) s0e1e0n5fD

1o(5) Some treatment
call geneout(1call media0(medarr)call geneout(2call pegs
npreci = 3
call cgcontrol(1,4,39dMRFKDDL(9) Some treatments after calculatio
if(iwatch.gt.0) call swatch(-88,iwatch
call plotxyz(99,0,0,0.D0,0.D0,0.D0,0.D0,0,0.D0,0.D0)
call cgcontrol(6,idum,idum
call counters_out(1write(6,*)'----- EGS5 sucessfully finished! ----en! General outpu! CG inpu! CG outpu(2) Initializatio
call counters_out(0)call block_set!

(7) Some treatment
call rluxinit
emaxe = 0.D0
! Autoset ecuts when emaxe=call hatch
call geneout(3call geneout(4call cgcontrol(4,idum,idum
call cgcontrol(5,idum,idum
if(iwatch.gt.0) call swatch(-99,iwatch!~~~~ (8) MAIN LOOP ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MAIN LOO
do i=1,ncasesg!

!

! incident particle charge! incident particle kinetic energ! source position of! source position of! source position of! cos(thetax), thetax=angle to x-axi! cos(thetay), thetay=angle to y-axi! cos(thetaz), thetaz=angle to z-axi! weigh

ELE&INP IRAYL=1 /EN
AL
AL
AENE
&INP AE=0.521,AP=0.010,UE=2.511,UP=2.0 /EN
PWL&INP /ENDEC&INP /END

ucnaicgv.data
RCC
RCC
EN
Z1
Z2Z3EN1120 000 0 -10+1+200-1-0 0 50 0 20 1