Initial commit from Polytechnique Montreal

9 files changed, 3158 insertions, 0 deletions

diff --git a/Dragon/data/proc_LZC/DimCell.c2m b/Dragon/data/proc_LZC/DimCell.c2m
new file mode 100644
index 0000000..e738481
--- /dev/null
+++ b/Dragon/data/proc_LZC/DimCell.c2m
@@ -0,0 +1,29 @@
+*DECK DimCell.c2m
+*----
+*  Name          : DimCell.c2m
+*  Type          : DRAGON procedure
+*  Use           : Dimensions for Gentilly-2 cell
+*  Author        : G. Marleau
+*  Date          : 2006/03/01
+*  Reference     : Gentilly-2 .
+*  
+*  Procedure called as: 
+*
+*  DimCell :: >>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 rcool     rpt     rgap     rct      LPitch     BLength   := 
+     5.168878  5.60320 6.44780  6.587482 28.575     49.53     ;
+:: <<rcool>> <<rpt>> <<rgap>> <<rct>>  <<LPitch>> <<BLength>> ;
+QUIT .
diff --git a/Dragon/data/proc_LZC/DimFuel.c2m b/Dragon/data/proc_LZC/DimFuel.c2m
new file mode 100644
index 0000000..d79c7a3
--- /dev/null
+++ b/Dragon/data/proc_LZC/DimFuel.c2m
@@ -0,0 +1,49 @@
+*DECK DimFuel.c2m
+*----
+*  Name          : DimFuel.c2m
+*  Type          : DRAGON procedure
+*  Use           : Fuel dimensions for Gentilly-2 cell
+*  Author        : G. Marleau
+*  Date          : 2006/03/01
+*  Reference     : Gentilly-2 .
+*  
+*  Description of the procedure: 
+*    DimFuel ::  <<Ring>> 
+*    >>npinp<< >>rpinp<< >>apinp<< >>rradf<< >>rrads<< ;
+*  Input data:
+*    Ring          : Ring number
+*  Output data:
+*    npinp         : number of pins in this ring
+*    rpinp         : ring radius 
+*    apinp         : first pin angle
+*    rradf         : outer fuel radius
+*    rrads         : outer sheat radius
+*----
+* Definition of the input and output data structures and 
+* of the input and output parameters for this procedure
+*
+MODULE    END:   ;
+INTEGER   Ring   ;
+  :: >>Ring<<  ;
+INTEGER      npinp1 npinp2 npinp3 npinp4 :=
+             1      6      12     18     ;
+REAL         rradf1   rrads1   rpinp1   apinp1 := 
+             0.6122   0.6540    0.0     0.0      ;    
+REAL         rradf2   rrads2   rpinp2   apinp2 :=
+             rradf1   rrads1   1.4885   0.0      ;
+REAL         rradf3   rrads3   rpinp3   apinp3 :=
+             rradf1   rrads1   2.8755   $Pi_R  npinp3 I_TO_R /  ;
+REAL         rradf4   rrads4   rpinp4   apinp4 :=
+             rradf1   rrads1   4.3305   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 DimFuel : " Ring ;
+ENDIF ;
+QUIT .
diff --git a/Dragon/data/proc_LZC/DimLZC.c2m b/Dragon/data/proc_LZC/DimLZC.c2m
new file mode 100644
index 0000000..af18d3d
--- /dev/null
+++ b/Dragon/data/proc_LZC/DimLZC.c2m
@@ -0,0 +1,52 @@
+*DECK DimLZC.c2m
+*----
+*  Name          : DimLZC.c2m
+*  Type          : DRAGON procedure
+*  Use           : Dimensions for Gentilly-2 Liquid Zone Controller
+*                  (cluster model)
+*  Author        : G. Marleau
+*  Date          : 2006/05/05
+  
+*  Description of the procedure: 
+*    DimLZC ::    
+*      >>RinnGT<<   >>RoutGT<< 
+*      >>RinnBa<< >>RoutBa<<  
+*      >>RinnFe<< >>RoutFe<<  
+*      >>RinnBu<< >>RoutBu<<  
+*      >>RinnSc<< >>RoutSc<<  
+*      >>Rcent<< ;
+*  Output data:
+*    RinnGT        : inner radius of guide tube
+*    RoutGT        : outer radius of guide tube
+*    RinnBa        : inner Balance tube radius
+*    RoutBa        : outer Balance tube radius
+*    RinnFe        : inner Feeder tube radius
+*    RoutFe        : outer Feeder tube radius
+*    RinnBu        : inner Bubbler tube radius
+*    RoutBu        : outer Bubbler tube radius
+*    RinnSc        : inner Scavanger tube radius
+*    RoutSc        : outer Scavanger tube radius
+*    Rcent         : annulus radius where the LZC are located
+*----
+* Definition of the input and output data structures and 
+* of the input and output parameters for this procedure
+*----
+REAL    RinnGT       RoutGT := 
+        11.430 2.0 / 11.684 2.0 / ;
+REAL    RinnBa       RoutBa  :=
+        1.814 2.0 /  1.966 2.0 / ;
+REAL    RinnFe       RoutFe  :=
+        2.972 2.0 /  3.124 2.0 / ;
+REAL    RinnBu       RoutBu  :=
+        0.483 2.0 /  0.635 2.0 / ;
+REAL    RinnSc       RoutSc  :=
+        1.608 2.0 /  1.760 2.0 / ;
+REAL    Rcent                :=
+        RinnGT 2.0 + 2.0 /   ;
+::    <<RinnGT>> <<RoutGT>> 
+      <<RinnBa>> <<RoutBa>>  
+      <<RinnFe>> <<RoutFe>>  
+      <<RinnBu>> <<RoutBu>>  
+      <<RinnSc>> <<RoutSc>>  
+      <<Rcent>> ;
+QUIT .
diff --git a/Dragon/data/proc_LZC/GeoCell.c2m b/Dragon/data/proc_LZC/GeoCell.c2m
new file mode 100755
index 0000000..97f14f9
--- /dev/null
+++ b/Dragon/data/proc_LZC/GeoCell.c2m
@@ -0,0 +1,306 @@
+*DECK GeoLZC.c2m
+*----------
+*  Name          : GeoLZC.c2m
+*  Mesh          : DRAGON procedure
+*  Use           : LZC geometry 
+*                  for G2 (cluster model)
+*  Author        : G. Marleau
+*  Date          : 2006/05/05
+*
+*  Procedure description: 
+*    Geometry := GeoLZC :: 
+*      <<PrtLev>> <<Type>> <<Mesh>> 
+*      <<NbAngles>> <<Density>> ;
+*
+*  Procedure parameters:
+*    Geometry      : XSM_FILE containing the LZC geometry 
+*                    for G2
+*    PrtLev        : Print Level
+*                    0  -> Minimum print
+*                    1  -> Standard print
+*                    10 -> Debug
+*    Type          : Geometry type
+*                    'XY' -> 2-D projection of geometry in XY plane
+*                    'XZ' -> 2-D projection of geometry in XZ plane
+*                    '3D' -> 3-D geometry (default value if Type 
+*                            is not known)
+*    Mesh          : Geometry Mesh
+*                    Coarse -> Corase mesh 3-D geometry
+*                    Fine   -> Fine mesh 3-D geometry (default) 
+*    NbAngles      : Suggested angular quadrature for tracking
+*    Density       : Suggested line density for tracking
+* 
+*----------
+* Define procedure parameters and options and read
+* the procedure input data
+*----
+PARAMETER Geometry :: EDIT 0 ::: XSM_FILE Geometry ; ;
+STRING    Type Mesh                ;
+INTEGER   PrtLev NbAngles          ;
+REAL      Density                  ;
+::  >>PrtLev<< >>Type<< >>Mesh<<  ;
+*----
+*  Local procedures and modules used in this procedure
+*----     
+PROCEDURE    DimFuel DimCell DimLZC      ;
+MODULE       GEO: END: ;
+*----
+*  Get overall cell dimensions using procedure DimCell
+*----
+REAL         rcool4 rpt rgap rct LPitch BLength ;
+DimCell  ::  
+  >>rcool4<< >>rpt<< >>rgap<< >>rct<< >>LPitch<< >>BLength<< ;
+IF PrtLev  9 > THEN
+  ECHO "Overall cell dimensions" ;
+  ECHO rcool4 rpt rgap rct LPitch BLength ;
+ENDIF ;
+*----
+*  Get ring by ring description of fuel using procedure DimFuel
+*----
+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  ;
+EVALUATE    Ring := 1 ;
+DimFuel  :: <<Ring>>
+  >>npinp1<< >>rpinp1<< >>apinp1<< >>rradf1<< >>rrads1<< ;
+EVALUATE    Ring := 2 ;
+DimFuel  :: <<Ring>>
+  >>npinp2<< >>rpinp2<< >>apinp2<< >>rradf2<< >>rrads2<< ;
+EVALUATE    Ring := 3 ;
+DimFuel  :: <<Ring>>
+  >>npinp3<< >>rpinp3<< >>apinp3<< >>rradf3<< >>rrads3<< ;
+EVALUATE    Ring := 4 ;
+DimFuel  :: <<Ring>> 
+  >>npinp4<< >>rpinp4<< >>apinp4<< >>rradf4<< >>rrads4<< ;
+IF PrtLev  9 > THEN
+  ECHO "Fuel rings 1 to 4" ;
+  ECHO  npinp1   npinp2   npinp3   npinp4  ;
+  ECHO  rradf1   rrads1   rpinp1   apinp1  ;    
+  ECHO  rradf2   rrads2   rpinp2   apinp2  ;
+  ECHO  rradf3   rrads3   rpinp3   apinp3  ;
+  ECHO  rradf4   rrads4   rpinp4   apinp4  ;
+ENDIF ;
+*----
+*  Get LZC using procedure DimLZC
+*----
+REAL    RinnGT       RoutGT 
+        RinnBa       RoutBa  
+        RinnFe       RoutFe  
+        RinnBu       RoutBu  
+        RinnSc       RoutSc 
+         Rcent              ;
+DimLZC ::    
+  >>RinnGT<<   >>RoutGT<< 
+  >>RinnBa<< >>RoutBa<<  
+  >>RinnFe<< >>RoutFe<<  
+  >>RinnBu<< >>RoutBu<< 
+  >>RinnSc<< >>RoutSc<<  
+  >>Rcent<< ;
+IF PrtLev  9 > THEN
+  ECHO "LZC dimensions" ;
+  ECHO "Guide tube      " RinnGT RoutGT ;
+  ECHO "Feeder tubes    " RinnBa RoutBa  RinnFe RoutFe ;
+  ECHO "Scavenger tubes " RinnBu RoutBu  RinnSc RoutSc ;
+  ECHO "LZC location    " Rcent ;
+ENDIF ;
+*----
+*  Evaluate derived dimensions for supercell description
+*----
+REAL HLPitch   := LPitch 2.0 / ;
+REAL HBLength  := BLength 2.0 / ;
+REAL LimCoCell := HLPitch rct - RoutGT + 2.0 / ;
+REAL IntFCell  := LPitch LimCoCell - ;
+REAL IntYCell  := HLPitch LimCoCell - ;
+REAL OffsetXL OffsetXR := LimCoCell 2.0 / LimCoCell -2.0 / ;
+* Left and right fuel cell X position
+REAL XLFmin        XLFint          XLFcent        XLFmax 
+     XRFmin        XRFint          XRFcent        XRFmax :=
+     LPitch -1.0 * IntFCell -1.0 * HLPitch -1.0 * LimCoCell -1.0 *
+     LimCoCell     IntFCell        HLPitch        LPitch  ;
+* Cell Y position
+REAL YFmin   YFintb  YFcent YFintt  YFmax   :=
+     XLFcent IntYCell -1.0 *  0.0    IntYCell XRFcent ;
+* Z position
+REAL ZFmin           ZFintf  ZFcent ZFintb  ZFmax   :=
+     HBLength -1.0 * XLFmax  0.0    XRFmin  HBLength ;
+IF PrtLev  9 > THEN
+  ECHO "Derived dimensions for supercell" ;
+  ECHO  HLPitch HBLength  LimCoCell IntFCell  OffsetXL OffsetXR ;
+  ECHO  XLFmin XLFint XLFcent XLFmax ;
+  ECHO  XRFmin XRFint XRFcent XRFmax ;    
+  ECHO  YFmin  YFintb YFcent  YFintt  YFmax ;
+  ECHO  ZFmin         ZFintf  ZFcent ZFintb  ZFmax   ;
+ENDIF ;
+*----
+*  Derived dimensions for fuel 
+*----
+REAL         eps1  := 0.1   ;
+REAL         rcool1   rcool2   rcool3          :=
+             rrads1   eps1  + rpinp2 rrads2 - + 2.0 /
+             rpinp2 rrads2  + rpinp3 rrads3 - + 2.0 /
+             rpinp3 rrads3  + rpinp4 rrads4 - + 2.0 / ;
+REAL         rmod  := HLPitch LimCoCell - eps1 - ;
+IF PrtLev  9 > THEN
+  ECHO "Derived dimensions for fuel" ;
+  ECHO  eps1 rcool1 rcool2 rcool3 rmod ;
+ENDIF ;
+*----
+*  Define mixtures
+*----
+INTEGER  mcL1 mpt mgap mct mmod :=   1  2  3  4  5 ;
+INTEGER  mfL1 mfL2 mfL3 mfL4    :=   6  7  8  9  ;
+INTEGER  mfR1 mfR2 mfR3 mfR4    :=  10 11 12 13 ;
+INTEGER  ms                     :=  14          ;
+INTEGER  mfill mgt              :=  15 16       ;
+INTEGER  ms11      ms12      ms13      ms14    :=  
+         17        18        19        20  ;
+INTEGER  ms21      ms22      ms23      ms24    :=  
+         ms11 4 +  ms12 4 +  ms13 4 +  ms14 4 + ;
+INTEGER  ms31      ms32      ms33      ms34    :=  
+         ms11 8 +  ms12 8 +  ms13 8 +  ms14 8 + ;
+INTEGER  mf11      mf12      mf13      mf14    :=  
+         ms11 12 + ms12 12 + ms13 12 + ms14 12 + ;
+INTEGER  mf21      mf22      mf23      mf24    :=  
+         ms11 16 + ms12 16 + ms13 16 + ms14 16 + ;
+REAL     PIO2 PIO6 := $Pi_R 2.0 / $Pi_R 6.0 /    ;
+REAL     ApinZs1        ApinZs2        ApinZs3 
+         ApinZf1        ApinZf2                    :=
+         0.0            PIO6 4.0 *     PIO6 8.0 *   
+         PIO6 2.0 *     PIO6 6.0 *                 ;
+REAL     ApinXs1        ApinXs2        ApinXs3 
+         ApinXf1        ApinXf2                    :=
+         PIO2 ApinZs1 - PIO2 ApinZs2 - PIO2 ApinZs3 -
+         PIO2 ApinZf1 - PIO2 ApinZf2 -             ;
+ECHO     ApinZs1 ApinZs2 ApinZs3 ApinZf1 ApinZf2 ;
+IF PrtLev  9 > THEN
+  ECHO "Mixtures for coolant, PT, gap, CT and moderator" ;
+  ECHO  mcL1 mpt mgap mct mmod ;
+  ECHO "Mixtures for left fuel in each ring and for sheathing" ;
+  ECHO  mfL1 mfL2 mfL3 mfL4 ms ;
+  ECHO "Mixtures for right fuel in each ring and for sheathing" ;
+  ECHO  mfR1 mfR2 mfR3 mfR4 ms ;
+  ECHO "Mixtures for scavenger tube 1" ;
+  ECHO  ms11 ms12 ms13 ms14 ;
+  ECHO "Mixtures for scavenger tube 2" ;
+  ECHO  ms21 ms22 ms23 ms24 ;
+  ECHO "Mixtures for scavenger tube 3" ;
+  ECHO  ms31 ms32 ms33 ms34 ;
+  ECHO "Mixtures for feeder tube 1" ;
+  ECHO  mf11 mf12 mf13 mf14 ;
+  ECHO "Mixtures for feeder tube 12" ;
+  ECHO  mf11 mf22 mf23 mf24 ;
+ENDIF ;
+*----
+*  Geometry
+*----
+  IF Mesh "Coarse" = THEN
+*----
+*  Coarse mesh geometry
+*----
+    Geometry := GEO: :: CAR2D 1 1 
+      X- REFL X+ SYME Y- REFL Y+ SYME 
+      CELL  FuelG
+*----
+*  Fuel on the left
+*----
+    ::: FuelG := GEO: CARCEL 7 1 1 
+      MESHX   <<YFmin>>  <<YFmax>>
+      MESHY   <<YFmin>>  <<YFmax>>
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> 
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBE 2 
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        RADIUS 0.00000 <<rradf1>> <<rrads1>>
+        MIX  <<mfL1>> <<ms>> ;
+      ::: ROD2 := GEO: TUBE 2 
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        RADIUS 0.00000 <<rradf2>> <<rrads2>> 
+        MIX  <<mfL2>> <<ms>>  ;
+      ::: ROD3 := GEO: TUBE 2 
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        RADIUS 0.00000 <<rradf3>> <<rrads3>>  
+        MIX  <<mfL3>> <<ms>>  ;
+      ::: ROD4 := GEO: TUBE 2 
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        RADIUS 0.00000 <<rradf4>> <<rrads4>> 
+        MIX  <<mfL4>> <<ms>>  ;
+    ;
+    ;
+    EVALUATE NbAngles Density := 15 30.0 ;
+  ELSE
+*----
+*  Fine mesh geometry (default)
+*----
+    Geometry := GEO: :: CAR2D 1 1 
+      X- REFL X+ REFL Y- REFL Y+ REFL 
+      CELL  FuelG
+*----
+*  Fuel on the left 
+*----
+    ::: FuelG := GEO: CARCEL 8 4 4
+      MESHX   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+      MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    <<rmod>>
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBE 2 
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        RADIUS 0.00000 <<rradf1>> <<rrads1>> SPLITR -3 1 
+        MIX  <<mfL1>> <<ms>> ;
+      ::: ROD2 := GEO: TUBE 2 
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        RADIUS 0.00000 <<rradf2>> <<rrads2>> SPLITR -3 1 
+        MIX  <<mfL2>> <<ms>> ;
+      ::: ROD3 := GEO: TUBE 2 
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        RADIUS 0.00000 <<rradf3>> <<rrads3>> SPLITR -3 1 
+        MIX  <<mfL3>> <<ms>> ;
+      ::: ROD4 := GEO: TUBE 2 
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        RADIUS 0.00000 <<rradf4>> <<rrads4>> SPLITR -3 1 
+        MIX  <<mfL4>> <<ms>> ;
+    ;
+    ;
+    EVALUATE NbAngles Density := 15 30.0 ;
+  ENDIF ;
+::  <<NbAngles>> <<Density>>          ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/proc_LZC/GeoLZC.c2m b/Dragon/data/proc_LZC/GeoLZC.c2m
new file mode 100755
index 0000000..45b89fe
--- /dev/null
+++ b/Dragon/data/proc_LZC/GeoLZC.c2m
@@ -0,0 +1,1020 @@
+*DECK GeoLZC.c2m
+*----------
+*  Name          : GeoLZC.c2m
+*  Mesh          : DRAGON procedure
+*  Use           : LZC geometry 
+*                  for G2 (cluster model)
+*  Author        : G. Marleau
+*  Date          : 2006/05/05
+*
+*  Procedure description: 
+*    Geometry := GeoLZC :: 
+*      <<PrtLev>> <<Type>> <<Mesh>> 
+*      <<NbAngles>> <<Density>> ;
+*
+*  Procedure parameters:
+*    Geometry      : XSM_FILE containing the LZC geometry 
+*                    for G2
+*    PrtLev        : Print Level
+*                    0  -> Minimum print
+*                    1  -> Standard print
+*                    10 -> Debug
+*    Type          : Geometry type
+*                    'XY' -> 2-D projection of geometry in XY plane
+*                    'XZ' -> 2-D projection of geometry in XZ plane
+*                    '3D' -> 3-D geometry (default value if Type 
+*                            is not known)
+*    Mesh          : Geometry Mesh
+*                    Coarse -> Corase mesh 3-D geometry
+*                    Fine   -> Fine mesh 3-D geometry (default) 
+*    NbAngles      : Suggested angular quadrature for tracking
+*    Density       : Suggested line density for tracking
+* 
+*----------
+* Define procedure parameters and options and read
+* the procedure input data
+*----
+PARAMETER Geometry :: EDIT 0 ::: XSM_FILE Geometry ; ;
+STRING    Type Mesh                ;
+INTEGER   PrtLev NbAngles          ;
+REAL      Density                  ;
+::  >>PrtLev<< >>Type<< >>Mesh<<  ;
+*----
+*  Local procedures and modules used in this procedure
+*----     
+PROCEDURE    DimFuel DimCell DimLZC      ;
+MODULE       GEO: END: ;
+*----
+*  Get overall cell dimensions using procedure DimCell
+*----
+REAL         rcool4 rpt rgap rct LPitch BLength ;
+DimCell  ::  
+  >>rcool4<< >>rpt<< >>rgap<< >>rct<< >>LPitch<< >>BLength<< ;
+IF PrtLev  9 > THEN
+  ECHO "Overall cell dimensions" ;
+  ECHO rcool4 rpt rgap rct LPitch BLength ;
+ENDIF ;
+*----
+*  Get ring by ring description of fuel using procedure DimFuel
+*----
+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  ;
+EVALUATE    Ring := 1 ;
+DimFuel  :: <<Ring>>
+  >>npinp1<< >>rpinp1<< >>apinp1<< >>rradf1<< >>rrads1<< ;
+EVALUATE    Ring := 2 ;
+DimFuel  :: <<Ring>>
+  >>npinp2<< >>rpinp2<< >>apinp2<< >>rradf2<< >>rrads2<< ;
+EVALUATE    Ring := 3 ;
+DimFuel  :: <<Ring>>
+  >>npinp3<< >>rpinp3<< >>apinp3<< >>rradf3<< >>rrads3<< ;
+EVALUATE    Ring := 4 ;
+DimFuel  :: <<Ring>> 
+  >>npinp4<< >>rpinp4<< >>apinp4<< >>rradf4<< >>rrads4<< ;
+IF PrtLev  9 > THEN
+  ECHO "Fuel rings 1 to 4" ;
+  ECHO  npinp1   npinp2   npinp3   npinp4  ;
+  ECHO  rradf1   rrads1   rpinp1   apinp1  ;    
+  ECHO  rradf2   rrads2   rpinp2   apinp2  ;
+  ECHO  rradf3   rrads3   rpinp3   apinp3  ;
+  ECHO  rradf4   rrads4   rpinp4   apinp4  ;
+ENDIF ;
+*----
+*  Get LZC using procedure DimLZC
+*----
+REAL    RinnGT       RoutGT 
+        RinnBa       RoutBa  
+        RinnFe       RoutFe  
+        RinnBu       RoutBu  
+        RinnSc       RoutSc 
+         Rcent              ;
+DimLZC ::    
+  >>RinnGT<<   >>RoutGT<< 
+  >>RinnBa<< >>RoutBa<<  
+  >>RinnFe<< >>RoutFe<<  
+  >>RinnBu<< >>RoutBu<< 
+  >>RinnSc<< >>RoutSc<<  
+  >>Rcent<< ;
+IF PrtLev  9 > THEN
+  ECHO "LZC dimensions" ;
+  ECHO "Guide tube      " RinnGT RoutGT ;
+  ECHO "Feeder tubes    " RinnBa RoutBa  RinnFe RoutFe ;
+  ECHO "Scavenger tubes " RinnBu RoutBu  RinnSc RoutSc ;
+  ECHO "LZC location    " Rcent ;
+ENDIF ;
+*----
+*  Evaluate derived dimensions for supercell description
+*----
+REAL HLPitch   := LPitch 2.0 / ;
+REAL HBLength  := BLength 2.0 / ;
+REAL LimCoCell := HLPitch rct - RoutGT + 2.0 / ;
+REAL IntFCell  := LPitch LimCoCell - ;
+REAL IntYCell  := HLPitch LimCoCell - ;
+REAL OffsetXL OffsetXR := LimCoCell 2.0 / LimCoCell -2.0 / ;
+* Left and right fuel cell X position
+REAL XLFmin        XLFint          XLFcent        XLFmax 
+     XRFmin        XRFint          XRFcent        XRFmax :=
+     LPitch -1.0 * IntFCell -1.0 * HLPitch -1.0 * LimCoCell -1.0 *
+     LimCoCell     IntFCell        HLPitch        LPitch  ;
+* Cell Y position
+REAL YFmin   YFintb  YFcent YFintt  YFmax   :=
+     XLFcent IntYCell -1.0 *  0.0    IntYCell XRFcent ;
+* Z position
+REAL ZFmin           ZFintf  ZFcent ZFintb  ZFmax   :=
+     HBLength -1.0 * XLFmax  0.0    XRFmin  HBLength ;
+IF PrtLev  9 > THEN
+  ECHO "Derived dimensions for supercell" ;
+  ECHO  HLPitch HBLength  LimCoCell IntFCell  OffsetXL OffsetXR ;
+  ECHO  XLFmin XLFint XLFcent XLFmax ;
+  ECHO  XRFmin XRFint XRFcent XRFmax ;    
+  ECHO  YFmin  YFintb YFcent  YFintt  YFmax ;
+  ECHO  ZFmin         ZFintf  ZFcent ZFintb  ZFmax   ;
+ENDIF ;
+*----
+*  Derived dimensions for fuel 
+*----
+REAL         eps1  := 0.1   ;
+REAL         rcool1   rcool2   rcool3          :=
+             rrads1   eps1  + rpinp2 rrads2 - + 2.0 /
+             rpinp2 rrads2  + rpinp3 rrads3 - + 2.0 /
+             rpinp3 rrads3  + rpinp4 rrads4 - + 2.0 / ;
+REAL         rmod  := HLPitch LimCoCell - eps1 - ;
+IF PrtLev  9 > THEN
+  ECHO "Derived dimensions for fuel" ;
+  ECHO  eps1 rcool1 rcool2 rcool3 rmod ;
+ENDIF ;
+*----
+*  Define mixtures
+*----
+INTEGER  mcL1 mpt mgap mct mmod :=   1  2  3  4  5 ;
+INTEGER  mfL1 mfL2 mfL3 mfL4    :=   6  7  8  9  ;
+INTEGER  mfR1 mfR2 mfR3 mfR4    :=  10 11 12 13 ;
+INTEGER  ms                     :=  14          ;
+INTEGER  mfill mgt              :=  15 16       ;
+INTEGER  ms11      ms12      ms13      ms14    :=  
+         17        18        19        20  ;
+INTEGER  ms21      ms22      ms23      ms24    :=  
+         ms11 4 +  ms12 4 +  ms13 4 +  ms14 4 + ;
+INTEGER  ms31      ms32      ms33      ms34    :=  
+         ms11 8 +  ms12 8 +  ms13 8 +  ms14 8 + ;
+INTEGER  mf11      mf12      mf13      mf14    :=  
+         ms11 12 + ms12 12 + ms13 12 + ms14 12 + ;
+INTEGER  mf21      mf22      mf23      mf24    :=  
+         ms11 16 + ms12 16 + ms13 16 + ms14 16 + ;
+REAL     PIO2 PIO6 := $Pi_R 2.0 / $Pi_R 6.0 /    ;
+REAL     ApinZs1        ApinZs2        ApinZs3 
+         ApinZf1        ApinZf2                    :=
+         0.0            PIO6 4.0 *     PIO6 8.0 *   
+         PIO6 2.0 *     PIO6 6.0 *                 ;
+REAL     ApinXs1        ApinXs2        ApinXs3 
+         ApinXf1        ApinXf2                    :=
+         PIO2 ApinZs1 - PIO2 ApinZs2 - PIO2 ApinZs3 -
+         PIO2 ApinZf1 - PIO2 ApinZf2 -             ;
+ECHO     ApinZs1 ApinZs2 ApinZs3 ApinZf1 ApinZf2 ;
+IF PrtLev  9 > THEN
+  ECHO "Mixtures for coolant, PT, gap, CT and moderator" ;
+  ECHO  mcL1 mpt mgap mct mmod ;
+  ECHO "Mixtures for left fuel in each ring and for sheathing" ;
+  ECHO  mfL1 mfL2 mfL3 mfL4 ms ;
+  ECHO "Mixtures for right fuel in each ring and for sheathing" ;
+  ECHO  mfR1 mfR2 mfR3 mfR4 ms ;
+  ECHO "Mixtures for scavenger tube 1" ;
+  ECHO  ms11 ms12 ms13 ms14 ;
+  ECHO "Mixtures for scavenger tube 2" ;
+  ECHO  ms21 ms22 ms23 ms24 ;
+  ECHO "Mixtures for scavenger tube 3" ;
+  ECHO  ms31 ms32 ms33 ms34 ;
+  ECHO "Mixtures for feeder tube 1" ;
+  ECHO  mf11 mf12 mf13 mf14 ;
+  ECHO "Mixtures for feeder tube 12" ;
+  ECHO  mf11 mf22 mf23 mf24 ;
+ENDIF ;
+*----
+*  Geometry
+*----
+IF Type "XY" = THEN
+  IF Mesh "Coarse" = THEN
+*----
+*  Coarse mesh geometry
+*----
+    Geometry := GEO: :: CAR2D 3 1 
+      X- REFL X+ REFL Y- REFL Y+ SYME 
+      CELL  FuelL LZC FuelR
+*----
+*  Fuel on the left
+*----
+    ::: FuelL := GEO: CARCEL 7 1 1 
+      MESHX   <<XLFmin>> <<XLFmax>>
+      MESHY   <<YFmin>>  <<YFmax>>
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    
+      OFFCENTER <<OffsetXL>> 0.0       0.0
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> 
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBE 2 
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        RADIUS 0.00000 <<rradf1>> <<rrads1>>
+        MIX  <<mfL1>> <<ms>> ;
+      ::: ROD2 := GEO: TUBE 2 
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        RADIUS 0.00000 <<rradf2>> <<rrads2>> 
+        MIX  <<mfL2>> <<ms>>  ;
+      ::: ROD3 := GEO: TUBE 2 
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        RADIUS 0.00000 <<rradf3>> <<rrads3>>  
+        MIX  <<mfL3>> <<ms>>  ;
+      ::: ROD4 := GEO: TUBE 2 
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        RADIUS 0.00000 <<rradf4>> <<rrads4>> 
+        MIX  <<mfL4>> <<ms>>  ;
+    ;
+*----
+*  LZC 
+*----
+    ::: LZC := GEO: CAR2D 1 1 
+      MESHX   <<XLFmax>> <<XRFmin>>
+      MESHY   <<YFmin>>  <<YFmax>>
+      MIX     <<mmod>>
+    ;
+*----
+*  Fuel on the right
+*----
+    ::: FuelR := GEO: CARCEL 7 1 1 
+      MESHX   <<XRFmin>> <<XRFmax>>
+      MESHY   <<YFmin>>  <<YFmax>>
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    
+      OFFCENTER <<OffsetXR>> 0.0       0.0
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> 
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBE 2 
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        RADIUS 0.00000 <<rradf1>> <<rrads1>>
+        MIX  <<mfR1>> <<ms>> ;
+      ::: ROD2 := GEO: TUBE 2 
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        RADIUS 0.00000 <<rradf2>> <<rrads2>> 
+        MIX  <<mfR2>> <<ms>>  ;
+      ::: ROD3 := GEO: TUBE 2 
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        RADIUS 0.00000 <<rradf3>> <<rrads3>>  
+        MIX  <<mfR3>> <<ms>>  ;
+      ::: ROD4 := GEO: TUBE 2 
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        RADIUS 0.00000 <<rradf4>> <<rrads4>> 
+        MIX  <<mfR4>> <<ms>>  ;
+    ;
+    ;
+    EVALUATE NbAngles Density := 15 30.0 ;
+  ELSE
+*----
+*  Fine mesh geometry (default)
+*----
+    Geometry := GEO: :: CAR2D 3 1 
+      X- REFL X+ REFL Y- REFL Y+ SYME 
+      CELL  FuelL LZC FuelR
+*----
+*  Fuel on the left 
+*----
+    ::: FuelL := GEO: CARCEL 8 3 4
+      MESHX   <<XLFmin>> <<XLFint>> <<XLFcent>> <<XLFmax>>
+      MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    <<rmod>>
+      OFFCENTER <<OffsetXL>> 0.0       0.0
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBE 2 
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        RADIUS 0.00000 <<rradf1>> <<rrads1>> SPLITR -3 1 
+        MIX  <<mfL1>> <<ms>> ;
+      ::: ROD2 := GEO: TUBE 2 
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        RADIUS 0.00000 <<rradf2>> <<rrads2>> SPLITR -3 1 
+        MIX  <<mfL2>> <<ms>> ;
+      ::: ROD3 := GEO: TUBE 2 
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        RADIUS 0.00000 <<rradf3>> <<rrads3>> SPLITR -3 1 
+        MIX  <<mfL3>> <<ms>> ;
+      ::: ROD4 := GEO: TUBE 2 
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        RADIUS 0.00000 <<rradf4>> <<rrads4>> SPLITR -3 1 
+        MIX  <<mfL4>> <<ms>> ;
+    ;
+*----
+*  LZC 
+*----
+    ::: LZC := GEO: CAR2D 2 4
+      MESHX   <<XLFmax>> 0.0        <<XRFmin>>
+      MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+      MIX     <<mmod>>   <<mmod>>
+              <<mmod>>   <<mmod>>
+              <<mmod>>   <<mmod>>
+              <<mmod>>   <<mmod>>
+    ;
+*----
+*  Fuel on the right 
+*----
+    ::: FuelR := GEO: CARCEL 8 3 4
+      MESHX   <<XRFmin>> <<XRFcent>> <<XRFint>> <<XRFmax>>
+      MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    <<rmod>>
+      OFFCENTER <<OffsetXR>> 0.0       0.0
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBE 2 
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        RADIUS 0.00000 <<rradf1>> <<rrads1>> SPLITR -3 1 
+        MIX  <<mfR1>> <<ms>> ;
+      ::: ROD2 := GEO: TUBE 2 
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        RADIUS 0.00000 <<rradf2>> <<rrads2>> SPLITR -3 1 
+        MIX  <<mfR2>> <<ms>> ;
+      ::: ROD3 := GEO: TUBE 2 
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        RADIUS 0.00000 <<rradf3>> <<rrads3>> SPLITR -3 1 
+        MIX  <<mfR3>> <<ms>> ;
+      ::: ROD4 := GEO: TUBE 2 
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        RADIUS 0.00000 <<rradf4>> <<rrads4>> SPLITR -3 1 
+        MIX  <<mfR4>> <<ms>> ;
+    ;
+    ;
+    EVALUATE NbAngles Density := 15 30.0 ;
+  ENDIF ;
+ELSEIF Type "XZ" = THEN
+  IF Mesh "Coarse" = THEN
+*----
+*  Coarse mesh geometry
+*----
+    Geometry := GEO: :: CAR2D 3 1
+      X- REFL X+ REFL Y- REFL Y+ REFL
+      CELL  FuelL LZC FuelR
+*----
+*  Fuel on the left
+*----
+    ::: FuelL := GEO: CAR2D 1 1
+      MESHX   <<XLFmin>> <<XLFmax>>
+      MESHY   <<ZFmin>>  <<ZFmax>>     
+      MIX     <<mmod>> 
+    ;
+*----
+*  LZC 
+*----
+    ::: LZC := GEO: CARCEL 2 1 1
+      MESHX   <<XLFmax>> <<XRFmin>>
+      MESHY   <<ZFmin>>  <<ZFmax>>     
+      RADIUS  0.00000    <<RinnGT>> <<RoutGT>>
+      MIX     <<mfill>>   <<mgt>>   <<mmod>>
+      CLUSTER S1 S2 S3 F1 F2 
+      ::: S1 := GEO: TUBE 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXs1>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>>
+        MIX <<ms11>> <<ms12>> <<ms13>> <<ms14>> ;
+      ::: S2 := GEO: TUBE 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXs2>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> 
+        MIX <<ms21>> <<ms22>> <<ms23>> <<ms24>> ;
+      ::: S3 := GEO: TUBE 4  
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXs3>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>>
+        MIX <<ms31>> <<ms32>> <<ms33>> <<ms34>> ;
+      ::: F1 := GEO: TUBE 4  
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXf1>>
+        RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>>
+        MIX <<mf11>> <<mf12>> <<mf13>> <<mf14>> ;
+      ::: F2 := GEO: TUBE 4  
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXf2>>
+        RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>>
+        MIX <<mf21>> <<mf22>> <<mf23>> <<mf24>> ;
+    ;
+*----
+*  Fuel on the right
+*----
+    ::: FuelR := GEO: CAR2D 1 1
+      MESHX   <<XRFmin>> <<XRFmax>>
+      MESHY   <<ZFmin>>  <<ZFmax>>     
+      MIX     <<mmod>> 
+    ;
+    ;
+    EVALUATE NbAngles Density := 15 30.0 ;
+  ELSE
+*----
+*  Fine mesh geometry (default)
+*----
+    Geometry := GEO: :: CAR2D 3 1
+      X- REFL X+ REFL Y- REFL Y+ REFL
+      CELL  FuelL LZC FuelR
+*----
+*  Fuel on the left 
+*----
+    ::: FuelL := GEO: CAR2D 3 4
+      MESHX   <<XLFmin>> <<XLFint>> <<XLFcent>> <<XLFmax>>
+      MESHY   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>     
+      MIX     <<mmod>> <<mmod>> <<mmod>> 
+              <<mmod>> <<mmod>> <<mmod>> 
+              <<mmod>> <<mmod>> <<mmod>> 
+              <<mmod>> <<mmod>> <<mmod>> 
+    ;
+*----
+*  LZC 
+*----
+    ::: LZC := GEO: CARCEL 2 2 4
+      MESHX   <<XLFmax>> 0.0        <<XRFmin>>
+      MESHY   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>     
+      RADIUS  0.00000    <<RinnGT>> <<RoutGT>>
+      MIX     <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+      CLUSTER S1 S2 S3 F1 F2 
+      ::: S1 := GEO: TUBE 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXs1>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>>
+        MIX <<ms11>> <<ms12>> <<ms13>> <<ms14>> ;
+      ::: S2 := GEO: TUBE 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXs2>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> 
+        MIX <<ms21>> <<ms22>> <<ms23>> <<ms24>> ;
+      ::: S3 := GEO: TUBE 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXs3>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>>
+        MIX <<ms31>> <<ms32>> <<ms33>> <<ms34>> ;
+      ::: F1 := GEO: TUBE 4  
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXf1>>
+        RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>>
+        MIX <<mf11>> <<mf12>> <<mf13>> <<mf14>> ;
+      ::: F2 := GEO: TUBE 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXf2>>
+        RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>>
+        MIX <<mf21>> <<mf22>> <<mf23>> <<mf24>> ;
+    ;
+*----
+*  Fuel on the right 
+*----
+    ::: FuelR := GEO: CAR2D 3 4
+      MESHX   <<XRFmin>> <<XRFcent>> <<XRFint>> <<XRFmax>>
+      MESHY   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>     
+      MIX     <<mmod>> <<mmod>> <<mmod>> 
+              <<mmod>> <<mmod>> <<mmod>> 
+              <<mmod>> <<mmod>> <<mmod>> 
+              <<mmod>> <<mmod>> <<mmod>> 
+    ;
+    ;
+    EVALUATE NbAngles Density := 15 30.0 ;
+  ENDIF ;
+ELSE
+  IF Mesh "Coarse" = THEN
+*----
+*  Coarse mesh geometry
+*----
+    Geometry := GEO: :: CAR3D 3 1 1
+      X- REFL X+ REFL Y- REFL Y+ SYME Z- REFL Z+ REFL
+      CELL  FuelL LZC FuelR
+*----
+*  Fuel on the left
+*----
+    ::: FuelL := GEO: CARCELZ 7 1 1 1
+      MESHX   <<XLFmin>> <<XLFmax>>
+      MESHY   <<YFmin>>  <<YFmax>>
+      MESHZ   <<ZFmin>>  <<ZFmax>>     
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    
+      OFFCENTER <<OffsetXL>> 0.0       0.0
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> 
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBEZ 2 1
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        MESHZ   <<ZFmin>> <<ZFmax>>     
+        RADIUS 0.00000 <<rradf1>> <<rrads1>>
+        MIX  <<mfL1>> <<ms>> ;
+      ::: ROD2 := GEO: TUBEZ 2 1
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        MESHZ   <<ZFmin>>  <<ZFmax>>     
+        RADIUS 0.00000 <<rradf2>> <<rrads2>> 
+        MIX  <<mfL2>> <<ms>>  ;
+      ::: ROD3 := GEO: TUBEZ 2 1
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        MESHZ   <<ZFmin>>  <<ZFmax>>     
+        RADIUS 0.00000 <<rradf3>> <<rrads3>>  
+        MIX  <<mfL3>> <<ms>>  ;
+      ::: ROD4 := GEO: TUBEZ 2 1
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        MESHZ   <<ZFmin>> <<ZFmax>>     
+        RADIUS 0.00000 <<rradf4>> <<rrads4>> 
+        MIX  <<mfL4>> <<ms>>  ;
+    ;
+*----
+*  LZC 
+*----
+    ::: LZC := GEO: CARCELY 2 1 1 1
+      MESHX   <<XLFmax>> <<XRFmin>>
+      MESHY   <<YFmin>>  <<YFmax>>
+      MESHZ   <<ZFmin>>  <<ZFmax>>     
+      RADIUS  0.00000    <<RinnGT>> <<RoutGT>>
+      MIX     <<mfill>>   <<mgt>>   <<mmod>>
+      CLUSTER S1 S2 S3 F1 F2 
+      ::: S1 := GEO: TUBEY 4 1 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZs1>>
+        MESHY <<YFmin>> <<YFmax>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>>
+        MIX <<ms11>> <<ms12>> <<ms13>> <<ms14>> ;
+      ::: S2 := GEO: TUBEY 4 1 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZs2>>
+        MESHY <<YFmin>> <<YFmax>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> 
+        MIX <<ms21>> <<ms22>> <<ms23>> <<ms24>> ;
+      ::: S3 := GEO: TUBEY 4 1 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZs3>>
+        MESHY <<YFmin>> <<YFmax>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>>
+        MIX <<ms31>> <<ms32>> <<ms33>> <<ms34>> ;
+      ::: F1 := GEO: TUBEY 4 1 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZf1>>
+        MESHY <<YFmin>> <<YFmax>>
+        RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>>
+        MIX <<mf11>> <<mf12>> <<mf13>> <<mf14>> ;
+      ::: F2 := GEO: TUBEY 4 1 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZf2>>
+        MESHY <<YFmin>> <<YFmax>>
+        RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>>
+        MIX <<mf21>> <<mf22>> <<mf23>> <<mf24>> ;
+    ;
+*----
+*  Fuel on the right
+*----
+    ::: FuelR := GEO: CARCELZ 7 1 1 1
+      MESHX   <<XRFmin>> <<XRFmax>>
+      MESHY   <<YFmin>>  <<YFmax>>
+      MESHZ   <<ZFmin>>  <<ZFmax>>     
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    
+      OFFCENTER <<OffsetXR>> 0.0       0.0
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> 
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBEZ 2 1
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        MESHZ   <<ZFmin>> <<ZFmax>>     
+        RADIUS 0.00000 <<rradf1>> <<rrads1>>
+        MIX  <<mfR1>> <<ms>> ;
+      ::: ROD2 := GEO: TUBEZ 2 1
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        MESHZ   <<ZFmin>>  <<ZFmax>>     
+        RADIUS 0.00000 <<rradf2>> <<rrads2>> 
+        MIX  <<mfR2>> <<ms>>  ;
+      ::: ROD3 := GEO: TUBEZ 2 1
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        MESHZ   <<ZFmin>>  <<ZFmax>>     
+        RADIUS 0.00000 <<rradf3>> <<rrads3>>  
+        MIX  <<mfR3>> <<ms>>  ;
+      ::: ROD4 := GEO: TUBEZ 2 1
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        MESHZ   <<ZFmin>> <<ZFmax>>     
+        RADIUS 0.00000 <<rradf4>> <<rrads4>> 
+        MIX  <<mfR4>> <<ms>>  ;
+    ;
+    ;
+    EVALUATE NbAngles Density := 8 100.0 ;
+  ELSE
+*----
+*  Fine mesh geometry (default)
+*----
+    Geometry := GEO: :: CAR3D 3 1 1
+      X- REFL X+ REFL Y- REFL Y+ SYME Z- REFL Z+ REFL
+      CELL  FuelL LZC FuelR
+*----
+*  Fuel on the left 
+*----
+    ::: FuelL := GEO: CARCELZ 8 3 4 4
+      MESHX   <<XLFmin>> <<XLFint>> <<XLFcent>> <<XLFmax>>
+      MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+      MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>     
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    <<rmod>>
+      OFFCENTER <<OffsetXL>> 0.0       0.0
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBEZ 2 4
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>     
+        RADIUS 0.00000 <<rradf1>> <<rrads1>> SPLITR -3 1 
+        MIX  <<mfL1>> <<ms>>
+             <<mfL1>> <<ms>>
+             <<mfL1>> <<ms>>
+             <<mfL1>> <<ms>>  ;
+      ::: ROD2 := GEO: TUBEZ 2 4
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>     
+        RADIUS 0.00000 <<rradf2>> <<rrads2>> SPLITR -3 1 
+        MIX  <<mfL2>> <<ms>>
+             <<mfL2>> <<ms>>
+             <<mfL2>> <<ms>>
+             <<mfL2>> <<ms>>  ;
+      ::: ROD3 := GEO: TUBEZ 2 4
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>     
+        RADIUS 0.00000 <<rradf3>> <<rrads3>> SPLITR -3 1 
+        MIX  <<mfL3>> <<ms>>
+             <<mfL3>> <<ms>>
+             <<mfL3>> <<ms>>
+             <<mfL3>> <<ms>>  ;
+      ::: ROD4 := GEO: TUBEZ 2 4
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>     
+        RADIUS 0.00000 <<rradf4>> <<rrads4>> SPLITR -3 1 
+        MIX  <<mfL4>> <<ms>>
+             <<mfL4>> <<ms>>
+             <<mfL4>> <<ms>>
+             <<mfL4>> <<ms>>  ;
+    ;
+*----
+*  LZC 
+*----
+    ::: LZC := GEO: CARCELY 2 2 4 4
+      MESHX   <<XLFmax>> 0.0        <<XRFmin>>
+      MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+      MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>     
+      RADIUS  0.00000    <<RinnGT>> <<RoutGT>>
+      MIX     <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+      CLUSTER S1 S2 S3 F1 F2 
+      ::: S1 := GEO: TUBEY 4 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZs1>>
+        MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>>
+        MIX <<ms11>> <<ms12>> <<ms13>> <<ms14>> 
+            <<ms11>> <<ms12>> <<ms13>> <<ms14>> 
+            <<ms11>> <<ms12>> <<ms13>> <<ms14>> 
+            <<ms11>> <<ms12>> <<ms13>> <<ms14>> ;
+      ::: S2 := GEO: TUBEY 4 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZs2>>
+        MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> 
+        MIX <<ms21>> <<ms22>> <<ms23>> <<ms24>> 
+            <<ms21>> <<ms22>> <<ms23>> <<ms24>> 
+            <<ms21>> <<ms22>> <<ms23>> <<ms24>> 
+            <<ms21>> <<ms22>> <<ms23>> <<ms24>> ;
+      ::: S3 := GEO: TUBEY 4 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZs3>>
+        MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>>
+        MIX <<ms31>> <<ms32>> <<ms33>> <<ms34>> 
+            <<ms31>> <<ms32>> <<ms33>> <<ms34>> 
+            <<ms31>> <<ms32>> <<ms33>> <<ms34>> 
+            <<ms31>> <<ms32>> <<ms33>> <<ms34>> ;
+      ::: F1 := GEO: TUBEY 4 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZf1>>
+        MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+        RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>>
+        MIX <<mf11>> <<mf12>> <<mf13>> <<mf14>> 
+            <<mf11>> <<mf12>> <<mf13>> <<mf14>> 
+            <<mf11>> <<mf12>> <<mf13>> <<mf14>> 
+            <<mf11>> <<mf12>> <<mf13>> <<mf14>> ;
+      ::: F2 := GEO: TUBEY 4 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZf2>>
+        MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+        RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>>
+        MIX <<mf21>> <<mf22>> <<mf23>> <<mf24>> 
+            <<mf21>> <<mf22>> <<mf23>> <<mf24>> 
+            <<mf21>> <<mf22>> <<mf23>> <<mf24>> 
+            <<mf21>> <<mf22>> <<mf23>> <<mf24>> ;
+    ;
+*----
+*  Fuel on the right 
+*----
+    ::: FuelR := GEO: CARCELZ 8 3 4 4
+      MESHX   <<XRFmin>> <<XRFcent>> <<XRFint>> <<XRFmax>>
+      MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+      MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>     
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    <<rmod>>
+      OFFCENTER <<OffsetXR>> 0.0       0.0
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBEZ 2 4
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>     
+        RADIUS 0.00000 <<rradf1>> <<rrads1>> SPLITR -3 1 
+        MIX  <<mfR1>> <<ms>>
+             <<mfR1>> <<ms>>
+             <<mfR1>> <<ms>>
+             <<mfR1>> <<ms>>  ;
+      ::: ROD2 := GEO: TUBEZ 2 4
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>     
+        RADIUS 0.00000 <<rradf2>> <<rrads2>> SPLITR -3 1 
+        MIX  <<mfR2>> <<ms>>
+             <<mfR2>> <<ms>>
+             <<mfR2>> <<ms>>
+             <<mfR2>> <<ms>>  ;
+      ::: ROD3 := GEO: TUBEZ 2 4
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>     
+        RADIUS 0.00000 <<rradf3>> <<rrads3>> SPLITR -3 1 
+        MIX  <<mfR3>> <<ms>>
+             <<mfR3>> <<ms>>
+             <<mfR3>> <<ms>>
+             <<mfR3>> <<ms>>  ;
+      ::: ROD4 := GEO: TUBEZ 2 4
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>     
+        RADIUS 0.00000 <<rradf4>> <<rrads4>> SPLITR -3 1 
+        MIX  <<mfR4>> <<ms>>
+             <<mfR4>> <<ms>>
+             <<mfR4>> <<ms>>
+             <<mfR4>> <<ms>>  ;
+    ;
+    ;
+    EVALUATE NbAngles Density := 16 1000.0 ;
+  ENDIF ;
+ENDIF ;
+::  <<NbAngles>> <<Density>>          ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/proc_LZC/GeoLZCM.c2m b/Dragon/data/proc_LZC/GeoLZCM.c2m
new file mode 100755
index 0000000..740d351
--- /dev/null
+++ b/Dragon/data/proc_LZC/GeoLZCM.c2m
@@ -0,0 +1,1020 @@
+*DECK GeoLZCM.c2m
+*----------
+*  Name          : GeoLZCM.c2m
+*  Mesh          : DRAGON procedure
+*  Use           : LZC geometry 
+*                  for G2 (cluster model)
+*  Author        : G. Marleau
+*  Date          : 2006/05/05
+*
+*  Procedure description: 
+*    Geometry := GeoLZCM :: 
+*      <<PrtLev>> <<Type>> <<Mesh>> 
+*      <<NbAngles>> <<Density>> ;
+*
+*  Procedure parameters:
+*    Geometry      : XSM_FILE containing the LZC geometry 
+*                    for G2
+*    PrtLev        : Print Level
+*                    0  -> Minimum print
+*                    1  -> Standard print
+*                    10 -> Debug
+*    Type          : Geometry type
+*                    'XY' -> 2-D projection of geometry in XY plane
+*                    'XZ' -> 2-D projection of geometry in XZ plane
+*                    '3D' -> 3-D geometry (default value if Type 
+*                            is not known)
+*    Mesh          : Geometry Mesh
+*                    Coarse -> Corase mesh 3-D geometry
+*                    Fine   -> Fine mesh 3-D geometry (default) 
+*    NbAngles      : Suggested angular quadrature for tracking
+*    Density       : Suggested line density for tracking
+* 
+*----------
+* Define procedure parameters and options and read
+* the procedure input data
+*----
+PARAMETER Geometry :: EDIT 0 ::: XSM_FILE Geometry ; ;
+STRING    Type Mesh                ;
+INTEGER   PrtLev NbAngles          ;
+REAL      Density                  ;
+::  >>PrtLev<< >>Type<< >>Mesh<<  ;
+*----
+*  Local procedures and modules used in this procedure
+*----     
+PROCEDURE    DimFuel DimCell DimLZC      ;
+MODULE       GEO: END: ;
+*----
+*  Get overall cell dimensions using procedure DimCell
+*----
+REAL         rcool4 rpt rgap rct LPitch BLength ;
+DimCell  ::  
+  >>rcool4<< >>rpt<< >>rgap<< >>rct<< >>LPitch<< >>BLength<< ;
+IF PrtLev  9 > THEN
+  ECHO "Overall cell dimensions" ;
+  ECHO rcool4 rpt rgap rct LPitch BLength ;
+ENDIF ;
+*----
+*  Get ring by ring description of fuel using procedure DimFuel
+*----
+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  ;
+EVALUATE    Ring := 1 ;
+DimFuel  :: <<Ring>>
+  >>npinp1<< >>rpinp1<< >>apinp1<< >>rradf1<< >>rrads1<< ;
+EVALUATE    Ring := 2 ;
+DimFuel  :: <<Ring>>
+  >>npinp2<< >>rpinp2<< >>apinp2<< >>rradf2<< >>rrads2<< ;
+EVALUATE    Ring := 3 ;
+DimFuel  :: <<Ring>>
+  >>npinp3<< >>rpinp3<< >>apinp3<< >>rradf3<< >>rrads3<< ;
+EVALUATE    Ring := 4 ;
+DimFuel  :: <<Ring>> 
+  >>npinp4<< >>rpinp4<< >>apinp4<< >>rradf4<< >>rrads4<< ;
+IF PrtLev  9 > THEN
+  ECHO "Fuel rings 1 to 4" ;
+  ECHO  npinp1   npinp2   npinp3   npinp4  ;
+  ECHO  rradf1   rrads1   rpinp1   apinp1  ;    
+  ECHO  rradf2   rrads2   rpinp2   apinp2  ;
+  ECHO  rradf3   rrads3   rpinp3   apinp3  ;
+  ECHO  rradf4   rrads4   rpinp4   apinp4  ;
+ENDIF ;
+*----
+*  Get LZC using procedure DimLZC
+*----
+REAL    RinnGT       RoutGT 
+        RinnBa       RoutBa  
+        RinnFe       RoutFe  
+        RinnBu       RoutBu  
+        RinnSc       RoutSc 
+         Rcent              ;
+DimLZC ::    
+  >>RinnGT<<   >>RoutGT<< 
+  >>RinnBa<< >>RoutBa<<  
+  >>RinnFe<< >>RoutFe<<  
+  >>RinnBu<< >>RoutBu<< 
+  >>RinnSc<< >>RoutSc<<  
+  >>Rcent<< ;
+IF PrtLev  9 > THEN
+  ECHO "LZC dimensions" ;
+  ECHO "Guide tube      " RinnGT RoutGT ;
+  ECHO "Feeder tubes    " RinnBa RoutBa  RinnFe RoutFe ;
+  ECHO "Scavenger tubes " RinnBu RoutBu  RinnSc RoutSc ;
+  ECHO "LZC location    " Rcent ;
+ENDIF ;
+*----
+*  Evaluate derived dimensions for supercell description
+*----
+REAL HLPitch   := LPitch 2.0 / ;
+REAL HBLength  := BLength 2.0 / ;
+REAL LimCoCell := HLPitch rct - RoutGT + 2.0 / ;
+REAL IntFCell  := LPitch LimCoCell - ;
+REAL IntYCell  := HLPitch LimCoCell - ;
+REAL OffsetXL OffsetXR := LimCoCell 2.0 / LimCoCell -2.0 / ;
+* Left and right fuel cell X position
+REAL XLFmin        XLFint          XLFcent        XLFmax 
+     XRFmin        XRFint          XRFcent        XRFmax :=
+     LPitch -1.0 * IntFCell -1.0 * HLPitch -1.0 * LimCoCell -1.0 *
+     LimCoCell     IntFCell        HLPitch        LPitch  ;
+* Cell Y position
+REAL YFmin   YFintb  YFcent YFintt  YFmax   :=
+     XLFcent IntYCell -1.0 *  0.0    IntYCell XRFcent ;
+* Z position
+REAL ZFmin           ZFintf  ZFcent ZFintb  ZFmax   :=
+     HBLength -1.0 * XLFmax  0.0    XRFmin  HBLength ;
+IF PrtLev  9 > THEN
+  ECHO "Derived dimensions for supercell" ;
+  ECHO  HLPitch HBLength  LimCoCell IntFCell  OffsetXL OffsetXR ;
+  ECHO  XLFmin XLFint XLFcent XLFmax ;
+  ECHO  XRFmin XRFint XRFcent XRFmax ;    
+  ECHO  YFmin  YFintb YFcent  YFintt  YFmax ;
+  ECHO  ZFmin         ZFintf  ZFcent ZFintb  ZFmax   ;
+ENDIF ;
+*----
+*  Derived dimensions for fuel 
+*----
+REAL         eps1  := 0.1   ;
+REAL         rcool1   rcool2   rcool3          :=
+             rrads1   eps1  + rpinp2 rrads2 - + 2.0 /
+             rpinp2 rrads2  + rpinp3 rrads3 - + 2.0 /
+             rpinp3 rrads3  + rpinp4 rrads4 - + 2.0 / ;
+REAL         rmod  := HLPitch LimCoCell - eps1 - ;
+IF PrtLev  9 > THEN
+  ECHO "Derived dimensions for fuel" ;
+  ECHO  eps1 rcool1 rcool2 rcool3 rmod ;
+ENDIF ;
+*----
+*  Define mixtures
+*----
+INTEGER  mcL1 mpt mgap mct mmod :=   1  2  3  4  5 ;
+INTEGER  mfL1 mfL2 mfL3 mfL4    :=   6  7  8  9  ;
+INTEGER  mfR1 mfR2 mfR3 mfR4    :=  10 11 12 13 ;
+INTEGER  ms                     :=  14          ;
+INTEGER  mfill mgt              :=  15 16       ;
+INTEGER  ms11      ms12      ms13      ms14    :=  
+         17        18        19        20  ;
+INTEGER  ms21      ms22      ms23      ms24    :=  
+         ms11 4 +  ms12 4 +  ms13 4 +  ms14 4 + ;
+INTEGER  ms31      ms32      ms33      ms34    :=  
+         ms11 8 +  ms12 8 +  ms13 8 +  ms14 8 + ;
+INTEGER  mf11      mf12      mf13      mf14    :=  
+         ms11 12 + ms12 12 + ms13 12 + ms14 12 + ;
+INTEGER  mf21      mf22      mf23      mf24    :=  
+         ms11 16 + ms12 16 + ms13 16 + ms14 16 + ;
+REAL     PIO2 PIO6 := $Pi_R 2.0 / $Pi_R 6.0 /    ;
+REAL     ApinZs1        ApinZs2        ApinZs3 
+         ApinZf1        ApinZf2                    :=
+         0.0            PIO6 4.0 *     PIO6 8.0 *   
+         PIO6 2.0 *     PIO6 6.0 *                 ;
+REAL     ApinXs1        ApinXs2        ApinXs3 
+         ApinXf1        ApinXf2                    :=
+         PIO2 ApinZs1 - PIO2 ApinZs2 - PIO2 ApinZs3 -
+         PIO2 ApinZf1 - PIO2 ApinZf2 -             ;
+ECHO     ApinZs1 ApinZs2 ApinZs3 ApinZf1 ApinZf2 ;
+IF PrtLev  9 > THEN
+  ECHO "Mixtures for coolant, PT, gap, CT and moderator" ;
+  ECHO  mcL1 mpt mgap mct mmod ;
+  ECHO "Mixtures for left fuel in each ring and for sheathing" ;
+  ECHO  mfL1 mfL2 mfL3 mfL4 ms ;
+  ECHO "Mixtures for right fuel in each ring and for sheathing" ;
+  ECHO  mfR1 mfR2 mfR3 mfR4 ms ;
+  ECHO "Mixtures for scavenger tube 1" ;
+  ECHO  ms11 ms12 ms13 ms14 ;
+  ECHO "Mixtures for scavenger tube 2" ;
+  ECHO  ms21 ms22 ms23 ms24 ;
+  ECHO "Mixtures for scavenger tube 3" ;
+  ECHO  ms31 ms32 ms33 ms34 ;
+  ECHO "Mixtures for feeder tube 1" ;
+  ECHO  mf11 mf12 mf13 mf14 ;
+  ECHO "Mixtures for feeder tube 12" ;
+  ECHO  mf11 mf22 mf23 mf24 ;
+ENDIF ;
+*----
+*  Geometry
+*----
+IF Type "XY" = THEN
+  IF Mesh "Coarse" = THEN
+*----
+*  Coarse mesh geometry
+*----
+    Geometry := GEO: :: CAR2D 3 1 
+      X- REFL X+ REFL Y- REFL Y+ SYME 
+      CELL  FuelL LZC FuelR
+*----
+*  Fuel on the left
+*----
+    ::: FuelL := GEO: CARCEL 7 1 1 
+      MESHX   <<XLFmin>> <<XLFmax>>
+      MESHY   <<YFmin>>  <<YFmax>>
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    
+      OFFCENTER <<OffsetXL>> 0.0       0.0
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> 
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBE 2 
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        RADIUS 0.00000 <<rradf1>> <<rrads1>>
+        MIX  <<mfL1>> <<ms>> ;
+      ::: ROD2 := GEO: TUBE 2 
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        RADIUS 0.00000 <<rradf2>> <<rrads2>> 
+        MIX  <<mfL2>> <<ms>>  ;
+      ::: ROD3 := GEO: TUBE 2 
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        RADIUS 0.00000 <<rradf3>> <<rrads3>>  
+        MIX  <<mfL3>> <<ms>>  ;
+      ::: ROD4 := GEO: TUBE 2 
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        RADIUS 0.00000 <<rradf4>> <<rrads4>> 
+        MIX  <<mfL4>> <<ms>>  ;
+    ;
+*----
+*  LZC 
+*----
+    ::: LZC := GEO: CAR2D 1 1 
+      MESHX   <<XLFmax>> <<XRFmin>>
+      MESHY   <<YFmin>>  <<YFmax>>
+      MIX     <<mmod>>
+    ;
+*----
+*  Fuel on the right
+*----
+    ::: FuelR := GEO: CARCEL 7 1 1 
+      MESHX   <<XRFmin>> <<XRFmax>>
+      MESHY   <<YFmin>>  <<YFmax>>
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    
+      OFFCENTER <<OffsetXR>> 0.0       0.0
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> 
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBE 2 
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        RADIUS 0.00000 <<rradf1>> <<rrads1>>
+        MIX  <<mfR1>> <<ms>> ;
+      ::: ROD2 := GEO: TUBE 2 
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        RADIUS 0.00000 <<rradf2>> <<rrads2>> 
+        MIX  <<mfR2>> <<ms>>  ;
+      ::: ROD3 := GEO: TUBE 2 
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        RADIUS 0.00000 <<rradf3>> <<rrads3>>  
+        MIX  <<mfR3>> <<ms>>  ;
+      ::: ROD4 := GEO: TUBE 2 
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        RADIUS 0.00000 <<rradf4>> <<rrads4>> 
+        MIX  <<mfR4>> <<ms>>  ;
+    ;
+    ;
+    EVALUATE NbAngles Density := 15 30.0 ;
+  ELSE
+*----
+*  Fine mesh geometry (default)
+*----
+    Geometry := GEO: :: CAR2D 3 1 
+      X- REFL X+ REFL Y- REFL Y+ SYME 
+      CELL  FuelL LZC FuelR
+*----
+*  Fuel on the left 
+*----
+    ::: FuelL := GEO: CARCEL 8 3 4
+      MESHX   <<XLFmin>> <<XLFint>> <<XLFcent>> <<XLFmax>>
+      MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    <<rmod>>
+      OFFCENTER <<OffsetXL>> 0.0       0.0
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBE 2 
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        RADIUS 0.00000 <<rradf1>> <<rrads1>> 
+        MIX  <<mfL1>> <<ms>> ;
+      ::: ROD2 := GEO: TUBE 2 
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        RADIUS 0.00000 <<rradf2>> <<rrads2>>  
+        MIX  <<mfL2>> <<ms>> ;
+      ::: ROD3 := GEO: TUBE 2 
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        RADIUS 0.00000 <<rradf3>> <<rrads3>>  
+        MIX  <<mfL3>> <<ms>> ;
+      ::: ROD4 := GEO: TUBE 2 
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        RADIUS 0.00000 <<rradf4>> <<rrads4>>  
+        MIX  <<mfL4>> <<ms>> ;
+    ;
+*----
+*  LZC 
+*----
+    ::: LZC := GEO: CAR2D 2 4
+      MESHX   <<XLFmax>> 0.0        <<XRFmin>>
+      MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+      MIX     <<mmod>>   <<mmod>>
+              <<mmod>>   <<mmod>>
+              <<mmod>>   <<mmod>>
+              <<mmod>>   <<mmod>>
+    ;
+*----
+*  Fuel on the right 
+*----
+    ::: FuelR := GEO: CARCEL 8 3 4
+      MESHX   <<XRFmin>> <<XRFcent>> <<XRFint>> <<XRFmax>>
+      MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    <<rmod>>
+      OFFCENTER <<OffsetXR>> 0.0       0.0
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBE 2 
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        RADIUS 0.00000 <<rradf1>> <<rrads1>> 
+        MIX  <<mfR1>> <<ms>> ;
+      ::: ROD2 := GEO: TUBE 2 
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        RADIUS 0.00000 <<rradf2>> <<rrads2>>  
+        MIX  <<mfR2>> <<ms>> ;
+      ::: ROD3 := GEO: TUBE 2 
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        RADIUS 0.00000 <<rradf3>> <<rrads3>> 
+        MIX  <<mfR3>> <<ms>> ;
+      ::: ROD4 := GEO: TUBE 2 
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        RADIUS 0.00000 <<rradf4>> <<rrads4>>  
+        MIX  <<mfR4>> <<ms>> ;
+    ;
+    ;
+    EVALUATE NbAngles Density := 15 30.0 ;
+  ENDIF ;
+ELSEIF Type "XZ" = THEN
+  IF Mesh "Coarse" = THEN
+*----
+*  Coarse mesh geometry
+*----
+    Geometry := GEO: :: CAR2D 3 1
+      X- REFL X+ REFL Y- REFL Y+ REFL
+      CELL  FuelL LZC FuelR
+*----
+*  Fuel on the left
+*----
+    ::: FuelL := GEO: CAR2D 1 1
+      MESHX   <<XLFmin>> <<XLFmax>>
+      MESHY   <<ZFmin>>  <<ZFmax>>     
+      MIX     <<mmod>> 
+    ;
+*----
+*  LZC 
+*----
+    ::: LZC := GEO: CARCEL 2 1 1
+      MESHX   <<XLFmax>> <<XRFmin>>
+      MESHY   <<ZFmin>>  <<ZFmax>>     
+      RADIUS  0.00000    <<RinnGT>> <<RoutGT>>
+      MIX     <<mfill>>   <<mgt>>   <<mmod>>
+      CLUSTER S1 S2 S3 F1 F2 
+      ::: S1 := GEO: TUBE 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXs1>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>>
+        MIX <<ms11>> <<ms12>> <<ms13>> <<ms14>> ;
+      ::: S2 := GEO: TUBE 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXs2>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> 
+        MIX <<ms21>> <<ms22>> <<ms23>> <<ms24>> ;
+      ::: S3 := GEO: TUBE 4  
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXs3>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>>
+        MIX <<ms31>> <<ms32>> <<ms33>> <<ms34>> ;
+      ::: F1 := GEO: TUBE 4  
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXf1>>
+        RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>>
+        MIX <<mf11>> <<mf12>> <<mf13>> <<mf14>> ;
+      ::: F2 := GEO: TUBE 4  
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXf2>>
+        RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>>
+        MIX <<mf21>> <<mf22>> <<mf23>> <<mf24>> ;
+    ;
+*----
+*  Fuel on the right
+*----
+    ::: FuelR := GEO: CAR2D 1 1
+      MESHX   <<XRFmin>> <<XRFmax>>
+      MESHY   <<ZFmin>>  <<ZFmax>>     
+      MIX     <<mmod>> 
+    ;
+    ;
+    EVALUATE NbAngles Density := 15 30.0 ;
+  ELSE
+*----
+*  Fine mesh geometry (default)
+*----
+    Geometry := GEO: :: CAR2D 3 1
+      X- REFL X+ REFL Y- REFL Y+ REFL
+      CELL  FuelL LZC FuelR
+*----
+*  Fuel on the left 
+*----
+    ::: FuelL := GEO: CAR2D 3 4
+      MESHX   <<XLFmin>> <<XLFint>> <<XLFcent>> <<XLFmax>>
+      MESHY   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>
+      MIX     <<mmod>> <<mmod>> <<mmod>> 
+              <<mmod>> <<mmod>> <<mmod>> 
+              <<mmod>> <<mmod>> <<mmod>> 
+              <<mmod>> <<mmod>> <<mmod>> 
+    ;
+*----
+*  LZC 
+*----
+    ::: LZC := GEO: CARCEL 2 2 4
+      MESHX   <<XLFmax>> 0.0        <<XRFmin>>
+      MESHY   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>
+      RADIUS  0.00000    <<RinnGT>> <<RoutGT>>
+      MIX     <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+      CLUSTER S1 S2 S3 F1 F2 
+      ::: S1 := GEO: TUBE 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXs1>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>>
+        MIX <<ms11>> <<ms12>> <<ms13>> <<ms14>> ;
+      ::: S2 := GEO: TUBE 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXs2>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> 
+        MIX <<ms21>> <<ms22>> <<ms23>> <<ms24>> ;
+      ::: S3 := GEO: TUBE 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXs3>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>>
+        MIX <<ms31>> <<ms32>> <<ms33>> <<ms34>> ;
+      ::: F1 := GEO: TUBE 4  
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXf1>>
+        RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>>
+        MIX <<mf11>> <<mf12>> <<mf13>> <<mf14>> ;
+      ::: F2 := GEO: TUBE 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinXf2>>
+        RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>>
+        MIX <<mf21>> <<mf22>> <<mf23>> <<mf24>> ;
+    ;
+*----
+*  Fuel on the right 
+*----
+    ::: FuelR := GEO: CAR2D 3 4
+      MESHX   <<XRFmin>> <<XRFcent>> <<XRFint>> <<XRFmax>>
+      MESHY   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>    
+      MIX     <<mmod>> <<mmod>> <<mmod>> 
+              <<mmod>> <<mmod>> <<mmod>> 
+              <<mmod>> <<mmod>> <<mmod>> 
+              <<mmod>> <<mmod>> <<mmod>> 
+    ;
+    ;
+    EVALUATE NbAngles Density := 15 30.0 ;
+  ENDIF ;
+ELSE
+  IF Mesh "Coarse" = THEN
+*----
+*  Coarse mesh geometry
+*----
+    Geometry := GEO: :: CAR3D 3 1 1
+      X- REFL X+ REFL Y- REFL Y+ SYME Z- REFL Z+ REFL
+      CELL  FuelL LZC FuelR
+*----
+*  Fuel on the left
+*----
+    ::: FuelL := GEO: CARCELZ 7 1 1 1
+      MESHX   <<XLFmin>> <<XLFmax>>
+      MESHY   <<YFmin>>  <<YFmax>>
+      MESHZ   <<ZFmin>>  <<ZFmax>>     
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    
+      OFFCENTER <<OffsetXL>> 0.0       0.0
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> 
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBEZ 2 1
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        MESHZ   <<ZFmin>> <<ZFmax>>     
+        RADIUS 0.00000 <<rradf1>> <<rrads1>>
+        MIX  <<mfL1>> <<ms>> ;
+      ::: ROD2 := GEO: TUBEZ 2 1
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        MESHZ   <<ZFmin>>  <<ZFmax>>     
+        RADIUS 0.00000 <<rradf2>> <<rrads2>> 
+        MIX  <<mfL2>> <<ms>>  ;
+      ::: ROD3 := GEO: TUBEZ 2 1
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        MESHZ   <<ZFmin>>  <<ZFmax>>     
+        RADIUS 0.00000 <<rradf3>> <<rrads3>>  
+        MIX  <<mfL3>> <<ms>>  ;
+      ::: ROD4 := GEO: TUBEZ 2 1
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        MESHZ   <<ZFmin>> <<ZFmax>>     
+        RADIUS 0.00000 <<rradf4>> <<rrads4>> 
+        MIX  <<mfL4>> <<ms>>  ;
+    ;
+*----
+*  LZC 
+*----
+    ::: LZC := GEO: CARCELY 2 1 1 1
+      MESHX   <<XLFmax>> <<XRFmin>>
+      MESHY   <<YFmin>>  <<YFmax>>
+      MESHZ   <<ZFmin>>  <<ZFmax>>     
+      RADIUS  0.00000    <<RinnGT>> <<RoutGT>>
+      MIX     <<mfill>>   <<mgt>>   <<mmod>>
+      CLUSTER S1 S2 S3 F1 F2 
+      ::: S1 := GEO: TUBEY 4 1 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZs1>>
+        MESHY <<YFmin>> <<YFmax>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>>
+        MIX <<ms11>> <<ms12>> <<ms13>> <<ms14>> ;
+      ::: S2 := GEO: TUBEY 4 1 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZs2>>
+        MESHY <<YFmin>> <<YFmax>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> 
+        MIX <<ms21>> <<ms22>> <<ms23>> <<ms24>> ;
+      ::: S3 := GEO: TUBEY 4 1 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZs3>>
+        MESHY <<YFmin>> <<YFmax>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>>
+        MIX <<ms31>> <<ms32>> <<ms33>> <<ms34>> ;
+      ::: F1 := GEO: TUBEY 4 1 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZf1>>
+        MESHY <<YFmin>> <<YFmax>>
+        RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>>
+        MIX <<mf11>> <<mf12>> <<mf13>> <<mf14>> ;
+      ::: F2 := GEO: TUBEY 4 1 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZf2>>
+        MESHY <<YFmin>> <<YFmax>>
+        RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>>
+        MIX <<mf21>> <<mf22>> <<mf23>> <<mf24>> ;
+    ;
+*----
+*  Fuel on the right
+*----
+    ::: FuelR := GEO: CARCELZ 7 1 1 1
+      MESHX   <<XRFmin>> <<XRFmax>>
+      MESHY   <<YFmin>>  <<YFmax>>
+      MESHZ   <<ZFmin>>  <<ZFmax>>     
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    
+      OFFCENTER <<OffsetXR>> 0.0       0.0
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> 
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBEZ 2 1
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        MESHZ   <<ZFmin>> <<ZFmax>>     
+        RADIUS 0.00000 <<rradf1>> <<rrads1>>
+        MIX  <<mfR1>> <<ms>> ;
+      ::: ROD2 := GEO: TUBEZ 2 1
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        MESHZ   <<ZFmin>>  <<ZFmax>>     
+        RADIUS 0.00000 <<rradf2>> <<rrads2>> 
+        MIX  <<mfR2>> <<ms>>  ;
+      ::: ROD3 := GEO: TUBEZ 2 1
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        MESHZ   <<ZFmin>>  <<ZFmax>>     
+        RADIUS 0.00000 <<rradf3>> <<rrads3>>  
+        MIX  <<mfR3>> <<ms>>  ;
+      ::: ROD4 := GEO: TUBEZ 2 1
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        MESHZ   <<ZFmin>> <<ZFmax>>     
+        RADIUS 0.00000 <<rradf4>> <<rrads4>> 
+        MIX  <<mfR4>> <<ms>>  ;
+    ;
+    ;
+    EVALUATE NbAngles Density := 8 50.0 ;
+  ELSE
+*----
+*  Fine mesh geometry (default)
+*----
+    Geometry := GEO: :: CAR3D 3 1 1
+      X- REFL X+ REFL Y- REFL Y+ SYME Z- REFL Z+ REFL
+      CELL  FuelL LZC FuelR
+*----
+*  Fuel on the left 
+*----
+    ::: FuelL := GEO: CARCELZ 8 3 4 4
+      MESHX   <<XLFmin>> <<XLFint>> <<XLFcent>> <<XLFmax>>
+      MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+      MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    <<rmod>>
+      OFFCENTER <<OffsetXL>> 0.0       0.0
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBEZ 2 4
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>
+        RADIUS 0.00000 <<rradf1>> <<rrads1>> 
+        MIX  <<mfL1>> <<ms>>
+             <<mfL1>> <<ms>>
+             <<mfL1>> <<ms>>
+             <<mfL1>> <<ms>>  ;
+      ::: ROD2 := GEO: TUBEZ 2 4
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>
+        RADIUS 0.00000 <<rradf2>> <<rrads2>> 
+        MIX  <<mfL2>> <<ms>>
+             <<mfL2>> <<ms>>
+             <<mfL2>> <<ms>>
+             <<mfL2>> <<ms>>  ;
+      ::: ROD3 := GEO: TUBEZ 2 4
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>
+        RADIUS 0.00000 <<rradf3>> <<rrads3>> 
+        MIX  <<mfL3>> <<ms>>
+             <<mfL3>> <<ms>>
+             <<mfL3>> <<ms>>
+             <<mfL3>> <<ms>>  ;
+      ::: ROD4 := GEO: TUBEZ 2 4
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>
+        RADIUS 0.00000 <<rradf4>> <<rrads4>>  
+        MIX  <<mfL4>> <<ms>>
+             <<mfL4>> <<ms>>
+             <<mfL4>> <<ms>>
+             <<mfL4>> <<ms>>  ;
+    ;
+*----
+*  LZC 
+*----
+    ::: LZC := GEO: CARCELY 2 2 4 4
+      MESHX   <<XLFmax>> 0.0        <<XRFmin>>
+      MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+      MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>
+      RADIUS  0.00000    <<RinnGT>> <<RoutGT>>
+      MIX     <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+              <<mfill>>   <<mgt>>   <<mmod>>
+      CLUSTER S1 S2 S3 F1 F2 
+      ::: S1 := GEO: TUBEY 4 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZs1>>
+        MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>>
+        MIX <<ms11>> <<ms12>> <<ms13>> <<ms14>> 
+            <<ms11>> <<ms12>> <<ms13>> <<ms14>> 
+            <<ms11>> <<ms12>> <<ms13>> <<ms14>> 
+            <<ms11>> <<ms12>> <<ms13>> <<ms14>> ;
+      ::: S2 := GEO: TUBEY 4 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZs2>>
+        MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> 
+        MIX <<ms21>> <<ms22>> <<ms23>> <<ms24>> 
+            <<ms21>> <<ms22>> <<ms23>> <<ms24>> 
+            <<ms21>> <<ms22>> <<ms23>> <<ms24>> 
+            <<ms21>> <<ms22>> <<ms23>> <<ms24>> ;
+      ::: S3 := GEO: TUBEY 4 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZs3>>
+        MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+        RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>>
+        MIX <<ms31>> <<ms32>> <<ms33>> <<ms34>> 
+            <<ms31>> <<ms32>> <<ms33>> <<ms34>> 
+            <<ms31>> <<ms32>> <<ms33>> <<ms34>> 
+            <<ms31>> <<ms32>> <<ms33>> <<ms34>> ;
+      ::: F1 := GEO: TUBEY 4 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZf1>>
+        MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+        RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>>
+        MIX <<mf11>> <<mf12>> <<mf13>> <<mf14>> 
+            <<mf11>> <<mf12>> <<mf13>> <<mf14>> 
+            <<mf11>> <<mf12>> <<mf13>> <<mf14>> 
+            <<mf11>> <<mf12>> <<mf13>> <<mf14>> ;
+      ::: F2 := GEO: TUBEY 4 4 
+        NPIN  1 RPIN <<Rcent>> APIN <<ApinZf2>>
+        MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+        RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>>
+        MIX <<mf21>> <<mf22>> <<mf23>> <<mf24>> 
+            <<mf21>> <<mf22>> <<mf23>> <<mf24>> 
+            <<mf21>> <<mf22>> <<mf23>> <<mf24>> 
+            <<mf21>> <<mf22>> <<mf23>> <<mf24>> ;
+    ;
+*----
+*  Fuel on the right 
+*----
+    ::: FuelR := GEO: CARCELZ 8 3 4 4
+      MESHX   <<XRFmin>> <<XRFcent>> <<XRFint>> <<XRFmax>>
+      MESHY   <<YFmin>>  <<YFintb>> <<YFcent>>  <<YFintt>> <<YFmax>>
+      MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>
+      RADIUS  0.00000    <<rcool1>> <<rcool2>>
+              <<rcool3>> <<rcool4>> <<rpt>>    
+              <<rgap>>   <<rct>>    <<rmod>>
+      OFFCENTER <<OffsetXR>> 0.0       0.0
+      MIX     <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>> 
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+              <<mcL1>>  <<mcL1>> <<mcL1>> <<mcL1>>   
+              <<mpt>>   <<mgap>> <<mct>>  <<mmod>> <<mmod>>
+      CLUSTER ROD1 ROD2 ROD3 ROD4
+      ::: ROD1 := GEO: TUBEZ 2 4
+        NPIN  <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>>
+        MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>
+        RADIUS 0.00000 <<rradf1>> <<rrads1>> 
+        MIX  <<mfR1>> <<ms>>
+             <<mfR1>> <<ms>>
+             <<mfR1>> <<ms>>
+             <<mfR1>> <<ms>>  ;
+      ::: ROD2 := GEO: TUBEZ 2 4
+        NPIN  <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>>
+        MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>
+        RADIUS 0.00000 <<rradf2>> <<rrads2>> 
+        MIX  <<mfR2>> <<ms>>
+             <<mfR2>> <<ms>>
+             <<mfR2>> <<ms>>
+             <<mfR2>> <<ms>>  ;
+      ::: ROD3 := GEO: TUBEZ 2 4
+        NPIN  <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>>
+        MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>
+        RADIUS 0.00000 <<rradf3>> <<rrads3>> 
+        MIX  <<mfR3>> <<ms>>
+             <<mfR3>> <<ms>>
+             <<mfR3>> <<ms>>
+             <<mfR3>> <<ms>>  ;
+      ::: ROD4 := GEO: TUBEZ 2 4
+        NPIN  <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>>
+        MESHZ   <<ZFmin>>  <<ZFintf>> <<ZFcent>>  <<ZFintb>> <<ZFmax>>
+        RADIUS 0.00000 <<rradf4>> <<rrads4>>  
+        MIX  <<mfR4>> <<ms>>
+             <<mfR4>> <<ms>>
+             <<mfR4>> <<ms>>
+             <<mfR4>> <<ms>>  ;
+    ;
+    ;
+    EVALUATE NbAngles Density := 16 100.0 ;
+  ENDIF ;
+ENDIF ;
+::  <<NbAngles>> <<Density>>          ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/proc_LZC/MicLZC.c2m b/Dragon/data/proc_LZC/MicLZC.c2m
new file mode 100644
index 0000000..81181b4
--- /dev/null
+++ b/Dragon/data/proc_LZC/MicLZC.c2m
@@ -0,0 +1,98 @@
+*DECK MicLZC.c2m
+*----------
+*  Name          : MicLZC.c2m
+*  Type          : DRAGON procedure
+*  Use           : Create Gentilly-2 library for LZC analysis.
+*                  (cluster model)
+*  Author        : G. Marleau
+*  Date          : 2006/03/01
+*  
+*  Procedure description: 
+*    MicroLib := MicLZC ::
+*      <<PrtLev>> <<LibSelect>> <<Position>> <<Type>> ;
+*  
+*    MicroLib      : XSM_FILE containing the microlib where:
+*                    Mix  1 = coolant at full density
+*                    Mix  2 = pressure tube
+*                    Mix  3 = helium for gap
+*                    Mix  4 = Zr for calandrial tube and adjuster
+*                    Mix  5 = moderator
+*                    Mix  6 = left central fuel
+*                    Mix  7 = left 6 pin fuel
+*                    Mix  8 = left 12 pin  fuel
+*                    Mix  9 = left outer fuel
+*                    Mix 10 = right central fuel
+*                    Mix 11 = right 6 pin fuel
+*                    Mix 12 = right 12 pin  fuel
+*                    Mix 13 = right outer fuel
+*                    Mix 14 = fuel sheating
+*                    Mix 15 = LZC zone to fill or empty
+*                    Mix 16 = Guide tube
+*                    Mix 17 = mixture in bubbler for scavenger 1
+*                    Mix 18 = mixture around bubbler for scavenger 1
+*                    Mix 19 = mixture in scavenger 1
+*                    Mix 20 = mixture around scavenger 1
+*                    Mix 21 = mixture in bubbler for scavenger 2
+*                    Mix 22 = mixture around bubbler for scavenger 2
+*                    Mix 23 = mixture in scavenger 2
+*                    Mix 24 = mixture around scavenger 2
+*                    Mix 25 = mixture in bubbler for scavenger 3
+*                    Mix 26 = mixture around bubbler for scavenger 3
+*                    Mix 27 = mixture in scavenger 3
+*                    Mix 28 = mixture around scavenger 3
+*                    Mix 29 = mixture in balance for feeder 1
+*                    Mix 30 = mixture around balance for feeder 1
+*                    Mix 31 = mixture in feeder 1
+*                    Mix 32 = mixture around feeder 1
+*                    Mix 33 = mixture in balance for feeder 2
+*                    Mix 34 = mixture around balance for feeder 2
+*                    Mix 35 = mixture in feeder 2
+*                    Mix 36 = mixture around feeder 2
+*    PrtLev        : Print Level
+*                    0  -> Minimum print
+*                    1  -> Standard print
+*                    10 -> Debug  
+*    LibSelect     : Selection of library with
+*                    WU  -> WLUP 69 groups WIMS-AECL library
+*                           (LibName = IAEA69).
+*                    WL  -> WIMFRITH 69 groups WIMS-AECL library
+*                           (LibName = WIMSLIB).
+*                    E5  -> ENDF/B-V 89 groups WIMS-AECL library
+*                           (LibName = E5WLIB).
+*                    E6  -> ENDF/B-VI 89 groups WIMS-AECL library
+*                           (LibName = E6MLIB).
+*    Position      : Cobalt adjusters position
+*                    NO  -> Adjuster and guide tubes absent
+*                    OUT -> Adjuster absent and guide tubes present
+*                    IN  -> Adjuster and guide tubes present
+*                           (default)
+*   Type           : 
+* 
+*----------
+* Define procedure parameters and options and read
+* the procedure input data
+*----
+PARAMETER MicroLib :: EDIT 0 ::: XSM_FILE MicroLib ; ;
+INTEGER   PrtLev ;
+STRING    Position LibSelect Type ;
+::  >>PrtLev<< >>LibSelect<< >>Position<< >>Type<< ;
+*----------
+*  Define procedures
+*----------
+PROCEDURE MicWUC MicWULZC  ;
+*PROCEDURE MicWU MicWLD  ;
+*PROCEDURE MicE5C MicE5D  ;
+*PROCEDURE MicE6C MicE6D  ;
+*----
+*  Get fuel XS
+*----
+IF LibSelect "WU" = THEN
+  MicroLib := MicWUC ::
+    <<PrtLev>> ;
+*----
+*  Get LZC cross sections
+*----
+  MicroLib := MicWULZC MicroLib ::
+    <<PrtLev>> <<Position>> <<Type>> ;
+ENDIF ;
+QUIT .
diff --git a/Dragon/data/proc_LZC/MicWUC.c2m b/Dragon/data/proc_LZC/MicWUC.c2m
new file mode 100644
index 0000000..d24bbec
--- /dev/null
+++ b/Dragon/data/proc_LZC/MicWUC.c2m
@@ -0,0 +1,189 @@
+*DECK MicWUC.c2m
+*----------
+*  Name          : MicWUC.c2m
+*  Type          : DRAGON procedure
+*  Use           : Create Gentilly-2 fuel microlib  
+*                  with cross sections from the WLUP based 69 
+*                  groups (LibName = IAEA69).
+*  Author        : G. Marleau
+*  Date          : 2006/03/01
+*
+*  Procedure description: 
+*    MicroLib := MicWUC ::
+*      <<PrtLev>> ;
+*  
+*  Procedure parameters:
+*    MicroLib      : XSM_FILE containing the new microlib where:
+*                    Mix  1 = coolant at full density
+*                    Mix  2 = pressure tube
+*                    Mix  3 = helium for gap and LZC
+*                    Mix  4 = Zr for calandrial tube and adjuster
+*                    Mix  5 = moderator
+*                    Mix  6 = left central fuel
+*                    Mix  7 = left 6 pin fuel
+*                    Mix  8 = left 12 pin  fuel
+*                    Mix  9 = left outer fuel
+*                    Mix 10 = right central fuel
+*                    Mix 11 = right 6 pin fuel
+*                    Mix 12 = right 12 pin  fuel
+*                    Mix 13 = right outer fuel
+*                    Mix 14 = fuel sheating
+*                    (see GeoLZC.c2m for mixture location in 
+*                     geometry)
+*    PrtLev        : Print Level
+*                    0  -> Minimum print
+*                    1  -> Standard print
+*                    10 -> Debug  
+* 
+*----------
+* Define procedure parameters and options and read
+* the procedure input data
+*----
+PARAMETER MicroLib :: EDIT 0 ::: XSM_FILE MicroLib ; ;
+INTEGER   PrtLev ;
+::  >>PrtLev<<   ;
+*----
+*  Local modules used in this procedure
+*----     
+MODULE  INFO: LIB:  DELETE:  ;
+*----------
+*  Define local fuel, coolant and moderator properties
+*  Fuel enrichment (%) =  0.71140;  T = 941.26 K; RHO =10.4375 g/cc;
+*  Coolant D2 AT (%)   = 99.222;    T = 560.66 K;
+*  Moderator D2 AT (%) = 99.911;    T = 345.66 K.
+*----------
+REAL
+  Ftemp   Fenri   Fdens 
+  Ctemp   Cprty 
+  Mtemp   Mprty            
+  Wprty          :=
+  941.26  0.7114  10.4375
+  560.66  99.30
+  345.66  99.92
+  0.0            ; 
+REAL
+  FwgtU8  FwgtU5  FwgtO16
+  CwgtH1  CwgtD2  CwgtO16  Cdens
+  MwgtH1  MwgtD2  MwgtO16  Mdens
+  WwgtH1  WwgtD2  WwgtO16  Wdens ; 
+REAL 
+  BNZrNb         BNZr         BNZrII         :=
+  0.00021        0.00031      0.00005962     ;  
+REAL 
+  B10ZrNb        B10Zr        B10ZrII 
+  B11ZrNb        B11Zr        B11ZrII        :=
+  0.199 BNZrNb * 0.199 BNZr * 0.199 BNZrII * 
+  0.801 BNZrNb * 0.801 BNZr * 0.801 BNZrII * ;  
+*----------
+*  Library and isotope names
+*----------
+STRING   LibName 
+         IsoU5    IsoU8    IsoU6    IsoPu9   IsoO16 
+         IsoH1    IsoD2    IsoB10   IsoB11   IsoBn
+         IsoZr    IsoFe    IsoFeA   IsoCr    IsoNi 
+         IsoNb    IsoC     IsoSi    IsoMn    IsoCu
+         IsoMo    IsoAl    IsoHe4   IsoN     IsoTi    
+         IsoCd    IsoCo59                          :=
+        "IAEA69"               
+        "2235"   "8238"   "236"    "6239"   "6016"
+        "3001"   "3002"   "10"     "11"     "1011"
+        "91"     "2056"   "2056"   "52"     "58"     
+        "93"     "2012"   "29"     "55"     "3063" 
+        "96"     "27"     "4"      "14"     "48"
+        "2113"   "1059"                            ; 
+*----------
+*  End region information
+*  endzrv   = Volume of Zr in endregion
+*  endcoolv = Volume of coolant in endregion
+*  endgasv  = Volume of gas (void) in endregion
+*  endtotv  = Total volume of endregion
+*  endzrf   = Volume fraction of Zr in endregion
+*  endcoolf = Volume fraction of coolant in endregion
+*  endgasf  = Volume fraction of gas (void) in endregion
+*---------- 
+REAL  endzrv endcoolv endgasv := 42.40  77.43    8.65 ;
+REAL  endtotv :=  endzrv endcoolv endgasv + +  ;          
+REAL  endzrf             endcoolf             endgasf  := 
+      endzrv endtotv /   endcoolv endtotv /   endgasv endtotv / ;
+*----------
+*  Guide tube perforation
+*---------- 
+REAL   PerfM1  PerfM2  PerfM3  :=   0.3566  0.02    0.005    ;
+REAL   PerfZ1        PerfZ2        PerfZ3        := 
+       1.0 PerfM1 -  1.0 PerfM2 -  1.0 PerfM3 -  ;
+*----------
+* COMPUTE WGT% FOR FUEL 
+* COMPUTE DENSITY AND WGT% FOR COOLANT AND MODERATOR 
+*----------         
+INFO: ::   
+  ENR: <<Fenri>> WGT%
+  LIB: WIMSD4 FIL: <<LibName>>
+  ISO: 3        <<IsoU5>>   <<IsoU8>>   <<IsoO16>> 
+  CALC WGT% UO2 >>FwgtU5<<  >>FwgtU8<<  >>FwgtO16<<   ;
+INFO: ::
+  TMP: <<Ctemp>> K 
+  PUR: <<Cprty>> WGT% 
+  CALC DENS WATER >>Cdens<<
+  LIB: WIMSD4 FIL: <<LibName>>
+  ISO: 3        <<IsoH1>>   <<IsoD2>>   <<IsoO16>> 
+  CALC WGT% D2O >>CwgtH1<<  >>CwgtD2<<  >>CwgtO16<<   ;
+INFO: ::
+  TMP: <<Mtemp>> K 
+  PUR: <<Mprty>> WGT% 
+  CALC DENS WATER >>Mdens<<
+  LIB: WIMSD4 FIL: <<LibName>>
+  ISO: 3        <<IsoH1>>   <<IsoD2>>   <<IsoO16>> 
+  CALC WGT% D2O >>MwgtH1<<  >>MwgtD2<<  >>MwgtO16<<   ;
+INFO: ::
+  TMP: <<Mtemp>> K 
+  PUR: <<Wprty>>  WGT% 
+  CALC DENS WATER >>Wdens<<
+  LIB: WIMSD4 FIL: <<LibName>>
+  ISO: 3        <<IsoH1>>   <<IsoD2>>   <<IsoO16>> 
+  CALC WGT% D2O >>WwgtH1<<  >>WwgtD2<<  >>WwgtO16<<   ;
+*----------
+*  Generate Microscopic cross section library
+*----------
+MicroLib := LIB: ::
+  EDIT <<PrtLev>>
+  NMIX 14  CTRA WIMS
+  DEPL LIB: WIMSD4 FIL: <<LibName>>
+  MIXS LIB: WIMSD4 FIL: <<LibName>>
+  MIX  1  <<Ctemp>>  <<Cdens>> 
+    H1  = <<IsoH1>>  <<CwgtH1>>     D2  = <<IsoD2>>  <<CwgtD2>>
+    O16 = <<IsoO16>> <<CwgtO16>>  
+  MIX  2  <<Ctemp>>  6.57
+    ZR  = <<IsoZr>>  97.5           NB  = <<IsoNb>>  2.5
+*    BN  = <<IsoBn>> <<BNZrNb>>   
+    B10 = <<IsoB10>>  <<B10ZrNb>>   B11 = <<IsoB11>>  <<B11ZrNb>>   
+  MIX  3  <<Mtemp>>  0.0014     
+    HE4 = <<IsoHe4>> 100.0        
+  MIX  4  <<Mtemp>>  6.44      
+    ZR  = <<IsoZr>>  99.71          FE  = <<IsoFe>>   0.16      
+    CR  = <<IsoCr>>  0.11           NI  = <<IsoNi>>   0.06    
+*    BN  = <<IsoBn>> <<BNZr>>   
+    B10 = <<IsoB10>> <<B10Zr>>      B11 = <<IsoB11>>  <<B11Zr>>   
+  MIX  5  <<Mtemp>>  <<Mdens>>  
+    H1  = <<IsoH1>>  <<MwgtH1>>     D2  = <<IsoD2>>  <<MwgtD2>> 
+    O16 = <<IsoO16>> <<MwgtO16>>  
+  MIX  6  <<Ftemp>>  <<Fdens>>    
+    O16 = <<IsoO16>> <<FwgtO16>>    U5  = <<IsoU5>>   <<FwgtU5>>  1
+    U8  = <<IsoU8>>  <<FwgtU8>>  1  U6  = <<IsoU6>>   0.0         1
+    PU9 = <<IsoPu9>> 0.0         1
+  MIX  7 COMB 6 100.0 
+  MIX  8 COMB 6 100.0 
+  MIX  9 COMB 6 100.0
+  MIX 10  <<Ftemp>>  <<Fdens>>    
+    O16 = <<IsoO16>> <<FwgtO16>>    U5  = <<IsoU5>>   <<FwgtU5>>  2
+    U8  = <<IsoU8>>  <<FwgtU8>>  2  U6  = <<IsoU6>>   0.0         2
+    PU9 = <<IsoPu9>> 0.0         2
+  MIX 11 COMB 10 100.0 
+  MIX 12 COMB 10 100.0
+  MIX 13 COMB 10 100.0
+  MIX 14  <<Ctemp>>  6.44       
+    ZR  = <<IsoZr>>  99.71          FE  = <<IsoFe>>  0.16      
+    CR  = <<IsoCr>>  0.11           NI  = <<IsoNi>>  0.06      
+*    BN   = <<IsoBn>> <<BNZr>>   
+    B10 = <<IsoB10>> <<B10Zr>>      B11 = <<IsoB11>>  <<B11Zr>>   
+  ;
+QUIT .
diff --git a/Dragon/data/proc_LZC/MicWULZC.c2m b/Dragon/data/proc_LZC/MicWULZC.c2m
new file mode 100644
index 0000000..0dcda38
--- /dev/null
+++ b/Dragon/data/proc_LZC/MicWULZC.c2m
@@ -0,0 +1,395 @@
+*DECK MicWULZC.c2m
+*----------
+*  Name          : MicWULZC.c2m
+*  Type          : DRAGON procedure
+*  Use           : Add LZC mixtures to Gentilly-2 fuel microlib 
+*                  with cross sections from the WLUP based 69 
+*                  groups  (LibName = IAEA69).
+*                  (cluster model)
+*  Author        : G. Marleau
+*  Date          : 2006/03/01
+*  
+*  Procedure description: 
+*    MicroLib := MicWULZC MicroLib ::
+*      <<PrtLev>> <<Position>> <<Type>> ;
+*  
+*    MicroLib      : XSM_FILE containing the updated microlib where:
+*                    Mix 15 = LZC zone to fill or empty
+*                    Mix 16 = Guide tube
+*                    Mix 17 = mixture in bubbler for scavenger 1
+*                    Mix 18 = mixture around bubbler for scavenger 1
+*                    Mix 19 = mixture in scavenger 1
+*                    Mix 20 = mixture around scavenger 1
+*                    Mix 21 = mixture in bubbler for scavenger 2
+*                    Mix 22 = mixture around bubbler for scavenger 2
+*                    Mix 23 = mixture in scavenger 2
+*                    Mix 24 = mixture around scavenger 2
+*                    Mix 25 = mixture in bubbler for scavenger 3
+*                    Mix 26 = mixture around bubbler for scavenger 3
+*                    Mix 27 = mixture in scavenger 3
+*                    Mix 28 = mixture around scavenger 3
+*                    Mix 29 = mixture in balance for feeder 1
+*                    Mix 30 = mixture around balance for feeder 1
+*                    Mix 31 = mixture in feeder 1
+*                    Mix 32 = mixture around feeder 1
+*                    Mix 33 = mixture in balance for feeder 2
+*                    Mix 34 = mixture around balance for feeder 2
+*                    Mix 35 = mixture in feeder 2
+*                    Mix 36 = mixture around feeder 2
+*    PrtLev        : Print Level
+*                    0  -> Minimum print
+*                    1  -> Standard print
+*                    10 -> Debug  
+*    Position      : LZC position
+*                    NO  -> LZC and guide tubes absent.
+*                    OUT -> LZC empty.
+*                    IN  -> LZC filled.
+*    Type          : LZC type
+*                    10  -> 1 scavenger tube 0 feeder tube (default)
+*                    21  -> 2 scavenger tubes 1 feeder tube
+*                    32  -> 3 scavenger tubes 2 feeder tubes
+*
+*
+* 
+*----------
+* Define procedure parameters and options and read
+* the procedure input data
+*----
+PARAMETER MicroLib :: EDIT 0 ::: XSM_FILE MicroLib ; ;
+INTEGER   PrtLev ;
+STRING    Position Type ;
+::  >>PrtLev<< >>Position<< >>Type<< ;
+*----------
+*  Define modules
+*----------
+MODULE  INFO: LIB:  DELETE:  ;
+*----------
+*  Library and isotope names
+*----------
+REAL    Mtemp Mprty := 345.66 99.92 ;
+REAL    MwgtH1  MwgtD2  MwgtO16  Mdens ;
+REAL    Wtemp Wprty := Mtemp  0.001  ;
+REAL    WwgtH1  WwgtD2  WwgtO16  Wdens ;
+REAL 
+  BNZrNb         BNZr         BNZrII         :=
+  0.00021        0.00031      0.00005962     ;  
+REAL 
+  B10ZrNb        B10Zr        B10ZrII 
+  B11ZrNb        B11Zr        B11ZrII        :=
+  0.199 BNZrNb * 0.199 BNZr * 0.199 BNZrII * 
+  0.801 BNZrNb * 0.801 BNZr * 0.801 BNZrII * ;  
+*----------
+*  Library and isotope names
+*----------
+STRING   LibName 
+         IsoU5    IsoU8    IsoU6    IsoPu9   IsoO16 
+         IsoH1    IsoD2    IsoB10   IsoB11   IsoBn
+         IsoZr    IsoFe    IsoFeA   IsoCr    IsoNi 
+         IsoNb    IsoC     IsoSi    IsoMn    IsoCu
+         IsoMo    IsoAl    IsoHe4   IsoN     IsoTi    
+         IsoCd    IsoCo59                          :=
+        "IAEA69"               
+        "2235"   "8238"   "236"    "6239"   "6016"
+        "3001"   "3002"   "10"     "11"     "1011"
+        "91"     "2056"   "2056"   "52"     "58"     
+        "93"     "2012"   "29"     "55"     "3063" 
+        "96"     "27"     "4"      "14"     "48"
+        "2113"   "1059"                            ; 
+INFO: ::
+  TMP: <<Mtemp>> K 
+  PUR: <<Mprty>> WGT% 
+  CALC DENS WATER >>Mdens<<
+  LIB: WIMSD4 FIL: <<LibName>>
+  ISO: 3        <<IsoH1>>   <<IsoD2>>   <<IsoO16>> 
+  CALC WGT% D2O >>MwgtH1<<  >>MwgtD2<<  >>MwgtO16<<   ;
+INFO: ::
+  TMP: <<Wtemp>> K 
+  PUR: <<Wprty>> WGT% 
+  CALC DENS WATER >>Wdens<<
+  LIB: WIMSD4 FIL: <<LibName>>
+  ISO: 3        <<IsoH1>>   <<IsoD2>>   <<IsoO16>> 
+  CALC WGT% D2O >>WwgtH1<<  >>WwgtD2<<  >>WwgtO16<<   ;
+*----------
+*  Complete Microscopic cross section library
+*----------
+IF Position "IN" = THEN
+  IF Type "32" = THEN
+    MicroLib   := LIB: MicroLib ::
+      EDIT <<PrtLev>>
+      NMIX 36  CTRA WIMS
+      MIXS LIB: WIMSD4 FIL: <<LibName>>
+      MIX  15  <<Wtemp>>  <<Wdens>>  
+        H1  = <<IsoH1>>  <<WwgtH1>>     D2  = <<IsoD2>>  <<WwgtD2>> 
+        O16 = <<IsoO16>> <<WwgtO16>>  
+      MIX 16 <<Mtemp>>   6.4003        
+        ZR  = <<IsoZr>>   98.2083        FE = <<IsoFe>>   0.135 
+        CR  = <<IsoCr>>   0.10           NI = <<IsoNi>>   0.055
+*        BN  = <<IsoBn>>   <<BNZrII>>
+        B10 = <<IsoB10>>  <<B10ZrII>>    B11 = <<IsoB11>>  <<B10ZrII>> 
+*Scavanger 1 
+      MIX  17 <<Mtemp>>  0.0014     
+        HE4 = <<IsoHe4>> 100.0  
+      MIX  18 COMB 16 100.0 
+      MIX  19 COMB 15 100.0  
+      MIX  20 COMB 16 100.0 
+*Scavanger 2 
+      MIX  21 COMB 17 100.0 
+      MIX  22 COMB 18 100.0 
+      MIX  23 COMB 19 100.0 
+      MIX  24 COMB 20 100.0 
+*Scavanger 3 
+      MIX  25 COMB 17 100.0 
+      MIX  26 COMB 18 100.0 
+      MIX  27 COMB 19 100.0 
+      MIX  28 COMB 20 100.0 
+*Feeder 1
+      MIX  29 COMB 17 100.0 
+      MIX  30 COMB 18 100.0 
+      MIX  31 COMB 19 100.0 
+      MIX  32 COMB 20 100.0 
+*Feeder 2
+      MIX  33 COMB 17 100.0 
+      MIX  34 COMB 18 100.0 
+      MIX  35 COMB 19 100.0 
+      MIX  36 COMB 20 100.0 
+      ;
+  ELSEIF Type "21" = THEN
+    MicroLib   := LIB: MicroLib ::
+      EDIT <<PrtLev>>
+      NMIX 36  CTRA WIMS
+      MIXS LIB: WIMSD4 FIL: <<LibName>>
+      MIX  15  <<Wtemp>>  <<Wdens>>  
+        H1  = <<IsoH1>>  <<WwgtH1>>     D2  = <<IsoD2>>  <<WwgtD2>> 
+        O16 = <<IsoO16>> <<WwgtO16>>  
+      MIX 16 <<Mtemp>>   6.4003        
+        ZR  = <<IsoZr>>   98.2083        FE = <<IsoFe>>   0.135 
+        CR  = <<IsoCr>>   0.10           NI = <<IsoNi>>   0.055
+*        BN  = <<IsoBn>>   <<BNZrII>>
+        B10 = <<IsoB10>>  <<B10ZrII>>    B11 = <<IsoB11>>  <<B10ZrII>> 
+*Scavanger 1 
+      MIX  17 <<Mtemp>>  0.0014     
+        HE4 = <<IsoHe4>> 100.0  
+      MIX  18 COMB 16 100.0 
+      MIX  19 COMB 15 100.0  
+      MIX  20 COMB 16 100.0 
+*Scavanger 2 
+      MIX  21 COMB 17 100.0 
+      MIX  22 COMB 18 100.0 
+      MIX  23 COMB 19 100.0 
+      MIX  24 COMB 20 100.0 
+*Scavanger 3 
+      MIX  25 COMB 15 100.0 
+      MIX  26 COMB 15 100.0 
+      MIX  27 COMB 15 100.0 
+      MIX  28 COMB 15 100.0 
+*Feeder 1
+      MIX  29 COMB 17 100.0 
+      MIX  30 COMB 18 100.0 
+      MIX  31 COMB 19 100.0 
+      MIX  32 COMB 20 100.0 
+*Feeder 2
+      MIX  33 COMB 15 100.0 
+      MIX  34 COMB 15 100.0 
+      MIX  35 COMB 15 100.0 
+      MIX  36 COMB 15 100.0 
+      ;
+  ELSE 
+    MicroLib   := LIB: MicroLib ::
+      EDIT <<PrtLev>>
+      NMIX 36  CTRA WIMS
+      MIXS LIB: WIMSD4 FIL: <<LibName>>
+      MIX  15  <<Wtemp>>  <<Wdens>>  
+        H1  = <<IsoH1>>  <<WwgtH1>>     D2  = <<IsoD2>>  <<WwgtD2>> 
+        O16 = <<IsoO16>> <<WwgtO16>>  
+      MIX 16 <<Mtemp>>   6.4003        
+        ZR  = <<IsoZr>>   98.2083        FE = <<IsoFe>>   0.135 
+        CR  = <<IsoCr>>   0.10           NI = <<IsoNi>>   0.055
+*        BN  = <<IsoBn>>   <<BNZrII>>
+        B10 = <<IsoB10>>  <<B10ZrII>>    B11 = <<IsoB11>>  <<B10ZrII>> 
+*Scavanger 1 
+      MIX  17 <<Mtemp>>  0.0014     
+        HE4 = <<IsoHe4>> 100.0  
+      MIX  18 COMB 16 100.0 
+      MIX  19 COMB 15 100.0  
+      MIX  20 COMB 16 100.0 
+*Scavanger 2 
+      MIX  21 COMB 15 100.0 
+      MIX  22 COMB 15 100.0 
+      MIX  23 COMB 15 100.0 
+      MIX  24 COMB 15 100.0 
+*Scavanger 3 
+      MIX  25 COMB 15 100.0 
+      MIX  26 COMB 15 100.0 
+      MIX  27 COMB 15 100.0 
+      MIX  28 COMB 15 100.0 
+*Feeder 1
+      MIX  29 COMB 15 100.0 
+      MIX  30 COMB 15 100.0 
+      MIX  31 COMB 15 100.0 
+      MIX  32 COMB 15 100.0 
+*Feeder 2
+      MIX  33 COMB 15 100.0 
+      MIX  34 COMB 15 100.0 
+      MIX  35 COMB 15 100.0 
+      MIX  36 COMB 15 100.0 
+      ;
+  ENDIF ;
+ELSEIF Position "OUT" = THEN
+  IF Type "32" = THEN
+    MicroLib   := LIB: MicroLib ::
+      EDIT <<PrtLev>>
+      NMIX 36  CTRA WIMS
+      MIXS LIB: WIMSD4 FIL: <<LibName>>
+      MIX  15  <<Mtemp>>  0.0014     
+        HE4 = <<IsoHe4>> 100.0        
+      MIX 16 <<Mtemp>>   6.4003        
+        ZR  = <<IsoZr>>   98.2083        FE = <<IsoFe>>   0.135 
+        CR  = <<IsoCr>>   0.10           NI = <<IsoNi>>   0.055
+*        BN  = <<IsoBn>>   <<BNZrII>>
+        B10 = <<IsoB10>>  <<B10ZrII>>    B11 = <<IsoB11>>  <<B10ZrII>>
+*Scavanger 1 
+      MIX  17 COMB 15 100.0 
+      MIX  18 COMB 16 100.0 
+      MIX  19  <<Wtemp>>  <<Wdens>>  
+        H1  = <<IsoH1>>  <<WwgtH1>>     D2  = <<IsoD2>>  <<WwgtD2>> 
+        O16 = <<IsoO16>> <<WwgtO16>>  
+      MIX  20 COMB 16 100.0 
+*Scavanger 2 
+      MIX  21 COMB 17 100.0 
+      MIX  22 COMB 18 100.0 
+      MIX  23 COMB 19 100.0 
+      MIX  24 COMB 20 100.0 
+*Scavanger 3 
+      MIX  25 COMB 17 100.0 
+      MIX  26 COMB 18 100.0 
+      MIX  27 COMB 19 100.0 
+      MIX  28 COMB 20 100.0 
+*Feeder 1
+      MIX  29 COMB 17 100.0 
+      MIX  30 COMB 18 100.0 
+      MIX  31 COMB 19 100.0 
+      MIX  32 COMB 20 100.0 
+*Feeder 2
+      MIX  33 COMB 17 100.0 
+      MIX  34 COMB 18 100.0 
+      MIX  35 COMB 19 100.0 
+      MIX  36 COMB 20 100.0 
+      ;
+  ELSEIF Type "21" = THEN
+    MicroLib   := LIB: MicroLib ::
+      EDIT <<PrtLev>>
+      NMIX 36  CTRA WIMS
+      MIXS LIB: WIMSD4 FIL: <<LibName>>
+      MIX  15  <<Mtemp>>  0.0014     
+        HE4 = <<IsoHe4>> 100.0        
+      MIX 16 <<Mtemp>>   6.4003        
+        ZR  = <<IsoZr>>   98.2083        FE = <<IsoFe>>   0.135 
+        CR  = <<IsoCr>>   0.10           NI = <<IsoNi>>   0.055
+*        BN  = <<IsoBn>>   <<BNZrII>>
+        B10 = <<IsoB10>>  <<B10ZrII>>    B11 = <<IsoB11>>  <<B10ZrII>> 
+*Scavanger 1 
+      MIX  17 COMB 15 100.0 
+      MIX  18 COMB 16 100.0 
+      MIX  19  <<Wtemp>>  <<Wdens>>  
+        H1  = <<IsoH1>>  <<WwgtH1>>     D2  = <<IsoD2>>  <<WwgtD2>> 
+        O16 = <<IsoO16>> <<WwgtO16>>  
+      MIX  20 COMB 16 100.0 
+*Scavanger 2 
+      MIX  21 COMB 17 100.0 
+      MIX  22 COMB 18 100.0 
+      MIX  23 COMB 19 100.0 
+      MIX  24 COMB 20 100.0 
+*Scavanger 3 
+      MIX  25 COMB 15 100.0 
+      MIX  26 COMB 15 100.0 
+      MIX  27 COMB 15 100.0 
+      MIX  28 COMB 15 100.0 
+*Feeder 1
+      MIX  29 COMB 17 100.0 
+      MIX  30 COMB 18 100.0 
+      MIX  31 COMB 19 100.0 
+      MIX  32 COMB 20 100.0 
+*Feeder 2
+      MIX  33 COMB 15 100.0 
+      MIX  34 COMB 15 100.0 
+      MIX  35 COMB 15 100.0 
+      MIX  36 COMB 15 100.0 
+      ;
+  ELSE 
+    MicroLib   := LIB: MicroLib ::
+      EDIT <<PrtLev>>
+      NMIX 36  CTRA WIMS
+      MIXS LIB: WIMSD4 FIL: <<LibName>>
+      MIX  15  <<Mtemp>>  0.0014     
+        HE4 = <<IsoHe4>> 100.0        
+      MIX 16 <<Mtemp>>   6.4003        
+        ZR  = <<IsoZr>>   98.2083        FE = <<IsoFe>>   0.135 
+        CR  = <<IsoCr>>   0.10           NI = <<IsoNi>>   0.055
+*        BN  = <<IsoBn>>   <<BNZrII>>
+        B10 = <<IsoB10>>  <<B10ZrII>>    B11 = <<IsoB11>>  <<B10ZrII>> 
+*Scavanger 1 
+      MIX  17 COMB 15 100.0 
+      MIX  18 COMB 16 100.0 
+      MIX  19  <<Wtemp>>  <<Wdens>>  
+        H1  = <<IsoH1>>  <<WwgtH1>>     D2  = <<IsoD2>>  <<WwgtD2>> 
+        O16 = <<IsoO16>> <<WwgtO16>>  
+      MIX  20 COMB 16 100.0 
+*Scavanger 2 
+      MIX  21 COMB 15 100.0 
+      MIX  22 COMB 15 100.0 
+      MIX  23 COMB 15 100.0 
+      MIX  24 COMB 15 100.0 
+*Scavanger 3 
+      MIX  25 COMB 15 100.0 
+      MIX  26 COMB 15 100.0 
+      MIX  27 COMB 15 100.0 
+      MIX  28 COMB 15 100.0 
+*Feeder 1
+      MIX  29 COMB 15 100.0 
+      MIX  30 COMB 15 100.0 
+      MIX  31 COMB 15 100.0 
+      MIX  32 COMB 15 100.0 
+*Feeder 2
+      MIX  33 COMB 15 100.0 
+      MIX  34 COMB 15 100.0 
+      MIX  35 COMB 15 100.0 
+      MIX  36 COMB 15 100.0 
+      ;
+  ENDIF ;
+ELSE
+  MicroLib   := LIB: MicroLib ::
+    EDIT <<PrtLev>>
+    NMIX 36  CTRA WIMS
+    MIXS LIB: WIMSD4 FIL: <<LibName>>
+    MIX  15  <<Mtemp>>  <<Mdens>>  
+      H1  = <<IsoH1>>  <<MwgtH1>>     D2  = <<IsoD2>>  <<MwgtD2>> 
+      O16 = <<IsoO16>> <<MwgtO16>>  
+    MIX  16 COMB 15 100.0 
+*Scavanger 1 
+    MIX  17 COMB 15 100.0 
+    MIX  18 COMB 15 100.0 
+    MIX  19 COMB 15 100.0 
+    MIX  20 COMB 15 100.0 
+*Scavanger 2 
+    MIX  21 COMB 15 100.0 
+    MIX  22 COMB 15 100.0 
+    MIX  23 COMB 15 100.0 
+    MIX  24 COMB 15 100.0 
+*Scavanger 3 
+    MIX  25 COMB 15 100.0 
+    MIX  26 COMB 15 100.0 
+    MIX  27 COMB 15 100.0 
+    MIX  28 COMB 15 100.0 
+*Feeder 1
+    MIX  29 COMB 15 100.0 
+    MIX  30 COMB 15 100.0 
+    MIX  31 COMB 15 100.0 
+    MIX  32 COMB 15 100.0 
+*Feeder 2
+    MIX  33 COMB 15 100.0 
+    MIX  34 COMB 15 100.0 
+    MIX  35 COMB 15 100.0 
+    MIX  36 COMB 15 100.0 
+    ;
+ENDIF ;
+QUIT .