summaryrefslogtreecommitdiff
path: root/Dragon/data/C2D20_proc
diff options
context:
space:
mode:
authorstainer_t <thomas.stainer@oecd-nea.org>2025-09-08 13:48:49 +0200
committerstainer_t <thomas.stainer@oecd-nea.org>2025-09-08 13:48:49 +0200
commit7dfcc480ba1e19bd3232349fc733caef94034292 (patch)
tree03ee104eb8846d5cc1a981d267687a729185d3f3 /Dragon/data/C2D20_proc
Initial commit from Polytechnique Montreal
Diffstat (limited to 'Dragon/data/C2D20_proc')
-rw-r--r--Dragon/data/C2D20_proc/DimACRCell.c2m35
-rw-r--r--Dragon/data/C2D20_proc/DimACRFuel.c2m73
-rw-r--r--Dragon/data/C2D20_proc/GeoACRC2DN.c2m1293
-rw-r--r--Dragon/data/C2D20_proc/MicACRE620.c2m295
-rw-r--r--Dragon/data/C2D20_proc/MicACRE6UH.c2m87
-rw-r--r--Dragon/data/C2D20_proc/MicACRE6UW.c2m114
6 files changed, 1897 insertions, 0 deletions
diff --git a/Dragon/data/C2D20_proc/DimACRCell.c2m b/Dragon/data/C2D20_proc/DimACRCell.c2m
new file mode 100644
index 0000000..0bbe226
--- /dev/null
+++ b/Dragon/data/C2D20_proc/DimACRCell.c2m
@@ -0,0 +1,35 @@
+*DECK DimACRCell.c2m
+*----
+* Name : DimACRCell.c2m
+* Type : DRAGON procedure
+* Use : Dimensions for ACR-1000 cell
+* Author : G. Marleau
+* Revisions : 2006/05/16 -> Replace ACR-700 by ACR-1000
+* Reference : fuel18n.pin and fuel20n.pin
+*
+* Procedure called as:
+*
+* DimACRCell ::
+* >>rcool<< >>rpt<< >>rgap<< >>rct<< >>LPitch<< >>BLength<< ;
+*
+* Output data:
+* rcool : outer coolant radius
+* rpt : outer radius of pressure tube
+* rgap : outer radius of gap
+* rct : outer radius of calandria tube
+* LPitch : Lattice pitch
+* BLength : Bundle length
+*----
+* Definition of the input and output data structures and
+* of the input and output parameters for this procedure
+*
+REAL rinpt thickpt rinct thickct :=
+ 5.20250 0.65 7.75 0.25 ;
+REAL rcool rpt
+ rgap rct
+ LPitch BLength :=
+ rinpt rinpt thickpt +
+ rinct rinct thickct +
+ 24.0 49.53 ;
+:: <<rcool>> <<rpt>> <<rgap>> <<rct>> <<LPitch>> <<BLength>> ;
+QUIT .
diff --git a/Dragon/data/C2D20_proc/DimACRFuel.c2m b/Dragon/data/C2D20_proc/DimACRFuel.c2m
new file mode 100644
index 0000000..c98b5ef
--- /dev/null
+++ b/Dragon/data/C2D20_proc/DimACRFuel.c2m
@@ -0,0 +1,73 @@
+*DECK DimACRFuel.c2m
+*----
+* Name : DimACRFuel.c2m
+* Type : DRAGON procedure
+* Use : Dimensions for ACR-1000 fuel
+* Author : G. Marleau
+* Revisions : 2006/05/16 -> Replace ACR-700 by ACR-1000
+* Reference : fuel18n.pin and fuel20n.pin
+*
+* Description of the procedure:
+* DimACRFuel ::
+* <<Ring>>
+* >>npinp<< >>rpinp<< >>apinp<< >>rradf<< >>rrads<< ;
+* Input data:
+* Ring : Ring number
+* Output data:
+* npinp : number of pins in this ring
+* rpinp : ring radius (cm)
+* apinp : first pin angle (radians)
+* rradf : outer fuel radius (cm)
+* rrads : outer sheat radius (cm)
+*----
+* Definition of the input and output data structures and
+* of the input and output parameters for this procedure
+*
+MODULE END: ;
+INTEGER Ring ;
+ :: >>Ring<< ;
+*----
+* Pin diameters are given in mm. Transform them
+* to radius in cm.
+*----
+REAL thikI thikA :=
+ 0.064 0.044 ;
+*----
+* Number of pins in each ring
+*----
+INTEGER npinp1 npinp2 npinp3 npinp4 :=
+ 1 7 14 21 ;
+*----
+* Radius of pins in each ring
+* Diameter is given in mm. Transform to radius in cm.
+*----
+REAL rradf1 rradf2 rradf3 rradf4 :=
+ 0.93025 0.53025 0.53025 0.53025 ;
+*----
+* Outer sheath radius in each ring
+*----
+REAL rrads1 rrads2 rrads3 rrads4 :=
+ 1.003383 0.573787 0.573787 0.573787 ;
+*----
+* Pins ring center location.
+* Diameter is given in mm. Transform to radius in cm.
+*----
+REAL rpinp1 rpinp2 rpinp3 rpinp4 :=
+ 0.0 1.725 2.953 4.384 ;
+*----
+* Pins angular location.
+*----
+REAL apinp1 apinp2 apinp3 apinp4 :=
+ 0.0 0.0 $Pi_R npinp3 I_TO_R / 0.0 ;
+IF Ring 4 = THEN
+ :: <<npinp4>> <<rpinp4>> <<apinp4>> <<rradf4>> <<rrads4>> ;
+ELSEIF Ring 3 = THEN
+ :: <<npinp3>> <<rpinp3>> <<apinp3>> <<rradf3>> <<rrads3>> ;
+ELSEIF Ring 2 = THEN
+ :: <<npinp2>> <<rpinp2>> <<apinp2>> <<rradf2>> <<rrads2>> ;
+ELSEIF Ring 1 = THEN
+ :: <<npinp1>> <<rpinp1>> <<apinp1>> <<rradf1>> <<rrads1>> ;
+ELSE
+ ECHO "WARNING: Invalid Ring number in DimACRFuel : " Ring ;
+ENDIF ;
+QUIT .
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 .
diff --git a/Dragon/data/C2D20_proc/MicACRE620.c2m b/Dragon/data/C2D20_proc/MicACRE620.c2m
new file mode 100644
index 0000000..d11cca5
--- /dev/null
+++ b/Dragon/data/C2D20_proc/MicACRE620.c2m
@@ -0,0 +1,295 @@
+*DECK MicACRE620.c2m
+*----
+* Name : MicACRE620.c2m
+* Type : DRAGON procedure
+* Use : To generate the ENDF/B-VI Microlib for the ACR
+* for 2.0 % enriched fuel.
+* Author : G. Marleau
+* Date : 2006/05/17 (Revision 0 by G. M.)
+* Reference : fuel20n.pin
+*
+* Procedure called as:
+*
+* MicLib := MicACRE620 ;
+*
+* Parameters of the procedure
+* MicLib : Microlib tabulated for the required parameters stored in an XSM file
+* The minimal set of mixtures included in this microlib is:
+* Mix 0 = void
+* Mix 1 = coolant at reference conditions
+* Mix 2 = pressure tube material
+* Mix 3 = gap material
+* Mix 4 = calendria tube material
+* Mix 5 = moderator
+* Mix 6 = fuel sheat
+* Mix 7 = central (first) Dy+Gd element
+* Mix 8 = second ring of fuel
+* Mix 9 = third ring of fuel
+* Mix 10 = last (fourth) ring of fuel
+* Mix 11 = coolant at 50 % density
+* Mix 12 = fuel at 1.8 % enrichment
+* Mix 13 = fuel at 2.0 % enrichment
+*
+* General properties of the mixtures in file fuel20n.pin
+* A) Coolant (light water): Mix 1
+* purety (D2O content) = 0.0 % (Cprty)
+* Boron = 0.0 % (CB)
+* Gd = 0.0 mg/kg (CGd)
+* Temperature = 573.16 K (Ctemp)
+* Density = computed g/cc (Cdens)
+* B) Pressure tube (ZrNb 2.65): Mix 2
+* Temperature = 573.66 K (PTtemp)
+* Density = 6.5018 g/cc (PTdens)
+* C) Gap (CO2): Mix 3
+* Temperature = 459.66 K (Gaptemp)
+* Density = 0.00128 g/cc (Gapdens)
+* D) Calandria tube (ZrII): Mix 4
+* Temperature = 346.16 K (CTtemp)
+* Density = 6.4635 g/cc (CTdens)
+* E) Moderator (heavy water): Mix 5
+* purety (D2O content) = 99.75 % (Mprty)
+* Boron = 0.0 mg/kg (MB)
+* Gd = 0.0 mg/kg (MGd)
+* Temperature = 346.16 K (Mtemp)
+* Density = computed g/cc (Mdens)
+* F) Fuel sheathing (ZrII): Mix 6
+* Temperature = 573.16 K (Shtemp)
+* Density = 7.028222 g/cc (Shdens)
+* G) Fuel
+* a) Central element (Dy+Gd mixture): Mix 7
+* Temperature = 960.16 K (Ftemp1)
+* Density = 6.09950 g/cc (Fdens1)
+* b) 7 pins ring (fuel 2): Mix 8
+* Enrichment = 2.0 %
+* Dy = 0.0 %
+* Temperature = 960.16 K (Ftemp2)
+* Density = 10.14937 g/cc (Fdens2)
+* c) 14 pins ring (fuel 3): Mix 9
+* Enrichment = 2.0 %
+* Dy = 0.0 %
+* Temperature = 960.16 K (Ftemp3)
+* Density = 10.14937 g/cc (Fdens3)
+* d) 21 pins ring (fuel 4): Mix 10
+* Enrichment = 2.0 %
+* Dy = 0.0 %
+* Temperature = 960.16 K (Ftemp4)
+* Density = 10.14937 g/cc (Fdens4)
+* H) Coolant at half density: Mix 11
+* purety (D2O content) = 0.0 %
+* Boron = 0.0 %
+* Temperature = 573.16 K
+* Density = (Mix 1)/2 g/cc
+*----
+* Definition of output data structures
+*----
+PARAMETER MicLib :: ::: XSM_FILE MicLib ; ;
+*----
+* Local modules used in this procedure
+*----
+MODULE INFO: LIB: DELETE: END: ;
+*----
+* Local parameters used in this procedure
+*----
+REAL Ctemp Cprty Cdens CB CGd := 573.16 0.0 0.0 0.0 0.0 ;
+REAL Mtemp Mprty Mdens MB MGd := 346.16 99.75 0.0 0.0 0.0 ;
+REAL PTtemp PTdens := Ctemp 6.5018 ;
+REAL Gaptemp Gapdens := 459.66 0.00128 ;
+REAL CTtemp CTdens := Mtemp 6.4635 ;
+REAL Shtemp Shdens := Ctemp 6.4635 1.132060 * ;
+REAL Ftemp1 Fdens1 := Ctemp 6.09950 ;
+REAL Ftemp2 Fdens2 := 960.16 10.14937 ;
+REAL Ftemp3 Fdens3 := 960.16 10.14937 ;
+REAL Ftemp4 Fdens4 := 960.16 10.14937 ;
+*----
+* Library and isotopes names
+*----
+STRING LibName := "E6MLIB" ;
+STRING IsoB10 IsoB11 :=
+ "B10" "B11" ;
+STRING IsoZr90 IsoZr91 IsoZr92 IsoZr94 IsoZr96 :=
+ "ZR90" "ZR91" "ZR92" "ZR94" "ZR96" ;
+STRING IsoZr90PT IsoZr91PT IsoZr92PT IsoZr94PT IsoZr96PT :=
+ "ZR90PT" "ZR91PT" "ZR92PT" "ZR94PT" "ZR96PT" ;
+STRING IsoZr90CT IsoZr91CT IsoZr92CT IsoZr94CT IsoZr96CT :=
+ "ZR90CT" "ZR91CT" "ZR92CT" "ZR94CT" "ZR96CT" ;
+STRING IsoZr90CL IsoZr91CL IsoZr92CL IsoZr94CL IsoZr96CL :=
+ "ZR90CL" "ZR91CL" "ZR92CL" "ZR94CL" "ZR96CL" ;
+STRING IsoCr50 IsoCr52 IsoCr53 IsoCr54 :=
+ "CR50" "CR52" "CR53" "CR54" ;
+STRING IsoFe54 IsoFe56 IsoFe57 IsoFe58 :=
+ "FE54" "FE56" "FE57" "FE58" ;
+STRING IsoNi58 IsoNi60 IsoNi61 IsoNi62 IsoNi64 :=
+ "NI58" "NI60" "NI61" "NI62" "NI64" ;
+STRING IsoDy160 IsoDy161 IsoDy162 IsoDy163 IsoDy164 :=
+ "DY160" "DY161" "DY162" "DY163" "DY164" ;
+STRING IsoGd152 IsoGd154 IsoGd155 IsoGd156 IsoGd157 IsoGd158
+ IsoGd160 :=
+ "GD152" "GD154" "GD155" "GD156" "GD157" "GD158"
+ "GD160" ;
+STRING IsoNb :=
+ "NB93" ;
+STRING IsoU4 IsoU5 IsoU8 IsoU6 IsoPu9 IsoO16
+ IsoH1 IsoD2 IsoC IsoMn55 :=
+ "U234" "U235" "U238" "U236" "PU239" "O16"
+ "H1H2O" "DD2O" "C" "MN55" ;
+*----
+* Coolant density and WGT%
+*----
+REAL CwgtH1 CwgtD2 CwgtO16 ;
+INFO: :: EDIT 0
+ TMP: <<Ctemp>> K
+ PUR: <<Cprty>> WGT%
+ CALC DENS WATER >>Cdens<<
+ LIB: WIMSAECL FIL: <<LibName>>
+ ISO: 3 <<IsoH1>> <<IsoD2>> <<IsoO16>>
+ CALC WGT% D2O >>CwgtH1<< >>CwgtD2<< >>CwgtO16<< ;
+*----
+* Moderator density and WGT%
+*----
+REAL MwgtH1 MwgtD2 MwgtO16 ;
+INFO: :: EDIT 0
+ TMP: <<Mtemp>> K
+ PUR: <<Mprty>> WGT%
+ CALC DENS WATER >>Mdens<<
+ LIB: WIMSAECL FIL: <<LibName>>
+ ISO: 3 <<IsoH1>> <<IsoD2>> <<IsoO16>>
+ CALC WGT% D2O >>MwgtH1<< >>MwgtD2<< >>MwgtO16<< ;
+*----
+* Create library
+*----
+MicLib := LIB: ::
+ EDIT 1
+ NMIX 13 CTRA WIMS
+ DEPL LIB: WIMSAECL FIL: <<LibName>>
+ MIXS LIB: WIMSAECL FIL: <<LibName>>
+*----
+* Coolant mixture
+ MIX 1 <<Ctemp>> <<Cdens>>
+ H1 = <<IsoH1>> <<CwgtH1>>
+ O16 = <<IsoO16>> <<CwgtO16>>
+*----
+* Pressure tube mixture
+ MIX 2 <<PTtemp>> <<PTdens>>
+ <<IsoZr90PT>> 49.172123 <<IsoZr91PT>> 10.842638
+ <<IsoZr92PT>> 16.755391 <<IsoZr94PT>> 17.349846
+ <<IsoZr96PT>> 2.854732
+ <<IsoNb>> 2.65
+ <<IsoFe54>> 0.007339 <<IsoFe56>> 0.119472
+ <<IsoFe57>> 0.002808 <<IsoFe58>> 0.000380
+ <<IsoCr50>> 0.000417 <<IsoCr52>> 0.008370
+ <<IsoCr53>> 0.000967 <<IsoCr54>> 0.000245
+ <<IsoNi58>> 0.004368 <<IsoNi60>> 0.001740
+ <<IsoNi61>> 0.000077 <<IsoNi62>> 0.000249
+ <<IsoNi64>> 0.000066
+ <<IsoMn55>> 0.005
+ <<IsoC>> 0.008
+ <<IsoO16>> 0.135
+ <<IsoB10>> 0.0000294
+*----
+* Gap mixture
+ MIX 3 <<Gaptemp>> <<Gapdens>>
+ <<IsoC>> 27.11 <<IsoO16>> 72.89
+*----
+* Calandria tube mixture
+ MIX 4 <<CTtemp>> 6.4635
+ <<IsoZr90CT>> 49.663374 <<IsoZr91CT>> 10.950961
+ <<IsoZr92CT>> 16.922785 <<IsoZr94CT>> 17.523178
+ <<IsoZr96CT>> 2.883252
+ <<IsoFe54>> 0.011856 <<IsoFe56>> 0.192993
+ <<IsoFe57>> 0.004537 <<IsoFe58>> 0.000614
+ <<IsoNi58>> 0.004704 <<IsoNi60>> 0.001874
+ <<IsoNi61>> 0.000083 <<IsoNi62>> 0.000268
+ <<IsoNi64>> 0.000071
+ <<IsoCr50>> 0.004174 <<IsoCr52>> 0.083699
+ <<IsoCr53>> 0.009674 <<IsoCr54>> 0.002453
+ <<IsoC>> 0.04
+ <<IsoO16>> 0.125
+ <<IsoNb>> 0.01
+ <<IsoMn55>> 0.005
+ <<IsoB10>> 0.0000698
+*----
+* Moderator mixture
+ MIX 5 <<Mtemp>> <<Mdens>>
+ H1 = <<IsoH1>> <<MwgtH1>> D2 = <<IsoD2>> <<MwgtD2>>
+ O16 = <<IsoO16>> <<MwgtO16>>
+*----
+* Sheathing mixture
+ MIX 6 <<Shtemp>> <<Shdens>>
+ <<IsoZr90CL>> 49.663374 <<IsoZr91CL>> 10.950961
+ <<IsoZr92CL>> 16.922785 <<IsoZr94CL>> 17.523178
+ <<IsoZr96CL>> 2.883252
+ <<IsoFe54>> 0.011856 <<IsoFe56>> 0.192993
+ <<IsoFe57>> 0.004537 <<IsoFe58>> 0.000614
+ <<IsoNi58>> 0.004704 <<IsoNi60>> 0.001874
+ <<IsoNi61>> 0.000083 <<IsoNi62>> 0.000268
+ <<IsoNi64>> 0.000071
+ <<IsoCr50>> 0.004174 <<IsoCr52>> 0.083699
+ <<IsoCr53>> 0.009674 <<IsoCr54>> 0.002453
+ <<IsoC>> 0.04
+ <<IsoO16>> 0.125
+ <<IsoNb>> 0.01
+ <<IsoMn55>> 0.005
+ <<IsoB10>> 0.0000698
+*----
+* Fuel mixture in central (first) position
+ MIX 7 <<Ftemp1>> <<Fdens1>>
+ O16 = <<IsoO16>> 22.28412
+ DY160 = <<IsoDy160>> 0.32158
+ DY161 = <<IsoDy161>> 2.44767
+ DY162 = <<IsoDy162>> 3.32248
+ DY163 = <<IsoDy163>> 3.26310
+ DY164 = <<IsoDy164>> 3.71562
+ Gd152 = <<IsoGd152>> 0.02348
+ Gd154 = <<IsoGd154>> 0.25925
+ Gd155 = <<IsoGd155>> 1.77150
+ Gd156 = <<IsoGd156>> 2.46598
+ Gd157 = <<IsoGd157>> 1.89744
+ ! Gd158 = <<IsoGd158>> 5.73194
+ Gd158 = <<IsoGd158>> 3.03084
+ Gd160 = <<IsoGd160>> 2.70110
+ Zr90 = <<IsoZr90CL>> 26.65326
+ Zr91 = <<IsoZr91CL>> 5.87699
+ zr92 = <<IsoZr92CL>> 9.08178
+ Zr94 = <<IsoZr94CL>> 9.40399
+ Zr96 = <<IsoZr96CL>> 1.54733
+*----
+* Fuel mixture in second ring
+ MIX 8 <<Ftemp2>> <<Fdens2>>
+ U5 = <<IsoU5>> 2.0 1
+ U8 = <<IsoU8>> 98.0 1
+ O16 = <<IsoO16>> 13.44172
+ U6 = <<IsoU6>> 0.0 1
+ PU9 = <<IsoPu9>> 0.0 1
+ PU240 0.0 1 AM241 0.0 1 PU241 0.0 1
+ AM242M 0.0 1 PU242 0.0 1 AM243 0.0 1
+*----
+* Fuel mixture in third ring identical to second ring
+ MIX 9 COMB 8 1.0
+*----
+* Fuel mixture in last (fourth) ring identical to second ring
+ MIX 10 COMB 8 1.0
+*----
+* Coolant at half density
+ MIX 11 COMB 1 0.5 0 0.5
+*----
+* Fuel at 1.8 % enrichment
+ MIX 12 <<Ftemp2>> <<Fdens2>>
+ U5 = <<IsoU5>> 1.8 1
+ U8 = <<IsoU8>> 98.2 1
+ O16 = <<IsoO16>> 13.44138
+ U6 = <<IsoU6>> 0.0 1
+ PU9 = <<IsoPu9>> 0.0 1
+ PU240 0.0 1 AM241 0.0 1 PU241 0.0 1
+ AM242M 0.0 1 PU242 0.0 1 AM243 0.0 1
+*----
+* Fuel at 2.0 % enrichment
+ MIX 13 <<Ftemp2>> <<Fdens2>>
+ U5 = <<IsoU5>> 2.0 1
+ U8 = <<IsoU8>> 98.0 1
+ O16 = <<IsoO16>> 13.44172
+ U6 = <<IsoU6>> 0.0 1
+ PU9 = <<IsoPu9>> 0.0 1
+ PU240 0.0 1 AM241 0.0 1 PU241 0.0 1
+ AM242M 0.0 1 PU242 0.0 1 AM243 0.0 1
+ ;
+QUIT .
diff --git a/Dragon/data/C2D20_proc/MicACRE6UH.c2m b/Dragon/data/C2D20_proc/MicACRE6UH.c2m
new file mode 100644
index 0000000..0e89b81
--- /dev/null
+++ b/Dragon/data/C2D20_proc/MicACRE6UH.c2m
@@ -0,0 +1,87 @@
+*DECK MicACRE6UH.c2m
+*----
+* Name : MicACRE6UH.c2m
+* Type : DRAGON procedure
+* Use : To update the ENDF/B-VI Microlib for the ACR
+* by adding the mixtures required
+* for homogenisation by mixture.
+* Author : G. Marleau, E. Varin
+* Date : 2005/07/19 (Revision 1 by G. M.)
+* Reference : MicACRE6.c2m
+*
+* Procedure called as:
+*
+* MicLib := MicACRE6UH MicLib ;
+*
+* Parameters of the procedure
+* MicLib : Microlib tabulated for the required parameters stored in an XSM file
+* At input, the minimal set of mixtures included in this microlib is:
+* Mix 0 = void
+* Mix 1 = coolant at reference conditions
+* Mix 2 = pressure tube material
+* Mix 3 = gap material
+* Mix 4 = calendria tube material
+* Mix 5 = moderator
+* Mix 6 = fuel sheat
+* Mix 7 = central (first) Dy+Gd element
+* Mix 8 = second ring of fuel
+* Mix 9 = third ring of fuel
+* Mix 10 = last (fourth) ring of fuel
+* Mix 11 = coolant at 50 % density
+* Mix 12 = Natural uranium
+* Mix 13 = Dy2O3
+* The mixtures added to MicLib by the current procedure are:
+* Mix 14 = Full density coolant for central (first) fuel
+* Mix 15 = Full density coolant for second ring of fuel
+* Mix 16 = Full density coolant for third ring of fuel
+* Mix 17 = Full density coolant for last (fourth) ring of fuel
+* Mix 18 = Sheat for central (first) fuel element
+* Mix 19 = Sheat for second ring of fuel
+* Mix 20 = Sheat for third ring of fuel
+* Mix 21 = Sheat for last (fourth) ring of fuel
+* Mix 22 = Half density coolant for central (first) fuel
+* Mix 23 = Half density coolant for second ring of fuel
+* Mix 24 = Half density coolant for third ring of fuel
+* Mix 25 = Half density coolant for last (fourth) ring of fuel
+* Note: 1- MIX 14, 15, 16 and 17 are identical to MIX 1
+* 2- MIX 18, 19, 20 and 21 are identical to MIX 6
+* 3- MIX 22, 23, 24 and 25 are identical to MIX 11
+*----
+* Definition of output data structures
+*----
+PARAMETER MicLib :: ::: XSM_FILE MicLib ; ;
+*----
+* Local modules used in this procedure
+*----
+MODULE LIB: DELETE: END: ;
+*----
+* Library name
+*----
+STRING LibName := "E6MLIB" ;
+*----
+* Update library
+*----
+MicLib := LIB: MicLib ::
+ EDIT 0
+ NMIX 25 CTRA WIMS
+ MIXS LIB: WIMSAECL FIL: <<LibName>>
+*----
+* Full density coolant mixtures
+ MIX 14 COMB 1 1.0
+ MIX 15 COMB 1 1.0
+ MIX 16 COMB 1 1.0
+ MIX 17 COMB 1 1.0
+*----
+* Sheat mixtures
+ MIX 18 COMB 6 1.0
+ MIX 19 COMB 6 1.0
+ MIX 20 COMB 6 1.0
+ MIX 21 COMB 6 1.0
+*----
+* Half density coolant mixtures
+ MIX 22 COMB 11 1.0
+ MIX 23 COMB 11 1.0
+ MIX 24 COMB 11 1.0
+ MIX 25 COMB 11 1.0
+ ;
+QUIT .
diff --git a/Dragon/data/C2D20_proc/MicACRE6UW.c2m b/Dragon/data/C2D20_proc/MicACRE6UW.c2m
new file mode 100644
index 0000000..fcd2780
--- /dev/null
+++ b/Dragon/data/C2D20_proc/MicACRE6UW.c2m
@@ -0,0 +1,114 @@
+*DECK MicACRE6UD.c2m
+*----
+* Name : MicACRE6UW.c2m
+* Type : DRAGON procedure
+* Use : To update the ENDF/B-VI Microlib for the ACR
+* by adding the mixtures required
+* for distributed self shielding with pin splitting.
+* Author : G. Marleau
+* Date : 2005/07/19 (Revision 1 by G. M.)
+* Reference : MicACRE6.c2m
+* MicACRE6UH.c2m
+*
+* Procedure called as:
+*
+* MicLib := MicACRE6UD MicLib ;
+*
+* Parameters of the procedure
+* MicLib : Microlib tabulated for the required parameters stored in an XSM file
+* At input, the minimal set of mixtures included in this microlib is:
+* Mix 0 = void
+* Mix 1 = coolant at reference conditions
+* Mix 2 = pressure tube material
+* Mix 3 = gap material
+* Mix 4 = calendria tube material
+* Mix 5 = moderator
+* Mix 6 = fuel sheat
+* Mix 7 = central (first) fuel element
+* Mix 8 = second ring of fuel
+* Mix 9 = third ring of fuel
+* Mix 10 = last (fourth) ring of fuel
+* Mix 11 = coolant at 50 % density
+* Mix 12 = Natural uranium
+* Mix 13 = Dy2O3
+* The mixtures 14-25 are added using MicACRE6UH.c2m
+* The mixtures added to MicLib by the current procedure are:
+* Mix 26 = 1-4 annulus for central (first) fuel
+* Mix 27 = 5-8 annulus for central (first) fuel
+* Mix 28 = 8-12 annulus for central (first) fuel
+* Mix 29 = first outward annulus for second ring of fuel
+* Mix 30 = first inward annulus for second ring of fuel
+* Mix 31 = second outward annulus for second ring of fuel
+* Mix 32 = second inward annulus for second ring of fuel
+* Mix 33 = third outward annulus for second ring of fuel
+* Mix 34 = third inward annulus for second ring of fuel
+* Mix 35 = first outward annulus for third ring of fuel
+* Mix 36 = first inward annulus for third ring of fuel
+* Mix 37 = second outward annulus for third ring of fuel
+* Mix 38 = second inward annulus for third ring of fuel
+* Mix 39 = third outward annulus for third ring of fuel
+* Mix 40 = third inward annulus for third ring of fuel
+* Mix 41 = first outward annulus for fourth ring of fuel
+* Mix 42 = first inward annulus for fourth ring of fuel
+* Mix 43 = second outward annulus for fourth ring of fuel
+* Mix 44 = second inward annulus for fourth ring of fuel
+* Mix 45 = third outward annulus for fourth ring of fuel
+* Mix 46 = third inward annulus for fourth ring of fuel
+* Mix 47 = fourth outward annulus for fourth ring of fuel
+* Mix 48 = fourth inward annulus for fourth ring of fuel
+* Note: 1- MIX 26, 27, 28 are identical to MIX 7
+* 2- MIX 29-34 are identical to MIX 8
+* 3- MIX 35-40 are identical to MIX 9
+* 4- MIX 41-48 are identical to MIX 10
+*----
+* Definition of output data structures
+*----
+PARAMETER MicLib :: ::: XSM_FILE MicLib ; ;
+*----
+* Local modules used in this procedure
+*----
+MODULE LIB: DELETE: END: ;
+*----
+* Library name
+*----
+STRING LibName := "E6MLIB" ;
+*----
+* Update library
+*----
+MicLib := LIB: MicLib ::
+ EDIT 1
+ NMIX 48 CTRA WIMS
+ MIXS LIB: WIMSAECL FIL: <<LibName>>
+*----
+* Central pin
+ MIX 26 COMB 7 1.0
+ MIX 27 COMB 7 1.0
+ MIX 28 COMB 7 1.0
+*----
+* Second ring of pins
+ MIX 29 COMB 8 1.0
+ MIX 30 COMB 8 1.0
+ MIX 31 COMB 8 1.0
+ MIX 32 COMB 8 1.0
+ MIX 33 COMB 8 1.0
+ MIX 34 COMB 8 1.0
+*----
+* Third ring of pins
+ MIX 35 COMB 9 1.0
+ MIX 36 COMB 9 1.0
+ MIX 37 COMB 9 1.0
+ MIX 38 COMB 9 1.0
+ MIX 39 COMB 9 1.0
+ MIX 40 COMB 9 1.0
+*----
+* Fourth ring of pins
+ MIX 41 COMB 10 1.0
+ MIX 42 COMB 10 1.0
+ MIX 43 COMB 10 1.0
+ MIX 44 COMB 10 1.0
+ MIX 45 COMB 10 1.0
+ MIX 46 COMB 10 1.0
+ MIX 47 COMB 10 1.0
+ MIX 48 COMB 10 1.0
+ ;
+QUIT .
generated by cgit v1.2.3 (git 2.46.0) at 2026-02-16 15:23:44 +0000