blob: f61d833fe93424997104ce3affc3ced4ed24d644 (
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
|
*----
* case ASSBLY_CASEA_1level_mpo
* 99-group eighth assembly in DRAGON -- 1-level scheme
* production of an MPO file for reflector feeding
* Author: A. Hebert
*----
HDF5_FILE MyMpo :: FILE './MyMpo.hdf' ;
LINKED_LIST GEOSS GEON2 TRACKSS TRACKN2 LIBRARY LIBRARY2 SYS FLUX2 EDIOBJ ;
MODULE SYBILT: G2S: SALT: MCCGT: USS: ASM: FLU: EDI: UTL: DELETE: LIB: MPO: END: ;
REAL 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 Lame ;
SEQ_ASCII ASSBLY_CASEA :: FILE './ASSBLY_CASEA_1level_g2s.txt' ;
SEQ_ASCII FIGASB_CASEA :: FILE './ASSBLY_CASEA_1level_g2s.eps' ;
SEQ_BINARY TF_EXC ;
INTEGER an2d := 8 ;
REAL densur := 20.0 ;
INTEGER istep maxstep istepNext maxautop iAutop := 0 77 0 8 1 ;
PROCEDURE Geo_SS Geo_N2 Mix_UOX ;
PROCEDURE assertS ;
************************************************************************
* BEGIN GEOMETRY *
************************************************************************
EVALUATE Rcomb4 := 0.4096 ;
EVALUATE Rcomb1 := 0.5 SQRT Rcomb4 * ;
EVALUATE Rcomb2 := 0.8 SQRT Rcomb4 * ;
EVALUATE Rcomb3 := 0.95 SQRT Rcomb4 * ;
EVALUATE R_int_TG := 0.5605 ;
EVALUATE R_ext_TG := 0.6225 ;
EVALUATE R_int_TI := 0.5725 ;
EVALUATE R_ext_TI := 0.6225 ;
EVALUATE R_int_G := 0.4180 ;
EVALUATE R_ext_G := 0.4750 ;
EVALUATE Cote := 1.26 ;
EVALUATE CoteL := 1.302 ;
EVALUATE Lame := CoteL Cote - ;
EVALUATE RmodeN1 := 0.670 ;
GEOSS := Geo_SS :: <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>>
<<R_int_TG>> <<R_ext_TG>> <<R_int_TI>>
<<R_ext_TI>> <<R_int_G>> <<R_ext_G>> <<Cote>>
<<CoteL>> ;
GEON2 := Geo_N2 :: <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>>
<<R_int_TG>> <<R_ext_TG>> <<R_int_TI>>
<<R_ext_TI>> <<R_int_G>> <<R_ext_G>> <<Cote>>
<<Lame>> ;
ASSBLY_CASEA FIGASB_CASEA := G2S: GEON2 ;
************************************************************************
* BEGIN TRACKING *
************************************************************************
! Level One
TRACKSS := SYBILT: GEOSS ::
EDIT 0
MAXR 500
MAXZ 1000000
TITLE 'TRACKING FOR ASSEMBLY SS'
QUA2 20 3
DP01 ;
! MOC tracking for flux calculation
TRACKN2 TF_EXC := SALT: ASSBLY_CASEA ::
EDIT 3
ALLG
TSPC EQW2 <<an2d>> <<densur>> REND
;
TRACKN2 := MCCGT: TRACKN2 TF_EXC ::
CACB 4 AAC 80 TMT EPSI 1E-5 MCU 10000 ;
************************************************************************
* END TRACKING *
************************************************************************
************************************************************************
* BEGIN LIBRARY READING *
************************************************************************
STRING Library := 'CLA99CEA93:CLA99CEA93_SS' ;
REAL dens_mod := 0.85 ;
REAL pbore := 975.0 ;
LIBRARY := Mix_UOX :: <<Library>> <<dens_mod>> <<pbore>> ;
************************************************************************
* END LIBRARY READING *
************************************************************************
************************************************************************
* BEGIN MyMpo INITIALIZATION *
************************************************************************
MyMpo := MPO: LIBRARY ::
EDIT 10
COMM '99-group feeding assembly'
REAC Total Absorption Fission FissionSpectrum NuFission CaptureEnergyCapture
FissionEnergyFission Diffusion Scattering
;
************************************************************************
* END MyMpo INITIALIZATION *
************************************************************************
************************************************************************
* RESONANCE SELF-SHIELDING *
************************************************************************
LIBRARY2 := USS: LIBRARY TRACKSS :: EDIT 1 PASS 2 ARM GRMIN 18
CALC REGI W1 U235 ALL
REGI W1 U238 3
REGI W2 U238 4
REGI W3 U238 5
REGI W4 U238 6
REGI W1 U238 9
REGI W2 U238 10
REGI W3 U238 11
REGI W4 U238 12
REGI W1 U238 94
REGI W2 U238 95
REGI W3 U238 96
REGI W4 U238 97
REGI W1 U238 114
REGI W2 U238 115
REGI W3 U238 116
REGI W4 U238 117
REGI W1 U238 122
REGI W2 U238 123
REGI W3 U238 124
REGI W4 U238 125
REGI W1 U238 126
REGI W2 U238 127
REGI W3 U238 128
REGI W4 U238 129
REGI W1 U238 154
REGI W2 U238 155
REGI W3 U238 156
REGI W4 U238 157
REGI W1 U238 158
REGI W2 U238 159
REGI W3 U238 160
REGI W4 U238 161
ENDC ;
************************************************************************
* BEGIN MOC FLUX CALCULATION *
************************************************************************
SYS := ASM: LIBRARY2 TRACKN2 TF_EXC :: ARM EDIT 1 ;
FLUX2 := FLU: LIBRARY2 SYS TRACKN2 TF_EXC :: EDIT 1 TYPE B B1 SIGS ;
EDIOBJ := EDI: FLUX2 LIBRARY2 TRACKN2 :: EDIT 5
MERG COMP
MICR RES
SAVE ON COND2
;
************************************************************************
* MyMpo FEEDING *
************************************************************************
MyMpo := MPO: MyMpo EDIOBJ :: EDIT 2 ;
assertS FLUX2 :: K-INFINITY 1 1.018207 ;
ECHO "test ASSBLY_CASEA_1level_mpo completed" ;
END: ;
|
