blob: 2580b1e2a1e96e94b1f0bf0ee96f0f12f60359ee (
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
|
*----
* TEST CASE TCM70
* FMAC-M generated homogeneous macrolib
* 1D slab sn photon testcase from ROZ6
*----
* Define STRUCTURES and MODULES used
*----
MODULE GEO: FMAC: MAC: SNT: ASM: FLU: OUT: DELETE: END: ;
PROCEDURE assertS ;
LINKED_LIST MACRO MACRO2 GEOM TRACK SYSTEM FLUX EEDITS1 GWT0 ;
SEQ_ASCII FMAC_M :: FILE './W_6MeV_ph_CSD.txt' ;
*----
* Set and track the geometry
*----
GEOM := GEO: :: CAR1D 2
X- VOID X+ VOID
MIX 1 2
MESHX 0.0 0.3 4.0
SPLITX 30 70
;
TRACK := SNT: GEOM ::
TITLE 'HOMOGENIZED W_6Mev_e_CSD BENCHMARK WITH FMAC-M'
EDIT 2 MAXR 1000 SCHM 1 DIAM 2 NLIVO SN 16 SCAT 16 MAXI 300 BTE ;
*----
* Compute the photon flux
*----
MACRO := FMAC: FMAC_M ::
EDIT 2
PARTICLE G
;
MACRO2 := MACRO ;
MACRO := MAC: MACRO MACRO2 ::
MIX 2 1 OLDL
;
MACRO2 := DELETE: MACRO2 ;
REAL qsour := 1.E10 2.0 / $Pi_R / ;
ECHO "qsour=" qsour ;
MACRO := MAC: MACRO ::
READ INPUT
MIX 1 FIXE
<<qsour>> 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.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 TRACK :: EDIT 2 ARM ;
FLUX := FLU: MACRO TRACK SYSTEM :: EDIT 1 TYPE S EXTE 100 ;
EEDITS1 := OUT: FLUX TRACK MACRO GEOM ::
EDIT 2 COND SOUR 0.0518134 INTG NONE ! 1/(\ro_{W} g/cm^{3})
;
*----
* Asserts
*----
GWT0 := EEDITS1 :: STEP UP GROUP STEP AT 1 ;
assertS GWT0 :: 'FLUX-INTG' 1 7.357735E-4 ;
assertS GWT0 :: 'FLUX-INTG' 5 8.507982E-4 ;
assertS GWT0 :: 'FLUX-INTG' 10 9.365114E-4 ;
assertS GWT0 :: 'FLUX-INTG' 15 9.739149E-4 ;
assertS GWT0 :: 'FLUX-INTG' 20 9.688310E-4 ;
assertS GWT0 :: 'FLUX-INTG' 25 9.216961E-4 ;
assertS GWT0 :: 'FLUX-INTG' 30 8.291884E-4 ;
*
ECHO "test TCM70 completed" ;
END: ;
|
