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

emaxe = 0.D0 ! dummy value to extract min(UE,UP+RM).write(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! ==========! ------------------------------! Close files (after HATCH call)! ------------------------------close(UNIT=KMPI)close(UNIT=KMPO)!-----------------------------------------------------------------------! Output medium and region information to file for calculation mode.!-----------------------------------------------------------------------write(6,*) ' Quantities associated with each media:'do j=1,nmedwrite(6,130) (media(i,j),i=1,24)130 FORMAT(/,1X,24A1)write(6,140) rhom(j),rlcm(j)140 FORMAT(5X,' Rho=',G15.7,' g/cm**3 RLC=',G15.7,' cm')write(6,150) ae(j),ue(j),ap(j),up(j)150 FORMAT(5X,' AE=',G15.7,' MeV UE=',G15.7,' MeV'/ 5X,' AP=',G* 15.7,' MeV UP=',G15.7,' MeV')end dowrite(6,160)160 FORMAT(/' Information of medium and cut-off for central region')i=imax/2+1j=jmax/2+1do k=1,kmaxirl=1+i+(j-1)*imax+(k-1)*ijmaxif (med(irl).eq.0) thenwrite(6,170) k,irl170 FORMAT(' Medium(',I3,'-th z bin, region:',I5,')= Vacuum')elsewrite(6,180) k,irl,(media(ii,med(irl)),ii=1,24),* ecut(irl),pcut(irl),rhor(irl)180 FORMAT(' Medium(',I3,'-th z bin, region:',I5,* ')=',24A1,/5X,'ECUT=',G10.5,' MeV, PCUT=',* G10.5, ' MeV, density=',F10.3)end ifend do2.6 Step 6: Initialization-for-howfar getvoxel 2.7 Step 7: Initialization-for-ausgaausgab (1) 5