cg Version

More documents

Recommendations

Info

include 'include/egs5 h.f' common include 2 include ! ----------------------! Auxiliary-code COMMONs! ----------------------include 'auxcommons/aux_h.f'include 'auxcommons/edata.f'include 'auxcommons/etaly1.f'include 'auxcommons/instuf.f'include 'auxcommons/lines.f'include 'auxcommons/watch.f'! Auxiliary-code "header" file! ------------------! cg related COMMONs! ------------------include 'auxcommons/geom_common.f' ! geom-common fileinteger irinn include CG CG common common/totals/* depe,deltae,spec(3,50),maxpictreal*8 depe,deltae,specinteger maxpict! Variables to scoreimplicit none open egs5 pegs5 pegs 7-26 close openpegs open 7-26 ! ----------! Open files! ----------open(6,FILE='egs5job.out',STATUS='unknown')open(4,FILE='egs5job.inp',STATUS='old')open(39,FILE='egs5job.pic',STATUS='unknown') 39 counters out(0) call 4.2 : Step 2 nmed Characteristic Dimensionpegs5 call medarr 24 24 24 Chracteristic Dimension (1 cm 1cm 1cm 1cm) nmed=4if(nmed.gt.MXMED) thenwrite(6,'(A,I4,A,I4,A/A)')* ' nmed (',nmed,') larger than MXMED (',MXMED,')',* ' MXMED in iclude/egs5_h.f must be increased.'stop2 EGS4 COMIN 7
end if! ==============call block_set! ==============! ---------------------------------! define media before calling PEGS5! ---------------------------------medarr(1)='NAImedarr(2)='ALmedarr(3)='QUARTZmedarr(4)='AIR-AT-NTPdo j=1,nmeddo i=1,24media(i,j)=medarr(j)(i:i)end doend do! Initialize some general variableschard(1) = 3.81d0 ! automatic step-size controlchard(2) = 0.1d0chard(3) = 0.5d0chard(4) = 5.0d0write(6,fmt="('chard =',5e12.5)") (chard(j),j=1,nmed)! -----------------------------------! Run KEK PEGS5 before calling HATCH! -----------------------------------write(6,100)100 FORMAT('PEGS5-call comes next'/)! ==========call pegs5! ==========4.3 : Step 3 npreci 3 CG CG data CG cg (6) geomgt call CG End of CG data ifto 39 PICT (CSTA-FREE CSTA) subroutine geomgt CG CG CEND cg nreg CG write(6,*) 'Read cg-related data'!-----------------------------------------------! initilize cg related parameter!-----------------------------------------------npreci=3 ! PICT data mode for CGView in free formatifti = 4ifto = 39''''! Input unit number for cg-data! Output unit number for PICTwrite(6,fmt="(' CG data')")call geomgt(ifti,6) ! Read in CG datawrite(6,fmt="(' End of CG data',/)")if(npreci.eq.3) write(ifto,fmt="('CSTA-FREE-TIME')")if(npreci.eq.2) write(ifto,fmt="('CSTA-TIME')")8
Page 1 and 2: KEK Internal 2011-5December 2011R/D
Page 3 and 4: Lecture Notes ofResponse calculatio
Page 5 and 6: 3 Details of user code 313.1 Main p
Page 7 and 8: 1 Combinatrial Geometry (CG)2 Combi
Page 9 and 10: ody 3 2 body + ; OR B23A1Figure 2
Page 11: Z4 +4Z5 +5 -3Z6 +6 -5END1 0 2 3 4 0
Page 15 and 16: ! pdf data for many sourcedeltae=0.
Page 17 and 18: ! ---------------------------------
Page 19 and 20: sigph=dsqrt((ph2s(ie)-avph*avph)/nc
Page 21 and 22: 6 6.1 NaI 1. Cs-137 (0.662MeV) 2.
Page 23 and 24: ekein=1.253ekein=1.333{ ! ---------
Page 25 and 26: ROuter vacuum (region 6)Air region
Page 27 and 28: nmed=2 medarr(1)='NAImedarr(2)='ALm
Page 31 and 32: EGS5 sample user code (ucnaicgv.f)R
Page 33 and 34: operator, it means that the region
Page 35 and 36: 2 Outlines of sample user code uccg
Page 37 and 38: include 'auxcommons/edata.f'include
Page 39 and 40: The material assignment is read in
Page 41 and 42: ! Write batch numberwrite(39,fmt="(
Page 43 and 44: ! ---------------------------------
Page 45 and 46: a distance to the boundary of regio
Page 47 and 48: { Add sampling routines for source
Page 51 and 52: to* xi0,yi0,zi0,avab,depes,depe2s,s
Page 55 and 56: Appendix: Full listings of ucnaicgv
Page 57 and 58: open(39,FILE='egs5job.pic',STATUS='
Page 59 and 60: write(6,140)140 FORMAT(/,' HATCH-ca
Page 61 and 62: ilines = ilines + 1write(6,280) eto
Page 63 and 64:
* 14X,' Positron')do ie=1,50elow=de
Page 65 and 66:
! ---------------------------------
Page 67 and 68:
ccyidd=y_np+tval0*v_npzidd=z_np+tva
show all

cg Version

Create successful ePaper yourself

Delete template?

Save as template?