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.
call randomset(rnnow)do ie=1,nsebinif(rnnow.le.escdf(ie)) go to 1000end do1000 ekein=esbin(ie)call randomset(rnnow)do ie=1,nsebinif(rnnow.le.escdf(ie)) go to 1000end do1000 if (ie.gt.nsebin) thenie=nsebinend ifsaspec(ie)=saspec(ie)+1.D0if (escdf(ie).eq.escdf(ie-1)) thenekein=esbin(ie-1)elseekein=esbin(ie-1)+(rnnow-escdf(ie-1))*(esbin(ie)-esbin(ie-1))/* (escdf(ie)-escdf(ie-1))end if !-----------------------------! Sampled source spectrum!-----------------------------!-----------------------------! Sampled source spectrum!-----------------------------do ie=2,nsebinsaspec(ie)=saspec(ie)/float(ncases)end dowrite(6,272)272 FORMAT(/' Comparison between sampled spectrum and pdf'* /23X,' Sampled pdf ',25X,' Sampled pdf '* )do ie=2,nsebin,2if(ie.eq.nsebin) thenwrite(6,274) esbin(ie),saspec(ie),escdf(ie)-escdf(ie-1)274 FORMAT(1X,G9.3,' MeV(upper)-- ',2G12.5)elsewrite(6,276) esbin(ie),saspec(ie),escdf(ie)-escdf(ie-1),* esbin(ie+1), saspec(ie+1),escdf(ie+1)-escdf(ie)276 FORMAT(1X,G9.3,' MeV(upper)-- ',2G12.5,3X, ' ',G9.3,* ' MeV(upper)-- ',2G12.5)end ifend do 280 FORMAT(/' Absorbed energy <strong>inside</strong> <strong>phantom</strong> for Co-60 photon'/280 FORMAT(/' Absorbed energy <strong>inside</strong> <strong>phantom</strong> for 100kV X-ray'/5. uc<strong>phantom</strong>cgv2.inp 21
Change language