blob: f0a579f7c19cf66097c4f26f924d76464528f687 (
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
|
*----
* TEST CASE TCM69
* High-order schemes for the continuous slowing-down term of the
* Boltzmann Fokker-Planck equation in Cartesian geometries.
*----
* Define STRUCTURES and MODULES used
*----
MODULE GEO: FMAC: MAC: SNT: ASM: FLU: PSOUR: OUT: DELETE: END: UTL:
HEAT: ABORT: ;
PROCEDURE assertS assertS_low ;
LINKED_LIST EMACRO EMACRO2 GEOM ETRACK SYSTEM
EFLUX ESOUR EEDITS DEPOS EFLX ;
SEQ_ASCII FMAC_M :: FILE './W_6MeV_e_ph_CSD.txt' ;
*----
* High-order energy schemes 1D Cartesian
*----
GEOM := GEO: :: CAR1D 1
X- VOID X+ VOID
MIX 1
MESHX 0.0 0.05
SPLITX 20
;
ETRACK := SNT: GEOM ::
EDIT 2
MAXR 1000
SCHM 1 DIAM 2
ESCHM 1 EDIAM 2
NLIVO
SN 8 SCAT 8
EPSI 1.E-5
MAXI 300
BFPG
;
EMACRO := FMAC: FMAC_M :: EDIT 2 PARTICLE BETA ;
EMACRO2 := EMACRO ;
EMACRO := MAC: EMACRO EMACRO2 :: MIX 1 1 OLDL ;
EMACRO2 := DELETE: EMACRO2 ;
ESOUR := PSOUR: EMACRO ETRACK GEOM ::
MONO 1
X-
INTG 1 1.0
DIR 1.0 0.0 0.0
;
SYSTEM := ASM: EMACRO ETRACK :: EDIT 2 ARM ;
EFLUX := FLU: EMACRO ETRACK SYSTEM ESOUR :: EDIT 1 TYPE S EXTE 100 ;
EEDITS := OUT: EFLUX ETRACK EMACRO GEOM :: EDIT 2 COND INTG NONE ;
*----
* Asserts in 1D Cartesian
*----
EFLX := EEDITS :: STEP UP GROUP STEP AT 1 ;
assertS EFLX :: 'FLUX-INTG' 1 3.990423E-3 ;
assertS EFLX :: 'FLUX-INTG' 5 6.039801E-3 ;
assertS EFLX :: 'FLUX-INTG' 15 5.647896E-3 ;
GEOM ETRACK EMACRO ESOUR SYSTEM EFLUX EEDITS EFLX := DELETE:
GEOM ETRACK EMACRO ESOUR SYSTEM EFLUX EEDITS EFLX ;
*----
* High-order energy schemes 2D Cartesian
*----
GEOM := GEO: :: CAR2D 1 1
X- VOID X+ VOID
Y- VOID Y+ VOID
MIX 1
MESHX 0.0 0.05
MESHY 0.0 0.05
SPLITX 5
SPLITY 5
;
ETRACK := SNT: GEOM ::
EDIT 2
MAXR 1000
SCHM 2 DIAM 1
ESCHM 2 EDIAM 1
NLIVO
SN 4 SCAT 6
EPSI 1.E-5
MAXI 300
BFPG
GQ 1
;
EMACRO := FMAC: FMAC_M :: EDIT 2 PARTICLE BETA ;
EMACRO2 := EMACRO ;
EMACRO := MAC: EMACRO EMACRO2 :: MIX 1 1 OLDL ;
EMACRO2 := DELETE: EMACRO2 ;
ESOUR := PSOUR: EMACRO ETRACK GEOM ::
ISO 1
INTG 1 1.0
XLIM 0.00 0.05
YLIM 0.00 0.02
;
SYSTEM := ASM: EMACRO ETRACK :: EDIT 2 ARM ;
EFLUX := FLU: EMACRO ETRACK SYSTEM ESOUR :: EDIT 1 TYPE S EXTE 100 ;
EEDITS := OUT: EFLUX ETRACK EMACRO GEOM :: EDIT 2 COND INTG NONE ;
*----
* Asserts in 2D Cartesian
*----
EFLX := EEDITS :: STEP UP GROUP STEP AT 1 ;
assertS EFLX :: 'FLUX-INTG' 1 1.930571E-6 ;
assertS EFLX :: 'FLUX-INTG' 5 1.930571E-6 ;
assertS EFLX :: 'FLUX-INTG' 8 2.869739E-6 ;
GEOM ETRACK EMACRO ESOUR SYSTEM EFLUX EEDITS EFLX := DELETE:
GEOM ETRACK EMACRO ESOUR SYSTEM EFLUX EEDITS EFLX ;
*----
* High-order energy schemes 3D Cartesian
*----
GEOM := GEO: :: CAR3D 1 1 1
X- VOID X+ VOID
Y- VOID Y+ VOID
Z- VOID Z+ VOID
MIX 1
MESHX 0.0 0.05
MESHY 0.0 0.05
MESHZ 0.0 0.05
SPLITX 2
SPLITY 2
SPLITZ 2
;
ETRACK := SNT: GEOM ::
EDIT 2
MAXR 1000
SCHM 2 DIAM 1
ESCHM 2 EDIAM 1
NLIVO
SN 4 SCAT 6
EPSI 1.E-5
MAXI 300
BFPG
GQ 2
;
EMACRO := FMAC: FMAC_M :: EDIT 2 PARTICLE BETA ;
EMACRO2 := EMACRO ;
EMACRO := MAC: EMACRO EMACRO2 :: MIX 1 1 OLDL ;
EMACRO2 := DELETE: EMACRO2 ;
ESOUR := PSOUR: EMACRO ETRACK GEOM ::
MONO 1
X-
INTG 1 1.0
DIR 1.0 0.0 0.0
YLIM 0.0 0.05
ZLIM 0.0 0.05
;
SYSTEM := ASM: EMACRO ETRACK :: EDIT 2 ARM ;
EFLUX := FLU: EMACRO ETRACK SYSTEM ESOUR :: EDIT 1 TYPE S EXTE 100 ;
EEDITS := OUT: EFLUX ETRACK EMACRO GEOM :: EDIT 2 COND INTG NONE ;
*----
* Asserts in 3D Cartesian
*----
EFLX := EEDITS :: STEP UP GROUP STEP AT 1 ;
assertS_low EFLX :: 'FLUX-INTG' 1 1.241492E-6 5.0E-3 ;
assertS_low EFLX :: 'FLUX-INTG' 2 5.177682E-7 5.0E-3 ;
assertS_low EFLX :: 'FLUX-INTG' 3 1.241492E-6 5.0E-3 ;
*----
* End of test
*----
ECHO "Test TCM69 completed" ;
END: ;
|
