summaryrefslogtreecommitdiff
path: root/Dragon/data/salmacro_proc/TDCM48.c2m
blob: ab2227aa7106066fefdf17983f9a39a40e4669c6 (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
155
156
157
158
159
160
161
162
163
164
165
*----
*  TEST CASE TDCM48
*  MACROSCOPIC CROSS SECTIONS
*  SIMPLE ASSEMBLY BENCHMARK WITH A WINDMILL DISCRETIZATION
*----
*  Define STRUCTURES and MODULES used
*----
LINKED_LIST GEON2 TRACK MACRO PIJ FLUX ;
SEQ_BINARY TF_EXC ;
SEQ_ASCII UOX_TBH :: FILE './GEOtst_g2s.tdt' ;
SEQ_ASCII GEOtstPS :: FILE './GEOtst2_g2s.ps' ;
MODULE GEO: G2S: MAC: ASM: FLU: END: ;
PROCEDURE geomProc ;
PROCEDURE assertS ;
*
*------------------------------------*                                   
* BEGIN GEOMETRY FOR FLUX CALCULATION*                                   
* -----------------------------------*                                   
REAL Rcomb1 := 0.288712 ;                                          
REAL Rcomb2 := 0.365195 ;                                          
REAL Rcomb3 := 0.397962 ;                                          
REAL Rcomb4 := 0.4083 ;                                            
REAL R_int_G := 0.4165 ;                                           
REAL R_ext_G := 0.4775 ;                                           
REAL Cote := 1.26 ;                                                
REAL Lame := 0.042 ;
REAL mesh1 := Cote ;
REAL mesh2 := mesh1 Cote + ;
REAL mesh3 := mesh2 Cote + ;
REAL mesh4 := mesh3 Lame + ;
*
GEON2 := GEO: :: CAR2D 4 4
  EDIT 0
  X- DIAG X+ REFL
  Y- SYME Y+ DIAG
  CELL
  C0101 C0201 C0301  Lame_V
        C0202 C0302  Lame_V
              C0303  Lame_V
                     Lame_C
  MESHX 0.0 <<mesh1>> <<mesh2>> <<mesh3>> <<mesh4>>
  MESHY 0.0 <<mesh1>> <<mesh2>> <<mesh3>> <<mesh4>>
  ::: Lame_C := GEO: CAR2D 1 1
    MESHX 0.0 <<Lame>> MESHY 0.0 <<Lame>>
    MIX 7 ;

  ::: Lame_V := GEO: CAR2D 1 1
    MESHX 0.0 <<Lame>> MESHY 0.0 <<Cote>>
    SPLITY 3 MIX 7 ;

  ::: C0101 := GEO: CARCEL 6
         MESHX 0.0 <<Cote>>
         MESHY 0.0 <<Cote>>
         SECT 4 6
         RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>>
                <<R_int_G>> <<R_ext_G>>
         MIX  1 2 3 4   5 6 7 7 7 7 7 7 7 7 7 7 7 7 ;

  ::: C0201 := GEO: C0101
         MIX  8 9 10 11   5 6 7 7 7 7 7 7 7 7 7 7 7 7 ;

  ::: C0202 := GEO: C0101
         MIX  12 13 14 15   5 6 7 7 7 7 7 7 7 7 7 7 7 7 ;

  ::: C0301 := GEO: C0101
         MIX  16 17 18 19   5 6 7 7 7 7 7 7 7 7 7 7 7 7 ;

  ::: C0302 := GEO: C0101
         MIX  20 21 22 23   5 6 7 7 7 7 7 7 7 7 7 7 7 7 ;

  ::: C0303 := GEO: C0101
         MIX  24 25 26 27   5 6 7 7 7 7 7 7 7 7 7 7 7 7 ;
;
UOX_TBH GEOtstPS := G2S: GEON2 :: DRAWMIX ;
                                                 
MACRO := MAC: ::
     EDIT 2 NGRO 1 NMIX 27 NIFI 1
     READ INPUT
     MIX 1 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 2 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 3 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 4 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 5 (*VOID*)
        TOTAL 0.001     SCAT 1 1 0.0005
     MIX 6 (*CLAD*)
        TOTAL 0.4029    SCAT 1 1 0.4000
     MIX 7 (*WATER*)
        TOTAL 0.3683    SCAT 1 1 0.35
     MIX 8 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 9 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 10 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 11 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 12 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 13 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 14 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 15 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 16 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 17 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 18 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 19 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 20 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 21 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 22 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 23 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 24 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 25 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 26 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
     MIX 27 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.07664  CHI 1.0
;
* 
* flux calculation
*
FLUX := geomProc UOX_TBH MACRO :: 8 0.01 ;
assertS FLUX :: 'K-EFFECTIVE' 1 1.031027 ;
ECHO "test TDCM48 completed" ;
END: ;
QUIT "LIST" .