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

if(iu.lt.0 .or. iu.ge.imax) iu=imaxif(jl.le.0) jl=1if(ju.le.0 .or. ju.ge.jmax) ju=jmaxif(kl.le.0) kl=1if(ku.le.0 .or. ku.ge.kmax) ku=kmax! Check medtmpif(medtmp.lt.0 .or. medtmp.gt.nmed) medtmp=1write(ifto,380) il,iu,jl,ju,kl,ku,medtmp,rhotmp380 FORMAT('+',3('(',I3,I4,')'), I4, F10.3)if (medtmp.ne.0) then! ---------------------------------------------------! Record 10a: ipeangsw,iedgeflsw,iaugersw,iraysw,! ipolarsw,incohrsw,iprofrsw,impacrsw! ---------------------------------------------------read(4,*) ipeangsw,iedgeflsw,iaugersw,iraysw,ipolarsw,* incohrsw,iprofrsw,impacrswwrite(ifto,390) ecutin,pcutin,ipeangsw,iedgeflsw,iaugersw,* iraysw,ipolarsw,incohrsw,iprofrsw,impacrsw390 FORMAT(' 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=il,iudo j=jl,judo k=kl,kuirl=1+i+(j-1)*imax+(k-1)*ijmaxmed(irl)=medtmpif (rhotmp.ne.0) rhor(irl)=rhotmpif (pcutin .gt. 0.) pcut(irl) = pcutinif (ecutin .gt. 0.) ecut(irl) = ecutin+ RMiphter(irl) = ipeangswiedgfl(irl) = iedgeflswiauger(irl) = iaugerswiraylr(irl) = irayswlpolar(irl) = ipolarswincohr(irl) = incohrswiprofr(irl) = iprofrswimpacr(irl) = impacrswend doend doend doelsedo i=il,iudo j=jl,judo k=kl,kuirl=1+i+(j-1)*imax+(k-1)*ijmaxmed(irl)=0end doend doend doend ifgo to 360400 continue! --------! Record 11 il,iu, jl,ju, kl,ku,izscan! ---------write(ifto,410)410 FORMAT(' 3 pairs defining lower,upper x,y,z indecses of dose',* 'regions'/' for which results are to be output'/* ' izscan non-zero : scan per page'/* ' One set of 6 per line, end with all zero')idgrp=0420 idgrp=idgrp+1write(ifto,430)430 FORMAT('$: ')read(4,*) idosl(idgrp),idosu(idgrp),jdosl(idgrp),jdosu(idgrp),* kdosl(idgrp),kdosu(idgrp),izscan(idgrp)if(idosl(idgrp).eq.0 .and. idosu(idgrp).eq.0) go to 460if(idosl(idgrp).le.0) idosl(idgrp)=1! End of define.64