summaryrefslogtreecommitdiff
path: root/Donjon/data/Reflector_beavrs_DF-ANM.x2m
blob: c99fd7026d4fca7938ab1c2d2915b81cca897476 (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
*---
*  Nodal Expansion Reflector DF-ANM equivalence prodedure
*---
LINKED_LIST Geom FVol Param MACRO ;
STRING nameFeed := "assbly_caseA_mpo_boron.hdf" ;
STRING nameRefl := "reflector_beavrs_DF-ANM_mpo.hdf" ;
STRING Palier := "BEAVRS" ;
STRING htype := "DF-ANM" ;
STRING LibType := "JEF22" ;
STRING NuclData := "DLIB_99" ;
STRING OptSph := "NOSP" ;
STRING OptNoal := "ALBE" ;
STRING OptSS := "TONE" ;
MODULE GEO: UTL: MCR: END: ;
HDF5_FILE MpoRefl :: FILE <<nameRefl>> ;
PROCEDURE REFL_MPO ;
!
Geom := GEO: ::
   CAR1D 12
   X- REFL X+ VOID
   MESHX -43.0 -21.5 -5.0 0.0 0.005 2.2225 3.9825 9.6975 10.1975 15.9125 21.5 21.505
         35.9125
   MIX     1  1  1     2        3         4         3          5          6          5
           !  fuel   gap       ss    ss+h2o        ss        h2o  ss_shield        h2o
                 7            8
           !   gap          h2o
   SPLITX  20 15 20  1         25        24        10          5          6          5
                 1            5
;

* Define volume fractions in 8 mixtures
* mix=1 : feeding assembly (first fraction is always set to 0)
* mix=2 : fuel-reflector gap
* mix=7 : right-most gap to compute albedo
* mix=8 : residual right reflector
INTEGER nbZones := 8 ;
REAL fvacier := 0.05 ;
REAL fveau := 1.0 fvacier - ;
FVol := UTL: ::
                  ! MIX  1  2  3  4           5  6  7  8
CREA 'H2O' <<nbZones>> = 0. 0. 0. <<fveau>>   1. 0. 1. 1.
CREA 'SS ' <<nbZones>> = 0. 1. 1. <<fvacier>> 0. 1. 0. 0.
;
UTL: FVol :: DUMP ;

* Define feeding assembly parameters (0 to 3 parameters allowed)
Param := UTL: :: CREA PARAMVALU 1 = 975.
                 STEP UP 'PARAMNAME' NEW LIST 1
                 CREA 1 2 = 'C-BO' 'RE  ' ;
UTL: Param :: DUMP ;

INTEGER nlf := 1 ;
MpoRefl := REFL_MPO Geom FVol Param :: <<nameFeed>> <<nameRefl>> <<Palier>> <<htype>>
                            <<nlf>> <<LibType>> <<NuclData>> <<OptSph>> <<OptNoal>>
                            <<OptSS>> ;

************************************************************************
* MPO extrapolation of the reflector at 350 ppm Boron
************************************************************************
REAL densB1 := 350.0 ;
MACRO :=  MCR:  MpoRefl :: EDIT 3 MACRO NMIX 1
   MPO MpoRefl output_0
   MIX  1 FROM 2 SET 'C-BORE' <<densB1>> SET 'Method' <<htype>> SET 'Palier' <<Palier>> ENDMIX
;
UTL: MACRO :: DIR DUMP ;
ECHO "test Reflector_beavrs_DF-ANM.x2m completed" ;
END: ;