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

if (iwatch .gt. 0) call swatch(iarg,iwatch)! ========================! ---------------------------------! Keep track of how deep stack gets! ---------------------------------if (np.gt.MXSTACK) thenwrite(6,100) np,MXSTACK100 FORMAT(//' In AUSGAB, np=',I3,' >= maximum stack',* ' allowed which is',I3/1X,79('*')//)stopend if! ------------------------! Set some local variables! ------------------------irl = ir(np)iql = iq(np)edepwt = edep*wt(np)! ----------------------------------------------------------------! Print out stack information (for limited number cases and lines)! ----------------------------------------------------------------if (ncount .le. nwrite .and. ilines .le. nlines) thenilines = ilines + 1write(6,101) e(np),x(np),y(np),z(np),u(np),v(np),w(np),* iql,irl,iarg101 FORMAT(7G15.7,3I5)end if! -----------------------------------------------------------! Keep track of energy deposition (for conservation purposes)! -----------------------------------------------------------if (iarg .gt. 5) returnesum(iql+2,irl,iarg+1) = esum(iql+2,irl,iarg+1) + edepwtnsum(iql+2,irl,iarg+1) = nsum(iql+2,irl,iarg+1) + 1i=mod(irl-1,imax)if (i.eq.0) i=imaxk=1+(irl-1-i)/ijmaxj=1+(irl-1-i-(k-1)*ijmax)/imaxif (irl.gt.1.and.edep.ne.0.D0) thendepe(i,j,k)=depe(i,j,k)+edepwtend if!--------------------------------! Check cross phantom surface!--------------------------------if(i.eq.imax/2+1.and.j.eq.jmax/2+1) then ! X-Y central regionif (abs(irl-irold).eq.ijmax.and.iq(np).eq.0) thenif ((w(np).gt.0.0.and.k.eq.2).or.* (w(np).le.0.0.and.k.eq.1)) thenif (dabs(w(np)).ge.0.0349) thencmod=dabs(w(np))elsecmod=0.01745end ifesing=e(np)dcon=encoea(esing)! PHOTX datafexps=fexps+e(np)*dcon*wt(np)/cmodif (w(np).lt.0.0) latch(np)=1if (w(np).gt.0.0.and.latch(np).eq.0) thenfaexp=faexp+e(np)*dcon*wt(np)/cmodend ifend ifend ifend if12