summaryrefslogtreecommitdiff
path: root/Dragon/data/salmacro_proc/TDCM33.c2m
blob: aac7a89d3021118f84d4a70465d7f03fa2226273 (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
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
*----
*  TEST CASE TDCM33
*  MACROSCOPIC CROSS SECTIONS
*  PULSE REACTOR WITH MOC
*----
*  Define STRUCTURES and MODULES used
*----
LINKED_LIST GEOM TRACK MACRO PIJ FLUX EDITION ;
MODULE GEO: G2S: SALT: MCCGT: MAC: ASM: FLU: EDI: GREP: END: ;
SEQ_BINARY TRACK_FIL ;
SEQ_ASCII Fic_SAL :: FILE './Fic_SAL.tdt' ;
SEQ_ASCII Fic_SAL.eps :: FILE './Fic_SAL.eps' ;
REAL Kinf ;
INTEGER an2d := 6 ;
REAL densur := 15.0 ;
PROCEDURE assertS ;
*
* evaluate pin positions in triangular assemblies
REAL x y := 12.6 13.91 ;
REAL rpin1 := x x * y y * + SQRT ;
REAL apin1 := y x / ARCTAN ;
EVALUATE y := y 1.26 - ;
REAL rpin2 := x x * y y * + SQRT ;
REAL apin2 := y x / ARCTAN ;
EVALUATE y := y 1.26 - ;
REAL rpin3 := x x * y y * + SQRT ;
REAL apin3 := y x / ARCTAN ;
EVALUATE y := y 1.26 - ;
REAL rpin4 := x x * y y * + SQRT ;
REAL apin4 := y x / ARCTAN ;
EVALUATE y := y 1.26 - ;
REAL rpin5 := x x * y y * + SQRT ;
REAL apin5 := y x / ARCTAN ;
EVALUATE y := y 1.26 - ;
REAL rpin6 := x x * y y * + SQRT ;
REAL apin6 := y x / ARCTAN ;
EVALUATE x y := 11.34 13.91 ;
REAL rpin7 := x x * y y * + SQRT ;
REAL apin7 := y x / ARCTAN ;
EVALUATE y := y 1.26 - ;
REAL rpin8 := x x * y y * + SQRT ;
REAL apin8 := y x / ARCTAN ;
EVALUATE y := y 1.26 - ;
REAL rpin9 := x x * y y * + SQRT ;
REAL apin9 := y x / ARCTAN ;
EVALUATE y := y 1.26 - ;
REAL rpin10 := x x * y y * + SQRT ;
REAL apin10 := y x / ARCTAN ;
EVALUATE y := y 1.26 - ;
REAL rpin11 := x x * y y * + SQRT ;
REAL apin11 := y x / ARCTAN ;
EVALUATE x y := 10.08 13.91 ;
REAL rpin12 := x x * y y * + SQRT ;
REAL apin12 := y x / ARCTAN ;
EVALUATE y := y 1.26 - ;
REAL rpin13 := x x * y y * + SQRT ;
REAL apin13 := y x / ARCTAN ;
EVALUATE y := y 1.26 - ;
REAL rpin14 := x x * y y * + SQRT ;
REAL apin14 := y x / ARCTAN ;
EVALUATE y := y 1.26 - ;
REAL rpin15 := x x * y y * + SQRT ;
REAL apin15 := y x / ARCTAN ;
EVALUATE x y := 8.82 13.91 ;
REAL rpin16 := x x * y y * + SQRT ;
REAL apin16 := y x / ARCTAN ;
EVALUATE y := y 1.26 - ;
REAL rpin17 := x x * y y * + SQRT ;
REAL apin17 := y x / ARCTAN ;
EVALUATE y := y 1.26 - ;
REAL rpin18 := x x * y y * + SQRT ;
REAL apin18 := y x / ARCTAN ;
EVALUATE x y := 7.56 13.91 ;
REAL rpin19 := x x * y y * + SQRT ;
REAL apin19 := y x / ARCTAN ;
EVALUATE y := y 1.26 - ;
REAL rpin20 := x x * y y * + SQRT ;
REAL apin20 := y x / ARCTAN ;
EVALUATE x y := 6.3 13.91 ;
REAL rpin21 := x x * y y * + SQRT ;
REAL apin21 := y x / ARCTAN ;
*
REAL apin22 := $Pi_R apin1 + CHS ;
REAL apin23 := $Pi_R apin2 + CHS ;
REAL apin24 := $Pi_R apin3 + CHS ;
REAL apin25 := $Pi_R apin4 + CHS ;
REAL apin26 := $Pi_R apin5 + CHS ;
REAL apin27 := $Pi_R apin6 + CHS ;
REAL apin28 := $Pi_R apin7 + CHS ;
REAL apin29 := $Pi_R apin8 + CHS ;
REAL apin30 := $Pi_R apin9 + CHS ;
REAL apin31 := $Pi_R apin10 + CHS ;
REAL apin32 := $Pi_R apin11 + CHS ;
REAL apin33 := $Pi_R apin12 + CHS ;
REAL apin34 := $Pi_R apin13 + CHS ;
REAL apin35 := $Pi_R apin14 + CHS ;
REAL apin36 := $Pi_R apin15 + CHS ;
REAL apin37 := $Pi_R apin16 + CHS ;
REAL apin38 := $Pi_R apin17 + CHS ;
REAL apin39 := $Pi_R apin18 + CHS ;
REAL apin40 := $Pi_R apin19 + CHS ;
REAL apin41 := $Pi_R apin20 + CHS ;
REAL apin42 := $Pi_R apin21 + CHS ;
*
REAL apin43 := $Pi_R apin1 + ;
REAL apin44 := $Pi_R apin2 + ;
REAL apin45 := $Pi_R apin3 + ;
REAL apin46 := $Pi_R apin4 + ;
REAL apin47 := $Pi_R apin5 + ;
REAL apin48 := $Pi_R apin6 + ;
REAL apin49 := $Pi_R apin7 + ;
REAL apin50 := $Pi_R apin8 + ;
REAL apin51 := $Pi_R apin9 + ;
REAL apin52 := $Pi_R apin10 + ;
REAL apin53 := $Pi_R apin11 + ;
REAL apin54 := $Pi_R apin12 + ;
REAL apin55 := $Pi_R apin13 + ;
REAL apin56 := $Pi_R apin14 + ;
REAL apin57 := $Pi_R apin15 + ;
REAL apin58 := $Pi_R apin16 + ;
REAL apin59 := $Pi_R apin17 + ;
REAL apin60 := $Pi_R apin18 + ;
REAL apin61 := $Pi_R apin19 + ;
REAL apin62 := $Pi_R apin20 + ;
REAL apin63 := $Pi_R apin21 + ;
*
REAL apin64 := apin1 CHS ;
REAL apin65 := apin2 CHS ;
REAL apin66 := apin3 CHS ;
REAL apin67 := apin4 CHS ;
REAL apin68 := apin5 CHS ;
REAL apin69 := apin6 CHS ;
REAL apin70 := apin7 CHS ;
REAL apin71 := apin8 CHS ;
REAL apin72 := apin9 CHS ;
REAL apin73 := apin10 CHS ;
REAL apin74 := apin11 CHS ;
REAL apin75 := apin12 CHS ;
REAL apin76 := apin13 CHS ;
REAL apin77 := apin14 CHS ;
REAL apin78 := apin15 CHS ;
REAL apin79 := apin16 CHS ;
REAL apin80 := apin17 CHS ;
REAL apin81 := apin18 CHS ;
REAL apin82 := apin19 CHS ;
REAL apin83 := apin20 CHS ;
REAL apin84 := apin21 CHS ;
*
REAL Rtest4 := 0.4083 ;
REAL Rtest1 := 0.5 SQRT Rtest4 * ;
REAL Rtest2 := 0.8 SQRT Rtest4 * ;
REAL Rtest3 := 0.95 SQRT Rtest4 * ;
REAL Rcomb4 := 0.46 ;
REAL Rcomb1 := 0.5 SQRT Rcomb4 * ;
REAL Rcomb2 := 0.8 SQRT Rcomb4 * ;
REAL Rcomb3 := 0.95 SQRT Rcomb4 * ;
*
GEOM := GEO: :: CARCEL 19
      X- VOID X+ VOID Y- VOID Y+ VOID
      RADIUS 0.0 <<Rtest1>> <<Rtest2>> <<Rtest3>> <<Rtest4>> 0.475 0.75
             1.55 2.4 3.92 4.1 6.1 7.7 8.47 12.5 12.9 13.75 14.9 16.2 17.55
      MESHX 0.0 26.46
      MESHY 0.0 29.08
      SPLITX 6 SPLITY 6
      MIX 7 8 9 10 11 12 13 14 15 16 17 18 19 7 20 7 7 7 7 7
      CLUSTER PIN1 PIN2 PIN3 PIN4 PIN5 PIN6 PIN7 PIN8 PIN9 PIN10 PIN11
              PIN12 PIN13 PIN14 PIN15 PIN16 PIN17 PIN18 PIN19 PIN20
              PIN21 PIN22 PIN23 PIN24 PIN25 PIN26 PIN27 PIN28 PIN29
              PIN30 PIN31 PIN32 PIN33 PIN34 PIN35 PIN36 PIN37 PIN38
              PIN39 PIN40 PIN41 PIN42 PIN43 PIN44 PIN45 PIN46 PIN47
              PIN48 PIN49 PIN50 PIN51 PIN52 PIN53 PIN54 PIN55 PIN56
              PIN57 PIN58 PIN59 PIN60 PIN61 PIN62 PIN63 PIN64 PIN65
              PIN66 PIN67 PIN68 PIN69 PIN70 PIN71 PIN72 PIN73 PIN74
              PIN75 PIN76 PIN77 PIN78 PIN79 PIN80 PIN81 PIN82 PIN83
              PIN84
      ::: PIN1 := GEO: TUBE 5
              RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> 0.5
              MIX 1 2 3 4 5
              NPIN 1 RPIN <<rpin1>> APIN <<apin1>>
      ;
      ::: PIN2 := GEO: PIN1
              NPIN 1 RPIN <<rpin2>> APIN <<apin2>>
      ;
      ::: PIN3 := GEO: PIN1
              NPIN 1 RPIN <<rpin3>> APIN <<apin3>>
      ;
      ::: PIN4 := GEO: PIN1
              NPIN 1 RPIN <<rpin4>> APIN <<apin4>>
      ;
      ::: PIN5 := GEO: PIN1
              NPIN 1 RPIN <<rpin5>> APIN <<apin5>>
      ;
      ::: PIN6 := GEO: PIN1
              NPIN 1 RPIN <<rpin6>> APIN <<apin6>>
      ;
      ::: PIN7 := GEO: PIN1
              NPIN 1 RPIN <<rpin7>> APIN <<apin7>>
      ;
      ::: PIN8 := GEO: TUBE 1
              RADIUS 0.0 0.5
              MIX 21
              NPIN 1 RPIN <<rpin8>> APIN <<apin8>>
      ;
      ::: PIN9 := GEO: PIN1
              NPIN 1 RPIN <<rpin9>> APIN <<apin9>>
      ;
      ::: PIN10 := GEO: PIN8
              NPIN 1 RPIN <<rpin10>> APIN <<apin10>>
      ;
      ::: PIN11 := GEO: PIN1
              NPIN 1 RPIN <<rpin11>> APIN <<apin11>>
      ;
      ::: PIN12 := GEO: PIN1
              NPIN 1 RPIN <<rpin12>> APIN <<apin12>>
      ;
      ::: PIN13 := GEO: PIN1
              NPIN 1 RPIN <<rpin13>> APIN <<apin13>>
      ;
      ::: PIN14 := GEO: PIN1
              NPIN 1 RPIN <<rpin14>> APIN <<apin14>>
      ;
      ::: PIN15 := GEO: PIN1
              NPIN 1 RPIN <<rpin15>> APIN <<apin15>>
      ;
      ::: PIN16 := GEO: PIN1
              NPIN 1 RPIN <<rpin16>> APIN <<apin16>>
      ;
      ::: PIN17 := GEO: PIN8
              NPIN 1 RPIN <<rpin17>> APIN <<apin17>>
      ;
      ::: PIN18 := GEO: PIN1
              NPIN 1 RPIN <<rpin18>> APIN <<apin18>>
      ;
      ::: PIN19 := GEO: PIN1
              NPIN 1 RPIN <<rpin19>> APIN <<apin19>>
      ;
      ::: PIN20 := GEO: PIN1
              NPIN 1 RPIN <<rpin20>> APIN <<apin20>>
      ;
      ::: PIN21 := GEO: PIN1
              NPIN 1 RPIN <<rpin21>> APIN <<apin21>>
      ;
*
      ::: PIN22 := GEO: PIN1
              NPIN 1 RPIN <<rpin1>> APIN <<apin22>>
      ;
      ::: PIN23 := GEO: PIN1
              NPIN 1 RPIN <<rpin2>> APIN <<apin23>>
      ;
      ::: PIN24 := GEO: PIN1
              NPIN 1 RPIN <<rpin3>> APIN <<apin24>>
      ;
      ::: PIN25 := GEO: PIN1
              NPIN 1 RPIN <<rpin4>> APIN <<apin25>>
      ;
      ::: PIN26 := GEO: PIN1
              NPIN 1 RPIN <<rpin5>> APIN <<apin26>>
      ;
      ::: PIN27 := GEO: PIN1
              NPIN 1 RPIN <<rpin6>> APIN <<apin27>>
      ;
      ::: PIN28 := GEO: PIN1
              NPIN 1 RPIN <<rpin7>> APIN <<apin28>>
      ;
      ::: PIN29 := GEO: PIN8
              NPIN 1 RPIN <<rpin8>> APIN <<apin29>>
      ;
      ::: PIN30 := GEO: PIN1
              NPIN 1 RPIN <<rpin9>> APIN <<apin30>>
      ;
      ::: PIN31 := GEO: PIN8
              NPIN 1 RPIN <<rpin10>> APIN <<apin31>>
      ;
      ::: PIN32 := GEO: PIN1
              NPIN 1 RPIN <<rpin11>> APIN <<apin32>>
      ;
      ::: PIN33 := GEO: PIN1
              NPIN 1 RPIN <<rpin12>> APIN <<apin33>>
      ;
      ::: PIN34 := GEO: PIN1
              NPIN 1 RPIN <<rpin13>> APIN <<apin34>>
      ;
      ::: PIN35 := GEO: PIN1
              NPIN 1 RPIN <<rpin14>> APIN <<apin35>>
      ;
      ::: PIN36 := GEO: PIN1
              NPIN 1 RPIN <<rpin15>> APIN <<apin36>>
      ;
      ::: PIN37 := GEO: PIN1
              NPIN 1 RPIN <<rpin16>> APIN <<apin37>>
      ;
      ::: PIN38 := GEO: PIN8
              NPIN 1 RPIN <<rpin17>> APIN <<apin38>>
      ;
      ::: PIN39 := GEO: PIN1
              NPIN 1 RPIN <<rpin18>> APIN <<apin39>>
      ;
      ::: PIN40 := GEO: PIN1
              NPIN 1 RPIN <<rpin19>> APIN <<apin40>>
      ;
      ::: PIN41 := GEO: PIN1
              NPIN 1 RPIN <<rpin20>> APIN <<apin41>>
      ;
      ::: PIN42 := GEO: PIN1
              NPIN 1 RPIN <<rpin21>> APIN <<apin42>>
      ;
*
      ::: PIN43 := GEO: PIN1
              NPIN 1 RPIN <<rpin1>> APIN <<apin43>>
      ;
      ::: PIN44 := GEO: PIN1
              NPIN 1 RPIN <<rpin2>> APIN <<apin44>>
      ;
      ::: PIN45 := GEO: PIN1
              NPIN 1 RPIN <<rpin3>> APIN <<apin45>>
      ;
      ::: PIN46 := GEO: PIN1
              NPIN 1 RPIN <<rpin4>> APIN <<apin46>>
      ;
      ::: PIN47 := GEO: PIN1
              NPIN 1 RPIN <<rpin5>> APIN <<apin47>>
      ;
      ::: PIN48 := GEO: PIN1
              NPIN 1 RPIN <<rpin6>> APIN <<apin48>>
      ;
      ::: PIN49 := GEO: PIN1
              NPIN 1 RPIN <<rpin7>> APIN <<apin49>>
      ;
      ::: PIN50 := GEO: PIN8
              NPIN 1 RPIN <<rpin8>> APIN <<apin50>>
      ;
      ::: PIN51 := GEO: PIN1
              NPIN 1 RPIN <<rpin9>> APIN <<apin51>>
      ;
      ::: PIN52 := GEO: PIN8
              NPIN 1 RPIN <<rpin10>> APIN <<apin52>>
      ;
      ::: PIN53 := GEO: PIN1
              NPIN 1 RPIN <<rpin11>> APIN <<apin53>>
      ;
      ::: PIN54 := GEO: PIN1
              NPIN 1 RPIN <<rpin12>> APIN <<apin54>>
      ;
      ::: PIN55 := GEO: PIN1
              NPIN 1 RPIN <<rpin13>> APIN <<apin55>>
      ;
      ::: PIN56 := GEO: PIN1
              NPIN 1 RPIN <<rpin14>> APIN <<apin56>>
      ;
      ::: PIN57 := GEO: PIN1
              NPIN 1 RPIN <<rpin15>> APIN <<apin57>>
      ;
      ::: PIN58 := GEO: PIN1
              NPIN 1 RPIN <<rpin16>> APIN <<apin58>>
      ;
      ::: PIN59 := GEO: PIN8
              NPIN 1 RPIN <<rpin17>> APIN <<apin59>>
      ;
      ::: PIN60 := GEO: PIN1
              NPIN 1 RPIN <<rpin18>> APIN <<apin60>>
      ;
      ::: PIN61 := GEO: PIN1
              NPIN 1 RPIN <<rpin19>> APIN <<apin61>>
      ;
      ::: PIN62 := GEO: PIN1
              NPIN 1 RPIN <<rpin20>> APIN <<apin62>>
      ;
      ::: PIN63 := GEO: PIN1
              NPIN 1 RPIN <<rpin21>> APIN <<apin63>>
      ;
*
      ::: PIN64 := GEO: PIN1
              NPIN 1 RPIN <<rpin1>> APIN <<apin64>>
      ;
      ::: PIN65 := GEO: PIN1
              NPIN 1 RPIN <<rpin2>> APIN <<apin65>>
      ;
      ::: PIN66 := GEO: PIN1
              NPIN 1 RPIN <<rpin3>> APIN <<apin66>>
      ;
      ::: PIN67 := GEO: PIN1
              NPIN 1 RPIN <<rpin4>> APIN <<apin67>>
      ;
      ::: PIN68 := GEO: PIN1
              NPIN 1 RPIN <<rpin5>> APIN <<apin68>>
      ;
      ::: PIN69 := GEO: PIN1
              NPIN 1 RPIN <<rpin6>> APIN <<apin69>>
      ;
      ::: PIN70 := GEO: PIN1
              NPIN 1 RPIN <<rpin7>> APIN <<apin70>>
      ;
      ::: PIN71 := GEO: PIN8
              NPIN 1 RPIN <<rpin8>> APIN <<apin71>>
      ;
      ::: PIN72 := GEO: PIN1
              NPIN 1 RPIN <<rpin9>> APIN <<apin72>>
      ;
      ::: PIN73 := GEO: PIN8
              NPIN 1 RPIN <<rpin10>> APIN <<apin73>>
      ;
      ::: PIN74 := GEO: PIN1
              NPIN 1 RPIN <<rpin11>> APIN <<apin74>>
      ;
      ::: PIN75 := GEO: PIN1
              NPIN 1 RPIN <<rpin12>> APIN <<apin75>>
      ;
      ::: PIN76 := GEO: PIN1
              NPIN 1 RPIN <<rpin13>> APIN <<apin76>>
      ;
      ::: PIN77 := GEO: PIN1
              NPIN 1 RPIN <<rpin14>> APIN <<apin77>>
      ;
      ::: PIN78 := GEO: PIN1
              NPIN 1 RPIN <<rpin15>> APIN <<apin78>>
      ;
      ::: PIN79 := GEO: PIN1
              NPIN 1 RPIN <<rpin16>> APIN <<apin79>>
      ;
      ::: PIN80 := GEO: PIN8
              NPIN 1 RPIN <<rpin17>> APIN <<apin80>>
      ;
      ::: PIN81 := GEO: PIN1
              NPIN 1 RPIN <<rpin18>> APIN <<apin81>>
      ;
      ::: PIN82 := GEO: PIN1
              NPIN 1 RPIN <<rpin19>> APIN <<apin82>>
      ;
      ::: PIN83 := GEO: PIN1
              NPIN 1 RPIN <<rpin20>> APIN <<apin83>>
      ;
      ::: PIN84 := GEO: PIN1
              NPIN 1 RPIN <<rpin21>> APIN <<apin84>>
      ;
;
Fic_SAL Fic_SAL.eps := G2S: GEOM ;

*----
* Macroscopic XS
*----
MACRO := MAC: ::
     EDIT 2 NGRO 1 NMIX 21 NIFI 1
     READ INPUT
     MIX 1 (*FUEL*)
        TOTAL 0.3652        SCAT 1 1 0.3234
        NUSIGF 0.089        CHI 1.0
     MIX 2 (*FUEL*)
        TOTAL 0.3652        SCAT 1 1 0.3234
        NUSIGF 0.089        CHI 1.0
     MIX 3 (*FUEL*)
        TOTAL 0.3652        SCAT 1 1 0.3234
        NUSIGF 0.089        CHI 1.0
     MIX 4 (*FUEL*)
        TOTAL 0.3652        SCAT 1 1 0.3234
        NUSIGF 0.089        CHI 1.0
     MIX 5 (*CLAD*)
        TOTAL 0.4029        SCAT 1 1 0.4000
     MIX 6 (*H2O*)
        TOTAL 0.3683        SCAT 1 1 0.3661
     MIX 7 (*FUEL*)
        TOTAL 0.3652        SCAT 1 1 0.3234
        NUSIGF 0.089        CHI 1.0
     MIX 8 (*FUEL*)
        TOTAL 0.3652        SCAT 1 1 0.3234
        NUSIGF 0.089        CHI 1.0
     MIX 9 (*FUEL*)
        TOTAL 0.3652        SCAT 1 1 0.3234
        NUSIGF 0.089        CHI 1.0
     MIX 10 (*FUEL*)
        TOTAL 0.3652        SCAT 1 1 0.3234
        NUSIGF 0.089        CHI 1.0
     MIX 11 (*CLAD*)
        TOTAL 0.4029        SCAT 1 1 0.4000
     MIX 12 (*H2O*)
        TOTAL 0.3683        SCAT 1 1 0.3661
     MIX 13 (*TUBE*)
        TOTAL 0.6029        SCAT 1 1 0.6000
     MIX 14 (*GAS*)
        TOTAL 0.0650        SCAT 1 1 0.0600
     MIX 15 (*PRESSURE TUBE*)
        TOTAL 0.4029        SCAT 1 1 0.4000
     MIX 16 (*H2O*)
        TOTAL 0.3683        SCAT 1 1 0.3661
     MIX 17 (*BETA TUBE*)
        TOTAL 0.4029        SCAT 1 1 0.4000
     MIX 18 (*VOID*)
        TOTAL 0.006         SCAT 1 1 0.0060
     MIX 19 (*ALPHA TUBE*)
        TOTAL 0.4029        SCAT 1 1 0.4000
     MIX 20 (*NICKEL TUBE*)
        TOTAL 0.4029        SCAT 1 1 0.4000
     MIX 21 (*BETA TUBE*)
        TOTAL 0.4029        SCAT 1 1 0.4000
;
*----
* Perform tracking
*----
TRACK TRACK_FIL := SALT: Fic_SAL ::
     EDIT 1
     TISO <<an2d>> <<densur>> LONG ;
TRACK := MCCGT: TRACK TRACK_FIL ::
  CACB 4 AAC 80 TMT EPSI 1E-5 MCU 100000 ;
*----
* Compute  flux
*----
PIJ := ASM: MACRO TRACK TRACK_FIL :: ARM ;
FLUX := FLU: PIJ MACRO TRACK TRACK_FIL ::
  TYPE K ;
GREP: FLUX :: GETVAL 'K-INFINITY  ' 1 1 1 >>Kinf<< ;
ECHO "@@@ Kinf=" Kinf ;
EDITION := EDI: FLUX MACRO TRACK Fic_SAL ::
  EDIT 3 MERG G2S 36
    RECT 0.0 4.41 0.0 4.846666667
    RECT 0.0 4.41 4.846666667 9.693333333
    RECT 0.0 4.41 9.693333333 14.54
    RECT 0.0 4.41 14.54 19.386666667
    RECT 0.0 4.41 19.386666667 24.233333333
    RECT 0.0 4.41 24.233333333 29.08

    RECT 4.41 8.82 0.0 4.846666667
    RECT 4.41 8.82 4.846666667 9.693333333
    RECT 4.41 8.82 9.693333333 14.54
    RECT 4.41 8.82 14.54 19.386666667
    RECT 4.41 8.82 19.386666667 24.233333333
    RECT 4.41 8.82 24.233333333 29.08

    RECT 8.82 13.23 0.0 4.846666667
    RECT 8.82 13.23 4.846666667 9.693333333
    RECT 8.82 13.23 9.693333333 14.54
    RECT 8.82 13.23 14.54 19.386666667
    RECT 8.82 13.23 19.386666667 24.233333333
    RECT 8.82 13.23 24.233333333 29.08

    RECT 13.23 17.64 0.0 4.846666667
    RECT 13.23 17.64 4.846666667 9.693333333
    RECT 13.23 17.64 9.693333333 14.54
    RECT 13.23 17.64 14.54 19.386666667
    RECT 13.23 17.64 19.386666667 24.233333333
    RECT 13.23 17.64 24.233333333 29.08

    RECT 17.64 22.05 0.0 4.846666667
    RECT 17.64 22.05 4.846666667 9.693333333
    RECT 17.64 22.05 9.693333333 14.54
    RECT 17.64 22.05 14.54 19.386666667
    RECT 17.64 22.05 19.386666667 24.233333333
    RECT 17.64 22.05 24.233333333 29.08

    RECT 22.05 26.46 0.0 4.846666667
    RECT 22.05 26.46 4.846666667 9.693333333
    RECT 22.05 26.46 9.693333333 14.54
    RECT 22.05 26.46 14.54 19.386666667
    RECT 22.05 26.46 19.386666667 24.233333333
    RECT 22.05 26.46 24.233333333 29.08
  SAVE ON TEST  ;
assertS FLUX :: 'K-EFFECTIVE' 1 1.036426 ;
ECHO "test TDCM33 completed" ;
END: ;
QUIT "LIST" .