diff options
Diffstat (limited to 'Dragon/data/C2D20_proc/GeoACRC2DN.c2m')
| -rw-r--r-- | Dragon/data/C2D20_proc/GeoACRC2DN.c2m | 1293 |
1 files changed, 1293 insertions, 0 deletions
diff --git a/Dragon/data/C2D20_proc/GeoACRC2DN.c2m b/Dragon/data/C2D20_proc/GeoACRC2DN.c2m new file mode 100644 index 0000000..145768f --- /dev/null +++ b/Dragon/data/C2D20_proc/GeoACRC2DN.c2m @@ -0,0 +1,1293 @@ +*DECK GeoACRC2D.c2m +*---- +* Name : GeoACRC2D.c2m +* Type : DRAGON procedure +* Use : Generate and track using EXCELT +* the ACR 2D cell models +* Author : G. Marleau +* Date : 2005/02/10 (Revision 0 by G. M.) +* 2005/07/19 (Revision 1 by G. M.) +* Reference : Report IGE-276 +* +* Procedure called as: +* +* Tracking Intline := GeoACRC2D :: +* <<Model>> <<Option>> <<State>> <<BCtrk>> <<DiscrOpt>> +* <<Fuel>> <<Cool>> <<PT>> <<CT>> <<Mode>> +* <<NbAngles>> <<TrkDens>> +* ; +* +* Parameters of the procedure +* Tracking : Tracking data structures +* Intline : Integration lines +* Input data: +* Model : Fuel 2D cell models (see report IGE-276) +* C2D1 -> 2D Cartesian cluster +* C2D2 -> 2D annular cluster +* C2D3 -> 2D Cartesian cell +* C2D4 -> 2D Cartesian cluster with independent pins +* C2D5 -> 2D Cartesian cluster for distributed self-shielding +* C2D6 -> 2D Cartesian cluster (WIMS discretization) +* Option : Geometry option: +* "Flux" -> Fine geometry for flux calculation +* "Shield" -> Coarse geometry for self-shielding +* calculation +* State : State of coolant +* "Present" -> Coolant present at reference density +* (coolant is mixture 1) +* "Part-void" -> Coolant has reduced density +* (coolant is mixture 11) +* "Absent" -> Coolant absent +* (coolant is mixture 0) +* BCtrk : Type of external boundary conditions +* used only for C2D1, C2D3 and C2D4 (TSPC or TISO) +* For Option = "Flux" the following mesh discretisation options +* must be provided +* DiscrOpt = Type of discretisation provided +* "Split" for integer split options +* "Thickness" for real region thickness +* For DiscrOpt = "Thickness" +* Fuel : real fuel radial Delta R +* negative value recommended as per SPLITR +* Cool : real coolant rings Delta R +* positive value recommended. +* PT : real pressure tube Delta R +* positive value recommended. +* CT : real calandria tube Delta R +* positive value recommended. +* Mode : real moderator Delta R +* positive value recommended. +* For DiscrOpt = "Split" +* Fuel : integer SPLITR for fuel +* negative value recommended as per SPLITR +* Cool : integer SPLITR for coolant +* positive value recommended. +* PT : integer SPLITR for pressure tube +* positive value recommended. +* CT : integer SPLITR for calandria tube +* positive value recommended. +* Mode : integer SPLITR for moderator +* positive value recommended. +* NbAngles : Number of angles for tracking. +* TrkDens : Tracking density. +*---- +* Definition of the input and output data structures and +* of the input and output parameters for this procedure +* +PARAMETER Tracking Intline :: + ::: LINKED_LIST Tracking ; + ::: SEQ_BINARY Intline ; ; +STRING Model Option State BCtrk DiscrOpt ; +:: >>Model<< >>Option<< >>State<< >>BCtrk<< ; +REAL FuelDR CoolDR PTDR CTDR ModeDR := + 0.0 0.0 0.0 0.0 0.0 ; +INTEGER FuelS CoolS PTS CTS ModeS := + 0 0 0 0 0 ; +INTEGER LocalIprt Isplit := 0 0 ; +*---- +* Validate procedure parameters +* Model and BCtrk +*---- +IF Model "C2D1" = Model "C2D3" = Model "C2D4" = + + THEN + IF LocalIprt 1 = THEN + ECHO "GeoACRC2D: Model is :" Model ; + ENDIF ; + IF BCtrk "TISO" = BCtrk "TSPC" = + THEN + ECHO "GeoACRC2D: BCtrk is :" BCtrk ; + ELSE + ECHO "WARNING: BCtrk in GeoACRC2D is invalid : " BCtrk ; + ENDIF ; +ELSEIF Model "C2D2" = THEN + IF LocalIprt 1 = THEN + ECHO "GeoACRC2D: Model is :" Model ; + ENDIF ; + IF BCtrk "TISO" = THEN + IF LocalIprt 1 = THEN + ECHO "GeoACRC2D: BCtrk is :" BCtrk ; + ENDIF ; + ELSE + ECHO "WARNING: BCtrk in GeoACRC2D is invalid : " BCtrk ; + ENDIF ; +ELSE + ECHO "WARNING: Model in GeoACRC2D is invalid : " Model ; +ENDIF ; +*---- +* Option and DiscrOpt +*---- +IF Option "Flux" = THEN + IF LocalIprt 1 = THEN + ECHO "GeoACRC2D: Option is :" Option ; + ENDIF ; + :: >>DiscrOpt<< ; + IF LocalIprt 1 = THEN + ECHO "GeoACRC2D: DiscrOpt is :" DiscrOpt ; + ENDIF ; + IF DiscrOpt "Thickness" = THEN + :: >>FuelDR<< >>CoolDR<< >>PTDR<< >>CTDR<< >>ModeDR<< ; + EVALUATE Isplit := 1 ; + ELSEIF DiscrOpt "Split" = THEN + :: >>FuelS<< >>CoolS<< >>PTS<< >>CTS<< >>ModeS<< ; + EVALUATE Isplit := -1 ; + ELSE + ECHO "WARNING: DiscrOpt in GeoACRC2D is invalid : " DiscrOpt ; + ENDIF ; +ELSEIF Option "Shield" = THEN + IF LocalIprt 1 = THEN + ECHO "GeoACRC2D: Option is :" Option ; + ENDIF ; +ELSE + ECHO "WARNING: Option in GeoACRC2D is invalid : " Option ; +ENDIF ; +*---- +* Get tracking parameters +*---- +INTEGER NbAngles ; +REAL TrkDens ; +:: >>NbAngles<< >>TrkDens<< ; +*---- +* State +*---- +IF State "Present" = State "Part-void" = + State "Absent" = + + THEN + IF LocalIprt 1 = THEN + ECHO "GeoACRC2D: State is :" State ; + ENDIF ; +ELSE + ECHO "WARNING: State in GeoACRC2D is invalid : " State ; +ENDIF ; +*---- +* Local modules used in this procedure +*---- +MODULE GEO: NXT: DELETE: ; +*---- +* Local procedures used in this procedure +*---- +PROCEDURE DimACRFuel DimACRCell ; +*---- +* Local data structures used in this procedure +*---- +LINKED_LIST LocalGeo ; +*---- +* Some parameters used in this procedure +*---- +STRING Title ; +STRING LocTrk := "TISO" ; +REAL eps1 eps2 := 0.05 0.05 ; +*---- +* ACR fuel dimensions +*---- +INTEGER Ring ; +INTEGER npinp1 npinp2 npinp3 npinp4 ; +REAL rradf1 rrads1 rpinp1 apinp1 ; +REAL rradf2 rrads2 rpinp2 apinp2 ; +REAL rradf3 rrads3 rpinp3 apinp3 ; +REAL rradf4 rrads4 rpinp4 apinp4 ; +*---- +* Central (first) fuel ring +*---- +EVALUATE Ring := 1 ; +DimACRFuel :: <<Ring>> + >>npinp1<< >>rpinp1<< >>apinp1<< >>rradf1<< >>rrads1<< ; +*---- +* Second fuel ring +*---- +EVALUATE Ring := 2 ; +DimACRFuel :: <<Ring>> + >>npinp2<< >>rpinp2<< >>apinp2<< >>rradf2<< >>rrads2<< ; +*---- +* Third fuel ring +*---- +EVALUATE Ring := 3 ; +DimACRFuel :: <<Ring>> + >>npinp3<< >>rpinp3<< >>apinp3<< >>rradf3<< >>rrads3<< ; +*---- +* Last (fourth) fuel ring +*---- +EVALUATE Ring := 4 ; +DimACRFuel :: <<Ring>> + >>npinp4<< >>rpinp4<< >>apinp4<< >>rradf4<< >>rrads4<< ; +*---- +* Overall cell dimensions +*---- +REAL rcool4 rpt rgap rct LPitch BLength ; +DimACRCell :: + >>rcool4<< >>rpt<< >>rgap<< >>rct<< >>LPitch<< >>BLength<< ; +*---- +* Derived dimensions +*---- +REAL frstann := rrads1 eps1 + ; +REAL rcool1 rcool2 rcool3 := + frstann rpinp2 rrads2 - + 2.0 / + rpinp2 rrads2 + rpinp3 rrads3 - + 2.0 / + rpinp3 rrads3 + rpinp4 rrads4 - + 2.0 / ; +REAL inttann := rcool1 eps1 + ; +REAL xymin xymax := LPitch 2.0 / -1.0 * LPitch 2.0 / ; +REAL equivrad := LPitch 2. ** $Pi_R / SQRT ; +REAL lastann := LPitch 2.0 / eps2 - ; +REAL Modratio := lastann rct - equivrad lastann - / ; +REAL xyimax := lastann 2.0 SQRT / eps2 + ; +REAL xyimin := xyimax -1.0 * ; +*---- +* For individual pins angular separation is apind(i) +*---- +REAL apind1 apind2 apind3 apind4 := + 0.0 + $Pi_R 2.0 * 7.0 / + $Pi_R 2.0 * 14.0 / + $Pi_R 2.0 * 21.0 / ; +REAL apin1a := apinp1 ; +REAL apin2a apin2b apin2c apin2d apin2e apin2f apin2g := + apinp2 apind2 0.0 * + + apinp2 apind2 1.0 * + + apinp2 apind2 2.0 * + + apinp2 apind2 3.0 * + + apinp2 apind2 4.0 * + + apinp2 apind2 5.0 * + + apinp2 apind2 6.0 * + ; +REAL apin3a apin3b apin3c apin3d apin3e apin3f apin3g + apin3h apin3i apin3j apin3k apin3l apin3m apin3n := + apinp3 apind3 0.0 * + + apinp3 apind3 1.0 * + + apinp3 apind3 2.0 * + + apinp3 apind3 3.0 * + + apinp3 apind3 4.0 * + + apinp3 apind3 5.0 * + + apinp3 apind3 6.0 * + + apinp3 apind3 7.0 * + + apinp3 apind3 8.0 * + + apinp3 apind3 9.0 * + + apinp3 apind3 10.0 * + + apinp3 apind3 11.0 * + + apinp3 apind3 12.0 * + + apinp3 apind3 13.0 * + ; +REAL apin4a apin4b apin4c apin4d apin4e apin4f apin4g + apin4h apin4i apin4j apin4k apin4l apin4m apin4n + apin4o apin4p apin4q apin4r apin4s apin4t apin4u := + apinp4 apind4 0.0 * + + apinp4 apind4 1.0 * + + apinp4 apind4 2.0 * + + apinp4 apind4 3.0 * + + apinp4 apind4 4.0 * + + apinp4 apind4 5.0 * + + apinp4 apind4 6.0 * + + apinp4 apind4 7.0 * + + apinp4 apind4 8.0 * + + apinp4 apind4 9.0 * + + apinp4 apind4 10.0 * + + apinp4 apind4 11.0 * + + apinp4 apind4 12.0 * + + apinp4 apind4 13.0 * + + apinp4 apind4 14.0 * + + apinp4 apind4 15.0 * + + apinp4 apind4 16.0 * + + apinp4 apind4 17.0 * + + apinp4 apind4 18.0 * + + apinp4 apind4 19.0 * + + apinp4 apind4 20.0 * + ; +*---- +* Spatial splitting +*---- +INTEGER CSplitR1 CSplitR2 CSplitR3 CSplitR4 := + 1 1 1 1 ; +REAL DeltaRC1 DeltaRC2 DeltaRC3 DeltaRC4 ; +INTEGER FSplitR1 FSplitR2 FSplitR3 FSplitR4 := + 1 1 1 1 ; +REAL DeltaRF1 DeltaRF2 DeltaRF3 DeltaRF4 ; +INTEGER PTSplitR CTSplitR := + 1 1 ; +REAL DeltaPT DeltaCT ; +INTEGER MSplitI MSplitO := + 1 1 ; +REAL DeltaMI DeltaMO ; +IF Isplit -1 = THEN +*---- +* Coolant +*---- + IF CoolS 0 > THEN + EVALUATE CSplitR1 CSplitR2 CSplitR3 CSplitR4 := + CoolS CoolS CoolS CoolS ; + ELSEIF CoolS 0 < THEN + EVALUATE CSplitR1 CSplitR2 CSplitR3 CSplitR4 := + CoolS CoolS CoolS CoolS ; + ENDIF ; +*---- +* Fuel +*---- + IF FuelS 0 > THEN + EVALUATE FSplitR1 FSplitR2 FSplitR3 FSplitR4 := + FuelS FuelS FuelS FuelS ; + ELSEIF FuelS 0 < THEN + EVALUATE FSplitR1 FSplitR2 FSplitR3 FSplitR4 := + FuelS FuelS FuelS FuelS ; + ENDIF ; +*---- +* PT and CT +*---- + IF PTS 0 > THEN + EVALUATE PTSplitR := PTS ; + ELSEIF PTS 0 < THEN + EVALUATE PTSplitR := PTS ; + ENDIF ; + IF CTS 0 > THEN + EVALUATE CTSplitR := CTS ; + ELSEIF CTS 0 < THEN + EVALUATE CTSplitR := CTS ; + ENDIF ; +*---- +* Moderator +*---- + IF ModeS 0 > THEN + EVALUATE MSplitI MSplitO := ModeS ModeS ; + ELSEIF ModeS 0 < THEN + EVALUATE MSplitI MSplitO := ModeS ModeS ; + ENDIF ; +ELSEIF Isplit 1 = THEN +*---- +* Coolant +*---- + IF CoolDR 0.0 > THEN + EVALUATE DeltaRC1 DeltaRC2 DeltaRC3 DeltaRC4 := + rcool1 frstann - CoolDR / + rcool2 rcool1 - CoolDR / + rcool3 rcool2 - CoolDR / + rcool4 rcool3 - CoolDR / ; + EVALUATE CSplitR1 CSplitR2 CSplitR3 CSplitR4 := + DeltaRC1 R_TO_I 1 + + DeltaRC2 R_TO_I 1 + + DeltaRC3 R_TO_I 1 + + DeltaRC4 R_TO_I 1 + ; + ELSEIF CoolDR 0.0 < THEN + EVALUATE DeltaRC1 DeltaRC2 DeltaRC3 DeltaRC4 := + rcool1 frstann - CoolDR / + rcool2 rcool1 - CoolDR / + rcool3 rcool2 - CoolDR / + rcool4 rcool3 - CoolDR / ; + EVALUATE CSplitR1 CSplitR2 CSplitR3 CSplitR4 := + DeltaRC1 R_TO_I 1 - + DeltaRC2 R_TO_I 1 - + DeltaRC3 R_TO_I 1 - + DeltaRC4 R_TO_I 1 - ; + ENDIF ; +*---- +* Fuel +*---- + IF FuelDR 0.0 > THEN + EVALUATE DeltaRF1 DeltaRF2 DeltaRF3 DeltaRF4 := + rradf1 FuelDR / + rradf2 FuelDR / + rradf3 FuelDR / + rradf4 FuelDR / ; + EVALUATE FSplitR1 FSplitR2 FSplitR3 FSplitR4 := + DeltaRF1 R_TO_I 1 + + DeltaRF2 R_TO_I 1 + + DeltaRF3 R_TO_I 1 + + DeltaRF4 R_TO_I 1 + ; + ELSEIF FuelDR 0.0 < THEN + EVALUATE DeltaRF1 DeltaRF2 DeltaRF3 DeltaRF4 := + rradf1 FuelDR / + rradf2 FuelDR / + rradf3 FuelDR / + rradf4 FuelDR / ; + EVALUATE FSplitR1 FSplitR2 FSplitR3 FSplitR4 := + DeltaRF1 R_TO_I 1 - + DeltaRF2 R_TO_I 1 - + DeltaRF3 R_TO_I 1 - + DeltaRF4 R_TO_I 1 - ; + ENDIF ; +*---- +* PT and CT +*---- + IF PTDR 0.0 > THEN + EVALUATE DeltaPT := rpt rcool4 - PTDR / ; + EVALUATE PTSplitR := DeltaPT R_TO_I 1 + ; + ELSEIF PTDR 0.0 < THEN + EVALUATE DeltaPT := rpt rcool4 - PTDR / ; + EVALUATE PTSplitR := DeltaPT R_TO_I 1 - ; + ENDIF ; + IF CTDR 0.0 > THEN + EVALUATE DeltaCT := rct rgap - CTDR / ; + EVALUATE CTSplitR := DeltaCT R_TO_I 1 + ; + ELSEIF CTDR 0.0 < THEN + EVALUATE DeltaCT := rct rgap - CTDR / ; + EVALUATE CTSplitR := DeltaCT R_TO_I 1 - ; + ENDIF ; +*---- +* Moderator splitting +*---- + IF ModeDR 0.0 > THEN + EVALUATE DeltaMI DeltaMO := + lastann rct - ModeDR / + equivrad lastann - ModeDR / ; + EVALUATE MSplitI MSplitO := + DeltaMI R_TO_I 1 + + DeltaMO R_TO_I 1 + ; + ELSEIF ModeDR 0.0 < THEN + EVALUATE DeltaMI DeltaMO := + lastann rct - ModeDR / + equivrad lastann - ModeDR / ; + EVALUATE MSplitI MSplitO := + DeltaMI R_TO_I 1 - + DeltaMO R_TO_I 1 - ; + ENDIF ; +ENDIF ; +*---- +* Select mixture number for cluster cell (C2D1, C2D2) +*---- +INTEGER mixPT mixGap mixCT mixM := 2 3 4 5 ; +INTEGER mixc mixc1 mixc2 mixc3 mixc4 := 1 14 15 16 17 ; +INTEGER mixs mixs1 mixs2 mixs3 mixs4 := 6 18 19 20 21 ; +INTEGER mixf1 mixf2 mixf3 mixf4 := 7 8 9 10 ; +IF State "Part-void" = THEN + EVALUATE mixc mixc1 mixc2 mixc3 mixc4 := 11 22 23 24 25 ; +ELSEIF State "Absent" = THEN + EVALUATE mixc mixc1 mixc2 mixc3 mixc4 := 0 0 0 0 0 ; +ENDIF ; +*---- +* Select mixture number for homogenized cell (C2D3) +*---- +INTEGER mixhPT mixhGap mixhCT mixhMI mixhMO := 1 2 3 4 5 ; +INTEGER mixhf1 mixhf2 mixhf3 mixhf4 := 6 7 8 9 ; +*---- +* Select mixture number for independent pins cluster cell (C2D4) +*---- +INTEGER mixr1a := mixf1 ; +INTEGER mixr2a mixr2b mixr2c mixr2d mixr2e mixr2f mixr2g := + mixf2 26 27 28 29 30 31 ; +INTEGER mixr3a mixr3b mixr3c mixr3d mixr3e mixr3f mixr3g + mixr3h mixr3i mixr3j mixr3k mixr3l mixr3m mixr3n := + mixf3 32 33 34 35 36 37 + 38 39 40 41 42 43 44 ; +INTEGER mixr4a mixr4b mixr4c mixr4d mixr4e mixr4f mixr4g + mixr4h mixr4i mixr4j mixr4k mixr4l mixr4m mixr4n + mixr4o mixr4p mixr4q mixr4r mixr4s mixr4t mixr4u := + mixf3 45 46 47 48 49 50 + 51 52 53 54 55 56 57 + 58 59 60 61 62 64 64 ; +*---- +* Select mixture number for distributed self-shielding / fuel discretization (C2D5) +*---- +!INTEGER mixf1a mixf1b mixf1c mixf1d := 7 26 30 34 ; +INTEGER mixf1a mixf1b mixf1c mixf1d mixf1e + mixf1f mixf1g mixf1h mixf1i mixf1j mixf1k := + 7 26 30 34 38 + 39 40 41 42 43 44 ; +INTEGER mixf2a mixf2b mixf2c mixf2d := 8 27 31 35 ; +INTEGER mixf3a mixf3b mixf3c mixf3d := 9 28 32 36 ; +INTEGER mixf4a mixf4b mixf4c mixf4d := 10 29 33 37 ; +REAL subf1 subf2 subf3 := 0.5 SQRT 0.8 SQRT 0.95 SQRT ; +REAL subp1 subp2 subp3 subp4 subp5 + subp6 subp7 subp8 subp9 subp10 := + 0.1 SQRT 0.2 SQRT 0.3 SQRT 0.4 SQRT 0.5 SQRT + 0.6 SQRT 0.7 SQRT 0.8 SQRT 0.9 SQRT 0.95 SQRT ; +!REAL rradf1a rradf1b rradf1c := +! rradf1 subf1 * rradf1 subf2 * rradf1 subf3 * ; +REAL rradf2a rradf2b rradf2c := + rradf2 subf1 * rradf2 subf2 * rradf2 subf3 * ; +REAL rradf3a rradf3b rradf3c := + rradf3 subf1 * rradf3 subf2 * rradf3 subf3 * ; +REAL rradf4a rradf4b rradf4c := + rradf4 subf1 * rradf4 subf2 * rradf4 subf3 * ; +REAL rradf1a rradf1b rradf1c rradf1d rradf1e + rradf1f rradf1g rradf1h rradf1i rradf1j := + rradf1 subp1 * rradf1 subp2 * rradf1 subp3 * + rradf1 subp4 * rradf1 subp5 * + rradf1 subp6 * rradf1 subp7 * rradf1 subp8 * + rradf1 subp9 * rradf1 subp10 * ; +*---- +REAL hradf2 hradf3 hradf4 := + 0.5 rradf2 * 0.5 rradf3 * 0.5 rradf4 * ; +*---- +* Select geometry +*---- +IF Model "C2D6" = THEN + EVALUATE LocTrk := BCtrk ; + IF Option "Shield" = THEN + LocalGeo := GEO: :: CARCEL 6 1 1 + EDIT 0 + X- REFL MESHX <<xymin>> <<xymax>> X+ REFL + Y- REFL MESHY <<xymin>> <<xymax>> Y+ REFL + RADIUS 0. <<frstann>> + <<rcool4>> <<rpt>> <<rgap>> <<rct>> <<lastann>> + MIX <<mixc>> + <<mixc>> <<mixPT>> <<mixGap>> <<mixCT>> <<mixM>> <<mixM>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 13 + MIX 26 26 26 26 27 27 27 27 28 28 28 28 <<mixs>> + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + RADIUS 0.00000 0.26847 0.37967 0.46500 0.53694 0.60031 + 0.65761 0.71030 0.75934 0.80540 0.84897 0.89041 + <<rradf1>> <<rrads1>> ; + ::: ROD2 := GEO: TUBE 4 1 2 + MESHY 0.0 <<hradf2>> <<rradf2>> + MESHX 0.0 <<rradf2>> + MIX 30 32 34 <<mixs>> + 29 31 33 <<mixs>> + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + RADIUS 0.00000 0.3050 0.4313 <<rradf2>> <<rrads2>> ; + ::: ROD3 := GEO: TUBE 4 1 2 + MESHY 0.0 <<hradf3>> <<rradf3>> + MESHX 0.0 <<rradf3>> + MIX 36 38 40 <<mixs>> + 35 37 39 <<mixs>> + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + RADIUS 0.00000 0.3050 0.4313 <<rradf3>> <<rrads3>> ; + ::: ROD4 := GEO: TUBE 5 1 2 + MESHY 0.0 <<hradf4>> <<rradf4>> + MESHX 0.0 <<rradf4>> + MIX 42 44 46 48 <<mixs>> + 41 43 45 47 <<mixs>> + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + RADIUS 0.00000 0.2641 0.3735 0.4575 <<rradf4>> <<rrads4>> ; + ; + ELSE + LocalGeo := GEO: :: CARCEL 9 1 1 + EDIT 0 + X- REFL MESHX <<xymin>> <<xymax>> X+ REFL + Y- REFL MESHY <<xymin>> <<xymax>> Y+ REFL + RADIUS 0. <<frstann>> <<rcool1>> <<rcool2>> <<rcool3>> + <<rcool4>> <<rpt>> <<rgap>> <<rct>> <<lastann>> + MIX <<mixc1>> <<mixc1>> <<mixc2>> <<mixc3>> + <<mixc4>> <<mixPT>> <<mixGap>> <<mixCT>> <<mixM>> <<mixM>> + SPLITR 1 <<CSplitR1>> <<CSplitR2>> <<CSplitR3>> + <<CSplitR4>> <<PTSplitR>> 1 <<CTSplitR>> <<MSplitI>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 13 + MIX 26 26 26 26 27 27 27 27 28 28 28 28 <<mixs1>> + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + RADIUS 0.00000 0.26847 0.37967 0.46500 0.53694 0.60031 + 0.65761 0.71030 0.75934 0.80540 0.84897 0.89041 + <<rradf1>> <<rrads1>> ; + ::: ROD2 := GEO: TUBE 4 1 2 + MESHY 0.0 <<hradf2>> <<rradf2>> + MESHX 0.0 <<rradf2>> + MIX 30 32 34 <<mixs2>> + 29 31 33 <<mixs2>> + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + RADIUS 0.00000 0.3050 0.4313 <<rradf2>> <<rrads2>> ; + ::: ROD3 := GEO: TUBE 4 1 2 + MESHY 0.0 <<hradf3>> <<rradf3>> + MESHX 0.0 <<rradf3>> + MIX 36 38 40 <<mixs3>> + 35 37 39 <<mixs3>> + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + RADIUS 0.00000 0.3050 0.4313 <<rradf3>> <<rrads3>> ; + ::: ROD4 := GEO: TUBE 5 1 2 + MESHY 0.0 <<hradf4>> <<rradf4>> + MESHX 0.0 <<rradf4>> + MIX 42 44 46 48 <<mixs4>> + 41 43 45 47 <<mixs4>> + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + RADIUS 0.00000 0.2641 0.3735 0.4575 <<rradf4>> <<rrads4>> ; + ; + ENDIF ; +ELSEIF Model "C2D5" = THEN + EVALUATE LocTrk := BCtrk ; + IF Option "Shield" = THEN + LocalGeo := GEO: :: CARCEL 6 1 1 + EDIT 0 + X- REFL MESHX <<xymin>> <<xymax>> X+ REFL + Y- REFL MESHY <<xymin>> <<xymax>> Y+ REFL + RADIUS 0. <<frstann>> + <<rcool4>> <<rpt>> <<rgap>> <<rct>> <<lastann>> + MIX <<mixc>> + <<mixc>> <<mixPT>> <<mixGap>> <<mixCT>> <<mixM>> <<mixM>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 12 + MIX <<mixf1a>> <<mixf1b>> <<mixf1c>> <<mixf1d>> <<mixf1e>> + <<mixf1f>> <<mixf1g>> <<mixf1h>> <<mixf1i>> <<mixf1j>> + <<mixf1k>> <<mixs>> + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + RADIUS 0.00000 <<rradf1a>> <<rradf1b>> + <<rradf1c>> <<rradf1d>> <<rradf1e>> + <<rradf1f>> <<rradf1g>> <<rradf1h>> + <<rradf1i>> <<rradf1j>> <<rradf1>> <<rrads1>> ; +! ::: ROD1 := GEO: TUBE 5 +! MIX <<mixf1a>> <<mixf1b>> <<mixf1c>> <<mixf1d>> <<mixs>> +! NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> +! RADIUS 0.00000 <<rradf1a>> <<rradf1b>> +! <<rradf1c>> <<rradf1>> <<rrads1>> ; + ::: ROD2 := GEO: TUBE 5 + MIX <<mixf2a>> <<mixf2b>> <<mixf2c>> <<mixf2d>> <<mixs>> + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + RADIUS 0.00000 <<rradf2a>> <<rradf2b>> + <<rradf2c>> <<rradf2>> <<rrads2>> ; + ::: ROD3 := GEO: TUBE 5 + MIX <<mixf3a>> <<mixf3b>> <<mixf3c>> <<mixf3d>> <<mixs>> + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + RADIUS 0.00000 <<rradf3a>> <<rradf3b>> + <<rradf3c>> <<rradf3>> <<rrads3>> ; + ::: ROD4 := GEO: TUBE 5 + MIX <<mixf4a>> <<mixf4b>> <<mixf4c>> <<mixf4d>> <<mixs>> + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + RADIUS 0.00000 <<rradf4a>> <<rradf4b>> + <<rradf4c>> <<rradf4>> <<rrads4>> ; + ; + ELSE + LocalGeo := GEO: :: CARCEL 9 1 1 + EDIT 0 + X- REFL MESHX <<xymin>> <<xymax>> X+ REFL + Y- REFL MESHY <<xymin>> <<xymax>> Y+ REFL + RADIUS 0. <<frstann>> <<rcool1>> <<rcool2>> <<rcool3>> + <<rcool4>> <<rpt>> <<rgap>> <<rct>> <<lastann>> + MIX <<mixc1>> <<mixc1>> <<mixc2>> <<mixc3>> + <<mixc4>> <<mixPT>> <<mixGap>> <<mixCT>> <<mixM>> <<mixM>> + SPLITR 1 <<CSplitR1>> <<CSplitR2>> <<CSplitR3>> + <<CSplitR4>> <<PTSplitR>> 1 <<CTSplitR>> <<MSplitI>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 12 + MIX <<mixf1a>> <<mixf1b>> <<mixf1c>> <<mixf1d>> <<mixf1e>> + <<mixf1f>> <<mixf1g>> <<mixf1h>> <<mixf1i>> <<mixf1j>> + <<mixf1k>> <<mixs1>> + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + RADIUS 0.00000 <<rradf1a>> <<rradf1b>> + <<rradf1c>> <<rradf1d>> <<rradf1e>> + <<rradf1f>> <<rradf1g>> <<rradf1h>> + <<rradf1i>> <<rradf1j>> <<rradf1>> <<rrads1>> ; +! ::: ROD1 := GEO: TUBE 5 +! MIX <<mixf1a>> <<mixf1b>> <<mixf1c>> <<mixf1d>> <<mixs1>> +! NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> +! RADIUS 0.00000 <<rradf1a>> <<rradf1b>> +! <<rradf1c>> <<rradf1>> <<rrads1>> ; + ::: ROD2 := GEO: TUBE 5 + MIX <<mixf2a>> <<mixf2b>> <<mixf2c>> <<mixf2d>> <<mixs2>> + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + RADIUS 0.00000 <<rradf2a>> <<rradf2b>> + <<rradf2c>> <<rradf2>> <<rrads2>> ; + ::: ROD3 := GEO: TUBE 5 + MIX <<mixf3a>> <<mixf3b>> <<mixf3c>> <<mixf3d>> <<mixs3>> + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + RADIUS 0.00000 <<rradf3a>> <<rradf3b>> + <<rradf3c>> <<rradf3>> <<rrads3>> ; + ::: ROD4 := GEO: TUBE 5 + MIX <<mixf4a>> <<mixf4b>> <<mixf4c>> <<mixf4d>> <<mixs4>> + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + RADIUS 0.00000 <<rradf4a>> <<rradf4b>> + <<rradf4c>> <<rradf4>> <<rrads4>> ; + ; + ENDIF ; +ELSEIF Model "C2D1" = THEN + EVALUATE LocTrk := BCtrk ; + IF Option "Shield" = THEN + LocalGeo := GEO: :: CARCEL 6 1 1 + EDIT 0 + X- REFL MESHX <<xymin>> <<xymax>> X+ REFL + Y- REFL MESHY <<xymin>> <<xymax>> Y+ REFL + RADIUS 0. <<frstann>> + <<rcool4>> <<rpt>> <<rgap>> <<rct>> <<lastann>> + MIX <<mixc>> + <<mixc>> <<mixPT>> <<mixGap>> <<mixCT>> <<mixM>> <<mixM>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 + MIX <<mixf1>> <<mixs>> NPIN <<npinp1>> + RPIN <<rpinp1>> APIN <<apinp1>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> ; + ::: ROD2 := GEO: TUBE 2 + MIX <<mixf2>> <<mixs>> NPIN <<npinp2>> + RPIN <<rpinp2>> APIN <<apinp2>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD3 := GEO: TUBE 2 + MIX <<mixf3>> <<mixs>> NPIN <<npinp3>> + RPIN <<rpinp3>> APIN <<apinp3>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD4 := GEO: TUBE 2 + MIX <<mixf4>> <<mixs>> NPIN <<npinp4>> + RPIN <<rpinp4>> APIN <<apinp4>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ; + ELSE + LocalGeo := GEO: :: CARCEL 9 1 1 + EDIT 0 + X- REFL MESHX <<xymin>> <<xymax>> X+ REFL + Y- REFL MESHY <<xymin>> <<xymax>> Y+ REFL + RADIUS 0. <<frstann>> <<rcool1>> <<rcool2>> <<rcool3>> + <<rcool4>> <<rpt>> <<rgap>> <<rct>> <<lastann>> + MIX <<mixc1>> <<mixc1>> <<mixc2>> <<mixc3>> + <<mixc4>> <<mixPT>> <<mixGap>> <<mixCT>> <<mixM>> <<mixM>> + SPLITR 1 <<CSplitR1>> <<CSplitR2>> <<CSplitR3>> + <<CSplitR4>> <<PTSplitR>> 1 <<CTSplitR>> <<MSplitI>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 + MIX <<mixf1>> <<mixs1>> NPIN <<npinp1>> + RPIN <<rpinp1>> APIN <<apinp1>> + SPLITR <<FSplitR1>> 1 + RADIUS 0.00000 <<rradf1>> <<rrads1>> ; + ::: ROD2 := GEO: TUBE 2 + MIX <<mixf2>> <<mixs2>> NPIN <<npinp2>> + RPIN <<rpinp2>> APIN <<apinp2>> + SPLITR <<FSplitR2>> 1 + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD3 := GEO: TUBE 2 + MIX <<mixf3>> <<mixs3>> NPIN <<npinp3>> + RPIN <<rpinp3>> APIN <<apinp3>> + SPLITR <<FSplitR3>> 1 + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD4 := GEO: TUBE 2 + MIX <<mixf4>> <<mixs4>> NPIN <<npinp4>> + RPIN <<rpinp4>> APIN <<apinp4>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ; + ENDIF ; +ELSEIF Model "C2D2" = THEN + IF Option "Shield" = THEN + LocalGeo := GEO: :: TUBE 7 + EDIT 0 + R+ REFL + RADIUS 0. <<frstann>> + <<rcool4>> <<rpt>> <<rgap>> <<rct>> <<lastann>> + <<equivrad>> + MIX <<mixc>> + <<mixc>> <<mixPT>> <<mixGap>> <<mixCT>> <<mixM>> <<mixM>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 + MIX <<mixf1>> <<mixs>> NPIN <<npinp1>> + RPIN <<rpinp1>> APIN <<apinp1>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> ; + ::: ROD2 := GEO: TUBE 2 + MIX <<mixf2>> <<mixs>> NPIN <<npinp2>> + RPIN <<rpinp2>> APIN <<apinp2>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD3 := GEO: TUBE 2 + MIX <<mixf3>> <<mixs>> NPIN <<npinp3>> + RPIN <<rpinp3>> APIN <<apinp3>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD4 := GEO: TUBE 2 + MIX <<mixf4>> <<mixs>> NPIN <<npinp4>> + RPIN <<rpinp4>> APIN <<apinp4>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ; + ELSE + LocalGeo := GEO: :: TUBE 10 + EDIT 0 + R+ REFL + RADIUS 0. <<frstann>> <<rcool1>> <<rcool2>> <<rcool3>> + <<rcool4>> <<rpt>> <<rgap>> <<rct>> <<lastann>> + <<equivrad>> + MIX <<mixc1>> <<mixc1>> <<mixc2>> <<mixc3>> + <<mixc4>> <<mixPT>> <<mixGap>> <<mixCT>> <<mixM>> <<mixM>> + SPLITR 1 <<CSplitR1>> <<CSplitR2>> <<CSplitR3>> + <<CSplitR4>> <<PTSplitR>> 1 <<CTSplitR>> <<MSplitI>> + <<MSplitO>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 + MIX <<mixf1>> <<mixs1>> NPIN <<npinp1>> + RPIN <<rpinp1>> APIN <<apinp1>> + SPLITR <<FSplitR1>> 1 + RADIUS 0.00000 <<rradf1>> <<rrads1>> ; + ::: ROD2 := GEO: TUBE 2 + MIX <<mixf2>> <<mixs2>> NPIN <<npinp2>> + RPIN <<rpinp2>> APIN <<apinp2>> + SPLITR <<FSplitR2>> 1 + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD3 := GEO: TUBE 2 + MIX <<mixf3>> <<mixs3>> NPIN <<npinp3>> + RPIN <<rpinp3>> APIN <<apinp3>> + SPLITR <<FSplitR3>> 1 + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD4 := GEO: TUBE 2 + MIX <<mixf4>> <<mixs4>> NPIN <<npinp4>> + RPIN <<rpinp4>> APIN <<apinp4>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ; + ENDIF ; +ELSEIF Model "C2D3" = THEN + EVALUATE LocTrk := BCtrk ; + IF Option "Shield" = THEN + LocalGeo := GEO: :: CARCEL 8 3 3 + EDIT 0 + X- REFL MESHX <<xymin>> <<xyimin>> <<xyimax>> <<xymax>> X+ REFL + Y- REFL MESHY <<xymin>> <<xyimin>> <<xyimax>> <<xymax>> Y+ REFL + RADIUS 0. <<rcool1>> <<rcool2>> <<rcool3>> + <<rcool4>> <<rpt>> <<rgap>> <<rct>> <<lastann>> + MIX <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + ; + ELSE + LocalGeo := GEO: :: CARCEL 8 3 3 + EDIT 0 + X- REFL MESHX <<xymin>> <<xyimin>> <<xyimax>> <<xymax>> X+ REFL + Y- REFL MESHY <<xymin>> <<xyimin>> <<xyimax>> <<xymax>> Y+ REFL + RADIUS 0. <<rcool1>> <<rcool2>> <<rcool3>> + <<rcool4>> <<rpt>> <<rgap>> <<rct>> <<lastann>> + MIX <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + <<mixhf1>> <<mixhf2>> <<mixhf3>> <<mixhf4>> + <<mixhPT>> <<mixhGap>> <<mixhCT>> <<mixhMI>> <<mixhMO>> + SPLITR <<FSplitR1>> <<CSplitR2>> <<CSplitR3>> + <<CSplitR4>> <<PTSplitR>> 1 <<CTSplitR>> <<MSplitI>> + SPLITX <<MSplitO>> 1 <<MSplitO>> + SPLITY <<MSplitO>> 1 <<MSplitO>> + ; + ENDIF ; +ELSEIF Model "C2D4" = THEN + EVALUATE LocTrk := BCtrk ; + IF Option "Shield" = THEN + LocalGeo := GEO: :: CARCEL 6 1 1 + EDIT 0 + X- REFL MESHX <<xymin>> <<xymax>> X+ REFL + Y- REFL MESHY <<xymin>> <<xymax>> Y+ REFL + RADIUS 0. <<frstann>> + <<rcool4>> <<rpt>> <<rgap>> <<rct>> <<lastann>> + MIX <<mixc>> + <<mixc>> <<mixPT>> <<mixGap>> <<mixCT>> <<mixM>> <<mixM>> + CLUSTER ROD1a + ROD2a ROD2b ROD2c ROD2d ROD2e ROD2f ROD2g + ROD3a ROD3b ROD3c ROD3d ROD3e ROD3f ROD3g + ROD3h ROD3i ROD3j ROD3k ROD3l ROD3m ROD3n + ROD4a ROD4b ROD4c ROD4d ROD4e ROD4f ROD4g + ROD4h ROD4i ROD4j ROD4k ROD4l ROD4m ROD4n + ROD4o ROD4p ROD4q ROD4r ROD4s ROD4t ROD4u + ::: ROD1a := GEO: TUBE 2 + MIX <<mixr1a>> <<mixs>> NPIN 1 + RPIN <<rpinp1>> APIN <<apin1a>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> ; + ::: ROD2a := GEO: TUBE 2 + MIX <<mixr2a>> <<mixs>> NPIN 1 + RPIN <<rpinp2>> APIN <<apin2a>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD2b := GEO: TUBE 2 + MIX <<mixr2b>> <<mixs>> NPIN 1 + RPIN <<rpinp2>> APIN <<apin2b>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD2c := GEO: TUBE 2 + MIX <<mixr2c>> <<mixs>> NPIN 1 + RPIN <<rpinp2>> APIN <<apin2c>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD2d := GEO: TUBE 2 + MIX <<mixr2d>> <<mixs>> NPIN 1 + RPIN <<rpinp2>> APIN <<apin2d>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD2e := GEO: TUBE 2 + MIX <<mixr2e>> <<mixs>> NPIN 1 + RPIN <<rpinp2>> APIN <<apin2e>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD2f := GEO: TUBE 2 + MIX <<mixr2f>> <<mixs>> NPIN 1 + RPIN <<rpinp2>> APIN <<apin2f>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD2g := GEO: TUBE 2 + MIX <<mixr2g>> <<mixs>> NPIN 1 + RPIN <<rpinp2>> APIN <<apin2g>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD3a := GEO: TUBE 2 + MIX <<mixr3a>> <<mixs>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3a>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3b := GEO: TUBE 2 + MIX <<mixr3b>> <<mixs>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3b>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3c := GEO: TUBE 2 + MIX <<mixr3c>> <<mixs>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3c>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3d := GEO: TUBE 2 + MIX <<mixr3d>> <<mixs>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3d>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3e := GEO: TUBE 2 + MIX <<mixr3e>> <<mixs>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3e>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3f := GEO: TUBE 2 + MIX <<mixr3f>> <<mixs>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3f>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3g := GEO: TUBE 2 + MIX <<mixr3g>> <<mixs>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3g>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3h := GEO: TUBE 2 + MIX <<mixr3h>> <<mixs>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3h>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3i := GEO: TUBE 2 + MIX <<mixr3i>> <<mixs>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3i>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3j := GEO: TUBE 2 + MIX <<mixr3j>> <<mixs>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3j>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3k := GEO: TUBE 2 + MIX <<mixr3k>> <<mixs>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3k>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3l := GEO: TUBE 2 + MIX <<mixr3l>> <<mixs>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3l>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3m := GEO: TUBE 2 + MIX <<mixr3m>> <<mixs>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3m>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3n := GEO: TUBE 2 + MIX <<mixr3n>> <<mixs>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3n>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD4a := GEO: TUBE 2 + MIX <<mixr4a>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4a>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4b := GEO: TUBE 2 + MIX <<mixr4b>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4b>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4c := GEO: TUBE 2 + MIX <<mixr4c>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4c>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4d := GEO: TUBE 2 + MIX <<mixr4d>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4d>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4e := GEO: TUBE 2 + MIX <<mixr4e>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4e>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4f := GEO: TUBE 2 + MIX <<mixr4f>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4f>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4g := GEO: TUBE 2 + MIX <<mixr4g>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4g>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4h := GEO: TUBE 2 + MIX <<mixr4h>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4h>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4i := GEO: TUBE 2 + MIX <<mixr4i>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4i>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4j := GEO: TUBE 2 + MIX <<mixr4j>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4j>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4k := GEO: TUBE 2 + MIX <<mixr4k>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4k>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4l := GEO: TUBE 2 + MIX <<mixr4l>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4l>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4m := GEO: TUBE 2 + MIX <<mixr4m>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4m>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4n := GEO: TUBE 2 + MIX <<mixr4n>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4n>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4o := GEO: TUBE 2 + MIX <<mixr4o>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4o>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4p := GEO: TUBE 2 + MIX <<mixr4p>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4p>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4q := GEO: TUBE 2 + MIX <<mixr4q>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4q>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4r := GEO: TUBE 2 + MIX <<mixr4r>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4r>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4s := GEO: TUBE 2 + MIX <<mixr4s>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4s>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4t := GEO: TUBE 2 + MIX <<mixr4t>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4t>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4u := GEO: TUBE 2 + MIX <<mixr4u>> <<mixs>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4u>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ; + ELSE + LocalGeo := GEO: :: CARCEL 9 1 1 + EDIT 0 + X- REFL MESHX <<xymin>> <<xymax>> X+ REFL + Y- REFL MESHY <<xymin>> <<xymax>> Y+ REFL + RADIUS 0. <<frstann>> <<rcool1>> <<rcool2>> <<rcool3>> + <<rcool4>> <<rpt>> <<rgap>> <<rct>> <<lastann>> + MIX <<mixc1>> <<mixc1>> <<mixc2>> <<mixc3>> + <<mixc4>> <<mixPT>> <<mixGap>> <<mixCT>> <<mixM>> <<mixM>> + SPLITR 1 <<CSplitR1>> <<CSplitR2>> <<CSplitR3>> + <<CSplitR4>> <<PTSplitR>> 1 <<CTSplitR>> <<MSplitI>> + CLUSTER ROD1a + ROD2a ROD2b ROD2c ROD2d ROD2e ROD2f ROD2g + ROD3a ROD3b ROD3c ROD3d ROD3e ROD3f ROD3g + ROD3h ROD3i ROD3j ROD3k ROD3l ROD3m ROD3n + ROD4a ROD4b ROD4c ROD4d ROD4e ROD4f ROD4g + ROD4h ROD4i ROD4j ROD4k ROD4l ROD4m ROD4n + ROD4o ROD4p ROD4q ROD4r ROD4s ROD4t ROD4u + ::: ROD1a := GEO: TUBE 2 + MIX <<mixr1a>> <<mixs1>> NPIN 1 + RPIN <<rpinp1>> APIN <<apin1a>> + SPLITR <<FSplitR1>> 1 + RADIUS 0.00000 <<rradf1>> <<rrads1>> ; + ::: ROD2a := GEO: TUBE 2 + MIX <<mixr2a>> <<mixs2>> NPIN 1 + RPIN <<rpinp2>> APIN <<apin2a>> + SPLITR <<FSplitR2>> 1 + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD2b := GEO: TUBE 2 + MIX <<mixr2b>> <<mixs2>> NPIN 1 + RPIN <<rpinp2>> APIN <<apin2b>> + SPLITR <<FSplitR2>> 1 + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD2c := GEO: TUBE 2 + MIX <<mixr2c>> <<mixs2>> NPIN 1 + RPIN <<rpinp2>> APIN <<apin2c>> + SPLITR <<FSplitR2>> 1 + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD2d := GEO: TUBE 2 + MIX <<mixr2d>> <<mixs2>> NPIN 1 + RPIN <<rpinp2>> APIN <<apin2d>> + SPLITR <<FSplitR2>> 1 + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD2e := GEO: TUBE 2 + MIX <<mixr2e>> <<mixs2>> NPIN 1 + RPIN <<rpinp2>> APIN <<apin2e>> + SPLITR <<FSplitR2>> 1 + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD2f := GEO: TUBE 2 + MIX <<mixr2f>> <<mixs2>> NPIN 1 + RPIN <<rpinp2>> APIN <<apin2f>> + SPLITR <<FSplitR2>> 1 + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD2g := GEO: TUBE 2 + MIX <<mixr2g>> <<mixs2>> NPIN 1 + RPIN <<rpinp2>> APIN <<apin2g>> + SPLITR <<FSplitR2>> 1 + RADIUS 0.00000 <<rradf2>> <<rrads2>> ; + ::: ROD3a := GEO: TUBE 2 + MIX <<mixr3a>> <<mixs3>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3a>> + SPLITR <<FSplitR3>> 1 + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3b := GEO: TUBE 2 + MIX <<mixr3b>> <<mixs3>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3b>> + SPLITR <<FSplitR3>> 1 + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3c := GEO: TUBE 2 + MIX <<mixr3c>> <<mixs3>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3c>> + SPLITR <<FSplitR3>> 1 + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3d := GEO: TUBE 2 + MIX <<mixr3d>> <<mixs3>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3d>> + SPLITR <<FSplitR3>> 1 + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3e := GEO: TUBE 2 + MIX <<mixr3e>> <<mixs3>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3e>> + SPLITR <<FSplitR3>> 1 + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3f := GEO: TUBE 2 + MIX <<mixr3f>> <<mixs3>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3f>> + SPLITR <<FSplitR3>> 1 + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3g := GEO: TUBE 2 + MIX <<mixr3g>> <<mixs3>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3g>> + SPLITR <<FSplitR3>> 1 + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3h := GEO: TUBE 2 + MIX <<mixr3h>> <<mixs3>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3h>> + SPLITR <<FSplitR3>> 1 + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3i := GEO: TUBE 2 + MIX <<mixr3i>> <<mixs3>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3i>> + SPLITR <<FSplitR3>> 1 + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3j := GEO: TUBE 2 + MIX <<mixr3j>> <<mixs3>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3j>> + SPLITR <<FSplitR3>> 1 + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3k := GEO: TUBE 2 + MIX <<mixr3k>> <<mixs3>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3k>> + SPLITR <<FSplitR3>> 1 + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3l := GEO: TUBE 2 + MIX <<mixr3l>> <<mixs3>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3l>> + SPLITR <<FSplitR3>> 1 + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3m := GEO: TUBE 2 + MIX <<mixr3m>> <<mixs3>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3m>> + SPLITR <<FSplitR3>> 1 + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD3n := GEO: TUBE 2 + MIX <<mixr3n>> <<mixs3>> NPIN 1 + RPIN <<rpinp3>> APIN <<apin3n>> + SPLITR <<FSplitR3>> 1 + RADIUS 0.00000 <<rradf3>> <<rrads3>> ; + ::: ROD4a := GEO: TUBE 2 + MIX <<mixr4a>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4a>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4b := GEO: TUBE 2 + MIX <<mixr4b>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4b>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4c := GEO: TUBE 2 + MIX <<mixr4c>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4c>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4d := GEO: TUBE 2 + MIX <<mixr4d>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4d>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4e := GEO: TUBE 2 + MIX <<mixr4e>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4e>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4f := GEO: TUBE 2 + MIX <<mixr4f>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4f>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4g := GEO: TUBE 2 + MIX <<mixr4g>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4g>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4h := GEO: TUBE 2 + MIX <<mixr4h>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4h>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4i := GEO: TUBE 2 + MIX <<mixr4i>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4i>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4j := GEO: TUBE 2 + MIX <<mixr4j>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4j>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4k := GEO: TUBE 2 + MIX <<mixr4k>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4k>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4l := GEO: TUBE 2 + MIX <<mixr4l>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4l>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4m := GEO: TUBE 2 + MIX <<mixr4m>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4m>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4n := GEO: TUBE 2 + MIX <<mixr4n>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4n>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4o := GEO: TUBE 2 + MIX <<mixr4o>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4o>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4p := GEO: TUBE 2 + MIX <<mixr4p>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4p>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4q := GEO: TUBE 2 + MIX <<mixr4q>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4q>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4r := GEO: TUBE 2 + MIX <<mixr4r>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4r>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4s := GEO: TUBE 2 + MIX <<mixr4s>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4s>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4t := GEO: TUBE 2 + MIX <<mixr4t>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4t>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ::: ROD4u := GEO: TUBE 2 + MIX <<mixr4u>> <<mixs4>> NPIN 1 + RPIN <<rpinp4>> APIN <<apin4u>> + SPLITR <<FSplitR4>> 1 + RADIUS 0.00000 <<rradf4>> <<rrads4>> ; + ; + ENDIF ; +ENDIF ; +*---- +* Process the geometry +*---- +EVALUATE Title := Model Option State LocTrk + + + ; +Tracking Intline := NXT: LocalGeo :: + EDIT 1 TITLE <<Title>> ANIS 2 + <<LocTrk>> <<NbAngles>> <<TrkDens>> ; + +*---- +* Clean up +*---- +LocalGeo := DELETE: LocalGeo ; +QUIT "LIST" . +*QUIT . |