blob: d6246b20affb3701b68ab9c2bea0fd6ebf7ff9dd (
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
|
*----
* TEST CASE TCM65
* FMAC-M generated homogeneous macrolib
* 1D slab sn testcase from ROZ6
* 1D slab sn neutron-photon coupled testcase from ROZ6
*----
* Define STRUCTURES and MODULES used
*----
MODULE GEO: FMAC: MAC: SNT: ASM: FLU: PSOUR: OUT: DELETE: END: ;
PROCEDURE assertS ;
LINKED_LIST MACRO GEOM GTRACK NTRACK SYSTEM FLUX GMACRO GSOUR GFLUX
GEDITS GWT0 ;
SEQ_ASCII FMAC_M :: FILE './caskmix_FMAC-M.txt' ;
*----
* Set and track the geometry
*----
GEOM := GEO: :: CAR1D 5
X- VOID X+ VOID
MIX 1 2 3 4 5
MESHX 0.0 10.0 20.0 30.0 40.0 50.0
SPLITX 5 5 5 5 5
;
GTRACK := SNT: GEOM ::
TITLE 'HOMOGENIZED CASKMIX BENCHMARK WITH FMAC-M'
EDIT 2 MAXR 1000 SN 16 SCAT 4 ;
NTRACK := GTRACK ;
*----
* Compute the neutron flux
*----
MACRO := FMAC: FMAC_M ::
EDIT 2
PARTICLE N
;
MACRO := MAC: MACRO ::
READ INPUT
MIX 1 FIXE
1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0 0.0 0.0 0.0
;
SYSTEM := ASM: MACRO NTRACK :: EDIT 1 ARM ;
FLUX := FLU: MACRO NTRACK SYSTEM :: EDIT 1 TYPE S EXTE 100 ;
SYSTEM := DELETE: SYSTEM ;
*----
* Set the gamma source
*----
GMACRO := FMAC: FMAC_M ::
EDIT 2
PARTICLE G
;
GSOUR := PSOUR: GMACRO GTRACK NTRACK FLUX ::
EDIT 1
PARTICLE N
;
END: ;
*----
* Compute the gamma flux
*----
SYSTEM := ASM: GMACRO GTRACK :: EDIT 1 ARM ;
GFLUX := FLU: GMACRO GTRACK SYSTEM GSOUR :: EDIT 1 TYPE S EXTE 100 ;
GEDITS := OUT: GFLUX GTRACK GMACRO GEOM ::
EDIT 2 INTG MIX
;
*----
* Asserts
*----
GWT0 := GEDITS :: STEP UP GROUP STEP AT 18 ;
assertS GWT0 :: 'FLUX-INTG' 1 7.89299488E-1 ;
assertS GWT0 :: 'FLUX-INTG' 2 8.66148663E+0 ;
assertS GWT0 :: 'FLUX-INTG' 3 5.81171036E-1 ;
assertS GWT0 :: 'FLUX-INTG' 4 1.09355543E-2 ;
assertS GWT0 :: 'FLUX-INTG' 5 3.78631288E-3 ;
*
ECHO "test TCM65 completed" ;
END: ;
|
