Initial commit from Polytechnique Montreal

1 files changed, 249 insertions, 0 deletions

diff --git a/Dragon/data/tmacro_proc/TCM47.c2m b/Dragon/data/tmacro_proc/TCM47.c2m
new file mode 100644
index 0000000..b5df858
--- /dev/null
+++ b/Dragon/data/tmacro_proc/TCM47.c2m
@@ -0,0 +1,249 @@
+*---
+* Regression test for MCCG:
+*---
+* TEST CASE TCM47
+* PRISMATIC 3D ANISOTROPIC SCATTERING BENCHMARK WITH VOID (SYMMETRIES)
+* MACROLIB-DEFINED CROSS SECTIONS
+*---
+
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOMRR GEOMRS GEOMSR GEOMSS GEOM2 
+            MACRO TRACK SYSTEM FLUX EDIT ;
+SEQ_BINARY TRK_FIL ;
+SEQ_ASCII OUT ;
+MODULE GREP: GEO: MAC: EXCELT: NXT: MCCGT: ASM: FLU: DELETE: END: ;
+PROCEDURE assertS ;
+REAL krr ksr krs kss k2d ;
+*---
+* Geometry Geom
+*---
+GEOMRR := GEO: :: CAR3D 5 5 2
+           X- REFL X+ VOID
+           Y- REFL Y+ VOID
+           Z- REFL Z+ REFL
+           CELL C1 C1 C2 C2 C3
+                C1 C1 C2 C2 C3
+                C2 C2 C2 C2 C3
+                C2 C2 C2 C3 C3
+                C3 C3 C3 C3 C0
+                C1 C1 C2 C2 C3
+                C1 C1 C2 C2 C3
+                C2 C2 C2 C2 C3
+                C2 C2 C2 C3 C3
+                C3 C3 C3 C3 C0
+        ::: C1 := GEO: CAR3D 1 1 1
+                MESHX 0.0 4.0 
+                MESHY 0.0 4.0 
+                MESHZ 0.0 1.0
+                MIX 1
+            ;
+        ::: C2 := GEO: C1
+                MIX 2
+            ;
+        ::: C3 := GEO: C1
+                MIX 3
+            ;
+        ::: C0 := GEO: C1
+                MIX 4
+            ;
+ ;
+GEOMSR := GEO: GEOMRR :: 
+           Z- SSYM Z+ REFL
+ ;
+GEOMRS := GEO: GEOMRR :: 
+           Z- REFL Z+ SSYM
+ ;
+GEOMSS := GEO: GEOMRR :: 
+           Z- SSYM Z+ SSYM
+ ;
+GEOM2 := GEO: :: CAR2D 5 5
+           X- REFL X+ VOID
+           Y- REFL Y+ VOID
+           CELL C1 C1 C2 C2 C3
+                C1 C1 C2 C2 C3
+                C2 C2 C2 C2 C3
+                C2 C2 C2 C3 C3
+                C3 C3 C3 C3 C0
+        ::: C1 := GEO: CAR2D 1 1
+                MESHX 0.0 4.0 
+                MESHY 0.0 4.0 
+                MIX 1
+            ;
+        ::: C2 := GEO: C1
+                MIX 2
+            ;
+        ::: C3 := GEO: C1
+                MIX 3
+            ;
+        ::: C0 := GEO: C1
+                MIX 4
+            ;
+ ;
+*---
+*  Macrolib
+*--
+MACRO := MAC: ::
+ EDIT 0 NGRO 1 NMIX 4 NIFI 1 ANIS 2
+ READ INPUT
+ MIX     1
+     TOTAL  0.25
+    NUSIGF  0.155
+       CHI  1.0
+      SCAT  1 1 0.13
+            1 1 0.07
+ MIX     2
+     TOTAL  0.25
+     SCAT   1 1 0.24
+            1 1 0.15
+ MIX     3
+     TOTAL  0.75
+     SCAT   1 1 0.0
+            1 1 0.0
+ MIX     4
+     TOTAL  1E-8
+     SCAT   1 1 0.0
+            1 1 0.0
+ ;
+
+******************************
+* NXT 3D PRISMATIC CALCULATION Z- REFL Z+ REFL
+******************************
+*----
+*  Tracking calculation 
+*----
+TRACK TRK_FIL := NXT: GEOMRR ::
+  TITLE '3D ANISOTROPIC BENCHMARK, 1 GROUP'
+  EDIT 0 ANIS 1 PRIZ 5.0 GAUS 4 TISO 4 5.0 ;
+*---
+*  Flux Calculation
+*---
+TRACK := MCCGT: TRACK TRK_FIL ::
+  EDIT 1
+  EPSI 1E-6 HDD 0.0
+  MAXI 10 KRYL 10 SCR 1 AAC 1 TMT ;
+
+SYSTEM := ASM: MACRO TRACK TRK_FIL ::
+  EDIT 0 ARM ;
+FLUX :=  FLU: MACRO TRACK TRK_FIL SYSTEM ::
+  EDIT 2 TYPE K EXTE 100 1E-5 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.032782E+00 ;
+GREP: FLUX :: GETVAL 'K-EFFECTIVE' 1 1 1 >>krr<< ; 
+TRACK TRK_FIL SYSTEM FLUX := DELETE:
+TRACK TRK_FIL SYSTEM FLUX ;
+******************************
+* NXT 3D PRISMATIC CALCULATION Z- REFL Z+ SSYM
+******************************
+*----
+*  Tracking calculation 
+*----
+TRACK TRK_FIL := NXT: GEOMRS ::
+  TITLE '3D ANISOTROPIC BENCHMARK, 1 GROUP'
+  EDIT 0 ANIS 1 PRIZ 5.0 GAUS 4 TISO 4 5.0 ;
+*---
+*  Flux Calculation
+*---
+TRACK := MCCGT: TRACK TRK_FIL ::
+  EDIT 1
+  EPSI 1E-6 HDD 0.0
+  MAXI 10 KRYL 10 SCR 1 AAC 1 TMT ;
+
+SYSTEM := ASM: MACRO TRACK TRK_FIL ::
+  EDIT 0 ARM ;
+FLUX :=  FLU: MACRO TRACK TRK_FIL SYSTEM ::
+  EDIT 2 TYPE K EXTE 100 1E-5 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.020030E+00 ;
+GREP: FLUX :: GETVAL 'K-EFFECTIVE' 1 1 1 >>krs<< ; 
+TRACK TRK_FIL SYSTEM FLUX := DELETE:
+TRACK TRK_FIL SYSTEM FLUX ;
+******************************
+* NXT 3D PRISMATIC CALCULATION Z- SSYM Z+ REFL
+******************************
+*----
+*  Tracking calculation 
+*----
+TRACK TRK_FIL := NXT: GEOMSR ::
+  TITLE '3D ANISOTROPIC BENCHMARK, 1 GROUP'
+  EDIT 0 ANIS 1 PRIZ 5.0 GAUS 4 TISO 4 5.0 ;
+*---
+*  Flux Calculation
+*---
+TRACK := MCCGT: TRACK TRK_FIL ::
+  EDIT 1
+  EPSI 1E-6 HDD 0.0
+  MAXI 10 KRYL 10 SCR 1 AAC 1 TMT ;
+
+SYSTEM := ASM: MACRO TRACK TRK_FIL ::
+  EDIT 0 ARM ;
+FLUX :=  FLU: MACRO TRACK TRK_FIL SYSTEM ::
+  EDIT 2 TYPE K EXTE 100 1E-5 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.020031E+00 ;
+GREP: FLUX :: GETVAL 'K-EFFECTIVE' 1 1 1 >>ksr<< ; 
+TRACK TRK_FIL SYSTEM FLUX := DELETE:
+TRACK TRK_FIL SYSTEM FLUX ;
+******************************
+* NXT 3D PRISMATIC CALCULATION Z+ SSYM Z- SSYM
+******************************
+*----
+*  Tracking calculation 
+*----
+TRACK TRK_FIL := NXT: GEOMSS ::
+  TITLE '3D ANISOTROPIC BENCHMARK, 1 GROUP'
+  EDIT 0 ANIS 1 PRIZ 5.0 GAUS 4 TISO 4 5.0 ;
+*---
+*  Flux Calculation
+*---
+TRACK := MCCGT: TRACK TRK_FIL ::
+  EDIT 1
+  EPSI 1E-6 HDD 0.0
+  MAXI 10 KRYL 10 SCR 1 AAC 1 TMT ;
+
+SYSTEM := ASM: MACRO TRACK TRK_FIL ::
+  EDIT 0 ARM ;
+FLUX :=  FLU: MACRO TRACK TRK_FIL SYSTEM ::
+  EDIT 2 TYPE K EXTE 100 1E-5 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 9.905481E-01 ;
+GREP: FLUX :: GETVAL 'K-EFFECTIVE' 1 1 1 >>kss<< ; 
+TRACK TRK_FIL SYSTEM FLUX := DELETE:
+TRACK TRK_FIL SYSTEM FLUX ;
+********************
+* NXT 2D CALCULATION
+********************
+*----
+*  Tracking calculation 
+*----
+TRACK TRK_FIL := NXT: GEOM2 ::
+  TITLE '3D ANISOTROPIC BENCHMARK, 1 GROUP'
+  EDIT 0 ANIS 1 TISO 4 5.0 ;
+*---
+*  Flux Calculation
+*---
+TRACK := MCCGT: TRACK TRK_FIL ::
+  EDIT 1 DGAU 4
+  EPSI 1E-6 HDD 0.0
+  MAXI 10 KRYL 10 SCR 1 AAC 1 TMT ;
+
+SYSTEM := ASM: MACRO TRACK TRK_FIL ::
+  EDIT 0 ARM ;
+FLUX :=  FLU: MACRO TRACK TRK_FIL SYSTEM ::
+  EDIT 2 TYPE K EXTE 100 1E-5 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 9.905450E-1 ;
+GREP: FLUX :: GETVAL 'K-EFFECTIVE' 1 1 1 >>k2d<< ; 
+TRACK TRK_FIL SYSTEM FLUX := DELETE:
+TRACK TRK_FIL SYSTEM FLUX ;
+
+ECHO "**********MCCG***********" ;
+ECHO "DELTA = K_RR - K_2D ="
+krr "-" k2d "=" krr k2d - 1E3 * "mk" ;
+ECHO "DELTA = K_RS - K_2D ="
+krs "-" k2d "=" krs k2d - 1E3 * "mk" ;
+ECHO "DELTA = K_SR - K_2D ="
+ksr "-" k2d "=" ksr k2d - 1E3 * "mk" ;
+ECHO "DELTA = K_SS - K_2D ="
+kss "-" k2d "=" kss k2d - 1E3 * "mk" ;
+ECHO "*************************" ;
+
+GEOMRR GEOM2 MACRO := DELETE: GEOMRR GEOM2 MACRO ;
+ECHO "test TCM47 completed" ;
+END: ;