summaryrefslogtreecommitdiff
path: root/Dragon/data/proc_LZC/MicWULZC.c2m
blob: 0dcda38511b10c81a09f553b915dee9fb4c5cd15 (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
*DECK MicWULZC.c2m
*----------
*  Name          : MicWULZC.c2m
*  Type          : DRAGON procedure
*  Use           : Add LZC mixtures to Gentilly-2 fuel microlib 
*                  with cross sections from the WLUP based 69 
*                  groups  (LibName = IAEA69).
*                  (cluster model)
*  Author        : G. Marleau
*  Date          : 2006/03/01
*  
*  Procedure description: 
*    MicroLib := MicWULZC MicroLib ::
*      <<PrtLev>> <<Position>> <<Type>> ;
*  
*    MicroLib      : XSM_FILE containing the updated microlib where:
*                    Mix 15 = LZC zone to fill or empty
*                    Mix 16 = Guide tube
*                    Mix 17 = mixture in bubbler for scavenger 1
*                    Mix 18 = mixture around bubbler for scavenger 1
*                    Mix 19 = mixture in scavenger 1
*                    Mix 20 = mixture around scavenger 1
*                    Mix 21 = mixture in bubbler for scavenger 2
*                    Mix 22 = mixture around bubbler for scavenger 2
*                    Mix 23 = mixture in scavenger 2
*                    Mix 24 = mixture around scavenger 2
*                    Mix 25 = mixture in bubbler for scavenger 3
*                    Mix 26 = mixture around bubbler for scavenger 3
*                    Mix 27 = mixture in scavenger 3
*                    Mix 28 = mixture around scavenger 3
*                    Mix 29 = mixture in balance for feeder 1
*                    Mix 30 = mixture around balance for feeder 1
*                    Mix 31 = mixture in feeder 1
*                    Mix 32 = mixture around feeder 1
*                    Mix 33 = mixture in balance for feeder 2
*                    Mix 34 = mixture around balance for feeder 2
*                    Mix 35 = mixture in feeder 2
*                    Mix 36 = mixture around feeder 2
*    PrtLev        : Print Level
*                    0  -> Minimum print
*                    1  -> Standard print
*                    10 -> Debug  
*    Position      : LZC position
*                    NO  -> LZC and guide tubes absent.
*                    OUT -> LZC empty.
*                    IN  -> LZC filled.
*    Type          : LZC type
*                    10  -> 1 scavenger tube 0 feeder tube (default)
*                    21  -> 2 scavenger tubes 1 feeder tube
*                    32  -> 3 scavenger tubes 2 feeder tubes
*
*
* 
*----------
* Define procedure parameters and options and read
* the procedure input data
*----
PARAMETER MicroLib :: EDIT 0 ::: XSM_FILE MicroLib ; ;
INTEGER   PrtLev ;
STRING    Position Type ;
::  >>PrtLev<< >>Position<< >>Type<< ;
*----------
*  Define modules
*----------
MODULE  INFO: LIB:  DELETE:  ;
*----------
*  Library and isotope names
*----------
REAL    Mtemp Mprty := 345.66 99.92 ;
REAL    MwgtH1  MwgtD2  MwgtO16  Mdens ;
REAL    Wtemp Wprty := Mtemp  0.001  ;
REAL    WwgtH1  WwgtD2  WwgtO16  Wdens ;
REAL 
  BNZrNb         BNZr         BNZrII         :=
  0.00021        0.00031      0.00005962     ;  
REAL 
  B10ZrNb        B10Zr        B10ZrII 
  B11ZrNb        B11Zr        B11ZrII        :=
  0.199 BNZrNb * 0.199 BNZr * 0.199 BNZrII * 
  0.801 BNZrNb * 0.801 BNZr * 0.801 BNZrII * ;  
*----------
*  Library and isotope names
*----------
STRING   LibName 
         IsoU5    IsoU8    IsoU6    IsoPu9   IsoO16 
         IsoH1    IsoD2    IsoB10   IsoB11   IsoBn
         IsoZr    IsoFe    IsoFeA   IsoCr    IsoNi 
         IsoNb    IsoC     IsoSi    IsoMn    IsoCu
         IsoMo    IsoAl    IsoHe4   IsoN     IsoTi    
         IsoCd    IsoCo59                          :=
        "IAEA69"               
        "2235"   "8238"   "236"    "6239"   "6016"
        "3001"   "3002"   "10"     "11"     "1011"
        "91"     "2056"   "2056"   "52"     "58"     
        "93"     "2012"   "29"     "55"     "3063" 
        "96"     "27"     "4"      "14"     "48"
        "2113"   "1059"                            ; 
INFO: ::
  TMP: <<Mtemp>> K 
  PUR: <<Mprty>> WGT% 
  CALC DENS WATER >>Mdens<<
  LIB: WIMSD4 FIL: <<LibName>>
  ISO: 3        <<IsoH1>>   <<IsoD2>>   <<IsoO16>> 
  CALC WGT% D2O >>MwgtH1<<  >>MwgtD2<<  >>MwgtO16<<   ;
INFO: ::
  TMP: <<Wtemp>> K 
  PUR: <<Wprty>> WGT% 
  CALC DENS WATER >>Wdens<<
  LIB: WIMSD4 FIL: <<LibName>>
  ISO: 3        <<IsoH1>>   <<IsoD2>>   <<IsoO16>> 
  CALC WGT% D2O >>WwgtH1<<  >>WwgtD2<<  >>WwgtO16<<   ;
*----------
*  Complete Microscopic cross section library
*----------
IF Position "IN" = THEN
  IF Type "32" = THEN
    MicroLib   := LIB: MicroLib ::
      EDIT <<PrtLev>>
      NMIX 36  CTRA WIMS
      MIXS LIB: WIMSD4 FIL: <<LibName>>
      MIX  15  <<Wtemp>>  <<Wdens>>  
        H1  = <<IsoH1>>  <<WwgtH1>>     D2  = <<IsoD2>>  <<WwgtD2>> 
        O16 = <<IsoO16>> <<WwgtO16>>  
      MIX 16 <<Mtemp>>   6.4003        
        ZR  = <<IsoZr>>   98.2083        FE = <<IsoFe>>   0.135 
        CR  = <<IsoCr>>   0.10           NI = <<IsoNi>>   0.055
*        BN  = <<IsoBn>>   <<BNZrII>>
        B10 = <<IsoB10>>  <<B10ZrII>>    B11 = <<IsoB11>>  <<B10ZrII>> 
*Scavanger 1 
      MIX  17 <<Mtemp>>  0.0014     
        HE4 = <<IsoHe4>> 100.0  
      MIX  18 COMB 16 100.0 
      MIX  19 COMB 15 100.0  
      MIX  20 COMB 16 100.0 
*Scavanger 2 
      MIX  21 COMB 17 100.0 
      MIX  22 COMB 18 100.0 
      MIX  23 COMB 19 100.0 
      MIX  24 COMB 20 100.0 
*Scavanger 3 
      MIX  25 COMB 17 100.0 
      MIX  26 COMB 18 100.0 
      MIX  27 COMB 19 100.0 
      MIX  28 COMB 20 100.0 
*Feeder 1
      MIX  29 COMB 17 100.0 
      MIX  30 COMB 18 100.0 
      MIX  31 COMB 19 100.0 
      MIX  32 COMB 20 100.0 
*Feeder 2
      MIX  33 COMB 17 100.0 
      MIX  34 COMB 18 100.0 
      MIX  35 COMB 19 100.0 
      MIX  36 COMB 20 100.0 
      ;
  ELSEIF Type "21" = THEN
    MicroLib   := LIB: MicroLib ::
      EDIT <<PrtLev>>
      NMIX 36  CTRA WIMS
      MIXS LIB: WIMSD4 FIL: <<LibName>>
      MIX  15  <<Wtemp>>  <<Wdens>>  
        H1  = <<IsoH1>>  <<WwgtH1>>     D2  = <<IsoD2>>  <<WwgtD2>> 
        O16 = <<IsoO16>> <<WwgtO16>>  
      MIX 16 <<Mtemp>>   6.4003        
        ZR  = <<IsoZr>>   98.2083        FE = <<IsoFe>>   0.135 
        CR  = <<IsoCr>>   0.10           NI = <<IsoNi>>   0.055
*        BN  = <<IsoBn>>   <<BNZrII>>
        B10 = <<IsoB10>>  <<B10ZrII>>    B11 = <<IsoB11>>  <<B10ZrII>> 
*Scavanger 1 
      MIX  17 <<Mtemp>>  0.0014     
        HE4 = <<IsoHe4>> 100.0  
      MIX  18 COMB 16 100.0 
      MIX  19 COMB 15 100.0  
      MIX  20 COMB 16 100.0 
*Scavanger 2 
      MIX  21 COMB 17 100.0 
      MIX  22 COMB 18 100.0 
      MIX  23 COMB 19 100.0 
      MIX  24 COMB 20 100.0 
*Scavanger 3 
      MIX  25 COMB 15 100.0 
      MIX  26 COMB 15 100.0 
      MIX  27 COMB 15 100.0 
      MIX  28 COMB 15 100.0 
*Feeder 1
      MIX  29 COMB 17 100.0 
      MIX  30 COMB 18 100.0 
      MIX  31 COMB 19 100.0 
      MIX  32 COMB 20 100.0 
*Feeder 2
      MIX  33 COMB 15 100.0 
      MIX  34 COMB 15 100.0 
      MIX  35 COMB 15 100.0 
      MIX  36 COMB 15 100.0 
      ;
  ELSE 
    MicroLib   := LIB: MicroLib ::
      EDIT <<PrtLev>>
      NMIX 36  CTRA WIMS
      MIXS LIB: WIMSD4 FIL: <<LibName>>
      MIX  15  <<Wtemp>>  <<Wdens>>  
        H1  = <<IsoH1>>  <<WwgtH1>>     D2  = <<IsoD2>>  <<WwgtD2>> 
        O16 = <<IsoO16>> <<WwgtO16>>  
      MIX 16 <<Mtemp>>   6.4003        
        ZR  = <<IsoZr>>   98.2083        FE = <<IsoFe>>   0.135 
        CR  = <<IsoCr>>   0.10           NI = <<IsoNi>>   0.055
*        BN  = <<IsoBn>>   <<BNZrII>>
        B10 = <<IsoB10>>  <<B10ZrII>>    B11 = <<IsoB11>>  <<B10ZrII>> 
*Scavanger 1 
      MIX  17 <<Mtemp>>  0.0014     
        HE4 = <<IsoHe4>> 100.0  
      MIX  18 COMB 16 100.0 
      MIX  19 COMB 15 100.0  
      MIX  20 COMB 16 100.0 
*Scavanger 2 
      MIX  21 COMB 15 100.0 
      MIX  22 COMB 15 100.0 
      MIX  23 COMB 15 100.0 
      MIX  24 COMB 15 100.0 
*Scavanger 3 
      MIX  25 COMB 15 100.0 
      MIX  26 COMB 15 100.0 
      MIX  27 COMB 15 100.0 
      MIX  28 COMB 15 100.0 
*Feeder 1
      MIX  29 COMB 15 100.0 
      MIX  30 COMB 15 100.0 
      MIX  31 COMB 15 100.0 
      MIX  32 COMB 15 100.0 
*Feeder 2
      MIX  33 COMB 15 100.0 
      MIX  34 COMB 15 100.0 
      MIX  35 COMB 15 100.0 
      MIX  36 COMB 15 100.0 
      ;
  ENDIF ;
ELSEIF Position "OUT" = THEN
  IF Type "32" = THEN
    MicroLib   := LIB: MicroLib ::
      EDIT <<PrtLev>>
      NMIX 36  CTRA WIMS
      MIXS LIB: WIMSD4 FIL: <<LibName>>
      MIX  15  <<Mtemp>>  0.0014     
        HE4 = <<IsoHe4>> 100.0        
      MIX 16 <<Mtemp>>   6.4003        
        ZR  = <<IsoZr>>   98.2083        FE = <<IsoFe>>   0.135 
        CR  = <<IsoCr>>   0.10           NI = <<IsoNi>>   0.055
*        BN  = <<IsoBn>>   <<BNZrII>>
        B10 = <<IsoB10>>  <<B10ZrII>>    B11 = <<IsoB11>>  <<B10ZrII>>
*Scavanger 1 
      MIX  17 COMB 15 100.0 
      MIX  18 COMB 16 100.0 
      MIX  19  <<Wtemp>>  <<Wdens>>  
        H1  = <<IsoH1>>  <<WwgtH1>>     D2  = <<IsoD2>>  <<WwgtD2>> 
        O16 = <<IsoO16>> <<WwgtO16>>  
      MIX  20 COMB 16 100.0 
*Scavanger 2 
      MIX  21 COMB 17 100.0 
      MIX  22 COMB 18 100.0 
      MIX  23 COMB 19 100.0 
      MIX  24 COMB 20 100.0 
*Scavanger 3 
      MIX  25 COMB 17 100.0 
      MIX  26 COMB 18 100.0 
      MIX  27 COMB 19 100.0 
      MIX  28 COMB 20 100.0 
*Feeder 1
      MIX  29 COMB 17 100.0 
      MIX  30 COMB 18 100.0 
      MIX  31 COMB 19 100.0 
      MIX  32 COMB 20 100.0 
*Feeder 2
      MIX  33 COMB 17 100.0 
      MIX  34 COMB 18 100.0 
      MIX  35 COMB 19 100.0 
      MIX  36 COMB 20 100.0 
      ;
  ELSEIF Type "21" = THEN
    MicroLib   := LIB: MicroLib ::
      EDIT <<PrtLev>>
      NMIX 36  CTRA WIMS
      MIXS LIB: WIMSD4 FIL: <<LibName>>
      MIX  15  <<Mtemp>>  0.0014     
        HE4 = <<IsoHe4>> 100.0        
      MIX 16 <<Mtemp>>   6.4003        
        ZR  = <<IsoZr>>   98.2083        FE = <<IsoFe>>   0.135 
        CR  = <<IsoCr>>   0.10           NI = <<IsoNi>>   0.055
*        BN  = <<IsoBn>>   <<BNZrII>>
        B10 = <<IsoB10>>  <<B10ZrII>>    B11 = <<IsoB11>>  <<B10ZrII>> 
*Scavanger 1 
      MIX  17 COMB 15 100.0 
      MIX  18 COMB 16 100.0 
      MIX  19  <<Wtemp>>  <<Wdens>>  
        H1  = <<IsoH1>>  <<WwgtH1>>     D2  = <<IsoD2>>  <<WwgtD2>> 
        O16 = <<IsoO16>> <<WwgtO16>>  
      MIX  20 COMB 16 100.0 
*Scavanger 2 
      MIX  21 COMB 17 100.0 
      MIX  22 COMB 18 100.0 
      MIX  23 COMB 19 100.0 
      MIX  24 COMB 20 100.0 
*Scavanger 3 
      MIX  25 COMB 15 100.0 
      MIX  26 COMB 15 100.0 
      MIX  27 COMB 15 100.0 
      MIX  28 COMB 15 100.0 
*Feeder 1
      MIX  29 COMB 17 100.0 
      MIX  30 COMB 18 100.0 
      MIX  31 COMB 19 100.0 
      MIX  32 COMB 20 100.0 
*Feeder 2
      MIX  33 COMB 15 100.0 
      MIX  34 COMB 15 100.0 
      MIX  35 COMB 15 100.0 
      MIX  36 COMB 15 100.0 
      ;
  ELSE 
    MicroLib   := LIB: MicroLib ::
      EDIT <<PrtLev>>
      NMIX 36  CTRA WIMS
      MIXS LIB: WIMSD4 FIL: <<LibName>>
      MIX  15  <<Mtemp>>  0.0014     
        HE4 = <<IsoHe4>> 100.0        
      MIX 16 <<Mtemp>>   6.4003        
        ZR  = <<IsoZr>>   98.2083        FE = <<IsoFe>>   0.135 
        CR  = <<IsoCr>>   0.10           NI = <<IsoNi>>   0.055
*        BN  = <<IsoBn>>   <<BNZrII>>
        B10 = <<IsoB10>>  <<B10ZrII>>    B11 = <<IsoB11>>  <<B10ZrII>> 
*Scavanger 1 
      MIX  17 COMB 15 100.0 
      MIX  18 COMB 16 100.0 
      MIX  19  <<Wtemp>>  <<Wdens>>  
        H1  = <<IsoH1>>  <<WwgtH1>>     D2  = <<IsoD2>>  <<WwgtD2>> 
        O16 = <<IsoO16>> <<WwgtO16>>  
      MIX  20 COMB 16 100.0 
*Scavanger 2 
      MIX  21 COMB 15 100.0 
      MIX  22 COMB 15 100.0 
      MIX  23 COMB 15 100.0 
      MIX  24 COMB 15 100.0 
*Scavanger 3 
      MIX  25 COMB 15 100.0 
      MIX  26 COMB 15 100.0 
      MIX  27 COMB 15 100.0 
      MIX  28 COMB 15 100.0 
*Feeder 1
      MIX  29 COMB 15 100.0 
      MIX  30 COMB 15 100.0 
      MIX  31 COMB 15 100.0 
      MIX  32 COMB 15 100.0 
*Feeder 2
      MIX  33 COMB 15 100.0 
      MIX  34 COMB 15 100.0 
      MIX  35 COMB 15 100.0 
      MIX  36 COMB 15 100.0 
      ;
  ENDIF ;
ELSE
  MicroLib   := LIB: MicroLib ::
    EDIT <<PrtLev>>
    NMIX 36  CTRA WIMS
    MIXS LIB: WIMSD4 FIL: <<LibName>>
    MIX  15  <<Mtemp>>  <<Mdens>>  
      H1  = <<IsoH1>>  <<MwgtH1>>     D2  = <<IsoD2>>  <<MwgtD2>> 
      O16 = <<IsoO16>> <<MwgtO16>>  
    MIX  16 COMB 15 100.0 
*Scavanger 1 
    MIX  17 COMB 15 100.0 
    MIX  18 COMB 15 100.0 
    MIX  19 COMB 15 100.0 
    MIX  20 COMB 15 100.0 
*Scavanger 2 
    MIX  21 COMB 15 100.0 
    MIX  22 COMB 15 100.0 
    MIX  23 COMB 15 100.0 
    MIX  24 COMB 15 100.0 
*Scavanger 3 
    MIX  25 COMB 15 100.0 
    MIX  26 COMB 15 100.0 
    MIX  27 COMB 15 100.0 
    MIX  28 COMB 15 100.0 
*Feeder 1
    MIX  29 COMB 15 100.0 
    MIX  30 COMB 15 100.0 
    MIX  31 COMB 15 100.0 
    MIX  32 COMB 15 100.0 
*Feeder 2
    MIX  33 COMB 15 100.0 
    MIX  34 COMB 15 100.0 
    MIX  35 COMB 15 100.0 
    MIX  36 COMB 15 100.0 
    ;
ENDIF ;
QUIT .