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

zbound(in+1)=zbound(in)+widthend doend ifif(nn.ne.nnn) thenwrite(ifto,280)end ifmaxz=maxz+nnnend doif (maxz.gt.LKMAX) thenwrite(ifto,'(A,I4,A,I4/A)') ' maxz(',maxz,* ') calculated in maxz < 0 mode is larger than',LKMAX,* ' You must increase LKMAX in auxcommon/aux_h.f 'stopend ifwrite(ifto,290) (zbound(i),i=1,maxz+1)end ifkmax=maxzijmax = imax*jmaxirmax = 1 + ijmax*kmaxnreg = irmaxwrite(ifto,320) imax,jmax,kmax,nreg320 FORMAT(' imax, jmxa, kmax, nreg =',4I8)! Check nregif (nreg .gt. MXREG) thenwrite(ifto,330) nreg330 FORMAT(' *** Stopped in getvoxel with nreg=',I5,' > MXREG')stopend ifwrite(ifto,340) nreg340 FORMAT(/,' number of region (nreg) =',I5,/,* ' nreg includes outside vacuum region (regin=1)')! Set all region except 1 set to medium=1.med(1)=0! ----------------------------------------------------------! Record 9: ipeangsw,iedgeflsw,iaugersw,iraysw,! ipolarsw,incohrsw,iprofrsw,impacrsw for medium 1! ----------------------------------------------------------read(4,*) ecutin,pcutin,ipeangsw,iedgeflsw,iaugersw,iraysw,* ipolarsw,incohrsw,iprofrsw,impacrswwrite(ifto,350) ecutin,pcutin,ipeangsw,iedgeflsw,iaugersw,* iraysw,ipolarsw,incohrsw,iprofrsw,impacrsw350 FORMAT(/' Medium 1'/' ecut =',G15.5,' and pcut =',G15.5/*' ipeangsw=',I3,',iedgeflsw=',I3,',iaugersw=',I3,'iraysw=',I3/*' ipolarsw=',I3,' ,incohrsw=',I3,',iprofrsw=',I3,',impacrsw=',I3)do i=2,irmaxmed(i)=1if (pcutin .gt. 0.)if (ecutin .gt. 0.)iphter(i) = ipeangswiedgfl(i) = iedgeflswiauger(i) = iaugerswiraylr(i) = irayswlpolar(i) = ipolarswincohr(i) = incohrswiprofr(i) = iprofrswimpacr(i) = impacrswend dopcut(i) = pcutinecut(i) = ecutin + RM! --------! Record 10 il,iu, jl,ju, kl,ku, medtmp, rhotmp, ecutin, pcutin! -------- (7I5,3F10.0) Line is repeated until a blank line found360 write(ifto,370)370 FORMAT(' Lower,upper i, j, k, medium, density')read(4,*) il,iu,jl,ju,kl,ku,medtmp,rhotmp,ecutin,pcutinif(il.eq.0 .and. iu.eq.0) go to 400! Check il etc.if(il.lt.0) il=163