summaryrefslogtreecommitdiff
path: root/Dragon/data/salmacro_proc/TDCM46.c2m
blob: cbcc2c8a8fed991c0611f6886205177323869f9d (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
*----
* TEST CASE TDCM46
* MACROSCOPIC CROSS SECTIONS
* ATBR 19 element ThO2 Cluster
*----
*  Define STRUCTURES and MODULES used
*----
LINKED_LIST MACRO ATBRTHS VOLMATF CP PIJ FLUX EDITION ;
SEQ_BINARY Track ;
MODULE MAC: GEO: G2S: SALT: ASM: FLU: GREP: EDI: END: ;
SEQ_ASCII THO2_DAT :: FILE './THO2_DAT_g2s.dat' ;
SEQ_ASCII Fil_SAL.ps :: FILE './tho2_fig.ps' ;
REAL Kinf ;
INTEGER an2d := 8 ;
REAL densur := 20.0 ;
PROCEDURE assertS ;
*----
* Macroscopic XS
*----
MACRO := MAC: ::
     EDIT 2 NGRO 1 NMIX 8 NIFI 1
     READ INPUT
     MIX 1 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.06564      CHI 1.0
     MIX 2 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.06564      CHI 1.0
     MIX 3 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.06564      CHI 1.0
     MIX 4 (*FUEL*)
        TOTAL 0.3652    SCAT 1 1 0.3234
        NUSIGF 0.06564      CHI 1.0
     MIX 5 (*GAINE*)
        TOTAL 0.4029        SCAT 1 1 0.4000
     MIX 6 (*H2O*)
        TOTAL 0.3683        SCAT 1 1 0.3661
     MIX 7 (*MODERATOR*)
        TOTAL 0.03683       SCAT 1 1 0.03661
     MIX 8 (*VOID*)
        TOTAL 0.001         SCAT 1 1 0.0
;
*----
* Geometry ATBRTHS : ThO2 Cluster
*----
ATBRTHS := GEO: :: HEX 1
  HBC COMPLETE ALBE 1.0
  CELL  C1 
  ::: C1 := GEO: HEXCEL 9
  SIDE 16.16581 
  RADIUS  0.00 0.50 0.54 7.00 7.15 8.80 9.35 10.20 10.35  10.50 
  MIX      6    5    5    5    6    5    8    5    7     7
  CLUSTER ROD1
  ::: ROD1 := GEO: TUBE 5 MIX 1 2 3 4 5
    NPIN 30 RPIN 8.0 APIN 0.0
    RADIUS 0.0 0.44548 0.56349 0.61405 0.63 0.7 ;
  ;
;
THO2_DAT Fil_SAL.ps := G2S: ATBRTHS :: DRAWMIX ;
VOLMATF Track := SALT: THO2_DAT ::                                       
     EDIT 3
     ALLG
     TISO <<an2d>> <<densur>>
;                                                                      
*----
* Compute  flux
*----
PIJ := ASM: MACRO VOLMATF Track ;
FLUX := FLU: PIJ MACRO VOLMATF ::
  TYPE K ;
GREP: FLUX :: GETVAL 'K-INFINITY  ' 1 1 1 >>Kinf<< ;
ECHO "@@@ Kinf=" Kinf ;
assertS FLUX :: 'K-EFFECTIVE' 1 1.040651 ;
EDITION := EDI: FLUX MACRO VOLMATF ::
  EDIT 3 MERG COMP SAVE ON TEST  ;
ECHO "test TDCM46 completed" ;
END: ;
QUIT "LIST" .