28.11.2014 Views

THE EGS5 CODE SYSTEM

THE EGS5 CODE SYSTEM

THE EGS5 CODE SYSTEM

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

iraylr(2)=1<br />

! Turn on rayleigh scattering in the slab<br />

iedgfl(2)=1<br />

! 1: Turn on fluorescence production in the slab<br />

! 0: Turn off fluorescence production in the slab<br />

! Note, above three parameters need to be set for all regions in which<br />

! there is particle transport - just region 2 in this case<br />

! --------------------------------------------------------<br />

! Random number seeds. Must be defined before call hatch<br />

! or defaults will be used. inseed (1- 2^31)<br />

! --------------------------------------------------------<br />

luxlev=1<br />

inseed=1<br />

write(6,120) inseed<br />

120 FORMAT(/,’ inseed=’,I12,5X,<br />

* ’ (seed for generating unique sequences of Ranlux)’)<br />

! =============<br />

call rluxinit<br />

! =============<br />

! Initialize the Ranlux random-number generator<br />

!-----------------------------------------------------------------------<br />

! Step 4: Determination-of-incident-particle-parameters<br />

!-----------------------------------------------------------------------<br />

! Define initial variables for 100 keV beam of photons normally incident<br />

! on the slab<br />

iqin=0<br />

! Incident photons<br />

! 100 keV<br />

ein=0.100<br />

xin=0.0<br />

yin=0.0<br />

zin=0.0<br />

! Incident at origin<br />

uin=0.0<br />

vin=0.0<br />

win=1.0<br />

! Moving along z axis<br />

irin=2<br />

! Starts in region 2, could be 1<br />

wtin=1.0<br />

! weight = 1 since no variance reduction used<br />

!-----------------------------------------------------------------------<br />

! Step 5: hatch-call<br />

!-----------------------------------------------------------------------<br />

! Maximum total energy of an electron for this problem must be<br />

! defined before hatch call<br />

emaxe = ein + RM<br />

180

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!