cg Version

! ----------------------------------------------------! Verify the normalization of source direction cosines! ----------------------------------------------------if(abs(uin*uin+vin*vin+win*win-1.0).gt.1.e-6) thenwrite(6,fmt="(' Following source direction cosines are not',1 ' normalized.',3e12.5)")uin,vin,winstopend if! ==========================================================call shower (iqin,etot,xin,yin,zin,uin,vin,win,irin,wtin)! ==========================================================! If some energy is deposited inside detector add pulse-height! and efficiency.if (depe .gt. 0.D0) thenie=depe/deltae + 1if (ie .gt. 50) ie = 50phs(ie)=phs(ie)+wtinph2s(ie)=ph2s(ie)+wtin*wtintefs=tefs + wtintef2s=tef2s + wtin*wtinif(depe .ge. ekein*0.999) thenpefs=pefs +wtinpef2s=pef2s +wtin*wtinend ifdepe = 0.D0end ifdo ntype=1,3do ie=1,50specs(ntype,ie)=specs(ntype,ie)+spec(ntype,ie)spec2s(ntype,ie)=spec2s(ntype,ie)+* spec(ntype,ie)*spec(ntype,ie)spec(ntype,ie)=0.D0end doend doncount = ncount + 1! ========================if(iwatch.gt.0) call swatch(-1,iwatch)! ========================end do! ========================if(iwatch.gt.0) call swatch(-88,iwatch)! ========================! Count total number of actual cases! -----------------------! End of CALL SHOWER loop! -----------------------call plotxyz(99,0,0,0.D0,0.D0,0.D0,0.D0,0,0.D0,0.D0)write(39,fmt="('9')")4.8.1 ! Set end of batch for CG Viewx ucnaicgv.f MCNP N xi i 12