Initial commit from Polytechnique Montreal

10 files changed, 1750 insertions, 0 deletions

diff --git a/Dragon/data/trowland_shem_proc/mox_6c_c1_sh.c2m b/Dragon/data/trowland_shem_proc/mox_6c_c1_sh.c2m
new file mode 100644
index 0000000..7ae2f0e
--- /dev/null
+++ b/Dragon/data/trowland_shem_proc/mox_6c_c1_sh.c2m
@@ -0,0 +1,205 @@
+*----
+*
+*  TEST CASE mox_6c_c1_sh
+*  FIRST MOX ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  281-GROUP JEF2.2 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK LIBRARY LIBRARY2 CP FLUX ;
+MODULE LIB: GEO: SYBILT: SHI: USS: ASM: FLU: DELETE: END: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 MODE2 := 
+        1 2 3 4 5 6 7 8 9 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+*
+EVALUATE RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 := 0.259307 
+0.343031 0.366715 0.388960 0.3996186 0.41 ;
+GEOM := GEO: :: TUBE 9
+      R+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.475 0.6045552 0.710879
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE2>>
+          ;
+TRACK := SYBILT: GEOM ::
+     TITLE 'FIRST MOX ROWLAND S BENCHMARK'
+     MAXR 20 QUA1 5 ;
+*----
+* SHI LEVEL 1
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ 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*)
+    ZR      = Zr0   3.88E-2 2
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 56 LEVEL 1 MXIT 3 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.219334 ;
+LIBRARY2 CP FLUX := DELETE: LIBRARY2 CP FLUX ;
+*----
+* SHI LEVEL 2
+*----
+LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 56 LEVEL 2 MXIT 3 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.219380 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* USS
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ SUBG      (*HELIOS TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ 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*)
+    ZR      = Zr0   3.88E-2 2
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 3 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.218508 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* Ribon extended
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ CALENDF 4 PTSL      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1 IRSET PTSL 1
+    U238     = U238  2.037E-2 1 IRSET PTSL 1
+    PU238    = Pu238 4.669E-5 1 IRSET PTSL 1
+    PU239    = Pu239 1.465E-3 1 IRSET PTSL 1
+    PU240    = Pu240 5.691E-4 1 IRSET PTSL 1
+    PU241    = Pu241 2.713E-4 1 IRSET PTSL 1
+    PU242    = Pu242 1.413E-4 1 IRSET PTSL 1
+    AM241    = Am241 3.028E-5 1 IRSET PTSL 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*)
+    ZR      = Zr0   3.88E-2 2 IRSET PTSL 1
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 3 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.218379 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* Ribon extended with correlation
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ CALENDF 4 PTSL      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1 IRSET PTSL 1
+    U238     = U238  2.037E-2 CORR 1 IRSET PTSL 1
+    PU238    = Pu238 4.669E-5 1 IRSET PTSL 1
+    PU239    = Pu239 1.465E-3 1 IRSET PTSL 1
+    PU240    = Pu240 5.691E-4 CORR 1 IRSET PTSL 1
+    PU241    = Pu241 2.713E-4 1 IRSET PTSL 1
+    PU242    = Pu242 1.413E-4 1 IRSET PTSL 1
+    AM241    = Am241 3.028E-5 1 IRSET PTSL 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*)
+    ZR      = Zr0   3.88E-2 2 IRSET PTSL 1
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 4 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.219567 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+ECHO "test mox_6c_c1_sh completed" ;
+END: ;
diff --git a/Dragon/data/trowland_shem_proc/mox_6c_c2_sh.c2m b/Dragon/data/trowland_shem_proc/mox_6c_c2_sh.c2m
new file mode 100644
index 0000000..9bc9f71
--- /dev/null
+++ b/Dragon/data/trowland_shem_proc/mox_6c_c2_sh.c2m
@@ -0,0 +1,205 @@
+*----
+*
+*  TEST CASE mox_6c_c2_sh
+*  SECOND MOX ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  281-GROUP JEF2.2 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK LIBRARY LIBRARY2 CP FLUX ;
+MODULE LIB: GEO: SYBILT: SHI: USS: ASM: FLU: DELETE: END: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 MODE2 := 
+        1 2 3 4 5 6 7 8 9 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+*
+EVALUATE RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 := 0.259307 
+0.343031 0.366715 0.388960 0.3996186 0.41 ;
+GEOM := GEO: :: TUBE 9
+      R+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.475 0.6045552 0.710879
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE2>>
+          ;
+TRACK := SYBILT: GEOM ::
+     TITLE 'SECOND MOX ROWLAND S BENCHMARK'
+     MAXR 20 QUA1 5 ;
+*----
+* SHI LEVEL 1
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 560.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*)
+    ZR      = Zr0  3.88E-2 2
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 56 LEVEL 1 MXIT 3 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.204594 ;
+LIBRARY2 CP FLUX := DELETE: LIBRARY2 CP FLUX ;
+*----
+* SHI LEVEL 2
+*----
+LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 56 LEVEL 2 MXIT 3 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.204738 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* USS
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ SUBG      (*HELIOS TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 560.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*)
+    ZR      = Zr0  3.88E-2 2
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 3 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.203734 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* Ribon extended
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ PTSL      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 560.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1 IRSET PTSL 1
+    U238     = U238  2.037E-2 1 IRSET PTSL 1
+    PU238    = Pu238 4.669E-5 1 IRSET PTSL 1
+    PU239    = Pu239 1.465E-3 1 IRSET PTSL 1
+    PU240    = Pu240 5.691E-4 1 IRSET PTSL 1
+    PU241    = Pu241 2.713E-4 1 IRSET PTSL 1
+    PU242    = Pu242 1.413E-4 1 IRSET PTSL 1
+    AM241    = Am241 3.028E-5 1 IRSET PTSL 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*)
+    ZR      = Zr0  3.88E-2 2 IRSET PTSL 1
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 3 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.203946 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* Ribon extended with correlation
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ PTSL      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 560.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1 IRSET PTSL 1
+    U238     = U238  2.037E-2 CORR 1 IRSET PTSL 1
+    PU238    = Pu238 4.669E-5 1 IRSET PTSL 1
+    PU239    = Pu239 1.465E-3 1 IRSET PTSL 1
+    PU240    = Pu240 5.691E-4 CORR 1 IRSET PTSL 1
+    PU241    = Pu241 2.713E-4 1 IRSET PTSL 1
+    PU242    = Pu242 1.413E-4 1 IRSET PTSL 1
+    AM241    = Am241 3.028E-5 1 IRSET PTSL 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*)
+    ZR      = Zr0  3.88E-2 2 IRSET PTSL 1
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 2 TRAN PASS 4 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.205297 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+ECHO "test mox_6c_c2_sh completed" ;
+END: ;
diff --git a/Dragon/data/trowland_shem_proc/mox_6c_c3_sh.c2m b/Dragon/data/trowland_shem_proc/mox_6c_c3_sh.c2m
new file mode 100644
index 0000000..5575615
--- /dev/null
+++ b/Dragon/data/trowland_shem_proc/mox_6c_c3_sh.c2m
@@ -0,0 +1,205 @@
+*----
+*
+*  TEST CASE mox_6c_c3_sh
+*  THIRD MOX ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  281-GROUP JEF2.2 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK LIBRARY LIBRARY2 CP FLUX ;
+MODULE LIB: GEO: SYBILT: SHI: USS: ASM: FLU: DELETE: END: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 MODE2 := 
+        1 2 3 4 5 6 7 8 9 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+*
+EVALUATE RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 := 0.259307 
+0.343031 0.366715 0.388960 0.3996186 0.41 ;
+GEOM := GEO: :: TUBE 9
+      R+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.475 0.6045552 0.710879
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE2>>
+          ;
+TRACK := SYBILT: GEOM ::
+     TITLE 'THIRD MOX ROWLAND S BENCHMARK'
+     MAXR 20 QUA1 5 ;
+*----
+* SHI LEVEL 1
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    U235     = U235  5.118E-5 1
+    U238     = U238  2.042E-2 1
+    PU239    = Pu239 1.972E-3 1
+    PU238    = Pu238 2.714E-5 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*)
+    ZR      = Zr0  3.88E-2 2
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 56 LEVEL 1 MXIT 3 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.261146 ;
+LIBRARY2 CP FLUX := DELETE: LIBRARY2 CP FLUX ;
+*----
+* SHI LEVEL 2
+*----
+LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 56 LEVEL 2 MXIT 3 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.261594 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* USS
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ SUBG      (*HELIOS TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    U235     = U235  5.118E-5 1
+    U238     = U238  2.042E-2 1
+    PU239    = Pu239 1.972E-3 1
+    PU238    = Pu238 2.714E-5 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*)
+    ZR      = Zr0  3.88E-2 2
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 3 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.260683 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* Ribon extended
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ PTSL      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    U235     = U235  5.118E-5 1 IRSET PTSL 1
+    U238     = U238  2.042E-2 1 IRSET PTSL 1
+    PU239    = Pu239 1.972E-3 1 IRSET PTSL 1
+    PU238    = Pu238 2.714E-5 1 IRSET PTSL 1
+    PU240    = Pu240 4.256E-4 1 IRSET PTSL 1
+    PU241    = Pu241 3.577E-5 1 IRSET PTSL 1
+    PU242    = Pu242 1.234E-5 1 IRSET PTSL 1
+    AM241    = Am241 1.406E-5 1 IRSET PTSL 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*)
+    ZR      = Zr0  3.88E-2 2 IRSET PTSL 1
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 3 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.260568 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* Ribon extended with correlation
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ PTSL      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 300.0 (*COMB0101*)
+    U235     = U235  5.118E-5 1 IRSET PTSL 1
+    U238     = U238  2.042E-2 CORR 1 IRSET PTSL 1
+    PU238    = Pu238 2.714E-5 1 IRSET PTSL 1
+    PU239    = Pu239 1.972E-3 1 IRSET PTSL 1
+    PU240    = Pu240 4.256E-4 CORR 1 IRSET PTSL 1
+    PU241    = Pu241 3.577E-5 1 IRSET PTSL 1
+    PU242    = Pu242 1.234E-5 1 IRSET PTSL 1
+    AM241    = Am241 1.406E-5 1 IRSET PTSL 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*)
+    ZR      = Zr0  3.88E-2 2 IRSET PTSL 1
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 4 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.261677 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+ECHO "test mox_6c_c3_sh completed" ;
+END: ;
diff --git a/Dragon/data/trowland_shem_proc/mox_6c_c4_sh.c2m b/Dragon/data/trowland_shem_proc/mox_6c_c4_sh.c2m
new file mode 100644
index 0000000..eaea149
--- /dev/null
+++ b/Dragon/data/trowland_shem_proc/mox_6c_c4_sh.c2m
@@ -0,0 +1,205 @@
+*----
+*
+*  TEST CASE mox_6c_c4_sh
+*  FOURTH MOX ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  281-GROUP JEF2.2 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK LIBRARY LIBRARY2 CP FLUX ;
+MODULE LIB: GEO: SYBILT: SHI: USS: ASM: FLU: DELETE: END: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 MODE2 := 
+        1 2 3 4 5 6 7 8 9 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+*
+EVALUATE RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 := 0.259307 
+0.343031 0.366715 0.388960 0.3996186 0.41 ;
+GEOM := GEO: :: TUBE 9
+      R+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.475 0.6045552 0.710879
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE2>>
+          ;
+TRACK := SYBILT: GEOM ::
+     TITLE 'FOURTH MOX ROWLAND S BENCHMARK'
+     MAXR 20 QUA1 5 ;
+*----
+* SHI LEVEL 1
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 560.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*)
+    ZR      = Zr0  3.88E-2 2
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 56 LEVEL 1 MXIT 3 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.247249 ;
+LIBRARY2 CP FLUX := DELETE: LIBRARY2 CP FLUX ;
+*----
+* SHI LEVEL 2
+*----
+LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 56 LEVEL 2 MXIT 3 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.247834 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* USS
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ SUBG      (*HELIOS TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 560.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*)
+    ZR      = Zr0  3.88E-2 2
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 3 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.246791 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* Ribon extended
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ PTSL      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 560.0 (*COMB0101*)
+    U235     = U235  5.118E-5 1 IRSET PTSL 1
+    U238     = U238  2.042E-2 1 IRSET PTSL 1
+    PU238    = Pu238 2.714E-5 1 IRSET PTSL 1
+    PU239    = Pu239 1.972E-3 1 IRSET PTSL 1
+    PU240    = Pu240 4.256E-4 1 IRSET PTSL 1
+    PU241    = Pu241 3.577E-5 1 IRSET PTSL 1
+    PU242    = Pu242 1.234E-5 1 IRSET PTSL 1
+    AM241    = Am241 1.406E-5 1 IRSET PTSL 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*)
+    ZR      = Zr0  3.88E-2 2 IRSET PTSL 1
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 3 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.247006 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* Ribon extended with correlation
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ PTSL      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 560.0 (*COMB0101*)
+    U235     = U235  5.118E-5 1 IRSET PTSL 1
+    U238     = U238  2.042E-2 CORR 1 IRSET PTSL 1
+    PU238    = Pu238 2.714E-5 1 IRSET PTSL 1
+    PU239    = Pu239 1.972E-3 1 IRSET PTSL 1
+    PU240    = Pu240 4.256E-4 CORR 1 IRSET PTSL 1
+    PU241    = Pu241 3.577E-5 1 IRSET PTSL 1
+    PU242    = Pu242 1.234E-5 1 IRSET PTSL 1
+    AM241    = Am241 1.406E-5 1 IRSET PTSL 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*)
+    ZR      = Zr0  3.88E-2 2 IRSET PTSL 1
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 4 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.248255 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+ECHO "test mox_6c_c4_sh completed" ;
+END: ;
diff --git a/Dragon/data/trowland_shem_proc/mox_6c_c5_sh.c2m b/Dragon/data/trowland_shem_proc/mox_6c_c5_sh.c2m
new file mode 100644
index 0000000..fe4e8cb
--- /dev/null
+++ b/Dragon/data/trowland_shem_proc/mox_6c_c5_sh.c2m
@@ -0,0 +1,119 @@
+*----
+*
+*  TEST CASE mox_6c_c5_sh
+*  FIFTH MOX ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  281-GROUP SHEM DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK LIBRARY LIBRARY2 CP CALC ;
+MODULE LIB: GEO: SYBILT: USS: ASM: FLU: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 MODE2 := 
+        1 2 3 4 5 6 7 8 9 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+*
+EVALUATE RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 := 0.259307 
+0.343031 0.366715 0.388960 0.3996186 0.41 ;
+GEOM := GEO: :: TUBE 9
+      R+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.475 0.6045552 0.710879
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE2>>
+          ;
+TRACK := SYBILT: GEOM ::
+     TITLE 'FIFTH MOX ROWLAND S BENCHMARK'
+     MAXR 20 QUA1 5 ;
+*----
+* USS
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ SUBG
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 900.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1 IRSET 0.0 NONE
+    U238     = U238  2.037E-2 1 IRSET 0.0 NONE
+    PU238    = Pu238 4.669E-5 1 IRSET 0.0 NONE
+    PU239    = Pu239 1.465E-3 1 IRSET 0.0 NONE
+    PU240    = Pu240 5.691E-4 1 IRSET 0.0 NONE
+    PU241    = Pu241 2.713E-4 1 IRSET 0.0 NONE
+    PU242    = Pu242 1.413E-4 1 IRSET 0.0 NONE
+    AM241    = Am241 3.028E-5 1 IRSET 0.0 NONE
+    O16      = O16   4.588E-2
+ MIX <<COMB0102>> 700.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1 IRSET 0.0 NONE
+    U238     = U238  2.037E-2 1 IRSET 0.0 NONE
+    PU238    = Pu238 4.669E-5 1 IRSET 0.0 NONE
+    PU239    = Pu239 1.465E-3 1 IRSET 0.0 NONE
+    PU240    = Pu240 5.691E-4 1 IRSET 0.0 NONE
+    PU241    = Pu241 2.713E-4 1 IRSET 0.0 NONE
+    PU242    = Pu242 1.413E-4 1 IRSET 0.0 NONE
+    AM241    = Am241 3.028E-5 1 IRSET 0.0 NONE
+    O16      = O16   4.588E-2
+ MIX <<COMB0103>> 500.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1 IRSET 0.0 NONE
+    U238     = U238  2.037E-2 1 IRSET 0.0 NONE
+    PU238    = Pu238 4.669E-5 1 IRSET 0.0 NONE
+    PU239    = Pu239 1.465E-3 1 IRSET 0.0 NONE
+    PU240    = Pu240 5.691E-4 1 IRSET 0.0 NONE
+    PU241    = Pu241 2.713E-4 1 IRSET 0.0 NONE
+    PU242    = Pu242 1.413E-4 1 IRSET 0.0 NONE
+    AM241    = Am241 3.028E-5 1 IRSET 0.0 NONE
+    O16      = O16   4.588E-2
+ MIX <<COMB0104>> 400.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1 IRSET 0.0 NONE
+    U238     = U238  2.037E-2 1 IRSET 0.0 NONE
+    PU238    = Pu238 4.669E-5 1 IRSET 0.0 NONE
+    PU239    = Pu239 1.465E-3 1 IRSET 0.0 NONE
+    PU240    = Pu240 5.691E-4 1 IRSET 0.0 NONE
+    PU241    = Pu241 2.713E-4 1 IRSET 0.0 NONE
+    PU242    = Pu242 1.413E-4 1 IRSET 0.0 NONE
+    AM241    = Am241 3.028E-5 1 IRSET 0.0 NONE
+    O16      = O16   4.588E-2
+ MIX <<COMB0105>> 350.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1 IRSET 0.0 NONE
+    U238     = U238  2.037E-2 1 IRSET 0.0 NONE
+    PU238    = Pu238 4.669E-5 1 IRSET 0.0 NONE
+    PU239    = Pu239 1.465E-3 1 IRSET 0.0 NONE
+    PU240    = Pu240 5.691E-4 1 IRSET 0.0 NONE
+    PU241    = Pu241 2.713E-4 1 IRSET 0.0 NONE
+    PU242    = Pu242 1.413E-4 1 IRSET 0.0 NONE
+    AM241    = Am241 3.028E-5 1 IRSET 0.0 NONE
+    O16      = O16   4.588E-2
+ MIX <<COMB0106>> 300.0 (*COMB0101*)
+    U235     = U235  5.105E-5 1 IRSET 0.0 NONE
+    U238     = U238  2.037E-2 1 IRSET 0.0 NONE
+    PU238    = Pu238 4.669E-5 1 IRSET 0.0 NONE
+    PU239    = Pu239 1.465E-3 1 IRSET 0.0 NONE
+    PU240    = Pu240 5.691E-4 1 IRSET 0.0 NONE
+    PU241    = Pu241 2.713E-4 1 IRSET 0.0 NONE
+    PU242    = Pu242 1.413E-4 1 IRSET 0.0 NONE
+    AM241    = Am241 3.028E-5 1 IRSET 0.0 NONE
+    O16      = O16   4.588E-2
+ MIX <<GAIN1>> 300.0 (*GAIN1*)
+    ZR      = Zr0   3.88E-2 2 IRSET 0.0 NONE
+ MIX <<MODE1>> 293.6 (*MODE1*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ MIX <<MODE2>> 293.6 (*MODE2*)
+    H1      = H1_H2O 4.744E-2
+    O16     = O16 2.372E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 5 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.200450 ;
+ECHO "test mox_6c_c5_sh completed" ;
+QUIT "LIST" .
diff --git a/Dragon/data/trowland_shem_proc/uo2_6c_c1_sh.c2m b/Dragon/data/trowland_shem_proc/uo2_6c_c1_sh.c2m
new file mode 100644
index 0000000..ddabbfd
--- /dev/null
+++ b/Dragon/data/trowland_shem_proc/uo2_6c_c1_sh.c2m
@@ -0,0 +1,182 @@
+*----
+*
+*  TEST CASE uo2_6c_c1_sh
+*  FIRST UO2 ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  281-GROUP JEF2.2 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK LIBRARY LIBRARY2 CP FLUX ;
+MODULE LIB: GEO: SYBILT: SHI: USS: ASM: FLU: DELETE: END: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 MODE2 := 
+        1 2 3 4 5 6 7 8 9 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+*
+EVALUATE RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 := 0.2529822 
+0.334664 0.3577709 0.3794733 0.3898718 0.40 ;
+GEOM := GEO: :: TUBE 9
+      R+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.45 0.5748331 0.6770275
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE2>>
+          ;
+*
+TRACK := SYBILT: GEOM ::
+  TITLE 'FIRST UO2 ROWLAND S BENCHMARK'
+  MAXR 20 QUA1 5 ;
+*----
+* SHI LEVEL 1
+*----
+LIBRARY := LIB: ::
+   EDIT 1
+   NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+   CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+   ANIS 2
+   ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 293.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>> 293.0 (*GAIN1*)
+    Zr0     = Zr0    4.3241E-2 2
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 6.6988E-2
+    O16     = O16    3.3494E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 6.6988E-2
+    O16     = O16    3.3494E-2
+;
+LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 56 LEVEL 1 MXIT 3 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.387970 ;
+LIBRARY2 CP FLUX := DELETE: LIBRARY2 CP FLUX ;
+*----
+* SHI LEVEL 2
+*----
+LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 56 LEVEL 2 MXIT 3 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.388679 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* USS
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ SUBG      (*HELIOS TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 293.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>> 293.0 (*GAIN1*)
+    Zr0     = Zr0    4.3241E-2 2
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 6.6988E-2
+    O16     = O16    3.3494E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 6.6988E-2
+    O16     = O16    3.3494E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 2 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.387089 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* Ribon extended
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ PTSL      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 293.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1 IRSET PTSL 1
+    U238    = U238   2.2604E-2 1 IRSET PTSL 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>> 293.0 (*GAIN1*)
+    Zr0     = Zr0    4.3241E-2 2 IRSET PTSL 1
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 6.6988E-2
+    O16     = O16    3.3494E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 6.6988E-2
+    O16     = O16    3.3494E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 2 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.387076 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* Ribon extended with correlation
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ PTSL      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 293.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U238    = U238   2.2604E-2 CORR 1 IRSET PTSL 1
+    U235    = U235   7.0803E-4 CORR 1 IRSET PTSL 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>> 293.0 (*GAIN1*)
+    Zr0     = Zr0    4.3241E-2 2 IRSET PTSL 1
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 6.6988E-2
+    O16     = O16    3.3494E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 6.6988E-2
+    O16     = O16    3.3494E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 4 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.386651 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+ECHO "test uo2_6c_c1_sh completed" ;
+END: ;
diff --git a/Dragon/data/trowland_shem_proc/uo2_6c_c2_sh.c2m b/Dragon/data/trowland_shem_proc/uo2_6c_c2_sh.c2m
new file mode 100644
index 0000000..6b635bc
--- /dev/null
+++ b/Dragon/data/trowland_shem_proc/uo2_6c_c2_sh.c2m
@@ -0,0 +1,182 @@
+*----
+*
+*  TEST CASE uo2_6c_c2_sh
+*  SECOND UO2 ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  281-GROUP JEF2.2 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK LIBRARY LIBRARY2 CP FLUX ;
+MODULE LIB: GEO: SYBILT: SHI: USS: ASM: FLU: DELETE: END: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 MODE2 := 
+        1 2 3 4 5 6 7 8 9 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+*
+EVALUATE RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 := 0.2529822 
+0.334664 0.3577709 0.3794733 0.3898718 0.40 ;
+GEOM := GEO: :: TUBE 9
+      R+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.45 0.5748331 0.6770275
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE2>>
+          ;
+*
+TRACK := SYBILT: GEOM ::
+  TITLE 'SECOND UO2 ROWLAND S BENCHMARK'
+  MAXR 20 QUA1 5 ;
+*----
+* SHI LEVEL 1
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 6 nutot nudel NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 293.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>> 293.0 (*GAIN1*)
+    Zr0     = Zr0    4.3241E-2 2
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+;
+LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 56 LEVEL 1 MXIT 3 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.335459 ;
+LIBRARY2 CP FLUX := DELETE: LIBRARY2 CP FLUX ;
+*----
+* SHI LEVEL 2
+*----
+LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 56 LEVEL 2 MXIT 3 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.336371 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* USS
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ SUBG      (*HELIOS TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 293.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>> 293.0 (*GAIN1*)
+    Zr0     = Zr0    4.3241E-2 2
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 2 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.334628 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* Ribon extended
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ PTSL      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 293.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1 IRSET PTSL 1
+    U238    = U238   2.2604E-2 1 IRSET PTSL 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>> 293.0 (*GAIN1*)
+    Zr0     = Zr0    4.3241E-2 2 IRSET PTSL 1
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 2 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.334246 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* Ribon extended with correlation
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ PTSL      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 293.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 CORR 1 IRSET PTSL 1
+    U238    = U238   2.2604E-2 CORR 1 IRSET PTSL 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>> 293.0 (*GAIN1*)
+    Zr0     = Zr0    4.3241E-2 2 IRSET PTSL 1
+ MIX <<MODE1>> 293.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ MIX <<MODE2>> 293.0 (*MODE2*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 4 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.333725 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+ECHO "test uo2_6c_c2_sh completed" ;
+END: ;
diff --git a/Dragon/data/trowland_shem_proc/uo2_6c_c3_sh.c2m b/Dragon/data/trowland_shem_proc/uo2_6c_c3_sh.c2m
new file mode 100644
index 0000000..14684be
--- /dev/null
+++ b/Dragon/data/trowland_shem_proc/uo2_6c_c3_sh.c2m
@@ -0,0 +1,182 @@
+*----
+*
+*  TEST CASE uo2_6c_c3_sh
+*  THIRD UO2 ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  281-GROUP JEF2.2 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK LIBRARY LIBRARY2 CP FLUX ;
+MODULE LIB: GEO: SYBILT: SHI: USS: ASM: FLU: DELETE: END: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 MODE2 := 
+        1 2 3 4 5 6 7 8 9 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+*
+EVALUATE RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 := 0.2529822 
+0.334664 0.3577709 0.3794733 0.3898718 0.40 ;
+GEOM := GEO: :: TUBE 9
+      R+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.45 0.5748331 0.6770275
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE2>>
+          ;
+*
+TRACK := SYBILT: GEOM ::
+  TITLE 'THIRD UO2 ROWLAND S BENCHMARK'
+  MAXR 20 QUA1 5 ;
+*----
+* SHI LEVEL 1
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ 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*)
+    Zr0     = Zr0    4.3241E-2 2
+ MIX <<MODE1>> 574.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ MIX <<MODE2>> 574.0 (*MODE2*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+;
+LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 56 LEVEL 1 MXIT 3 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.303292 ;
+LIBRARY2 CP FLUX := DELETE: LIBRARY2 CP FLUX ;
+*----
+* SHI LEVEL 2
+*----
+LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 56 LEVEL 2 MXIT 3 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.304779 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* USS
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ SUBG      (*HELIOS TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ 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*)
+    Zr0     = Zr0    4.3241E-2 2
+ MIX <<MODE1>> 574.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ MIX <<MODE2>> 574.0 (*MODE2*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 2 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.302482 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* Ribon extended
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ PTSL      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 900.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1 IRSET PTSL 1
+    U238    = U238   2.2604E-2 1 IRSET PTSL 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*)
+    Zr0     = Zr0    4.3241E-2 2 IRSET PTSL 1
+ MIX <<MODE1>> 574.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ MIX <<MODE2>> 574.0 (*MODE2*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 2 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.303266 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* Ribon extended with correlation
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ PTSL      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 900.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 CORR 1 IRSET PTSL 1
+    U238    = U238   2.2604E-2 CORR 1 IRSET PTSL 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*)
+    Zr0     = Zr0    4.3241E-2 2 IRSET PTSL 1
+ MIX <<MODE1>> 574.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ MIX <<MODE2>> 574.0 (*MODE2*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 4 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.303032 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+ECHO "test uo2_6c_c3_sh completed" ;
+END: ;
diff --git a/Dragon/data/trowland_shem_proc/uo2_6c_c4_sh.c2m b/Dragon/data/trowland_shem_proc/uo2_6c_c4_sh.c2m
new file mode 100644
index 0000000..98a3bd0
--- /dev/null
+++ b/Dragon/data/trowland_shem_proc/uo2_6c_c4_sh.c2m
@@ -0,0 +1,182 @@
+*----
+*
+*  TEST CASE uo2_6c_c4_sh
+*  THIRD UO2 ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  281-GROUP JEF2.2 DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK LIBRARY LIBRARY2 CP FLUX ;
+MODULE LIB: GEO: SYBILT: SHI: USS: ASM: FLU: DELETE: END: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 MODE2 := 
+        1 2 3 4 5 6 7 8 9 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+*
+EVALUATE RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 := 0.2529822 
+0.334664 0.3577709 0.3794733 0.3898718 0.40 ;
+GEOM := GEO: :: TUBE 9
+      R+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.45 0.5748331 0.6770275
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE2>>
+          ;
+*
+TRACK := SYBILT: GEOM ::
+  TITLE 'FOURTH UO2 ROWLAND S BENCHMARK'
+  MAXR 20 QUA1 5 ;
+*----
+* SHI LEVEL 1
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 574.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>> 574.0 (*GAIN1*)
+    Zr0     = Zr0    4.3241E-2 2
+ MIX <<MODE1>> 574.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ MIX <<MODE2>> 574.0 (*MODE2*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+;
+LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 56 LEVEL 1 MXIT 3 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.316859 ;
+LIBRARY2 CP FLUX := DELETE: LIBRARY2 CP FLUX ;
+*----
+* SHI LEVEL 2
+*----
+LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 56 LEVEL 2 MXIT 3 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.318105 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* USS
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ SUBG      (*HELIOS TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 574.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>> 574.0 (*GAIN1*)
+    Zr0     = Zr0    4.3241E-2 2
+ MIX <<MODE1>> 574.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ MIX <<MODE2>> 574.0 (*MODE2*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 2 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.315926 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* Ribon extended
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ PTSL      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 574.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1 IRSET PTSL 1
+    U238    = U238   2.2604E-2 1 IRSET PTSL 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>> 574.0 (*GAIN1*)
+    Zr0     = Zr0    4.3241E-2 2 IRSET PTSL 1
+ MIX <<MODE1>> 574.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ MIX <<MODE2>> 574.0 (*MODE2*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 2 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.316123 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+*----
+* Ribon extended with correlation
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ PTSL      (*CALENDF TYPE PROBABILITY TABLES*)
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 574.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 CORR 1 IRSET PTSL 1
+    U238    = U238   2.2604E-2 CORR 1 IRSET PTSL 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>> 574.0 (*GAIN1*)
+    Zr0     = Zr0    4.3241E-2 2 IRSET PTSL 1
+ MIX <<MODE1>> 574.0 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ MIX <<MODE2>> 574.0 (*MODE2*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 4 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: PIJ ;
+FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ;
+assertS FLUX :: K-EFFECTIVE 1 1.315763 ;
+LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ;
+ECHO "test uo2_6c_c4_sh completed" ;
+END: ;
diff --git a/Dragon/data/trowland_shem_proc/uo2_6c_c5_sh.c2m b/Dragon/data/trowland_shem_proc/uo2_6c_c5_sh.c2m
new file mode 100644
index 0000000..50aa4b7
--- /dev/null
+++ b/Dragon/data/trowland_shem_proc/uo2_6c_c5_sh.c2m
@@ -0,0 +1,83 @@
+*----
+*
+*  TEST CASE uo2_6c_c5_sh
+*  FIFTH UO2 ROWLAND'S BENCHMARK
+*  DISTRIBUTED SELF-SHIELDING
+*  281-GROUP SHEM DRAGLIB
+*
+*  Author: A. Hebert
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM TRACK LIBRARY LIBRARY2 CP CALC ;
+MODULE LIB: GEO: SYBILT: USS: ASM: FLU: ;
+INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1
+        MODE1 MODE2 := 
+        1 2 3 4 5 6 7 8 9 ; (* DISTRIBUTED SELF-SHIELDING *)
+REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ;
+PROCEDURE assertS ;
+*
+EVALUATE RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 := 0.2529822 
+0.334664 0.3577709 0.3794733 0.3898718 0.40 ;
+GEOM := GEO: :: TUBE 9
+      R+ REFL
+      RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>>
+                 <<RCOMB5>> <<RCOMB6>> 0.45 0.5748331 0.6770275
+      MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>>
+          <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE2>>
+          ;
+TRACK := SYBILT: GEOM ::
+     TITLE 'FIFTH UO2 ROWLAND S BENCHMARK'
+     MAXR 20 QUA1 5 ;
+*----
+* USS
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 9    (*MAXIMUM OF MATERIAL MIXTURES*)
+ CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*)
+ ANIS 2
+ ADED 4 NELAS N4N N2N N3N
+ SUBG
+*
+ MIXS LIB: DRAGON FIL: DLIB_SH
+ MIX <<COMB0101>> 1200.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1 IRSET 0.0 NONE
+    U238    = U238   2.2604E-2 1 IRSET 0.0 NONE
+ MIX <<COMB0102>> 1000.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1 IRSET 0.0 NONE
+    U238    = U238   2.2604E-2 1 IRSET 0.0 NONE
+ MIX <<COMB0103>> 800.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1 IRSET 0.0 NONE
+    U238    = U238   2.2604E-2 1 IRSET 0.0 NONE
+ MIX <<COMB0104>> 700.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1 IRSET 0.0 NONE
+    U238    = U238   2.2604E-2 1 IRSET 0.0 NONE
+ MIX <<COMB0105>> 650.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1 IRSET 0.0 NONE
+    U238    = U238   2.2604E-2 1 IRSET 0.0 NONE
+ MIX <<COMB0106>> 600.0 (*COMB0101*)
+    O16     = O16    4.6624E-2
+    U235    = U235   7.0803E-4 1 IRSET 0.0 NONE
+    U238    = U238   2.2604E-2 1 IRSET 0.0 NONE
+ MIX <<GAIN1>> 600.0 (*GAIN1*)
+    Zr0     = Zr0    4.3241E-2 2 IRSET 0.0 NONE
+ MIX <<MODE1>> 573.6 (*MODE1*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ MIX <<MODE2>> 573.6 (*MODE2*)
+    H1      = H1_H2O 4.6892E-2
+    O16     = O16    2.3446E-2
+ ;
+LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 5 GRMIN 56 ;
+CP := ASM: LIBRARY2 TRACK :: EDIT 1 PIJ ;
+CALC := FLU: CP LIBRARY2 TRACK :: EDIT 1 TYPE K ;
+assertS CALC :: K-EFFECTIVE 1 1.301500 ;
+ECHO "test uo2_6c_c5_sh completed" ;
+QUIT "LIST" .