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

if(idosu(idgrp).le.0 .or. idosu(idgrp).ge.imax) idosu(idgrp)=imaxif(jdosl(idgrp).le.0) jdosl(idgrp)=1if(jdosu(idgrp).le.0 .or. jdosu(idgrp).ge.jmax) jdosu(idgrp)=jmaxif(kdosl(idgrp).le.0) kdosl(idgrp)=1if(kdosu(idgrp).le.0 .or. kdosu(idgrp).ge.kmax) kdosu(idgrp)=kmaxwrite(ifto,450) idosl(idgrp),idosu(idgrp),jdosl(idgrp),* jdosu(idgrp),kdosl(idgrp),kdosu(idgrp),izscan(idgrp)450 FORMAT('+',T5,3(I6,I4),I6)go to 420460 continueidgrp=idgrp-1if(idgrp.gt.LMXDOS) thenwrite(ifto,470) idgrp,LMXDOS470 FORMAT(' idgrp(=',I5,') must be less than LMXDOS(=',I5,')'/* ' Or you must chnage LMXDOS in xyzdose_h.f')end if! -----------------! Record 12: iwatch! -----------------read(4,*) iwatchwrite(ifto,790) iwatch790 FORMAT(//,' SWATCH tracking switch: iwatch=',I2,* ' (0=off, 1=each interaction, 2=each step)')! --------------------------------------! Record 13: ibrdst,iprdst,ibrspl,nbrspl! --------------------------------------read(4,*) ibrdst,iprdst,ibrspl,nbrsplwrite(ifto,800) ibrdst,iprdst,ibrspl,nbrspl800 FORMAT(/,' IBRDST=',I2,/,' IPRDST=',I2,/,' IBRSPL=',I2,' (NBRSPL='*,I5,')')if (ibrspl .gt. 0) thenif (nbrspl .gt. 0) thenfbrspl = 1.0/float(nbrspl)elsewrite(ifto,810) ibrspl,nbrspl810 FORMAT(//,' Stopped in Getvoxel with IBRSPL=',* I5,' and NBRSPL=',I5)stopend ifend ifreturnend! --------------! Return to MAIN! --------------!-------------------------last line of getvoxel.f------------------------!-------------------------------ausgab.f--------------------------------! Version: 031108-1300! 060801-1000! Reference: SLAC-265 (p.19-20, Appendix 2)!-----------------------------------------------------------------------!23456789|123456789|123456789|123456789|123456789|123456789|123456789|12! ----------------------------------------------------------------------! Required subroutine for use with the EGS5 Code System! ----------------------------------------------------------------------! A simple AUSGAB to:!! 1) Score energy deposition! 2) Print out stack information! 3) Print out particle transport information (if switch is turned on)!! ----------------------------------------------------------------------subroutine ausgab(iarg)implicit none65