cg Version
cg Version
cg Version
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
{ Add sampling routines for source photon energy sampling.Just after! ----------------------! Select incident energy! ----------------------add followings.call randomset(rnnow)if(rnnow.le.0.5) thenekein=1.173elseekein=1.333end if{ Modify output statement concerning the source energy.Changewrite(6,340) ekin340 FORMAT(' Results for ',G15.5,'MeV photon'/)towrite(6,340)340 FORMAT(' Results for Co-60 gamma-ray (1.173 and 1.333 MeV)'/) Run ucnai<strong>cg</strong>v2.f by egs5run.{ In the case of Linux or CygwinEnter ucnai<strong>cg</strong>v2 as the user code.Simply enter "return" as the le name for unit 4 and 25.Enter 1 for "Does this user code read from the terminal?".{ In the case of DOSegs5run ucnai<strong>cg</strong>v2 Compare the calculated results with nai.out. Conrm that there are 2peaks in the response corresponding to the source energies.3. Isotropic source cp ucnai<strong>cg</strong>v.f ucnai<strong>cg</strong>v3.f cp ucnai<strong>cg</strong>v.data ucnai<strong>cg</strong>v3.data cp ucnai<strong>cg</strong>v.inp ucnai<strong>cg</strong>v3.inp Modify ucnai<strong>cg</strong>v3.f as follows:{ Add a source direction sampling routine.Change! ----------------------! Select incident angle! ----------------------to! ----------------------! Select incident angle! ----------------------275 call randomset(rnnow)zi0=rnnowcall randomset(rnnow)xi0=2.0*rnnow-1.0call randomset(rnnow)yi0=2.0*rnnow-1.0rr0=dsqrt(xi0*xi0+yi0*yi0+zi0*zi0)if(rr0.gt.1.0) go to 275win = zi0/rr0uin = xi0/rr0vin = yi0/rr042