blob: b1dd7f5e0d9a387da2751316d324ac6c41ebac38 (
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
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
|
*----
* TEST CASE TCM53
* MACROSCOPIC CROSS SECTIONS
* MONTE-CARLO MULTIGROUP METHOD
* CANDU-NG UNIT CELL WITH NXT:
*
*----
* Define STRUCTURES and MODULES used
*----
MODULE
NXT: MAC: DELETE: MC: GREP: ABORT: SYBILT: ASM: FLU: GEO: END:
;
REAL VALUE REFVALUE DELTA ;
PROCEDURE
CHB_NXT_MC
;
LINKED_LIST
GEOMF
MACRO
TRAC2
OUTMC
MACRO2 GEOM2 TRACK2 SYS FLUX
;
SEQ_ASCII
plot.ps
;
STRING
state SolveurSS SolveurFL polquad
Library
;
INTEGER
splitc splitm nban2d nban3d nban2dss
;
REAL
density densityss
;
INTEGER
CF11 CF12 CF13 CF14
CF21i CF22i CF23i CF24i
CF31i CF32i CF33i CF34i
CF41i CF42i CF43i CF44i
CF21e CF22e CF23e CF24e
CF31e CF32e CF33e CF34e
CF41e CF42e CF43e CF44e
VF11 VF12 VF13 VF14
VF21i VF22i VF23i VF24i
VF31i VF32i VF33i VF34i
VF41i VF42i VF43i VF44i
VF21e VF22e VF23e VF24e
VF31e VF32e VF33e VF34e
VF41e VF42e VF43e VF44e
COOL PT HEG CT
MOD SH VOID :=
1 1 1 1
2 2 2 2
3 3 3 3
4 4 4 4
2 2 2 2
3 3 3 3
4 4 4 4
1 1 1 1
2 2 2 2
3 3 3 3
4 4 4 4
2 2 2 2
3 3 3 3
4 4 4 4
6 6 5 5
6 5 5
;
*
* geometrical parameters
*
EVALUATE state := "Voided" ; ! Cooled/Voided/Chboard
EVALUATE splitc := 1 ;
EVALUATE splitm := 1 ;
EVALUATE nban2dss := 8 ;
EVALUATE densityss := 20.0 ;
EVALUATE nban2d := 12 ;
EVALUATE density := 50.0 ;
EVALUATE polquad := "LCMD" ;
EVALUATE nban3d := 2 ;
*
* geometries
*
GEOMF := CHB_NXT_MC ::
<<state>>
<<splitc>> <<splitm>>
<<CF11>> <<CF12>> <<CF13>> <<CF14>>
<<CF21i>> <<CF22i>> <<CF23i>> <<CF24i>>
<<CF31i>> <<CF32i>> <<CF33i>> <<CF34i>>
<<CF41i>> <<CF42i>> <<CF43i>> <<CF44i>>
<<CF21e>> <<CF22e>> <<CF23e>> <<CF24e>>
<<CF31e>> <<CF32e>> <<CF33e>> <<CF34e>>
<<CF41e>> <<CF42e>> <<CF43e>> <<CF44e>>
<<VF11>> <<VF12>> <<VF13>> <<VF14>>
<<VF21i>> <<VF22i>> <<VF23i>> <<VF24i>>
<<VF31i>> <<VF32i>> <<VF33i>> <<VF34i>>
<<VF41i>> <<VF42i>> <<VF43i>> <<VF44i>>
<<VF21e>> <<VF22e>> <<VF23e>> <<VF24e>>
<<VF31e>> <<VF32e>> <<VF33e>> <<VF34e>>
<<VF41e>> <<VF42e>> <<VF43e>> <<VF44e>>
<<COOL>> <<PT>> <<HEG>> <<CT>>
<<MOD>> <<SH>> <<VOID>>
;
*
* geometry tracking
*
TRAC2 := NXT: GEOMF ::
EDIT 1 ANIS 1 MC PSPC
;
*
* macrolib
*
MACRO := MAC: ::
EDIT 0 NGRO 2 NMIX 6 NIFI 1
READ INPUT
MIX 1
TOTAL 2.222E-01 8.3333E-01
NUSIGF 0.000E+00 2.3600E-01
CHI 1.000E+00 0.0000E+00
SCAT 1 1 0.1921 2 2 0.7533 0.02
MIX 2
TOTAL 2.222E-01 8.3333E-01
NUSIGF 0.000E+00 2.3600E-01
CHI 1.000E+00 0.0000E+00
SCAT 1 1 0.1921 2 2 0.7483 0.02
MIX 3
TOTAL 2.222E-01 7.5333E-01
NUSIGF 0.000E+00 2.36000E-01
CHI 1.000E+00 0.00000E+00
SCAT 1 1 0.1921 2 2 0.7033 0.02
MIX 4
TOTAL 2.322E-01 7.3333E-01
NUSIGF 0.000E+00 2.36000E-01
CHI 1.000E+00 0.00000E+00
SCAT 1 1 0.2121 2 2 0.7033 0.02
MIX 5
TOTAL 0.0 0.0
SCAT 1 1 0.0 2 2 0.0 0.0
MIX 6
TOTAL 0.1677 1.1211
SCAT 1 1 0.1265 2 2 1.1011 0.04
;
*
* Monte-Carlo calculation
*
OUTMC := MC: TRAC2 MACRO ::
KCODE 300 400 6000
SEED 23412
EDIT 1
TALLY
MERG COMP
COND NONE
ENDT
;
MACRO2 := OUTMC :: STEP UP 'MACROLIB' ;
GEOM2 := GEO: :: CARCEL 0
X- REFL X+ REFL Y- REFL Y+ REFL
MIX 1
MESHX 0.0 1.0
MESHY 0.0 1.0
;
TRACK2 := SYBILT: GEOM2 ::
TITLE 'INFINITE MEDIUM'
EDIT 1 MAXR 5 QUA1 5 ;
SYS := ASM: MACRO2 TRACK2 ;
FLUX := FLU: SYS MACRO2 TRACK2 :: TYPE K ;
*
* Non regression assert
*
EVALUATE REFVALUE := 1.0088 ;
*
GREP: FLUX :: GETVAL K-EFFECTIVE 1 >>VALUE<< ;
EVALUATE DELTA := VALUE REFVALUE - REFVALUE / ABS ;
IF DELTA 5.0E-3 < THEN
PRINT "TEST SUCCESSFUL; DELTA=" DELTA ;
ELSE
PRINT "------------" ;
PRINT "TEST FAILURE" ;
PRINT "------------" ;
PRINT "REFERENCE=" REFVALUE " CALCULATED=" VALUE ;
ABORT: ;
ENDIF ;
ECHO "test TCM53 completed" ;
END: ;
QUIT "LIST" .
|
