Initial commit from Polytechnique Montreal

11 files changed, 1976 insertions, 0 deletions

diff --git a/Dragon/data/trowland_shem295_jeff3.1_proc/gd_295_kec1.c2m b/Dragon/data/trowland_shem295_jeff3.1_proc/gd_295_kec1.c2m
new file mode 100644
index 0000000..ef84375
--- /dev/null
+++ b/Dragon/data/trowland_shem295_jeff3.1_proc/gd_295_kec1.c2m
@@ -0,0 +1,91 @@
+*----
+*
+*  TEST CASE gduo2_295_kec1
+*  DISTRIBUTED SELF-SHIELDING
+*  295-GROUP JEF3.1.1 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK_SS TRACK LIBRARY LIBRARY2 CP CALC ;
+SEQ_ASCII UOX_TBH ;
+SEQ_BINARY TF_EXC ;
+MODULE LIB: GEO: SYBILT: G2S: SALT: MAC: USS: TONE: ASM: FLU: DELETE: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 := 
+        1 2 3 4 5 6 7 8 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+REAL P_Zr90     P_Zr91     P_Zr92    P_Zr94    P_Zr96    :=
+     51.45E-02  11.22E-02  17.15E-02 17.38E-02 2.80E-02  ;
+REAL N_Zr0 := 3.88E-2 ;
+REAL N_Zr90       N_Zr91       N_Zr92       N_Zr94       N_Zr96    :=
+     N_Zr0 P_Zr90 *  N_Zr0 P_Zr91 *  N_Zr0 P_Zr92 *  N_Zr0 P_Zr94 *
+     N_Zr0 P_Zr96 *  ;
+INTEGER an2d := 12 ;
+REAL densur := 20.0 ;
+*
+EVALUATE RCOMB6 := 0.41 ;
+EVALUATE RCOMB1 := 0.4 SQRT RCOMB6 * ;
+EVALUATE RCOMB2 := 0.7 SQRT RCOMB6 * ;
+EVALUATE RCOMB3 := 0.8 SQRT RCOMB6 * ;
+EVALUATE RCOMB4 := 0.9 SQRT RCOMB6 * ;
+EVALUATE RCOMB5 := 0.95 SQRT RCOMB6 * ;
+GEOM := GEO: :: CARCEL 8
+      X- REFL X+ REFL Y- REFL Y+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.475 0.6045552
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE1>>
+      MESHX 0.0 1.26 MESHY 0.0 1.26
+          ;
+TRACK_SS := SYBILT: GEOM ::
+     TITLE 'UO2 ROWLAND S BENCHMARK'
+     MAXR 20 QUA2 20 3 DP01 MAXZ 200000 ;
+UOX_TBH := G2S: GEOM ;
+TRACK TF_EXC := SALT: UOX_TBH ::
+  EDIT 3
+  ALLG
+  TSPC EQW2 <<an2d>> <<densur>> REND
+;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ RSE       (*RESONANCE SPECTRUM EXPANSION METHOD*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 294.0 (*COMB0101*)
+    O16     = O16    4.59036E-2
+    U235    = U235   7.90006E-4 1
+    U238    = U238   2.08874E-2 CORR 1
+    Gd154   = Gd154  3.78440E-5 1
+    Gd155   = Gd155  2.55262E-4 1
+    Gd156   = Gd156  3.50791E-4 1
+    Gd157   = Gd157  2.66479E-4 CORR 1
+    Gd158   = Gd158  4.20283E-4 1
+    Gd160   = Gd160  3.65230E-4 1
+MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 294.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 294.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16    2.372E-2
+;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 0.3288226 ;
+QUIT "LIST" .
diff --git a/Dragon/data/trowland_shem295_jeff3.1_proc/mox_295_kec1.c2m b/Dragon/data/trowland_shem295_jeff3.1_proc/mox_295_kec1.c2m
new file mode 100644
index 0000000..06c9f2b
--- /dev/null
+++ b/Dragon/data/trowland_shem295_jeff3.1_proc/mox_295_kec1.c2m
@@ -0,0 +1,176 @@
+*----
+*
+*  TEST CASE mox_295_kec1
+*  MOX ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  295-GROUP JEF3.1.1 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK_SS TRACK LIBRARY LIBRARY2 CP CALC ;
+SEQ_ASCII UOX_TBH ;
+SEQ_BINARY TF_EXC ;
+MODULE LIB: GEO: SYBILT: G2S: SALT: MAC: USS: TONE: ASM: FLU: DELETE: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 := 
+        1 2 3 4 5 6 7 8 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+REAL P_Zr90     P_Zr91     P_Zr92    P_Zr94    P_Zr96    :=
+     51.45E-02  11.22E-02  17.15E-02 17.38E-02 2.80E-02  ;
+REAL N_Zr0 := 3.88E-2 ;
+REAL N_Zr90       N_Zr91       N_Zr92       N_Zr94       N_Zr96    :=
+     N_Zr0 P_Zr90 *  N_Zr0 P_Zr91 *  N_Zr0 P_Zr92 *  N_Zr0 P_Zr94 *
+     N_Zr0 P_Zr96 *  ;
+INTEGER an2d := 12 ;
+REAL densur := 20.0 ;
+*
+EVALUATE RCOMB6 := 0.41 ;
+EVALUATE RCOMB1 := 0.4 SQRT RCOMB6 * ;
+EVALUATE RCOMB2 := 0.7 SQRT RCOMB6 * ;
+EVALUATE RCOMB3 := 0.8 SQRT RCOMB6 * ;
+EVALUATE RCOMB4 := 0.9 SQRT RCOMB6 * ;
+EVALUATE RCOMB5 := 0.95 SQRT RCOMB6 * ;
+GEOM := GEO: :: CARCEL 8
+      X- REFL X+ REFL Y- REFL Y+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.475 0.6045552
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE1>>
+      MESHX 0.0 1.26 MESHY 0.0 1.26
+          ;
+TRACK_SS := SYBILT: GEOM ::
+     TITLE 'UO2 ROWLAND S BENCHMARK'
+     MAXR 20 QUA2 20 3 DP01 MAXZ 200000 ;
+UOX_TBH := G2S: GEOM ;
+TRACK TF_EXC := SALT: UOX_TBH ::
+  EDIT 3
+  ALLG
+  TSPC <<an2d>> <<densur>> REND
+;
+*
+* Subgroup projection method
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ CALENDF 3 PT      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.228971 ;
+*
+* Tone's method
+LIBRARY CP CALC := DELETE: LIBRARY CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY := TONE: LIBRARY TRACK_SS :: EDIT 1 TRAN MXIT 3 SPH GRMIN 52 ;
+CP := ASM: LIBRARY TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.227709 ;
+*
+* Resonance spectrum expansion method
+LIBRARY LIBRARY2 CP CALC := DELETE: LIBRARY LIBRARY2 CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ RSE       (*RSE TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.229507 ;
+ECHO "test mox_295_kec1 completed" ;
+QUIT "LIST" .
diff --git a/Dragon/data/trowland_shem295_jeff3.1_proc/mox_295_kec2.c2m b/Dragon/data/trowland_shem295_jeff3.1_proc/mox_295_kec2.c2m
new file mode 100644
index 0000000..ae2cd66
--- /dev/null
+++ b/Dragon/data/trowland_shem295_jeff3.1_proc/mox_295_kec2.c2m
@@ -0,0 +1,176 @@
+*----
+*
+*  TEST CASE mox_295_kec2
+*  MOX ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  295-GROUP JEF3.1.1 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK_SS TRACK LIBRARY LIBRARY2 CP CALC ;
+SEQ_ASCII UOX_TBH ;
+SEQ_BINARY TF_EXC ;
+MODULE LIB: GEO: SYBILT: G2S: SALT: MAC: USS: TONE: ASM: FLU: DELETE: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 := 
+        1 2 3 4 5 6 7 8 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+REAL P_Zr90     P_Zr91     P_Zr92    P_Zr94    P_Zr96    :=
+     51.45E-02  11.22E-02  17.15E-02 17.38E-02 2.80E-02  ;
+REAL N_Zr0 := 3.88E-2 ;
+REAL N_Zr90       N_Zr91       N_Zr92       N_Zr94       N_Zr96    :=
+     N_Zr0 P_Zr90 *  N_Zr0 P_Zr91 *  N_Zr0 P_Zr92 *  N_Zr0 P_Zr94 *
+     N_Zr0 P_Zr96 *  ;
+INTEGER an2d := 12 ;
+REAL densur := 20.0 ;
+*
+EVALUATE RCOMB6 := 0.41 ;
+EVALUATE RCOMB1 := 0.4 SQRT RCOMB6 * ;
+EVALUATE RCOMB2 := 0.7 SQRT RCOMB6 * ;
+EVALUATE RCOMB3 := 0.8 SQRT RCOMB6 * ;
+EVALUATE RCOMB4 := 0.9 SQRT RCOMB6 * ;
+EVALUATE RCOMB5 := 0.95 SQRT RCOMB6 * ;
+GEOM := GEO: :: CARCEL 8
+      X- REFL X+ REFL Y- REFL Y+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.475 0.6045552
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE1>>
+      MESHX 0.0 1.26 MESHY 0.0 1.26
+          ;
+TRACK_SS := SYBILT: GEOM ::
+     TITLE 'UO2 ROWLAND S BENCHMARK'
+     MAXR 20 QUA2 20 3 DP01 MAXZ 200000 ;
+UOX_TBH := G2S: GEOM ;
+TRACK TF_EXC := SALT: UOX_TBH ::
+  EDIT 3
+  ALLG
+  TSPC <<an2d>> <<densur>> REND
+;
+*
+* Subgroup projection method
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ CALENDF 3 PT      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 600.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.212714 ;
+*
+* Tone's method
+LIBRARY CP CALC := DELETE: LIBRARY CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 600.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY := TONE: LIBRARY TRACK_SS :: EDIT 1 TRAN MXIT 3 SPH GRMIN 52 ;
+CP := ASM: LIBRARY TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.211497 ;
+*
+** Resonance spectrum expansion method
+LIBRARY LIBRARY2 CP CALC := DELETE: LIBRARY LIBRARY2 CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ RSE       (*RSE TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 600.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.212991 ;
+ECHO "test mox_295_kec2 completed" ;
+QUIT "LIST" .
diff --git a/Dragon/data/trowland_shem295_jeff3.1_proc/mox_295_kec3.c2m b/Dragon/data/trowland_shem295_jeff3.1_proc/mox_295_kec3.c2m
new file mode 100644
index 0000000..fd4ea3f
--- /dev/null
+++ b/Dragon/data/trowland_shem295_jeff3.1_proc/mox_295_kec3.c2m
@@ -0,0 +1,176 @@
+*----
+*
+*  TEST CASE mox_295_kec3
+*  MOX ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  295-GROUP JEF3.1.1 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK_SS TRACK LIBRARY LIBRARY2 CP CALC ;
+SEQ_ASCII UOX_TBH ;
+SEQ_BINARY TF_EXC ;
+MODULE LIB: GEO: SYBILT: G2S: SALT: MAC: USS: TONE: ASM: FLU: DELETE: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 := 
+        1 2 3 4 5 6 7 8 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+REAL P_Zr90     P_Zr91     P_Zr92    P_Zr94    P_Zr96    :=
+     51.45E-02  11.22E-02  17.15E-02 17.38E-02 2.80E-02  ;
+REAL N_Zr0 := 3.88E-2 ;
+REAL N_Zr90       N_Zr91       N_Zr92       N_Zr94       N_Zr96    :=
+     N_Zr0 P_Zr90 *  N_Zr0 P_Zr91 *  N_Zr0 P_Zr92 *  N_Zr0 P_Zr94 *
+     N_Zr0 P_Zr96 *  ;
+INTEGER an2d := 12 ;
+REAL densur := 20.0 ;
+*
+EVALUATE RCOMB6 := 0.41 ;
+EVALUATE RCOMB1 := 0.4 SQRT RCOMB6 * ;
+EVALUATE RCOMB2 := 0.7 SQRT RCOMB6 * ;
+EVALUATE RCOMB3 := 0.8 SQRT RCOMB6 * ;
+EVALUATE RCOMB4 := 0.9 SQRT RCOMB6 * ;
+EVALUATE RCOMB5 := 0.95 SQRT RCOMB6 * ;
+GEOM := GEO: :: CARCEL 8
+      X- REFL X+ REFL Y- REFL Y+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.475 0.6045552
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE1>>
+      MESHX 0.0 1.26 MESHY 0.0 1.26
+          ;
+TRACK_SS := SYBILT: GEOM ::
+     TITLE 'UO2 ROWLAND S BENCHMARK'
+     MAXR 20 QUA2 20 3 DP01 MAXZ 200000 ;
+UOX_TBH := G2S: GEOM ;
+TRACK TF_EXC := SALT: UOX_TBH ::
+  EDIT 3
+  ALLG
+  TSPC <<an2d>> <<densur>> REND
+;
+*
+* Subgroup projection method
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ CALENDF 3 PT      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    U235     = U235  5.118E-5 1
+    U238     = U238  2.042E-2 1
+    PU238    = Pu238 2.714E-5 1
+    PU239    = Pu239 1.972E-3 1
+    PU240    = Pu240 4.256E-4 1
+    PU241    = Pu241 3.577E-5 1
+    PU242    = Pu242 1.234E-5 1
+    AM241    = Am241 1.406E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.274984 ;
+*
+* Tone's method
+LIBRARY CP CALC := DELETE: LIBRARY CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    U235     = U235  5.118E-5 1
+    U238     = U238  2.042E-2 1
+    PU238    = Pu238 2.714E-5 1
+    PU239    = Pu239 1.972E-3 1
+    PU240    = Pu240 4.256E-4 1
+    PU241    = Pu241 3.577E-5 1
+    PU242    = Pu242 1.234E-5 1
+    AM241    = Am241 1.406E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY := TONE: LIBRARY TRACK_SS :: EDIT 1 TRAN MXIT 3 SPH GRMIN 52 ;
+CP := ASM: LIBRARY TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.273706 ;
+*
+* Resonance spectrum expansion method
+LIBRARY LIBRARY2 CP CALC := DELETE: LIBRARY LIBRARY2 CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ RSE       (*RSE TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    U235     = U235  5.118E-5 1
+    U238     = U238  2.042E-2 1
+    PU238    = Pu238 2.714E-5 1
+    PU239    = Pu239 1.972E-3 1
+    PU240    = Pu240 4.256E-4 1
+    PU241    = Pu241 3.577E-5 1
+    PU242    = Pu242 1.234E-5 1
+    AM241    = Am241 1.406E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.275405 ;
+ECHO "test mox_295_kec3 completed" ;
+QUIT "LIST" .
diff --git a/Dragon/data/trowland_shem295_jeff3.1_proc/mox_295_kec4.c2m b/Dragon/data/trowland_shem295_jeff3.1_proc/mox_295_kec4.c2m
new file mode 100644
index 0000000..e51c7d0
--- /dev/null
+++ b/Dragon/data/trowland_shem295_jeff3.1_proc/mox_295_kec4.c2m
@@ -0,0 +1,176 @@
+*----
+*
+*  TEST CASE mox_295_kec4
+*  MOX ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  295-GROUP JEF3.1.1 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK_SS TRACK LIBRARY LIBRARY2 CP CALC ;
+SEQ_ASCII UOX_TBH ;
+SEQ_BINARY TF_EXC ;
+MODULE LIB: GEO: SYBILT: G2S: SALT: MAC: USS: TONE: ASM: FLU: DELETE: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 := 
+        1 2 3 4 5 6 7 8 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+REAL P_Zr90     P_Zr91     P_Zr92    P_Zr94    P_Zr96    :=
+     51.45E-02  11.22E-02  17.15E-02 17.38E-02 2.80E-02  ;
+REAL N_Zr0 := 3.88E-2 ;
+REAL N_Zr90       N_Zr91       N_Zr92       N_Zr94       N_Zr96    :=
+     N_Zr0 P_Zr90 *  N_Zr0 P_Zr91 *  N_Zr0 P_Zr92 *  N_Zr0 P_Zr94 *
+     N_Zr0 P_Zr96 *  ;
+INTEGER an2d := 12 ;
+REAL densur := 20.0 ;
+*
+EVALUATE RCOMB6 := 0.41 ;
+EVALUATE RCOMB1 := 0.4 SQRT RCOMB6 * ;
+EVALUATE RCOMB2 := 0.7 SQRT RCOMB6 * ;
+EVALUATE RCOMB3 := 0.8 SQRT RCOMB6 * ;
+EVALUATE RCOMB4 := 0.9 SQRT RCOMB6 * ;
+EVALUATE RCOMB5 := 0.95 SQRT RCOMB6 * ;
+GEOM := GEO: :: CARCEL 8
+      X- REFL X+ REFL Y- REFL Y+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.475 0.6045552
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE1>>
+      MESHX 0.0 1.26 MESHY 0.0 1.26
+          ;
+TRACK_SS := SYBILT: GEOM ::
+     TITLE 'UO2 ROWLAND S BENCHMARK'
+     MAXR 20 QUA2 20 3 DP01 MAXZ 200000 ;
+UOX_TBH := G2S: GEOM ;
+TRACK TF_EXC := SALT: UOX_TBH ::
+  EDIT 3
+  ALLG
+  TSPC <<an2d>> <<densur>> REND
+;
+*
+* Subgroup projection method
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ CALENDF 3 PT      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 600.0 (*COMB0101*)
+    U235     = U235  5.118E-5 1
+    U238     = U238  2.042E-2 1
+    PU238    = Pu238 2.714E-5 1
+    PU239    = Pu239 1.972E-3 1
+    PU240    = Pu240 4.256E-4 1
+    PU241    = Pu241 3.577E-5 1
+    PU242    = Pu242 1.234E-5 1
+    AM241    = Am241 1.406E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.259702 ;
+*
+* Tone's method
+LIBRARY CP CALC := DELETE: LIBRARY CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 600.0 (*COMB0101*)
+    U235     = U235  5.118E-5 1
+    U238     = U238  2.042E-2 1
+    PU238    = Pu238 2.714E-5 1
+    PU239    = Pu239 1.972E-3 1
+    PU240    = Pu240 4.256E-4 1
+    PU241    = Pu241 3.577E-5 1
+    PU242    = Pu242 1.234E-5 1
+    AM241    = Am241 1.406E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY := TONE: LIBRARY TRACK_SS :: EDIT 1 TRAN MXIT 3 SPH GRMIN 52 ;
+CP := ASM: LIBRARY TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.258466 ;
+*
+* Resonance spectrum expansion method
+LIBRARY LIBRARY2 CP CALC := DELETE: LIBRARY LIBRARY2 CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ RSE       (*RSE TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 600.0 (*COMB0101*)
+    U235     = U235  5.118E-5 1
+    U238     = U238  2.042E-2 1
+    PU238    = Pu238 2.714E-5 1
+    PU239    = Pu239 1.972E-3 1
+    PU240    = Pu240 4.256E-4 1
+    PU241    = Pu241 3.577E-5 1
+    PU242    = Pu242 1.234E-5 1
+    AM241    = Am241 1.406E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.259876 ;
+ECHO "test mox_295_kec4 completed" ;
+QUIT "LIST" .
diff --git a/Dragon/data/trowland_shem295_jeff3.1_proc/mox_295_kec5.c2m b/Dragon/data/trowland_shem295_jeff3.1_proc/mox_295_kec5.c2m
new file mode 100644
index 0000000..af6d2a1
--- /dev/null
+++ b/Dragon/data/trowland_shem295_jeff3.1_proc/mox_295_kec5.c2m
@@ -0,0 +1,311 @@
+*----
+*
+*  TEST CASE mox_295_kec5
+*  MOX ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  295-GROUP JEF3.1.1 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK_SS TRACK LIBRARY LIBRARY2 CP CALC ;
+SEQ_ASCII UOX_TBH ;
+SEQ_BINARY TF_EXC ;
+MODULE LIB: GEO: SYBILT: G2S: SALT: MAC: USS: TONE: ASM: FLU: DELETE: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 := 
+        1 2 3 4 5 6 7 8 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+REAL P_Zr90     P_Zr91     P_Zr92    P_Zr94    P_Zr96    :=
+     51.45E-02  11.22E-02  17.15E-02 17.38E-02 2.80E-02  ;
+REAL N_Zr0 := 3.88E-2 ;
+REAL N_Zr90       N_Zr91       N_Zr92       N_Zr94       N_Zr96    :=
+     N_Zr0 P_Zr90 *  N_Zr0 P_Zr91 *  N_Zr0 P_Zr92 *  N_Zr0 P_Zr94 *
+     N_Zr0 P_Zr96 *  ;
+INTEGER an2d := 12 ;
+REAL densur := 20.0 ;
+*
+EVALUATE RCOMB6 := 0.41 ;
+EVALUATE RCOMB1 := 0.4 SQRT RCOMB6 * ;
+EVALUATE RCOMB2 := 0.7 SQRT RCOMB6 * ;
+EVALUATE RCOMB3 := 0.8 SQRT RCOMB6 * ;
+EVALUATE RCOMB4 := 0.9 SQRT RCOMB6 * ;
+EVALUATE RCOMB5 := 0.95 SQRT RCOMB6 * ;
+GEOM := GEO: :: CARCEL 8
+      X- REFL X+ REFL Y- REFL Y+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.475 0.6045552
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE1>>
+      MESHX 0.0 1.26 MESHY 0.0 1.26
+          ;
+TRACK_SS := SYBILT: GEOM ::
+     TITLE 'UO2 ROWLAND S BENCHMARK'
+     MAXR 20 QUA2 20 3 DP01 MAXZ 200000 ;
+UOX_TBH := G2S: GEOM ;
+TRACK TF_EXC := SALT: UOX_TBH ::
+  EDIT 3
+  ALLG
+  TSPC <<an2d>> <<densur>> REND
+;
+*
+* Subgroup projection method (ECCO full correlation)
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ CALENDF 3 PT      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 900.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0102>> 700.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0103>> 500.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0104>> 400.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0105>> 350.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0106>> 300.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 NOCO GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.210817 ;
+*
+* Tone's method
+LIBRARY CP CALC := DELETE: LIBRARY CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 900.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0102>> 700.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0103>> 500.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0104>> 400.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0105>> 350.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0106>> 300.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY := TONE: LIBRARY TRACK_SS :: EDIT 1 TRAN MXIT 3 SPH GRMIN 52 ;
+CP := ASM: LIBRARY TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.208753 ;
+*
+* Resonance spectrum expansion method (ECCO full correlation)
+LIBRARY LIBRARY2 CP CALC := DELETE: LIBRARY LIBRARY2 CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ RSE       (*RSE TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 900.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0102>> 700.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0103>> 500.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0104>> 400.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0105>> 350.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<COMB0106>> 300.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1
+    U238     = U238  2.037E-2 1
+    PU238    = Pu238 4.669E-5 1
+    PU239    = Pu239 1.465E-3 1
+    PU240    = Pu240 5.691E-4 1
+    PU241    = Pu241 2.713E-4 1
+    PU242    = Pu242 1.413E-4 1
+    AM241    = Am241 3.028E-5 1
+    O16      = O16   4.588E-2
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 NOCO GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.209550 ;
+ECHO "test mox_295_kec5 completed" ;
+QUIT "LIST" .
diff --git a/Dragon/data/trowland_shem295_jeff3.1_proc/uo2_295_kec1.c2m b/Dragon/data/trowland_shem295_jeff3.1_proc/uo2_295_kec1.c2m
new file mode 100644
index 0000000..60fe9d6
--- /dev/null
+++ b/Dragon/data/trowland_shem295_jeff3.1_proc/uo2_295_kec1.c2m
@@ -0,0 +1,193 @@
+*----
+*
+*  TEST CASE uo2_295_kec1
+*  UO2 ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  295-GROUP JEF3.1.1 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK_SS TRACK LIBRARY LIBRARY2 CP CALC ;
+SEQ_ASCII UOX_TBH ;
+SEQ_BINARY TF_EXC ;
+MODULE LIB: GEO: SYBILT: G2S: SALT: MAC: USS: TONE: AUTO: ASM: FLU: DELETE: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 := 
+        1 2 3 4 5 6 7 8 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+REAL P_Zr90     P_Zr91     P_Zr92    P_Zr94    P_Zr96    :=
+     51.45E-02  11.22E-02  17.15E-02 17.38E-02 2.80E-02  ;
+REAL N_Zr0 := 4.3241E-2 ;
+REAL N_Zr90       N_Zr91       N_Zr92       N_Zr94       N_Zr96    :=
+     N_Zr0 P_Zr90 *  N_Zr0 P_Zr91 *  N_Zr0 P_Zr92 *  N_Zr0 P_Zr94 *
+     N_Zr0 P_Zr96 *  ;
+INTEGER an2d := 12 ;
+REAL densur := 20.0 ;
+*
+EVALUATE RCOMB6 := 0.4 ;
+EVALUATE RCOMB1 := 0.4 SQRT RCOMB6 * ;
+EVALUATE RCOMB2 := 0.7 SQRT RCOMB6 * ;
+EVALUATE RCOMB3 := 0.8 SQRT RCOMB6 * ;
+EVALUATE RCOMB4 := 0.9 SQRT RCOMB6 * ;
+EVALUATE RCOMB5 := 0.95 SQRT RCOMB6 * ;
+GEOM := GEO: :: CARCEL 8
+      X- REFL X+ REFL Y- REFL Y+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.45 0.5748331
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE1>>
+      MESHX 0.0 1.2 MESHY 0.0 1.2
+          ;
+TRACK_SS := SYBILT: GEOM ::
+     TITLE 'UO2 ROWLAND S BENCHMARK'
+     MAXR 20 QUA2 20 3 DP01 MAXZ 200000 ;
+UOX_TBH := G2S: GEOM ;
+TRACK TF_EXC := SALT: UOX_TBH ::
+  EDIT 3
+  ALLG
+  TSPC <<an2d>> <<densur>> REND
+;
+*
+* Subgroup projection method
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ CALENDF 3 PT (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 6.6988E-2
+    O16     = O16    3.3494E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.388247 ;
+*
+* Tone's method
+LIBRARY CP CALC := DELETE: LIBRARY CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 6.6988E-2
+    O16     = O16    3.3494E-2
+ ;
+LIBRARY := TONE: LIBRARY TRACK_SS :: EDIT 1 TRAN MXIT 3 SPH GRMIN 52 ;
+CP := ASM: LIBRARY TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.387176 ;
+*
+* Resonance spectrum expansion method
+LIBRARY LIBRARY2 CP CALC := DELETE: LIBRARY LIBRARY2 CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ RSE       (*RSE TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 6.6988E-2
+    O16     = O16    3.3494E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.388969 ;
+*
+* Autosecol UFG method
+LIBRARY LIBRARY2 CP CALC := DELETE: LIBRARY LIBRARY2 CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8     (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL  (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ SUBG       (*PHYSICAL TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 1
+    Zr91    = Zr91   <<N_Zr91>> 1
+    Zr92    = Zr92   <<N_Zr92>> 1
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 6.6988E-2
+    O16     = O16    3.3494E-2
+ ;
+LIBRARY2 := AUTO: LIBRARY TRACK_SS :: EDIT 3 GRMIN 52 MAXT 60000 SEED 123456 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.387987 ;
+ECHO "test uo2_295_kec1 completed" ;
+QUIT "LIST" .
diff --git a/Dragon/data/trowland_shem295_jeff3.1_proc/uo2_295_kec2.c2m b/Dragon/data/trowland_shem295_jeff3.1_proc/uo2_295_kec2.c2m
new file mode 100644
index 0000000..d5b91d6
--- /dev/null
+++ b/Dragon/data/trowland_shem295_jeff3.1_proc/uo2_295_kec2.c2m
@@ -0,0 +1,158 @@
+*----
+*
+*  TEST CASE uo2_295_kec2
+*  UO2 ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  295-GROUP JEF3.1.1 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK_SS TRACK LIBRARY LIBRARY2 CP CALC ;
+SEQ_ASCII UOX_TBH ;
+SEQ_BINARY TF_EXC ;
+MODULE LIB: GEO: SYBILT: G2S: SALT: MAC: USS: TONE: ASM: FLU: DELETE: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 := 
+        1 2 3 4 5 6 7 8 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+REAL P_Zr90     P_Zr91     P_Zr92    P_Zr94    P_Zr96    :=
+     51.45E-02  11.22E-02  17.15E-02 17.38E-02 2.80E-02  ;
+REAL N_Zr0 := 4.3241E-2 ;
+REAL N_Zr90       N_Zr91       N_Zr92       N_Zr94       N_Zr96    :=
+     N_Zr0 P_Zr90 *  N_Zr0 P_Zr91 *  N_Zr0 P_Zr92 *  N_Zr0 P_Zr94 *
+     N_Zr0 P_Zr96 *  ;
+INTEGER an2d := 12 ;
+REAL densur := 20.0 ;
+*
+EVALUATE RCOMB6 := 0.4 ;
+EVALUATE RCOMB1 := 0.4 SQRT RCOMB6 * ;
+EVALUATE RCOMB2 := 0.7 SQRT RCOMB6 * ;
+EVALUATE RCOMB3 := 0.8 SQRT RCOMB6 * ;
+EVALUATE RCOMB4 := 0.9 SQRT RCOMB6 * ;
+EVALUATE RCOMB5 := 0.95 SQRT RCOMB6 * ;
+GEOM := GEO: :: CARCEL 8
+      X- REFL X+ REFL Y- REFL Y+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.45 0.5748331
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE1>>
+      MESHX 0.0 1.2 MESHY 0.0 1.2
+          ;
+TRACK_SS := SYBILT: GEOM ::
+     TITLE 'UO2 ROWLAND S BENCHMARK'
+     MAXR 20 QUA2 20 3 DP01 MAXZ 200000 ;
+UOX_TBH := G2S: GEOM ;
+TRACK TF_EXC := SALT: UOX_TBH ::
+  EDIT 3
+  ALLG
+  TSPC <<an2d>> <<densur>> REND
+;
+*
+* Subgroup projection method
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ CALENDF 3 PT (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.335245 ;
+*
+* Tone's method
+LIBRARY CP CALC := DELETE: LIBRARY CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY := TONE: LIBRARY TRACK_SS :: EDIT 1 TRAN MXIT 3 SPH GRMIN 52 ;
+CP := ASM: LIBRARY TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.334218 ;
+*
+* Resonance spectrum expansion method
+LIBRARY LIBRARY2 CP CALC := DELETE: LIBRARY LIBRARY2 CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ RSE       (*RSE TYPE PROBABILITY TABLES*)
+**
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 300.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.335770 ;
+ECHO "test uo2_295_kec2 completed" ;
+QUIT "LIST" .
diff --git a/Dragon/data/trowland_shem295_jeff3.1_proc/uo2_295_kec3.c2m b/Dragon/data/trowland_shem295_jeff3.1_proc/uo2_295_kec3.c2m
new file mode 100644
index 0000000..363c0fa
--- /dev/null
+++ b/Dragon/data/trowland_shem295_jeff3.1_proc/uo2_295_kec3.c2m
@@ -0,0 +1,158 @@
+*----
+*
+*  TEST CASE uo2_295_kec3
+*  UO2 ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  295-GROUP JEF3.1.1 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK_SS TRACK LIBRARY LIBRARY2 CP CALC ;
+SEQ_ASCII UOX_TBH ;
+SEQ_BINARY TF_EXC ;
+MODULE LIB: GEO: SYBILT: G2S: SALT: MAC: USS: TONE: ASM: FLU: DELETE: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 := 
+        1 2 3 4 5 6 7 8 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+REAL P_Zr90     P_Zr91     P_Zr92    P_Zr94    P_Zr96    :=
+     51.45E-02  11.22E-02  17.15E-02 17.38E-02 2.80E-02  ;
+REAL N_Zr0 := 4.3241E-2 ;
+REAL N_Zr90       N_Zr91       N_Zr92       N_Zr94       N_Zr96    :=
+     N_Zr0 P_Zr90 *  N_Zr0 P_Zr91 *  N_Zr0 P_Zr92 *  N_Zr0 P_Zr94 *
+     N_Zr0 P_Zr96 *  ;
+INTEGER an2d := 12 ;
+REAL densur := 20.0 ;
+*
+EVALUATE RCOMB6 := 0.4 ;
+EVALUATE RCOMB1 := 0.4 SQRT RCOMB6 * ;
+EVALUATE RCOMB2 := 0.7 SQRT RCOMB6 * ;
+EVALUATE RCOMB3 := 0.8 SQRT RCOMB6 * ;
+EVALUATE RCOMB4 := 0.9 SQRT RCOMB6 * ;
+EVALUATE RCOMB5 := 0.95 SQRT RCOMB6 * ;
+GEOM := GEO: :: CARCEL 8
+      X- REFL X+ REFL Y- REFL Y+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.45 0.5748331
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE1>>
+      MESHX 0.0 1.2 MESHY 0.0 1.2
+          ;
+TRACK_SS := SYBILT: GEOM ::
+     TITLE 'UO2 ROWLAND S BENCHMARK'
+     MAXR 20 QUA2 20 3 DP01 MAXZ 200000 ;
+UOX_TBH := G2S: GEOM ;
+TRACK TF_EXC := SALT: UOX_TBH ::
+  EDIT 3
+  ALLG
+  TSPC <<an2d>> <<densur>> REND
+;
+*
+* Subgroup projection method
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ CALENDF 3 PT      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 900.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 600.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 600.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.302678 ;
+*
+* Tone's method
+LIBRARY CP CALC := DELETE: LIBRARY CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 900.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 600.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 600.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY := TONE: LIBRARY TRACK_SS :: EDIT 1 TRAN MXIT 3 SPH GRMIN 52 ;
+CP := ASM: LIBRARY TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.300903 ;
+*
+* Resonance spectrum expansion method
+LIBRARY LIBRARY2 CP CALC := DELETE: LIBRARY LIBRARY2 CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ RSE       (*RSE TYPE PROBABILITY TABLES*)
+**
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 900.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 600.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 600.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.303022 ;
+ECHO "test uo2_295_kec3 completed" ;
+QUIT "LIST" .
diff --git a/Dragon/data/trowland_shem295_jeff3.1_proc/uo2_295_kec4.c2m b/Dragon/data/trowland_shem295_jeff3.1_proc/uo2_295_kec4.c2m
new file mode 100644
index 0000000..45a349d
--- /dev/null
+++ b/Dragon/data/trowland_shem295_jeff3.1_proc/uo2_295_kec4.c2m
@@ -0,0 +1,158 @@
+*----
+*
+*  TEST CASE uo2_295_kec4
+*  UO2 ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  295-GROUP JEF3.1.1 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK_SS TRACK LIBRARY LIBRARY2 CP CALC ;
+SEQ_ASCII UOX_TBH ;
+SEQ_BINARY TF_EXC ;
+MODULE LIB: GEO: SYBILT: G2S: SALT: MAC: USS: TONE: ASM: FLU: DELETE: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 := 
+        1 2 3 4 5 6 7 8 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+REAL P_Zr90     P_Zr91     P_Zr92    P_Zr94    P_Zr96    :=
+     51.45E-02  11.22E-02  17.15E-02 17.38E-02 2.80E-02  ;
+REAL N_Zr0 := 4.3241E-2 ;
+REAL N_Zr90       N_Zr91       N_Zr92       N_Zr94       N_Zr96    :=
+     N_Zr0 P_Zr90 *  N_Zr0 P_Zr91 *  N_Zr0 P_Zr92 *  N_Zr0 P_Zr94 *
+     N_Zr0 P_Zr96 *  ;
+INTEGER an2d := 12 ;
+REAL densur := 20.0 ;
+*
+EVALUATE RCOMB6 := 0.4 ;
+EVALUATE RCOMB1 := 0.4 SQRT RCOMB6 * ;
+EVALUATE RCOMB2 := 0.7 SQRT RCOMB6 * ;
+EVALUATE RCOMB3 := 0.8 SQRT RCOMB6 * ;
+EVALUATE RCOMB4 := 0.9 SQRT RCOMB6 * ;
+EVALUATE RCOMB5 := 0.95 SQRT RCOMB6 * ;
+GEOM := GEO: :: CARCEL 8
+      X- REFL X+ REFL Y- REFL Y+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.45 0.5748331
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE1>>
+      MESHX 0.0 1.2 MESHY 0.0 1.2
+          ;
+TRACK_SS := SYBILT: GEOM ::
+     TITLE 'UO2 ROWLAND S BENCHMARK'
+     MAXR 20 QUA2 20 3 DP01 MAXZ 200000 ;
+UOX_TBH := G2S: GEOM ;
+TRACK TF_EXC := SALT: UOX_TBH ::
+  EDIT 3
+  ALLG
+  TSPC <<an2d>> <<densur>> REND
+;
+*
+* Subgroup projection method
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ CALENDF 3 PT      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 600.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 600.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 600.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.314612 ;
+*
+* Tone's method
+LIBRARY CP CALC := DELETE: LIBRARY CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 600.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 600.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 600.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY := TONE: LIBRARY TRACK_SS :: EDIT 1 TRAN MXIT 3 SPH GRMIN 52 ;
+CP := ASM: LIBRARY TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.313378 ;
+*
+* Resonance spectrum expansion method
+LIBRARY LIBRARY2 CP CALC := DELETE: LIBRARY LIBRARY2 CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ RSE       (*RSE TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 600.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*)
+ MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*)
+ MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*)
+ MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*)
+ MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*)
+ MIX <<GAIN1>> 600.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 600.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 3 GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.314974 ;
+ECHO "test uo2_295_kec4 completed" ;
+QUIT "LIST" .
diff --git a/Dragon/data/trowland_shem295_jeff3.1_proc/uo2_295_kec5.c2m b/Dragon/data/trowland_shem295_jeff3.1_proc/uo2_295_kec5.c2m
new file mode 100644
index 0000000..31e3401
--- /dev/null
+++ b/Dragon/data/trowland_shem295_jeff3.1_proc/uo2_295_kec5.c2m
@@ -0,0 +1,203 @@
+*----
+*
+*  TEST CASE uo2_295_kec5
+*  UO2 ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  295-GROUP JEF3.1.1 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK_SS TRACK LIBRARY LIBRARY2 CP CALC ;
+SEQ_ASCII UOX_TBH ;
+SEQ_BINARY TF_EXC ;
+MODULE LIB: GEO: SYBILT: G2S: SALT: MAC: USS: TONE: ASM: FLU: DELETE: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 := 
+        1 2 3 4 5 6 7 8 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+REAL P_Zr90     P_Zr91     P_Zr92    P_Zr94    P_Zr96    :=
+     51.45E-02  11.22E-02  17.15E-02 17.38E-02 2.80E-02  ;
+REAL N_Zr0 := 4.3241E-2 ;
+REAL N_Zr90       N_Zr91       N_Zr92       N_Zr94       N_Zr96    :=
+     N_Zr0 P_Zr90 *  N_Zr0 P_Zr91 *  N_Zr0 P_Zr92 *  N_Zr0 P_Zr94 *
+     N_Zr0 P_Zr96 *  ;
+INTEGER an2d := 12 ;
+REAL densur := 20.0 ;
+*
+EVALUATE RCOMB6 := 0.4 ;
+EVALUATE RCOMB1 := 0.4 SQRT RCOMB6 * ;
+EVALUATE RCOMB2 := 0.7 SQRT RCOMB6 * ;
+EVALUATE RCOMB3 := 0.8 SQRT RCOMB6 * ;
+EVALUATE RCOMB4 := 0.9 SQRT RCOMB6 * ;
+EVALUATE RCOMB5 := 0.95 SQRT RCOMB6 * ;
+GEOM := GEO: :: CARCEL 8
+      X- REFL X+ REFL Y- REFL Y+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.45 0.5748331
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE1>>
+      MESHX 0.0 1.2 MESHY 0.0 1.2
+          ;
+TRACK_SS := SYBILT: GEOM ::
+     TITLE 'UO2 ROWLAND S BENCHMARK'
+     MAXR 20 QUA2 20 3 DP01 MAXZ 200000 ;
+UOX_TBH := G2S: GEOM ;
+TRACK TF_EXC := SALT: UOX_TBH ::
+  EDIT 3
+  ALLG
+  TSPC <<an2d>> <<densur>> REND
+;
+*
+* Subgroup projection method (ECCO full correlation)
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORTECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ CALENDF 3 PT      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 1200.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0102>> 1000.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0103>> 800.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0104>> 700.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0105>> 650.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0106>> 600.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<GAIN1>> 600.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 600.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 4 NOCO GRMIN 52 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.301416 ;
+*
+* Tone's method
+LIBRARY CP CALC := DELETE: LIBRARY CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 1200.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0102>> 1000.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0103>> 800.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0104>> 700.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0105>> 650.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0106>> 600.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<GAIN1>> 600.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 600.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY := TONE: LIBRARY TRACK_SS :: EDIT 1 TRAN MXIT 3 SPH GRMIN 52 ;
+CP := ASM: LIBRARY TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.298397 ;
+*
+* Resonance spectrum expansion method (ECCO full correlation)
+LIBRARY LIBRARY2 CP CALC := DELETE: LIBRARY LIBRARY2 CP CALC ;
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 8    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORTECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ RSE       (*RSE TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_295
+ MIX <<COMB0101>> 1200.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0102>> 1000.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0103>> 800.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0104>> 700.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0105>> 650.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<COMB0106>> 600.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1
+    U238    = U238   2.2604E-2 1
+ MIX <<GAIN1>> 600.0 (*GAIN1*)
+    Zr90    = Zr90   <<N_Zr90>> 2
+    Zr91    = Zr91   <<N_Zr91>> 2
+    Zr92    = Zr92   <<N_Zr92>> 2
+    Zr94    = Zr94   <<N_Zr94>>
+    Zr96    = Zr96   <<N_Zr96>>
+ MIX <<MODE1>> 600.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK_SS :: EDIT 1 TRAN PASS 4 NOCO GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK TF_EXC :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.299549 ;
+ECHO "test uo2_295_kec5 completed" ;
+QUIT "LIST" .