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

310 FORMAT(/' Exposure in free air (using mu_en) =', G15.5,'+-',G15.* 5,' Gy/incident'/ ' Exposure at phantom surface (using mu_en) ='* , G15.5,'+-',G15.5,'Gy/incident'/ ' Backscattering factor =',G15* .5,'+-',G15.5)elsewrite(6,320) faexpa,faexrr,fexpsa,fexerr320 FORMAT(/' Exposure in free air (using mu_en) =', G15.5,'+-',G15.* 5,' Gy/incident'/ ' Exposure at phantom surface (using mu_en) ='* , G15.5,'+-',G15.5,'Gy/incident')end ifgetvoxel (Z-scan X-scan)2.10 subroutine getvoxel(ifto)getvoxel 4 4 1. Record 1 : (80 )2. Record 2 : (nmed)3. Record 3 : 24 pegs 4. Record 4 : characteristic dimension5. Record 5 : X-, Y-, Z- (maxx,maxy,maxz)6. Record 6 : X- maxx >0 maxx + 1 X-1 1 maxx x- LIMAX 22 LIMAX maxx auxcommon/au h.f LIMAX maxx < 0 X-abs(maxx) 1 1 maxx LIMAX LIMAX 22 LIMAX maxx auxcommon/au h.f LIMAX 7. Record 7 : Y-8. Record 8 : Z-9. Record 9 : ecut, pcut (0: o, 1:on)10