Lecture Notes of Dose distribution calculation inside phantom with ...

Output the material data and parameters of each region to the result le (unit 6).! Define possible maximum total energy of electron before hatchemaxe = ekein + RMwrite(6,110)110 format(/' Call hatch to get cross-section data')! ------------------------------! Open files (before HATCH call)! ------------------------------open(UNIT=KMPI,FILE='pgs5job.pegs5dat',STATUS='old')open(UNIT=KMPO,FILE='egs5job.dummy',STATUS='unknown')write(6,120)120 FORMAT(/,' HATCH-call comes next',/)! ==========call hatch! ==========2.6 Step 6: Initialization-for-howfarDene various parameters used for the geometry denition in this step. In this user code, thispart is done at subroutine getvoxcel.2.7 Step 7: Initialization-for-ausgabInitialize various variables to be calculated and set the number of detectors and a number ofhistories to calculate (ncases) and to store trajectory data (maxpict).ncount = 0ilines = 0nwrite = 10nlines = 25idin = -1totke = 0.wtsum = 0.!---------------------! Clear variables!---------------------! Zero the dosedo k=1,kmaxdo j=1,jmaxdo i=1,imaxdepe(i,j,k)=0.D0depeh(i,j,k)=0.D0depeh2(i,j,k)=0.D0end doend doend dofaexp=0.D0faexps=0.D0faexp2s=0.D0fexps=0.D0fexpss=0.D0fexps2s=0.D0! =========================call ecnsv1(0,nreg,totke)call ntally(0,nreg)! =========================!--------------------------! History number!--------------------------! History numberncases=10000031