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

end ifend ifelse if (u(np) .lt. 0.0) thendist = -(x(np) - xbound(irx))/u(np)if (dist .lt. ustep) thenustep = distif (irx .ne. 1) thenirnew=irl-1elseirnew = 1end ifend ifend if! -----------------! Check Y-direction! -----------------dnearl=min(dnearl,(ybound(iry+1)-y(np)),(y(np)-ybound(iry)))if (v(np) .gt. 0.0) thendist = (ybound(iry+1)-y(np))/v(np)if (dist .lt. ustep) thenustep=distif (iry .ne. jmax) thenirnew=irl+imaxelseirnew=1end ifend ifelse if (v(np) .lt. 0.0) thendist = -(y(np) - ybound(iry))/v(np)if (dist .lt. ustep) thenustep = distif (iry .ne. 1) thenirnew=irl-imaxelseirnew = 1end ifend ifend ifdnear(np)=dnearlreturnend! -----------------------! Return to ELECTR/PHOTON! -----------------------!--------------------------last line of howfar.f------------------------!-------------------------------encoea.f--------------------------------! Version: 030831-1300! Reference: SLAC-265 (p.19-20, Appendix 2)!-----------------------------------------------------------------------!23456789|123456789|123456789|123456789|123456789|123456789|123456789|12!! double precision encoea(energy)! Function to evaluate the energy absorption coefficient of air.! (Tables and Graphs oh photon mass attenuation coefficients and! energy-absorption coefficients for photon energies 1 keV to! 20 MeV for elements Z=1 to 92 and some dosimetric materials,! S. M. Seltzer and J. H. Hubbell 1995, Japanese Society of! Radiological Technology)!----------------------------------------------------------------------double precision function encoea(energy)real*8 hnu(38)/0.001,0.0015,0.002,0.003,0.0032029,0.0032029,* 0.004,0.005,0.006,0.008,0.01,0.015,0.02,0.03,0.04,* 0.05,0.06,0.08,0.10,0.15,0.2,0.3,0.4,0.5,0.6,0.8,1.0,* 1.25,1.5,2.0,3.0,4.0,5.0,6.0,8.0,10.0,15.0,20.0/real*8 enmu(38)/3599., 1188., 526.2, 161.4, 133.0, 146.0,* 76.36, 39.31, 22.70, 9.446, 4.742, 1.334, 0.5389,* 0.1537,0.06833,0.04098,0.03041,0.02407,0.02325,0.02496,* 0.02672,0.02872,0.02949,0.02966,0.02953,0.02882,0.02789,* 0.02666,0.02547,0.02345,0.02057,0.01870,0.01740,0.01647,* 0.01525,0.01450,0.01353,0.01311/real*8 energy,enm1,hnu1,ene0,slopeinteger i69