*DECK Geo_N1
*----
*  Name          : Geo_N1.c2m
*  Type          : DRAGON procedure
*  Use           : Geometry generation for 1st Level Flux Calculation
*                  of the CASEA assembly
*  Author        : A. Canbakan
*
*  Procedure called as: 
*
*GEON1 := Geo_N1 :: <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>>
*                   <<R_int_TG>> <<R_ext_TG>> <<R_int_TI>> 
*                   <<R_ext_TI>> <<R_int_G>> <<R_ext_G>> <<Cote>>
*                   <<CoteL>> <<RmodeN1>> ;
*
*  Input data   :
*    Rcomb1     :  50/100 of outer radius of fuel (cm)
*    Rcomb2     :  80/100 of outer radius of fuel (cm)
*    Rcomb3     :  95/100 of outer radius of fuel (cm)
*    Rcomb4     : 100/100 of outer radius of fuel (cm)
*    R_int_TG   : Inner radius of cladding of guide tube (cm)
*    R_ext_TG   : Outer radius of cladding of guide tube (cm)
*    R_int_TI   : Inner radius of cladding of guide tube (cm)
*    R_ext_TI   : Outer radius of cladding of guide tube (cm)
*    R_int_G    : Inner radius of cladding of fuel tube (cm)
*    R_ext_G    : Outer radius of cladding of fuel tube (cm)
*    Cote       : Lattice pitch (cm)
*    CoteL      : Lattice pitch + Water space (cm)
*    RmodeN1    : Parameter for flux calculation in Level 1 (cm)
*
*  Output data  :
*    GEON1      : Geometry for Self-Shielding Calculation

PARAMETER  GEON1  ::  
       EDIT 0 
           ::: LINKED_LIST GEON1  ; 
   ;
*----
*  Modules used in this procedure
*----
MODULE  GEO: END: ;

*----
*  Input data recovery
*----
*                                                                      ;
REAL Rcomb1       Rcomb2       Rcomb3       Rcomb4     ;
:: >>Rcomb1<<   >>Rcomb2<<   >>Rcomb3<<   >>Rcomb4<<   ;
REAL R_int_TG     R_ext_TG     R_int_TI     R_ext_TI   ;
:: >>R_int_TG<< >>R_ext_TG<< >>R_int_TI<< >>R_ext_TI<< ;
REAL R_int_G      R_ext_G      Cote         CoteL      ;
:: >>R_int_G<<  >>R_ext_G<<  >>Cote<<     >>CoteL<<    ;
REAL RmodeN1   ;
:: >>RmodeN1<< ;
REAL Lame := CoteL Cote - ;

GEON1 := GEO: :: CAR2D 10 10
  EDIT 0
  X- DIAG X+ REFL
  Y- SYME Y+ DIAG
  CELL TI C1 C1 TG C1 C1 TG C4 C6 C9
          C2 C2 C1 C2 C2 C1 C2 C6 C9
             C2 C1 C2 C2 C1 C2 C6 C9
                TG C1 C1 TG C4 C6 C9
                   C2 C1 C1 C2 C6 C9
                      TG C1 C3 C6 C9
                         C2 C3 C6 C9
                            C5 C7 C9
                               C8 C9
                                  C10

  MERGE 11 1  1 10  1  1  9  4  6 12
           2  2  1  2  2  1  2  6 12
              2  1  2  2  1  2  6 12
                10  1  1  9  4  6 12
                    2  1  1  2  6 12
                       9  1  3  6 12
                          2  3  6 12
                             5  7 12
                                8 12
                                  13 

  TURN  A  A  E  A  A  E  A  A  A  A
           A  E  F  A  E  D  A  A  A
              C  B  G  C  H  G  A  A
                 A  G  C  A  G  A  A
                    A  B  D  A  A  A
                       A  A  A  A  A
                          A  G  A  A
                             A  A  A
                                A  A
                                   A
  ::: C1 := GEO: CARCEL 7
         MESHX 0.0 <<Cote>>
         MESHY 0.0 <<Cote>>
         RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>>
                <<R_int_G>> <<R_ext_G>> <<RmodeN1>>
         MIX  3 4 5 6   7 8 1 1
  ;
  ::: C2 := GEO: C1
         MIX 9 10 11 12   7 8 1 1
  ;
  ::: C3 := GEO: C1
         MIX 114 115 116 117   7 8 1 1
  ;
  ::: C4 := GEO: C1
         MIX 94 95 96 97   7 8 1 1
  ;
  ::: C5 := GEO: C1
         MIX 122 123 124 125   7 8 1 1
  ;
  ::: C6 := GEO: CARCEL 7
         MESHX 0.0 <<Cote>>
         MESHY 0.0 <<Cote>>
         RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>>
                <<R_int_G>> <<R_ext_G>> <<RmodeN1>>
         MIX 126 127 128 129   7 8 1 1
  ;
  ::: C7 := GEO: C6
         MIX 154 155 156 157   7 8 1 1
  ;
  ::: C8 := GEO: CARCEL 7
         MESHX 0.0 <<Cote>>
         MESHY 0.0 <<Cote>>
         RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>>
                <<R_int_G>> <<R_ext_G>> <<RmodeN1>>
         MIX 158 159 160 161   7 8 1 1
  ;
  ::: TG := GEO: CARCEL 2
         MESHX 0.0 <<Cote>>
         MESHY 0.0 <<Cote>>
         RADIUS 0.0 <<R_int_TG>> <<R_ext_TG>>
         MIX 1 25 1
  ;
  ::: TI := GEO: CARCEL 2
         MESHX 0.0 <<Cote>>
         MESHY 0.0 <<Cote>>
         RADIUS 0.0 <<R_int_TI>> <<R_ext_TI>>
         MIX 1 2 1
  ;
  ::: C9 := GEO: CARCEL 0
         MESHX 0.0 <<Lame>>
         MESHY 0.0 <<Cote>>
         MIX 162
  ;
  ::: C10 := GEO: CARCEL 0
         MESHX 0.0 <<Lame>>
         MESHY 0.0 <<Lame>>
         MIX 162
  ;
;

END: ;
QUIT .