Lecture Notes of Dose distribution calculation inside phantom with ...
Lecture Notes of Dose distribution calculation inside phantom with ...
Lecture Notes of Dose distribution calculation inside phantom with ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
do k=kdosl(idose),kdosu(idose)do j=jdosl(idose),jdosu(idose),4jl=jju=min(j+3,jdosu(idose))write(6,400) zbound(k),zbound(k+1),k400 FORMAT(/T15,'for z=',F10.3,' to',F10.3,5X,'k=',I3)write(6,410) (ybound(jj),jj=jl,ju+1)410 FORMAT(' Ybounds:',F7.3,F12.3,3F17.3)write(6,420) (jj, jj=jl,ju)420 FORMAT(T10,'j=',T17,5(I4,13X))write(6,430) xbound(idosl(idose))430 FORMAT(' Xbounds (',F10.3,')')do i=idosl(idose),idosu(idose)write(6,440) xbound(i+1),i,(dose(i,jj,k),* min(99.9, 100.*doseun(i,jj,k)/dose(i,jj,k)),* jj=jl,ju)440 FORMAT(F8.3,I4,4(1PE11.3,'-',0PF4.1,'%') )if (mod(i,5).eq.0) thenwrite(6,fmt="(' ')")end ifend doif(mod(ipage,nperpg).eq.0.and.(ju.ne.jdosu(idose).* or.k.ne.kdosu(idose))) thenwrite(6,330) Titleipage=1elseipage=ipage+1end ifend do! end j-loopend do! end k-loopend if! end <strong>of</strong> x scan per page outputend do! end <strong>of</strong> idose loop! =============================call ecnsv1(1,nreg,totke)! =============================! ====================call counters_out(1)! ====================stopend!-------------------------last line <strong>of</strong> main code------------------------!----------------------------getvoxel.f---------------------------------! Version: 060907-1300! Read material and geomery data from egs5job.inp! 121119-1300 Fixed degug errors!-----------------------------------------------------------------------!23456789|123456789|123456789|123456789|123456789|123456789|123456789|12! ----------------------------------------------------------------------! Auxiliary subroutine for use <strong>with</strong> the EGS5 Code System! ----------------------------------------------------------------------! This is a data-entry subprogram for use <strong>with</strong> a general-purpose! egs5 user code to do cartesian coordinate dose deposition studies.! Every voxel (volume element) can have different materials and/or! varying densities (for use <strong>with</strong> CT data).! Basic pats <strong>of</strong> this subroutine related <strong>with</strong> geometry taken from! xyzdos.mor.! ----------------------------------------------------------------------!! voxels are labeled by indicies (i,j,k) and defined by:! xbound(i)
Change language