summaryrefslogtreecommitdiff
path: root/Dragon/data/OSC_CASEA_2level_rse_proc/Geo_N1.c2m
blob: 38b61462abb5c027a97f99abe2deb3f076fdd0ec (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
*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 .