summaryrefslogtreecommitdiff
path: root/Dragon/data/twimsE_proc
diff options
context:
space:
mode:
Diffstat (limited to 'Dragon/data/twimsE_proc')
-rw-r--r--Dragon/data/twimsE_proc/TCWE01.c2m129
-rw-r--r--Dragon/data/twimsE_proc/TCWE02.c2m179
-rw-r--r--Dragon/data/twimsE_proc/TCWE03.c2m172
-rw-r--r--Dragon/data/twimsE_proc/TCWE04.c2m151
-rw-r--r--Dragon/data/twimsE_proc/TCWE05.c2m192
-rw-r--r--Dragon/data/twimsE_proc/TCWE05Lib.c2m68
-rw-r--r--Dragon/data/twimsE_proc/TCWE06.c2m155
-rw-r--r--Dragon/data/twimsE_proc/TCWE07.c2m142
-rw-r--r--Dragon/data/twimsE_proc/TCWE08.c2m121
-rw-r--r--Dragon/data/twimsE_proc/TCWE09.c2m74
-rw-r--r--Dragon/data/twimsE_proc/TCWE10.c2m104
-rw-r--r--Dragon/data/twimsE_proc/TCWE11.c2m171
-rw-r--r--Dragon/data/twimsE_proc/TCWE12.c2m107
-rw-r--r--Dragon/data/twimsE_proc/TCWE13.c2m146
-rw-r--r--Dragon/data/twimsE_proc/TCWE14.c2m121
-rw-r--r--Dragon/data/twimsE_proc/TCWE15.c2m159
-rw-r--r--Dragon/data/twimsE_proc/TCWE17.c2m337
-rw-r--r--Dragon/data/twimsE_proc/TCWE17Lib.c2m97
-rw-r--r--Dragon/data/twimsE_proc/TCWE31.c2m172
-rw-r--r--Dragon/data/twimsE_proc/TCWE32.c2m97
20 files changed, 2894 insertions, 0 deletions
diff --git a/Dragon/data/twimsE_proc/TCWE01.c2m b/Dragon/data/twimsE_proc/TCWE01.c2m
new file mode 100644
index 0000000..a18e4da
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE01.c2m
@@ -0,0 +1,129 @@
+*----
+* TEST CASE TCWE01
+* MOSTELLER BENCHMARK: 1-D ANNULAR CELL AND 2-D CARTESIAN CELL
+* WIMSE 69 GROUPS LIBRARY FILE WIMS-E
+*
+* REF: R. Mosteller et al. Nucl. Sci. Eng. 107, 265 (1991)
+*
+*----
+* Define STRUCTURES and MODULES used
+*----
+LINKED_LIST
+ MOSTELA MOSTELC DISCR LIBRARY CP CALC OUT ;
+SEQ_BINARY
+ TRKSPC ;
+MODULE
+ LIB: GEO: SYBILT: EXCELT: SHI: ASM: FLU: EDI: DELETE: END: ;
+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.83243E-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 * ;
+
+REAL P_B10 P_B11 := 0.199 0.801 ;
+REAL N_B0 := 1.02133E-5 ;
+REAL N_B10 N_B11 := N_B0 P_B10 * N_B0 P_B11 * ;
+*----
+* Microscopic cross sections from file WIMS-E format WIMSE
+*----
+LIBRARY := LIB: ::
+ NMIX 3 CTRA WIMS
+ MIXS LIB: WIMSE FIL: WIMS-E
+ MIX 1 600.0 O16 = '8016' 4.61309E-2
+ U235 = '92235' 1.66078E-4 1
+ U238 = '92238' 2.28994E-2 1
+ MIX 2 600.0
+ Zr90 = '40090' <<N_Zr90>>
+ Zr91 = '40091' <<N_Zr91>> 2
+ Zr92 = '40092' <<N_Zr92>>
+ Zr94 = '40094' <<N_Zr94>>
+ Zr96 = '40096' <<N_Zr96>>
+ MIX 3 600.0
+ H1H2O = '1001' 4.42326E-2 O16H2O = '8016' 2.21163E-2
+ B10 = '5010' <<N_B10>> B11 = '5011' <<N_B11>>
+ ;
+*----
+* Geometry MOSTELA : annular 3 region geometry
+* MOSTELC : Cartesian 3 region geometry
+*----
+MOSTELA := GEO: :: TUBE 3
+ R+ REFL RADIUS 0.0 0.39306 0.45802 0.71206 SPLITR 2 1 1
+ MIX 1 2 3 ;
+MOSTELC := GEO: :: CARCEL 2
+ X- REFL X+ REFL MESHX 0.0 1.26209
+ Y- REFL Y+ REFL MESHY 0.0 1.26209
+ RADIUS 0.0 0.39306 0.45802 SPLITR 2 1
+ MIX 1 2 3 ;
+*----
+* Case 1 -- annular
+* Self-Shielding calculation SYBIL
+* Transport calculation SYBIL
+* Flux calculation for K no leakage
+*----
+DISCR := SYBILT: MOSTELA ::
+ TITLE 'TCWE01: MOSTELLER BENCHMARK (SYBIL / SYBIL)'
+ MAXR 4 QUA1 5 ;
+LIBRARY := SHI: LIBRARY DISCR :: EDIT 0 NOLJ ;
+CP := ASM: LIBRARY DISCR ;
+CALC := FLU: CP LIBRARY DISCR ::
+ TYPE K ;
+assertS CALC :: 'K-EFFECTIVE' 1 0.8289803 ;
+OUT := EDI: LIBRARY DISCR CALC ::
+ EDIT 4 MERG MIX 1 2 3 COND 4.0 SAVE ;
+DISCR CP := DELETE: DISCR CP ;
+*----
+* Case 2 -- Cartesian
+* Self-Shielding calculation SYBIL
+* Transport calculation SYBIL
+* Flux calculation for K no leakage
+*----
+DISCR := SYBILT: MOSTELC ::
+ TITLE 'TCWE01: MOSTELLER BENCHMARK (SYBIL / SYBIL)'
+ MAXR 4 QUA1 5 QUA2 6 5 ;
+LIBRARY := SHI: LIBRARY DISCR :: EDIT 0 NOLJ ;
+CP := ASM: LIBRARY DISCR ;
+CALC := FLU: CALC CP LIBRARY DISCR ::
+ TYPE K ;
+assertS CALC :: 'K-EFFECTIVE' 1 0.8290998 ;
+OUT := EDI: OUT LIBRARY DISCR CALC ::
+ EDIT 1 MERG MIX 1 2 3 COND 4.0 STAT ALL REFE 1 ;
+CALC DISCR CP := DELETE: CALC DISCR CP ;
+*----
+* Case 3 -- annular
+* Self-Shielding calculation EXCEL-ISO
+* Transport calculation EXCEL-ISO
+* Flux calculation for K no leakage
+*----
+DISCR TRKSPC := EXCELT: MOSTELC ::
+ TITLE 'TCWE01: MOSTELLER BENCHMARK (EXCELL)'
+ MAXR 4 TRAK TISO 12 20.0 ;
+LIBRARY := SHI: LIBRARY DISCR TRKSPC :: EDIT 0 NOLJ ;
+CP := ASM: LIBRARY DISCR TRKSPC ;
+CALC := FLU: CP LIBRARY DISCR ::
+ TYPE K ;
+assertS CALC :: 'K-EFFECTIVE' 1 0.8302385 ;
+OUT := EDI: OUT LIBRARY DISCR CALC ::
+ EDIT 1 MERG MIX 1 2 3 COND 4.0 STAT ALL REFE 1 ;
+DISCR TRKSPC CP := DELETE: DISCR TRKSPC CP ;
+*----
+* Case 4 -- Cartesian
+* Self-Shielding calculation EXCEL-SPC
+* Transport calculation EXCEL-SPC
+* Flux calculation for K no leakage
+*----
+DISCR TRKSPC := EXCELT: MOSTELC ::
+ TITLE 'TCWE01: MOSTELLER BENCHMARK (EXCELL)'
+ MAXR 4 TRAK TSPC 12 20.0 ;
+LIBRARY := SHI: LIBRARY DISCR TRKSPC :: EDIT 0 NOLJ ;
+CP := ASM: LIBRARY DISCR TRKSPC ;
+CALC := FLU: CALC CP LIBRARY DISCR ::
+ TYPE K ;
+assertS CALC :: 'K-EFFECTIVE' 1 0.8325990 ;
+OUT := EDI: OUT LIBRARY DISCR CALC ::
+ EDIT 1 MERG MIX 1 2 3 COND 4.0 STAT ALL REFE 1 ;
+DISCR TRKSPC CP := DELETE: DISCR TRKSPC CP ;
+ECHO "test TCWE01 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE02.c2m b/Dragon/data/twimsE_proc/TCWE02.c2m
new file mode 100644
index 0000000..afb01e4
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE02.c2m
@@ -0,0 +1,179 @@
+*----
+* TEST CASE TCWE02
+* 17 X 17 PWR ASSEMBLY WITHOUT POISON
+* WIMSE 69 GROUPS LIBRARY FILE WIMS-E
+*
+* REF: none
+*
+*----
+* Define STRUCTURES and MODULES used
+*----
+LINKED_LIST
+ ASSMB DISCR LIBRARY CP CALC OUT DATABASE ISOT SPHGEOM MTRACK ;
+SEQ_ASCII
+ res ;
+MODULE
+ LIB: GEO: SYBILT: BIVACT: SHI: ASM: FLU: EDI: COMPO: SPH: DELETE:
+ END: ;
+PROCEDURE assertS ;
+*----
+* Microscopic cross sections from file WIMS-E format WIMSE
+*----
+LIBRARY := LIB: ::
+ NMIX 8 CTRA WIMS
+ MIXS LIB: WIMSE FIL: WIMS-E
+ MIX 1 579.9
+ H1H2O = '1001' 4.76690E-2 O16H2O = '8016' 2.38345E-2
+ B10 = '5010' 4.73825E-06 B11 = '5011' 1.90721E-05
+ MIX 2 579.9 O16 = '8016' 3.06711E-4
+ Cr52 = '24052' 7.54987E-5 Fe56 = '26056' 1.47624E-4
+ Zr90 = '40090' 2.15381E-02 Zr91 = '40091' 4.69693E-03 2
+ Zr92 = '40092' 7.17935E-03 Zr94 = '40094' 7.27563E-03
+ Zr96 = '40096' 1.17214E-03
+ MIX 3 579.9
+ H1H2O = '1001' 4.65292E-2 O16H2O = '8016' 2.32646E-2
+ Cr52 = '24052' 4.79927E-5
+ Fe56 = '26056' 4.45845E-5 Ni58 = '28058' 1.13521E-4
+ Mo92 = '42092' 5.91461E-07 Mo95 = '42095' 3.70930E-07
+ Mo96 = '42096' 6.40880E-07 Mo97 = '42097' 3.86878E-07
+ Mo98 = '42098' 9.80802E-07 Mo100 = '42100' 3.93419E-07
+ Al27 = '13027' 2.35231E-6
+ Mn55 = '25055' 4.15901E-7
+ B10 = '5010' 4.63194E-06 B11 = '5011' 1.86442E-05
+ Zr90 = '40090' 4.59154E-04 Zr91 = '40091' 1.00130E-04
+ Zr92 = '40092' 1.53051E-04 Zr94 = '40094' 1.55104E-04
+ Zr96 = '40096' 2.49880E-05
+ MIX 4 579.9 O16 = '8016' 2.87335E-4
+ Cr52 = '24052' 7.07291E-5 Fe56 = '26056' 1.38298E-4
+ Zr90 = '40090' 2.01774E-02 Zr91 = '40091' 4.40020E-03 2
+ Zr92 = '40092' 6.72580E-03 Zr94 = '40094' 6.81600E-03
+ Zr96 = '40096' 1.09809E-03
+ MIX 5 579.9
+ H1H2O = '1001' 4.71346E-2 O16H2O = '8016' 2.35673E-2
+ Cr52 = '24052' 2.24991E-5
+ Fe56 = '26056' 2.09013E-5 Ni58 = '28058' 5.32188E-5
+ Mo92 = '42092' 2.77278E-07 Mo95 = '42095' 1.73892E-07
+ Mo96 = '42096' 3.00446E-07 Mo97 = '42097' 1.81369E-07
+ Mo98 = '42098' 4.59801E-07 Mo100 = '42100' 1.84435E-07
+ Al27 = '13027' 1.10277E-6
+ Mn55 = '25055' 1.94976E-7
+ B10 = '5010' 4.68840E-06 B11 = '5011' 1.88714E-05
+ Zr90 = '40090' 2.15252E-04 Zr91 = '40091' 4.69413E-05
+ Zr92 = '40092' 7.17508E-05 Zr94 = '40094' 7.27131E-05
+ Zr96 = '40096' 1.17144E-05
+ MIX 6 579.9
+ H1H2O = '1001' 4.71676E-2 O16H2O = '8016' 2.35838E-2
+ Cr52 = '24052' 2.11122E-5
+ Fe56 = '26056' 1.96130E-5 Ni58 = '28058' 4.99383E-5
+ Mo92 = '42092' 2.60187E-07 Mo95 = '42095' 1.63174E-07
+ Mo96 = '42096' 2.81927E-07 Mo97 = '42097' 1.70190E-07
+ Mo98 = '42098' 4.31460E-07 Mo100 = '42100' 1.73067E-07
+ Al27 = '13027' 1.03479E-6
+ Mn55 = '25055' 1.82957E-7
+ B10 = '5010' 4.69148E-06 B11 = '5011' 1.88838E-05
+ Zr90 = '40090' 2.01984E-04 Zr91 = '40091' 4.40478E-05
+ Zr92 = '40092' 6.73280E-05 Zr94 = '40094' 6.82309E-05
+ Zr96 = '40096' 1.09923E-05
+ MIX 7 579.9
+ H1H2O = '1001' 4.72020E-2 O16H2O = '8016' 2.36010E-2
+ Cr52 = '24052' 1.96591E-5
+ Fe56 = '26056' 1.82630E-5 Ni58 = '28058' 4.65011E-5
+ Mo92 = '42092' 2.42278E-07 Mo95 = '42095' 1.51943E-07
+ Mo96 = '42096' 2.62522E-07 Mo97 = '42097' 1.58476E-07
+ Mo98 = '42098' 4.01763E-07 Mo100 = '42100' 1.61155E-07
+ Al27 = '13027' 9.63569E-7
+ Mn55 = '25055' 1.70365E-7
+ B10 = '5010' 4.69469E-06 B11 = '5011' 1.88967E-05
+ Zr90 = '40090' 1.88082E-04 Zr91 = '40091' 4.10161E-05
+ Zr92 = '40092' 6.26939E-05 Zr94 = '40094' 6.35347E-05
+ Zr96 = '40096' 1.02357E-05
+ MIX 8 933.6 O16 = '8016' 4.49355E-2
+ U235 = '92235' 7.39237E-4 1
+ U238 = '92238' 2.17285E-2 1
+ ;
+*----
+* Geometry ASSMB : a 17 X 17 normal PWR assembly
+* contains C1 : cell without fuel
+* C2 : normal fuel cell
+* C3 : peripheral cell
+* C4 : corner cell
+*----
+ASSMB := GEO: :: CAR2D 9 9
+ X- DIAG X+ REFL Y- SYME Y+ DIAG
+ CELL C1 C2 C2 C1 C2 C2 C1 C2 C3
+ C2 C2 C2 C2 C2 C2 C2 C3
+ C2 C2 C2 C2 C2 C2 C3
+ C1 C2 C2 C1 C2 C3
+ C2 C2 C2 C2 C3
+ C1 C2 C2 C3
+ C2 C2 C3
+ C2 C3
+ C4
+ MERGE 1 3 12 11 12 12 11 12 15
+ 4 6 5 6 6 5 6 8
+ 13 5 6 6 5 6 8
+ 2 5 5 10 5 8
+ 13 5 5 6 8
+ 2 5 7 8
+ 13 7 8
+ 14 8
+ 9
+ TURN H H B H H B H H A
+ H G G H G G H A
+ A E E F E E A
+ H H F H H A
+ H E G H A
+ H H A A
+ H A A
+ A A
+ A
+ ::: C1 := GEO: CARCEL 2
+ MESHX 0.0 1.26472 MESHY 0.0 1.26472
+ RADIUS 0.0 0.572435 0.613142 MIX 1 2 3 ;
+ ::: C2 := GEO: C1 RADIUS 0.0 0.412660 0.474364 MIX 8 4 5 ;
+ ::: C3 := GEO: C2 MESHX 0.0 1.31472 MIX 8 4 6 ;
+ ::: C4 := GEO: C3 MESHY 0.0 1.31472 MIX 8 4 7 ;
+ ;
+*----
+* Self-Shielding calculation SYBIL
+* Transport calculation SYBIL
+* Flux calculation for B1 homogeneous leakage
+* Editing using SPH model for transport-diffusion
+*----
+DISCR := SYBILT: ASSMB ::
+ TITLE 'TCWE02: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON'
+ MAXR 400 QUA2 6 3 ;
+LIBRARY := SHI: LIBRARY DISCR :: EDIT 0 NOLJ ;
+CP := ASM: LIBRARY DISCR ;
+CALC := FLU: CP LIBRARY DISCR ::
+ TYPE B B1 PNL ;
+assertS CALC :: 'K-INFINITY' 1 1.258811 ;
+
+OUT := EDI: LIBRARY DISCR CALC ASSMB ::
+ EDIT 3 UPS SAVE MICR RES MERGE CELL COND 4.0
+ ;
+SPHGEOM := OUT :: STEP UP 'MACRO-GEOM' ;
+MTRACK := BIVACT: SPHGEOM :: PRIM 1 2 EDIT 0 ;
+OUT := SPH: OUT MTRACK ;
+SPHGEOM MTRACK := DELETE: SPHGEOM MTRACK ;
+
+DATABASE := COMPO: ::
+ EDIT 5
+ COMM 'Multi-parameter reactor database' ENDC
+ INIT
+ ;
+DATABASE := COMPO: DATABASE OUT ::
+ EDIT 3
+ ;
+res := DATABASE ;
+
+ISOT := DATABASE :: STEP UP default
+ STEP UP MIXTURES STEP AT 5
+ STEP UP CALCULATIONS STEP AT 1
+ STEP UP ISOTOPESLIST STEP AT 1 ;
+assertS ISOT :: 'NWT0' 1 5.182448E-01 ;
+assertS ISOT :: 'NWT0' 2 1.297516E-01 ;
+
+ECHO "test TCWE02 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE03.c2m b/Dragon/data/twimsE_proc/TCWE03.c2m
new file mode 100644
index 0000000..53d1dae
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE03.c2m
@@ -0,0 +1,172 @@
+*----
+* TEST CASE TCWE03
+* MULTICELL HEXAGONAL ASSEMBLY WITH POISON
+* WIMS-E Library
+*
+* REF: none
+*
+*----
+* Define STRUCTURES and MODULES used
+*----
+LINKED_LIST
+ ASSMBH DISCR LIBRARY CP CALC OUT DATABASE ISOT SPHGEOM MTRACK ;
+SEQ_ASCII
+ res ;
+MODULE
+ GEO: SYBILT: BIVACT: LIB: SHI: ASM: FLU: EDI: COMPO: SPH: DELETE:
+ END: ;
+PROCEDURE assertS ;
+*----
+* Microscopic cross sections from file WIMS-E format WIMSE
+*----
+LIBRARY := LIB: ::
+ NMIX 11 CTRA WIMS
+ MIXS LIB: WIMSE FIL: WIMS-E
+ MIX 1 579.9
+ H1H2O = '1001' 4.76690E-2 O16H2O = '8016' 2.38345E-2
+ B10 = '5010' 4.73825E-06 B11 = '5011' 1.90721E-05
+ MIX 2 579.9 O16 = '8016' 3.06711E-4
+ Cr52 = '24052' 7.54987E-5 Fe56 = '26056' 1.47624E-4
+ Zr90 = '40090' 2.15381E-02 Zr91 = '40091' 4.69693E-03 2
+ Zr92 = '40092' 7.17935E-03 Zr94 = '40094' 7.27563E-03
+ Zr96 = '40096' 1.17214E-03
+ MIX 3 579.9
+ H1H2O = '1001' 4.65292E-2 O16H2O = '8016' 2.32646E-2
+ Cr52 = '24052' 4.79927E-5
+ Fe56 = '26056' 4.45845E-5 Ni58 = '28058' 1.13521E-4
+ Mo92 = '42092' 5.91461E-07 Mo95 = '42095' 3.70930E-07
+ Mo96 = '42096' 6.40880E-07 Mo97 = '42097' 3.86878E-07
+ Mo98 = '42098' 9.80802E-07 Mo100 = '42100' 3.93419E-07
+ Mn55 = '25055' 4.15901E-7
+ B10 = '5010' 4.63194E-06 B11 = '5011' 1.86442E-05
+ Zr90 = '40090' 4.59154E-04 Zr91 = '40091' 1.00130E-04
+ Zr92 = '40092' 1.53051E-04 Zr94 = '40094' 1.55104E-04
+ Zr96 = '40096' 2.49880E-05
+ MIX 4 933.6 O16 = '8016' 4.49355E-2
+ U235 = '92235' 7.39237E-4 1
+ U238 = '92238' 2.17285E-2 1
+ MIX 5 579.9 In115 = '49115' 7.57464E-3
+ Cd113 = '48113' 2.62493E-3
+ Ag109 = '47109' 4.49188E-2
+ MIX 6 579.9 Cr52 = '24052' 1.52702E-2
+ Fe56 = '26056' 5.57670E-2 Ni58 = '28058' 7.51418E-3
+ Mn55 = '25055' 8.02943E-4
+ MIX 7 579.9
+ H1H2O = '1001' 3.06466E-2 O16H2O = '8016' 1.53233E-2
+ Fe56 = '26056' 5.27485E-5 Cr52 = '24052' 2.69769E-5
+ B10 = '5010' 3.04623E-06 B11 = '5011' 1.22615E-05
+ Zr90 = '40090' 7.69589E-03 Zr91 = '40091' 1.67829E-03 2
+ Zr92 = '40092' 2.56530E-03 Zr94 = '40094' 2.59970E-03
+ Zr96 = '40096' 4.18824E-04
+ MIX 8 579.9
+ H1H2O = '1001' 4.65292E-2 O16H2O = '8016' 2.32646E-2
+ Cr52 = '24052' 4.79927E-5
+ Fe56 = '26056' 4.45845E-5 Ni58 = '28058' 1.13521E-4
+ Mo92 = '42092' 5.91461E-07 Mo95 = '42095' 3.70930E-07
+ Mo96 = '42096' 6.40880E-07 Mo97 = '42097' 3.86878E-07
+ Mo98 = '42098' 9.80802E-07 Mo100 = '42100' 3.93419E-07
+ Al27 = '13027' 2.35231E-6
+ Mn55 = '25055' 4.15901E-7
+ B10 = '5010' 4.63194E-06 B11 = '5011' 1.86442E-05
+ Zr90 = '40090' 4.59154E-04 Zr91 = '40091' 1.00130E-04
+ Zr92 = '40092' 1.53051E-04 Zr94 = '40094' 1.55104E-04
+ Zr96 = '40096' 2.49880E-05
+ MIX 9 579.9 O16 = '8016' 2.87335E-4
+ Cr52 = '24052' 7.07291E-5 Fe56 = '26056' 1.38298E-4
+ Zr90 = '40090' 2.01774E-02 Zr91 = '40091' 4.40020E-03 2
+ Zr92 = '40092' 6.72580E-03 Zr94 = '40094' 6.81600E-03
+ Zr96 = '40096' 1.09809E-03
+ MIX 10 579.9
+ H1H2O = '1001' 4.71346E-2 O16H2O = '8016' 2.35673E-2
+ Cr52 = '24052' 2.24991E-5
+ Fe56 = '26056' 2.09013E-5 Ni58 = '28058' 5.32188E-5
+ Mo92 = '42092' 2.77278E-07 Mo95 = '42095' 1.73892E-07
+ Mo96 = '42096' 3.00446E-07 Mo97 = '42097' 1.81369E-07
+ Mo98 = '42098' 4.59801E-07 Mo100 = '42100' 1.84435E-07
+ Al27 = '13027' 1.10277E-6
+ Mn55 = '25055' 1.94976E-7
+ B10 = '5010' 4.68840E-06 B11 = '5011' 1.88714E-05
+ Zr90 = '40090' 2.15252E-04 Zr91 = '40091' 4.69413E-05
+ Zr92 = '40092' 7.17508E-05 Zr94 = '40094' 7.27131E-05
+ Zr96 = '40096' 1.17144E-05
+ MIX 11 579.9
+ H1H2O = '1001' 4.71676E-2 O16H2O = '8016' 2.35838E-2
+ Cr52 = '24052' 2.11122E-5
+ Fe56 = '26056' 1.96130E-5 Ni58 = '28058' 4.99383E-5
+ Mo92 = '42092' 2.60187E-07 Mo95 = '42095' 1.63174E-07
+ Mo96 = '42096' 2.81927E-07 Mo97 = '42097' 1.70190E-07
+ Mo98 = '42098' 4.31460E-07 Mo100 = '42100' 1.73067E-07
+ Al27 = '13027' 1.03479E-6
+ Mn55 = '25055' 1.82957E-7
+ B10 = '5010' 4.69148E-06 B11 = '5011' 1.88838E-05
+ Zr90 = '40090' 2.01984E-04 Zr91 = '40091' 4.40478E-05
+ Zr92 = '40092' 6.73280E-05 Zr94 = '40094' 6.82309E-05
+ Zr96 = '40096' 1.09923E-05
+ ;
+*----
+* Geometry ASSMBH : hexagonal assembly with poison
+* contains C1 : cell without fuel
+* C2 : poison cell
+* C3 : normal fuel cell
+* C4 : peripheral cell
+*----
+ASSMBH := GEO: :: HEX 36
+ HBC S30 REFL
+ CELL C1 C3 C3 C3 C3 C3 C2 C3 C3 C3 C2 C3 C3 C3 C3 C3 C3 C2
+ C3 C3 C2 C3 C3 C3 C3 C3 C3 C3 C3 C3 C4 C4 C4 C4 C4 C4
+ TURN A A A A A A A A B D A I A C F J B A
+ F A A E E A A E A A A A A A A A A A
+ MERGE 1 2 3 4 5 4 6 7 8 7 9 8 10 7 7 4 7 11
+ 12 13 14 15 12 16 17 12 16 18 18 19 20 21 21 22 22 23
+ ::: C1 := GEO: HEXCEL 2
+ SIDE 0.707297 RADIUS 0.0 0.412282 0.475917
+ MIX 1 2 3 ;
+ ::: C2 := GEO: HEXCEL 5
+ SIDE 0.707297 RADIUS 0.0 0.25057 0.354359 0.436 0.486 0.6125
+ MIX 5 5 5 6 7 8 ;
+ ::: C3 := GEO: C1 MIX 4 9 10 ;
+ ::: C4 := GEO: C3 MIX 4 9 11 ;
+ ;
+*----
+* Self-Shielding calculation JPM
+* Transport calculation SYBIL
+* Flux calculation for B1 homogeneous leakage
+* Editing using SPH model for transport-diffusion
+*----
+DISCR := SYBILT: ASSMBH ::
+ TITLE 'TCWE03: MULTICELL HEXAGONAL ASSEMBLY WITH POISON'
+ MAXR 400 MAXZ 15000 QUA2 6 3 ;
+LIBRARY := SHI: LIBRARY DISCR :: EDIT 0 NOLJ ;
+CP := ASM: LIBRARY DISCR ;
+CALC := FLU: CP LIBRARY DISCR ::
+ TYPE B B1 PNL ;
+assertS CALC :: 'K-INFINITY' 1 0.6724899 ;
+OUT := EDI: LIBRARY DISCR CALC ASSMBH ::
+ EDIT 3 UPS SAVE MICR RES MERGE CELL COND 4.0
+ ;
+SPHGEOM := OUT :: STEP UP 'MACRO-GEOM' ;
+SPHGEOM := GEO: SPHGEOM :: SPLITH 0 ;
+MTRACK := BIVACT: SPHGEOM :: PRIM 1 2 ;
+OUT := SPH: OUT MTRACK ;
+SPHGEOM MTRACK := DELETE: SPHGEOM MTRACK ;
+
+DATABASE := COMPO: ::
+ EDIT 5
+ COMM 'Multi-parameter reactor database' ENDC
+ INIT
+ ;
+DATABASE := COMPO: DATABASE OUT ::
+ EDIT 3
+ ;
+res := DATABASE ;
+
+ISOT := DATABASE :: STEP UP default
+ STEP UP MIXTURES STEP AT 5
+ STEP UP CALCULATIONS STEP AT 1
+ STEP UP ISOTOPESLIST STEP AT 1 ;
+assertS ISOT :: 'NWT0' 1 1.483070E+00 ;
+assertS ISOT :: 'NWT0' 2 1.400084E-01 ;
+
+ECHO "test TCWE03 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE04.c2m b/Dragon/data/twimsE_proc/TCWE04.c2m
new file mode 100644
index 0000000..e9ab59f
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE04.c2m
@@ -0,0 +1,151 @@
+*----
+* TEST CASE TCWE04
+* WIMS-E Library
+* ANNULAR MOSTELLER BENCHMARK WITH BURNUP
+*
+* REF: R. Mosteller et al. Nucl. Sci. Eng. 107, 265 (1991)
+*
+*----
+* Define variables
+*----
+INTEGER
+ istep := 1 ;
+REAL
+ evobeg evoend ;
+REAL
+ step2 step3 step4 step5 :=
+ 1.0 27.1739 67.9348 135.8696 ;
+*----
+* Define STRUCTURES and MODULES used
+*----
+PARAMETER res :: ::: SEQ_ASCII res ; ;
+LINKED_LIST
+ LIBRARY MOSTELAS MOSTELA TRACKS TRACK SYS FLUX BURNUP
+ EDITION DATABASE ISOT ;
+MODULE
+ GEO: SYBILT: LIB: SHI: ASM: FLU: EVO: EDI: COMPO: DELETE: END: ;
+PROCEDURE assertS ;
+*----
+* Depletion data from file WIMS-E format WIMSE
+* Microscopic cross sections from file WIMS-E format WIMSE
+*----
+LIBRARY := LIB: ::
+ NMIX 3 CTRA WIMS
+ DEPL LIB: WIMSE FIL: WIMS-E
+ MIXS LIB: WIMSE FIL: WIMS-E
+ MIX 1 600.0 O16 = '8016' 4.61309E-2
+ U235 = '92235' 1.66078E-4 1
+ U238 = '92238' 2.28994E-2 1
+ U236 = '92236' 0.0 1
+ Pu239 = '94239' 0.0 1
+ MIX 2 600.0
+ Zr90 = '40090' 1.97179E-02 Zr91 = '40091' 4.29999E-03 2
+ Zr92 = '40092' 6.57262E-03 Zr94 = '40094' 6.66076E-03
+ Zr96 = '40096' 1.07308E-03
+ MIX 3 600.0
+ H1H2O = '1001' 4.42326E-2 O16H2O = '8016' 2.21163E-2
+ B10 = '5010' 2.03245E-06 B11 = '5011' 8.18085E-06
+ ;
+*----
+* Geometry MOSTELAS : 3 regions annular cell for self-shielding
+* MOSTELA : 4 regions annular cell for transport
+*----
+MOSTELAS := GEO: :: TUBE 3
+ R+ REFL RADIUS 0.0 0.39306 0.45802 0.71206
+ MIX 1 2 3 ;
+MOSTELA := GEO: MOSTELAS ::
+ SPLITR 2 1 1 ;
+*----
+* Create the reactor database
+*----
+DATABASE := COMPO: ::
+ EDIT 5
+ COMM 'Multi-parameter reactor database' ENDC
+ PARA 'BURN' IRRA
+ PARA 'FLUB' FLUB
+ INIT
+ ;
+*----
+* Self-Shielding calculation SYBIL
+* Transport calculation SYBIL
+* Flux calculation for keff with imposed buckling
+* using B1 homogeneous leakage model
+*----
+TRACKS := SYBILT: MOSTELAS ::
+ TITLE 'TCWE04: MOSTELLER BENCHMARK WITH BURNUP'
+ EDIT 1 MAXR 3 ;
+LIBRARY := SHI: LIBRARY TRACKS :: EDIT 0 NOLJ ;
+TRACK := SYBILT: MOSTELA ::
+ TITLE 'TCWE04: MOSTELLER BENCHMARK WITH BURNUP'
+ EDIT 1 MAXR 4 ;
+SYS := ASM: LIBRARY TRACK ;
+FLUX := FLU: SYS LIBRARY TRACK ::
+ TYPE K B1 PNL BUCK 0.2948E-2 ;
+EDITION := EDI: LIBRARY TRACK FLUX ::
+ EDIT 3 MICR RES MERG COMP COND 4.0 SAVE ;
+*----
+* Burnup loop: for first step BURNUP is created
+* while for other steps it is modified
+* two burnup per step:
+* 1) get a first approximation of final composition followed
+* by a transport calculation
+* 2) use approximation for final flux distribution to get a
+* better approximation for final composition
+*----
+EVALUATE evoend := 0.0 ;
+WHILE evoend step2 < DO
+ EVALUATE evobeg := evoend ;
+ EVALUATE evoend := step2 ;
+ IF istep 1 = THEN
+ BURNUP LIBRARY := EVO: LIBRARY FLUX TRACK ::
+ SAVE <<evobeg>> DAY POWR 36.8
+ DEPL <<evobeg>> <<evoend>> DAY POWR 36.8
+ SET <<evoend>> DAY ;
+ ELSE
+ BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX TRACK ::
+ SAVE <<evobeg>> DAY POWR 36.8
+ NOEX DEPL <<evobeg>> <<evoend>> DAY POWR 36.8
+ SET <<evoend>> DAY ;
+ ENDIF ;
+ LIBRARY := SHI: LIBRARY TRACKS :: EDIT 0 NOLJ ;
+ SYS := DELETE: SYS ;
+ SYS := ASM: LIBRARY TRACK ;
+ FLUX := FLU: FLUX SYS LIBRARY TRACK ::
+ TYPE K B1 PNL BUCK 0.2948E-2 ;
+
+ BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX TRACK ::
+ SAVE <<evoend>> DAY POWR 36.8
+ DEPL <<evobeg>> <<evoend>> DAY POWR 36.8
+ SET <<evoend>> DAY ;
+ LIBRARY := SHI: LIBRARY TRACKS :: EDIT 0 NOLJ ;
+ SYS := DELETE: SYS ;
+ SYS := ASM: LIBRARY TRACK ;
+ FLUX := FLU: FLUX SYS LIBRARY TRACK ::
+ TYPE K B1 PNL BUCK 0.2948E-2 ;
+
+ EDITION := EDI: EDITION LIBRARY TRACK FLUX ::
+ EDIT 3 SAVE ;
+ BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX TRACK ::
+ SAVE <<evoend>> DAY POWR 36.8 ;
+ DATABASE := COMPO: DATABASE EDITION BURNUP ::
+ EDIT 3
+ SET <<evoend>> DAY
+ ;
+ EVALUATE step2 step3 step4 step5 :=
+ step3 step4 step5 step2 ;
+ EVALUATE istep := istep 1 + ;
+ENDWHILE ;
+assertS FLUX :: 'K-EFFECTIVE' 1 0.7397701 ;
+
+res := DATABASE ;
+
+ISOT := DATABASE :: STEP UP default
+ STEP UP MIXTURES STEP AT 1
+ STEP UP CALCULATIONS STEP AT 2
+ STEP UP ISOTOPESLIST STEP AT 1 ;
+assertS ISOT :: 'NWT0' 1 3.840256E+01 ;
+assertS ISOT :: 'NWT0' 2 1.686708E+01 ;
+
+ECHO "test TCWE04 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE05.c2m b/Dragon/data/twimsE_proc/TCWE05.c2m
new file mode 100644
index 0000000..93d8f95
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE05.c2m
@@ -0,0 +1,192 @@
+*----
+* TEST CASE TCWE05
+* CANDU-6 ANNULAR CELL
+* WIMS-E Library
+* POWER (KW) = 615.00000
+* BURN POWER (KW/KG) = 31.97130
+* URANIUM MASS = 19.23600
+* UO2 REAL DENSITY = 10.59300
+* UO2 EFF DENSITY = 10.43750
+* UO2 TEMPERATURE = 941.28998
+* ENRICHMENT = 0.71140
+* COOLANT D2 AT % = 99.222
+* MODERATOR D2 AT % = 99.911
+* NUMBER OF DAYS = 50
+*
+*----
+* Define variables and initialize
+* Burnup paremeters
+* a) Power
+* = 31.9713 kw/kg for 0.0 to 300.0 days
+* b) Burnup time interval Delt
+* = 1 day for 0 to 1 day
+* = 4 days for 1 to 5 days
+* = 5 days for 5 to 10 days
+* = 10 days for 10 to 50 days
+* = 20 days for 50 to 150 days
+* = 50 days for 150 to 300 days
+* c) Days with burnup interval changes
+* = 1.0, 5.0, 10.0, 50.0, 150.0 and 300.0 days
+* d) Burnup control time variables Timei, Timef
+* Timei = initial time
+* Timef = final time
+*----
+REAL
+ Power Delt Timec Timei Timef :=
+ 31.9713 1.0 1.0 0.0 0.0 ;
+*----
+* Define STRUCTURES and MODULES used
+*----
+LINKED_LIST
+ LIBRARY CANDU6S CANDU6F VOLMATS VOLMATF PIJ FLUX BURNUP EDITION
+ DATABASE ISOT ;
+SEQ_BINARY
+ INTLINS INTLINF ;
+SEQ_ASCII
+ database ;
+MODULE
+ GEO: EXCELT: SHI: ASM: FLU: EVO: EDI: COMPO: DELETE: END: ;
+PROCEDURE assertS ;
+*----
+* Depletion data from file WIMS-E format WIMSE
+* Microscopic cross sections from file WIMS-E format WIMSE
+*----
+PROCEDURE TCWE05Lib ;
+INTEGER iedit := 1 ;
+LIBRARY := TCWE05Lib :: <<iedit>> ;
+*----
+* Geometry CANDU6S : 13 regions annular cluster for self-shielding
+* CANDU6F : 31 regions annular cluster for transport
+*----
+CANDU6S := GEO: :: TUBE 5
+ R+ REFL RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 16.12171
+ MIX 1 2 3 4 5
+ CLUSTER ROD1 ROD2 ROD3 ROD4
+ ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000
+ RADIUS 0.00000 0.6122 0.6540 ;
+ ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ;
+ ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ;
+ ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ;
+ ;
+CANDU6F := GEO: CANDU6S :: SPLITR 6 1 1 1 10
+ ::: ROD1 := GEO: ROD1 SPLITR 2 1 ;
+ ::: ROD2 := GEO: ROD2 SPLITR 2 1 ;
+ ::: ROD3 := GEO: ROD3 SPLITR 2 1 ;
+ ::: ROD4 := GEO: ROD4 SPLITR 2 1 ;
+ ;
+*----
+* Create the reactor database
+*----
+DATABASE := COMPO: ::
+ EDIT 5
+ STEP UP 'moderator'
+ COMM 'Multi-parameter reactor database for moderator' ENDC
+ INIT
+ STEP UP 'fuel'
+ COMM 'Multi-parameter reactor database for fuel' ENDC
+ PARA 'BURN' IRRA
+ PARA 'FLUB' FLUB
+ INIT
+ ;
+*----
+* Self-Shielding calculation EXCEL
+* Transport calculation EXCEL
+* Flux calculation for keff
+*----
+VOLMATS INTLINS := EXCELT: CANDU6S ::
+ TITLE 'TCWE05: CANDU-6 ANNULAR POWER= 31.971 FUEL TEMP= 941.29'
+ EDIT 0 MAXR 13
+ ALLG BATCH 100
+ TRAK TISO 5 10.0 SYMM 12 ;
+LIBRARY := SHI: LIBRARY VOLMATS INTLINS :: EDIT 0 NOLJ ;
+VOLMATF INTLINF := EXCELT: CANDU6F ::
+ TITLE 'TCWE05: CANDU-6 ANNULAR POWER= 31.971 FUEL TEMP= 941.29'
+ EDIT 0 MAXR 31
+ ALLG BATCH 100
+ TRAK TISO 5 10.0 SYMM 12 ;
+PIJ := ASM: LIBRARY VOLMATF INTLINF ;
+FLUX := FLU: PIJ LIBRARY VOLMATF ::
+ TYPE K ;
+assertS FLUX :: 'K-INFINITY' 1 1.124888 ;
+EDITION := EDI: LIBRARY VOLMATF FLUX ::
+ COND 4.0 MERGE MIX 0 0 0 0 1 0 0 0 0 0 MICR RES SAVE ON 'moderator' ;
+DATABASE := COMPO: DATABASE EDITION ::
+ EDIT 3
+ STEP UP *
+ ;
+EDITION := EDI: EDITION LIBRARY VOLMATF FLUX ::
+ COND 4.0 MERGE COMP MICR 1 Xe135 SAVE ON 'fuel' ;
+*----
+* Burnup loop: for first step BURNUP is created
+* while for other steps it is modified
+*----
+WHILE Timei Timec < DO
+ EVALUATE Timef := Timei Delt + ;
+ IF Timei 0.0 = THEN
+ BURNUP LIBRARY := EVO: LIBRARY FLUX VOLMATF ::
+ DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ;
+ ELSE
+ BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX VOLMATF ::
+ NOEX DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ;
+ ENDIF ;
+ LIBRARY := SHI: LIBRARY VOLMATS INTLINS :: EDIT 0 NOLJ ;
+ PIJ := DELETE: PIJ ;
+ PIJ := ASM: LIBRARY VOLMATF INTLINF ;
+ FLUX := FLU: FLUX PIJ LIBRARY VOLMATF ::
+ TYPE K ;
+ EDITION := EDI: EDITION LIBRARY VOLMATF FLUX ::
+ SAVE ON 'fuel' ;
+ BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX VOLMATF ::
+ SAVE <<Timef>> DAY POWR <<Power>> ;
+ DATABASE := COMPO: DATABASE EDITION BURNUP LIBRARY ::
+ EDIT 3
+ STEP UP *
+ SET <<Timef>> DAY
+ ;
+*----
+* change delta t for burnup and final time if required
+*----
+ IF Timef Timec = THEN
+ IF Timec 150.0 = THEN
+ EVALUATE Delt Timec := 50.0 300.0 ;
+ ENDIF ;
+ IF Timec 50.0 = THEN
+ EVALUATE Delt Timec := 20.0 150.0 ;
+ ENDIF ;
+ IF Timec 10.0 = THEN
+ EVALUATE Delt Timec := 10.0 50.0 ;
+ ENDIF ;
+ IF Timec 5.0 = THEN
+ EVALUATE Delt Timec := 5.0 10.0 ;
+ ENDIF ;
+ IF Timec 1.0 = THEN
+ EVALUATE Delt Timec := 4.0 5.0 ;
+ ENDIF ;
+ ENDIF ;
+ EVALUATE Timei := Timef ;
+ENDWHILE ;
+assertS FLUX :: 'K-INFINITY' 1 0.9690224 ;
+*----
+* Export and access the database
+*----
+database := DATABASE ;
+
+ISOT := DATABASE :: STEP UP fuel
+ STEP UP MIXTURES STEP AT 1
+ STEP UP CALCULATIONS STEP AT 5
+ STEP UP ISOTOPESLIST STEP AT 1 ;
+assertS ISOT :: 'NWT0' 1 9.663608E0 ;
+assertS ISOT :: 'NWT0' 2 2.579294E1 ;
+ISOT := DELETE: ISOT ;
+
+ISOT := DATABASE :: STEP UP moderator
+ STEP UP MIXTURES STEP AT 1
+ STEP UP CALCULATIONS STEP AT 1
+ STEP UP ISOTOPESLIST STEP AT 1 ;
+assertS ISOT :: 'NWT0' 1 9.295185E-2 ;
+assertS ISOT :: 'NWT0' 2 3.190355E-1 ;
+
+INTLINF INTLINS := DELETE: INTLINF INTLINS ;
+ECHO "test TCWE05 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE05Lib.c2m b/Dragon/data/twimsE_proc/TCWE05Lib.c2m
new file mode 100644
index 0000000..9114411
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE05Lib.c2m
@@ -0,0 +1,68 @@
+*----
+* Procedure TCWE05Lib
+* Create Library for test CASE TCWE05
+* Calling :
+* LIBRARY := TCWE05Lib :: iedit ;
+* with :
+* LIBRARY = Linked list containing the result of LIB: for
+* TCWE05
+* iprint = print level for LIB: module
+*----
+* Define PARAMETERS,STRUCTURES and MODULES used
+*----
+PARAMETER LIBRARY :: ::: LINKED_LIST LIBRARY ; ;
+MODULE LIB: END: ;
+*----
+* Define and read LIB: EDIT option
+INTEGER iedit ;
+:: >>iedit<< ;
+*----
+* Depletion data from file WIMS-E format WIMSE
+* Microscopic cross sections from file WIMS-E format WIMSE
+*----
+LIBRARY := LIB: ::
+ EDIT <<iedit>>
+ NMIX 17 CTRA WIMS
+ DEPL LIB: WIMSE FIL: WIMS-E
+ MIXS LIB: WIMSE FIL: WIMS-E
+ MIX 1 560.66 0.81212 O16 = '8016' 7.99449E-1
+ D2D2O = '1002' 1.99768E-1 H1H2O = '1001' 7.83774E-4
+ MIX 2 560.66 6.57 Nb93 = '41093' 2.50000
+ B10 = '5010' 4.17900E-05 B11 = '5011' 1.68210E-04
+ Zr90 = '40090' 5.01638E+01 Zr91 = '40091' 1.09395E+01 2
+ Zr92 = '40092' 1.67213E+01 Zr94 = '40094' 1.69455E+01
+ Zr96 = '40096' 2.73000E+00
+ MIX 3 345.66 0.0014 He4 = '2004' 1.00000E+2
+ MIX 4 345.66 6.44 Fe56 = '26056' 1.60000E-1
+ Ni58 = '28058' 6.00000E-2 Cr52 = '24052' 1.10000E-1
+ B10 = '5010' 6.16900E-05 B11 = '5011' 2.48310E-04
+ Zr90 = '40090' 5.13008E+01 Zr91 = '40091' 1.11875E+01 2
+ Zr92 = '40092' 1.71003E+01 Zr94 = '40094' 1.73296E+01
+ Zr96 = '40096' 2.79188E+00
+ MIX 5 345.66 1.082885 O16 = '8016' 7.98895E-1
+ D2D2O = '1002' 2.01016E-1 H1H2O = '1001' 8.96000E-5
+ MIX 6 941.29 10.4375010 O16 = '8016' 1.18473E+1
+ Xe135 = '54135' 0.0
+ U235 = '92235' 6.27118E-1 1
+ U238 = '92238' 8.75256E+1 1
+ U236 = '92236' 0.0 1
+ Pu239 = '94239' 0.0 1
+ MIX 7 COMB 6 1.0
+ MIX 8 COMB 6 1.0
+ MIX 9 COMB 6 1.0
+ MIX 10 560.66 6.44 Fe56 = '26056' 1.60000E-1
+ Ni58 = '28058' 6.00000E-2 Cr52 = '24052' 1.10000E-1
+ B10 = '5010' 6.16900E-05 B11 = '5011' 2.48310E-04
+ Zr90 = '40090' 5.13008E+01 Zr91 = '40091' 1.11875E+01 2
+ Zr92 = '40092' 1.71003E+01 Zr94 = '40094' 1.73296E+01
+ Zr96 = '40096' 2.79188E+00
+ MIX 11 COMB 10 1.0
+ MIX 12 COMB 10 1.0
+ MIX 13 COMB 10 1.0
+ MIX 14 COMB 1 1.0
+ MIX 15 COMB 1 1.0
+ MIX 16 COMB 1 1.0
+ MIX 17 COMB 1 1.0
+ ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE06.c2m b/Dragon/data/twimsE_proc/TCWE06.c2m
new file mode 100644
index 0000000..287c539
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE06.c2m
@@ -0,0 +1,155 @@
+*----
+* TEST CASE TCWE06
+* CANDU-6 CARTESIAN CELL
+* WIMS-E Library
+* STAINLESS STELL RODS IN 3D SUPERCELL
+*
+* REF: R. Roy et al. Ann. Nucl. Energy 21, 115 (1994)
+*
+*----
+* Define STRUCTURES and MODULES used
+*----
+LINKED_LIST
+ LIBRARY LIBRARY2 CANDU6F CANDU6S TRACK SYS FLUX EDITION BCO BCI
+ SSRODS ;
+SEQ_BINARY
+ INTLIN ;
+MODULE
+ GEO: EXCELT: LIB: SHI: ASM: FLU: EDI: DELETE: END: ;
+PROCEDURE assertS ;
+*----
+* Depletion data from file WIMS-E format WIMSE
+* Microscopic cross sections from file WIMS-E format WIMSE
+*----
+PROCEDURE TCWE05Lib ;
+INTEGER iedit := 1 ;
+LIBRARY := TCWE05Lib :: <<iedit>> ;
+*----
+* CELL CALCULATION
+* Geometry CANDU6S : 14 regions Cartesian cluster for self-shielding
+* CANDU6F : 32 regions Cartesian cluster for transport
+* BCO : 48 regions 3D Cartesian geometry
+* BCI : 48 regions 3D Cartesian geometry
+*----
+CANDU6S := GEO: :: CARCEL 5
+ X+ REFL X- REFL MESHX -14.2875 14.2875
+ Y+ REFL Y- REFL MESHY -14.2875 14.2875
+ RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 14.00
+ MIX 1 2 3 4 5 5
+ CLUSTER ROD1 ROD2 ROD3 ROD4
+ ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000
+ RADIUS 0.00000 0.6122 0.6540 ;
+ ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ;
+ ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ;
+ ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ;
+ ;
+CANDU6F := GEO: CANDU6S :: SPLITR 6 1 1 1 10
+ ::: ROD1 := GEO: ROD1 SPLITR 2 1 ;
+ ::: ROD2 := GEO: ROD2 SPLITR 2 1 ;
+ ::: ROD3 := GEO: ROD3 SPLITR 2 1 ;
+ ::: ROD4 := GEO: ROD4 SPLITR 2 1 ;
+ ;
+*----
+* Self-Shielding calculation EXCEL
+* Transport calculation EXCEL
+* Flux calculation for keff
+*----
+TRACK INTLIN := EXCELT: CANDU6S ::
+ TITLE 'TCWE06: CANDU-6 CARTESIAN FUEL TEMP= 941.29'
+ EDIT 0 MAXR 14
+ ALLG BATCH 100
+ TRAK TISO 29 20.0 SYMM 4 ;
+LIBRARY := SHI: LIBRARY TRACK INTLIN :: EDIT 0 NOLJ ;
+TRACK INTLIN := DELETE: TRACK INTLIN ;
+TRACK INTLIN := EXCELT: CANDU6F ::
+ TITLE 'TCWE06: CANDU-6 CARTESIAN FUEL TEMP= 941.29'
+ EDIT 0 MAXR 32
+ ALLG BATCH 100
+ TRAK TISO 29 20.0 SYMM 4 ;
+SYS := ASM: LIBRARY TRACK INTLIN ::
+ EDIT 0 ;
+FLUX := FLU: SYS LIBRARY TRACK ::
+ TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.124544 ;
+*----
+* Microscopic cross sections from WIMSLIB for stainless steel
+* MIX 1,2,3 from EDI: : fuel, structure material and moderator
+* MIX 4 stainless steel rods
+*----
+EDITION := EDI: LIBRARY TRACK FLUX ::
+ EDIT 0 MERGE MIX 1 2 2 2 3 1 1 1 1 1 SAVE ON SSRODS ;
+SSRODS := EDITION :: STEP UP SSRODS STEP UP MACROLIB ;
+LIBRARY2 := LIB: SSRODS ::
+ EDIT 0 NMIX 4 CTRA WIMS
+ MIXS LIB: WIMSE FIL: WIMS-E
+ MIX 4 345.66 Fe56 = '26056' 6.19027E-2
+ Cr52 = '24052' 1.56659E-3 Ni58 = '28058' 6.83337E-3
+ Si29 = '14029' 7.79072E-4 C12 = '6000' 1.46552E-4
+ Mn55 = '25055' 1.25431E-3
+ ;
+EDITION TRACK INTLIN SYS FLUX := DELETE:
+ EDITION TRACK INTLIN SYS FLUX ;
+*----
+* SUPERCELL CALCULATION
+* Geometry BCO : 27 regions 3D Cartesian geometry with rods out
+* BCI : 27 regions 3D Cartesian geometry with rods in
+*----
+BCO := GEO: :: CAR3D 3 2 2
+ X- REFL X+ SYME Y- REFL Y+ SYME Z- REFL Z+ SYME
+ CELL M MX MX MX FXY MXY M MX BX MX FXY BXY
+ TURN A A A F A A A A A F A A
+ ::: M := GEO: CAR3D 1 1 1 MIX 3
+ MESHX 0.0 7.14375
+ MESHY 0.0 7.14375
+ MESHZ -8.25500 8.25500 SPLITZ 2 ;
+ ::: MX := GEO: M MESHX -7.14375 +7.14375 SPLITX 2 ;
+ ::: MXY := GEO: MX MESHY -7.14375 +7.14375 SPLITY 2 ;
+ ::: BX := GEO: CARCELY 2 1 MIX 3 3 3
+ MESHX -7.14375 7.14375 SPLITX 2
+ MESHY 0.0 7.14375
+ MESHZ -8.25500 8.25500 SPLITZ 2
+ RADIUS 0.0 3.5100 3.8100 ;
+ ::: BXY := GEO: BX MESHY -7.14375 +7.14375 SPLITY 2 ;
+ ::: FXY := GEO: CARCELZ 2 1 MIX 1 2 3
+ MESHX -7.14375 7.14375 SPLITX 2
+ MESHY -7.14375 7.14375 SPLITY 2
+ MESHZ -8.25500 8.25500 SPLITZ 2
+ RADIUS 0.0 5.16890 6.58750 ;
+ ;
+BCI := GEO: BCO ::
+ ::: BX := GEO: BX MIX 3 4 3 ;
+ ::: BXY := GEO: BXY MIX 3 4 3 ;
+ ;
+*----
+* Transport calculation EXCEL
+* Flux calculation for keff
+* Homogenized properties for rod out
+*----
+TRACK INTLIN := EXCELT: BCO ::
+ EDIT 0 MAXR 40 TRAK TISO 2 1.0 ;
+SYS := ASM: LIBRARY2 TRACK INTLIN ::
+ EDIT 0 ;
+FLUX := FLU: SYS LIBRARY2 TRACK ::
+ TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.131641 ;
+EDITION := EDI: LIBRARY2 TRACK FLUX ::
+ EDIT 2 MERG COMP COND 4.0 SAVE ON 'NOBC' ;
+SYS TRACK INTLIN := DELETE: SYS TRACK INTLIN ;
+*----
+* Transport calculation EXCEL
+* Flux calculation for keff
+* Homogenized properties for rod in
+*----
+TRACK INTLIN := EXCELT: BCI ::
+ EDIT 0 MAXR 40 TRAK TISO 2 1.0 ;
+SYS := ASM: LIBRARY2 TRACK INTLIN ::
+ EDIT 0 ;
+FLUX := FLU: FLUX SYS LIBRARY2 TRACK ::
+ TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.046070 ;
+EDITION := EDI: EDITION LIBRARY2 TRACK FLUX ::
+ EDIT 2 MERG COMP COND 4.0 STAT DELS REFE 'NOBC' ;
+TRACK INTLIN SYS := DELETE: TRACK INTLIN SYS ;
+ECHO "test TCWE06 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE07.c2m b/Dragon/data/twimsE_proc/TCWE07.c2m
new file mode 100644
index 0000000..afbf9c3
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE07.c2m
@@ -0,0 +1,142 @@
+*----
+* TEST CASE TCWE07
+* CANDU-6 CARTESIAN CELL
+* WIMS-E Library
+* TEST VARIOUS LEAKAGE OPTIONS
+*
+*----
+* Define STRUCTURES and MODULES used
+*----
+LINKED_LIST
+ LIBRARY CANDU6S CANDU6T CANDU6SV CANDU6TV TRACK
+ SYS FLUX EDITION ;
+MODULE
+ GEO: EXCELT: LIB: SHI: ASM: FLU: EDI: DELETE: END: ;
+SEQ_BINARY
+ INTLIN ;
+PROCEDURE assertS ;
+*----
+* Depletion data from file WIMS-E format WIMSE
+* Microscopic cross sections from file WIMS-E format WIMSE
+*----
+PROCEDURE TCWE05Lib ;
+INTEGER iedit := 1 ;
+LIBRARY := TCWE05Lib :: <<iedit>> ;
+*----
+* Geometry CANDU6S : GEOMETRY FOR SELF-SHIELDING (NO VOID)
+* CANDU6F : GEOMETRY FOR TRANSPORT (NO VOID)
+* CANDU6FV: GEOMETRY FOR TRANSPORT (COOLANT VOID)
+* CANDU6FV: GEOMETRY FOR TRANSPORT (COOLANT VOID)
+*----
+CANDU6S := GEO: :: CARCEL 5
+ X+ REFL X- REFL MESHX -14.2875 14.2875
+ Y+ REFL Y- REFL MESHY -14.2875 14.2875
+ RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 14.00
+ MIX 1 2 3 4 5 5
+ CLUSTER ROD1 ROD2 ROD3 ROD4
+ ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000
+ RADIUS 0.00000 0.6122 0.6540 ;
+ ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ;
+ ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ;
+ ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ;
+ ;
+CANDU6T := GEO: CANDU6S :: SPLITR 6 1 1 1 10
+ ::: ROD1 := GEO: ROD1 SPLITR 2 1 ;
+ ::: ROD2 := GEO: ROD2 SPLITR 2 1 ;
+ ::: ROD3 := GEO: ROD3 SPLITR 2 1 ;
+ ::: ROD4 := GEO: ROD4 SPLITR 2 1 ;
+ ;
+CANDU6SV := GEO: CANDU6S :: MIX 0 2 3 4 5 5 ;
+CANDU6TV := GEO: CANDU6SV :: SPLITR 6 1 1 1 10
+ ::: ROD1 := GEO: ROD1 SPLITR 2 1 ;
+ ::: ROD2 := GEO: ROD2 SPLITR 2 1 ;
+ ::: ROD3 := GEO: ROD3 SPLITR 2 1 ;
+ ::: ROD4 := GEO: ROD4 SPLITR 2 1 ;
+ ;
+*----
+* CASE WITH NO VOID
+* Self-Shielding calculation EXCEL
+* Transport calculation EXCEL
+* Flux TYPE K AND B WITH VARIOUS LEAKAGE OPTIONS
+*----
+TRACK INTLIN := EXCELT: CANDU6S ::
+ TITLE 'TCWE07: CANDU-6 CARTESIAN FUEL TEMP= 941.29'
+ EDIT 0 MAXR 14
+ ALLG BATCH 100
+ TRAK TISO 7 20.0 SYMM 4 ;
+LIBRARY := SHI: LIBRARY TRACK INTLIN :: EDIT 0 NOLJ ;
+TRACK INTLIN := DELETE: TRACK INTLIN ;
+TRACK INTLIN := EXCELT: CANDU6T ::
+ TITLE 'TCWE07: CANDU-6 CARTESIAN FUEL TEMP= 941.29'
+ EDIT 0 MAXR 32 ANIS 2
+ ALLG BATCH 100
+ TRAK TISO 7 20.0 SYMM 4 ;
+SYS := ASM: LIBRARY TRACK INTLIN ::
+ EDIT 0 PIJK ;
+FLUX := FLU: SYS LIBRARY TRACK ::
+ TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.124517 ;
+EDITION := EDI: LIBRARY TRACK FLUX ::
+ EDIT 3 SAVE COND 4.0 TAKE REGI 1 4 7 10 16 24 ;
+FLUX := FLU: FLUX SYS LIBRARY TRACK ::
+ TYPE B B1 PNL ;
+assertS FLUX :: 'K-INFINITY' 1 1.116083 ;
+EDITION := EDI: EDITION LIBRARY TRACK FLUX ::
+ EDIT 3 SAVE COND 4.0 TAKE REGI 1 4 7 10 16 24 ;
+FLUX := FLU: FLUX SYS LIBRARY TRACK ::
+ TYPE B B1 TIBERE ;
+assertS FLUX :: 'K-INFINITY' 1 1.116019 ;
+EDITION := EDI: EDITION LIBRARY TRACK FLUX ::
+ EDIT 3 SAVE COND 4.0 TAKE REGI 1 4 7 10 16 24 ;
+FLUX SYS := DELETE: FLUX SYS ;
+SYS := ASM: LIBRARY TRACK INTLIN ::
+ EDIT 0 PIJ ECCO ;
+FLUX := FLU: SYS LIBRARY TRACK ::
+ TYPE B B1 ECCO ;
+assertS FLUX :: 'K-INFINITY' 1 1.116024 ;
+EDITION := EDI: EDITION LIBRARY TRACK FLUX ::
+ EDIT 3 SAVE COND 4.0 TAKE REGI 1 4 7 10 16 24 ;
+FLUX SYS TRACK INTLIN := DELETE: FLUX SYS TRACK INTLIN ;
+*----
+* CASE WITH COOLANT VOIDED
+* Self-Shielding calculation EXCEL
+* Transport calculation EXCEL
+* Flux TYPE K AND B WITH VARIOUS LEAKAGE OPTIONS
+*----
+TRACK INTLIN := EXCELT: CANDU6SV ::
+ TITLE 'TCWE07: CANDU-6 CARTESIAN FUEL TEMP= 941.29'
+ EDIT 0 MAXR 14 TRAK TISO 7 20.0 SYMM 4 ;
+LIBRARY := SHI: LIBRARY TRACK INTLIN :: EDIT 0 NOLJ ;
+TRACK INTLIN := DELETE: TRACK INTLIN ;
+TRACK INTLIN := EXCELT: CANDU6TV ::
+ TITLE 'TCWE07: CANDU-6 CARTESIAN FUEL TEMP= 941.29'
+ EDIT 0 MAXR 32 ANIS 2 TRAK TISO 7 20.0 SYMM 4 ;
+SYS := ASM: LIBRARY TRACK INTLIN ::
+ EDIT 0 PIJK ;
+FLUX := FLU: SYS LIBRARY TRACK ::
+ TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.144657 ;
+EDITION := EDI: EDITION LIBRARY TRACK FLUX ::
+ EDIT 3 SAVE COND 4.0 TAKE REGI 1 4 7 10 16 24 ;
+FLUX := FLU: FLUX SYS LIBRARY TRACK ::
+ TYPE B B1 PNL ;
+assertS FLUX :: 'K-INFINITY' 1 1.136610 ;
+EDITION := EDI: EDITION LIBRARY TRACK FLUX ::
+ EDIT 3 SAVE COND 4.0 TAKE REGI 1 4 7 10 16 24 ;
+FLUX := FLU: FLUX SYS LIBRARY TRACK ::
+ TYPE B B1 TIBERE ;
+assertS FLUX :: 'K-INFINITY' 1 1.136823 ;
+EDITION := EDI: EDITION LIBRARY TRACK FLUX ::
+ EDIT 3 SAVE COND 4.0 TAKE REGI 1 4 7 10 16 24 ;
+FLUX SYS := DELETE: FLUX SYS ;
+SYS := ASM: LIBRARY TRACK INTLIN ::
+ EDIT 0 PIJ ECCO ;
+FLUX := FLU: SYS LIBRARY TRACK ::
+ TYPE B B1 ECCO ;
+assertS FLUX :: 'K-INFINITY' 1 1.136833 ;
+EDITION := EDI: EDITION LIBRARY TRACK FLUX ::
+ EDIT 3 SAVE COND 4.0 TAKE REGI 1 4 7 10 16 24 ;
+FLUX SYS TRACK INTLIN := DELETE: FLUX SYS TRACK INTLIN ;
+ECHO "test TCWE07 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE08.c2m b/Dragon/data/twimsE_proc/TCWE08.c2m
new file mode 100644
index 0000000..cb329ad
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE08.c2m
@@ -0,0 +1,121 @@
+*----
+* TEST CASE TCWW08
+* HOMOGENEOUS DEPLETION CASE
+* WIMS-E Library
+*
+* REF: None
+*
+*----
+* Define variables
+* Burnup paremeters
+* a) Power
+* = 600.0 kw/kg for 0.0 to 1000.0 days
+* = 0.0 kw/kg for 1000.0 to 2000.0 days
+* b) Burnup time interval Delt
+* = 10 days for 0 to 50 days
+* = 50 days for 50 to 500 days
+* = 100 days for 500 to 1000 days
+* = 1000 days for 1000 to 2000 days
+* c) Editing time Timec
+* = 0.0, 50.0, 500.0, 1000.0 and 2000.0 days
+* d) Burnup control time variables Timei, Timef, TotalTime
+* Timei = initial time
+* Timef = final time
+* TotalTime = Final time reached
+* d) Print variable Iprint
+* = 1 reduced print
+* = 3 full print
+*----
+REAL
+ Power Delt Timec Timei Timef TotalTime :=
+ 600.0 10.0 50.0 0.0 0.0 2000.0 ;
+INTEGER
+ Iprint := 1 ;
+*----
+* Define STRUCTURES and MODULES used
+*----
+LINKED_LIST
+ LIBRARY HOM TRACK PIJ FLUX BURNUP EDITION ;
+MODULE
+ GEO: SYBILT: LIB: SHI: ASM: FLU: EVO: EDI: DELETE: END: ;
+PROCEDURE assertS ;
+*----
+* Depletion data from file WIMS-E format WIMSE
+* Microscopic cross sections from file WIMS-E format WIMSE
+*----
+LIBRARY := LIB: ::
+ NMIX 1 CTRA WIMS
+ DEPL LIB: WIMSE FIL: WIMS-E
+ MIXS LIB: WIMSE FIL: WIMS-E
+ MIX 1 300.0 H1 = '1001' 2.00000E+1
+ U235 = '92235' 1.0 1
+ U236 = '92238' 0.0 1
+ ;
+*----
+* Geometry HOM : Homogeneous geometry
+*----
+HOM := GEO: :: HOMOGE
+ MIX 1 ;
+*----
+* Self-Shielding calculation SYBIL
+* Transport calculation SYBIL
+* Flux calculation for keff
+*----
+TRACK := SYBILT: HOM ::
+ TITLE 'TCWW08: HOMOGENEOUS BENCHMARK WITH BURNUP' ;
+LIBRARY := SHI: LIBRARY TRACK :: EDIT 0 NOLJ ;
+PIJ := ASM: LIBRARY TRACK ;
+FLUX := FLU: PIJ LIBRARY TRACK ::
+ TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.843108 ;
+EDITION := EDI: LIBRARY TRACK FLUX ::
+ COND 4.0 MERGE COMP SAVE ;
+*----
+* Burnup loop: for first step BURNUP is created
+* while for other steps it is modified
+*----
+WHILE Timei TotalTime < DO
+ EVALUATE Timef := Timei Delt + ;
+ IF Timef Timec = THEN
+ EVALUATE Iprint := 3 ;
+ ELSE
+ EVALUATE Iprint := 1 ;
+ ENDIF ;
+ IF Timei 0.0 = THEN
+ BURNUP LIBRARY := EVO: LIBRARY FLUX TRACK ::
+ EDIT <<Iprint>> DEPL <<Timei>> <<Timef>> DAY
+ POWR <<Power>> ;
+ ELSE
+ BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX TRACK ::
+ EDIT <<Iprint>> NOEX DEPL <<Timei>> <<Timef>> DAY
+ POWR <<Power>> ;
+ ENDIF ;
+ LIBRARY := SHI: LIBRARY TRACK :: EDIT 0 NOLJ ;
+ PIJ := DELETE: PIJ ;
+ PIJ := ASM: LIBRARY TRACK ;
+ FLUX := FLU: FLUX PIJ LIBRARY TRACK ::
+ TYPE K ;
+ IF Iprint 3 = THEN
+ EDITION := EDI: EDITION LIBRARY TRACK FLUX ::
+ SAVE ;
+ ENDIF ;
+*----
+* change delta t for burnup, final time and power if required
+*----
+ IF Timef Timec = THEN
+ IF Timec 1000.0 = THEN
+ EVALUATE Power Delt Timec := 0.0 1000.0 2000.0 ;
+ ENDIF ;
+ IF Timec 500.0 = THEN
+ EVALUATE Delt Timec := 100.0 1000.0 ;
+ ENDIF ;
+ IF Timec 50.0 = THEN
+ EVALUATE Delt Timec := 50.0 500.0 ;
+ ENDIF ;
+ ENDIF ;
+ EVALUATE Timei := Timef ;
+ENDWHILE ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.322846 ;
+ECHO "test TCWE08 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE09.c2m b/Dragon/data/twimsE_proc/TCWE09.c2m
new file mode 100644
index 0000000..a50e6d1
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE09.c2m
@@ -0,0 +1,74 @@
+*----
+* TEST CASE TCWE09
+* MOSTELLER BENCHMARK FOR 1-D ANNULAR CELL
+* WIMS-E Library
+* REFLECTIVE AND VOID BC
+*
+* REF: None
+*
+*----
+* Define STRUCTURES and MODULES used
+*----
+LINKED_LIST
+ MOSTELA MOSTELV VOLMAT LIBRARY PIJ FLUX OUT ;
+MODULE
+ LIB: GEO: SYBILT: SHI: ASM: FLU: EDI: DELETE: END: ;
+PROCEDURE assertS ;
+*----
+* Microscopic cross section from file WIMS-E format WIMSE
+*----
+LIBRARY := LIB: ::
+ NMIX 3 CTRA WIMS
+ MIXS LIB: WIMSE FIL: WIMS-E
+ MIX 1 600.0 O16 = '8016' 4.61309E-2
+ U235 = '92235' 1.66078E-4 1
+ U238 = '92238' 2.28994E-2 1
+ MIX 2 600.0
+ Zr90 = '40090' 1.97179E-02 Zr91 = '40091' 4.29999E-03 2
+ Zr92 = '40092' 6.57262E-03 Zr94 = '40094' 6.66076E-03
+ Zr96 = '40096' 1.07308E-03
+ MIX 3 600.0
+ H1H2O = '1001' 4.42326E-2 O16H2O = '8016' 2.21163E-2
+ B10 = '5010' 2.03245E-06 B11 = '5011' 8.18085E-06
+ ;
+*----
+* Geometry MOSTELA : Annular cell with reflective BC
+* MOSTELV : Annular cell with void BC
+*----
+MOSTELA := GEO: :: TUBE 3
+ RADIUS 0.0 0.39306 0.45802 0.71206
+ SPLITR 2 1 1
+ MIX 1 2 3
+ R+ REFL ;
+MOSTELV := GEO: MOSTELA ::
+ R+ VOID ;
+*----
+* Self-Shielding calculation SYBIL
+* Transport calculation SYBIL
+* Flux calculation for keff
+*----
+VOLMAT := SYBILT: MOSTELA ::
+ TITLE 'TCWE09: SYBIL TRACK MOSTELLER BENCHMARK REFLECTIVE BC '
+ MAXR 4 QUA1 5 ;
+LIBRARY := SHI: LIBRARY VOLMAT :: EDIT 0 NOLJ ;
+PIJ := ASM: LIBRARY VOLMAT ;
+FLUX := FLU: PIJ LIBRARY VOLMAT ::
+ TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 0.8289741 ;
+OUT := EDI: LIBRARY VOLMAT FLUX ::
+ EDIT 4 MERG MIX 1 2 3 COND 4.0 SAVE ;
+PIJ VOLMAT := DELETE: PIJ VOLMAT ;
+VOLMAT := SYBILT: MOSTELV ::
+ TITLE 'TCWE09: SYBIL TRACK MOSTELLER BENCHMARK VOID BC '
+ MAXR 4 QUA1 5 ;
+PIJ := ASM: LIBRARY VOLMAT ;
+FLUX := FLU: FLUX PIJ LIBRARY VOLMAT ::
+ TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 9.918292E-03 ;
+OUT := EDI: OUT LIBRARY VOLMAT FLUX ::
+ EDIT 4 MERG MIX 1 2 3 COND 4.0 SAVE ;
+OUT FLUX PIJ LIBRARY VOLMAT := DELETE:
+ OUT FLUX PIJ LIBRARY VOLMAT ;
+ECHO "test TCWE09 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE10.c2m b/Dragon/data/twimsE_proc/TCWE10.c2m
new file mode 100644
index 0000000..0005914
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE10.c2m
@@ -0,0 +1,104 @@
+*----
+* TEST CASE TCWE10
+* MOSTELLER BENCHMARK FOR 1-D ANNULAR CELL
+* WIMS-E Library
+* FIXED SOURCE PROBLEM IN MULTIPLICATIVE MEDIA
+*
+* REF: None
+*
+*----
+* Define STRUCTURES and MODULES used
+*----
+LINKED_LIST
+ MOSTELA VOLMAT LIBRARY PIJ FLUX OUT ;
+MODULE
+ LIB: MAC: GEO: SYBILT: SHI: ASM: FLU: EDI: DELETE: END: ;
+PROCEDURE assertS assertV ;
+*----
+* Microscopic cross section from file WIMS-E format WIMSE
+* Fixed source of 1.0E5 in group 6
+*----
+LIBRARY := LIB: ::
+ EDIT 0 NMIX 3 CTRA WIMS
+ MIXS LIB: WIMSE FIL: WIMS-E
+ MIX 1 600.0 O16 = '8016' 4.61309E-2
+ U235 = '92235' 1.66078E-4 1
+ U238 = '92238' 2.28994E-2 1
+ MIX 2 600.0
+ Zr90 = '40090' 1.97179E-02 Zr91 = '40091' 4.29999E-03 2
+ Zr92 = '40092' 6.57262E-03 Zr94 = '40094' 6.66076E-03
+ Zr96 = '40096' 1.07308E-03
+ MIX 3 600.0
+ H1H2O = '1001' 4.42326E-2 O16H2O = '8016' 2.21163E-2
+ B10 = '5010' 2.03245E-06 B11 = '5011' 8.18085E-06
+ ;
+LIBRARY := MAC: LIBRARY ::
+ EDIT 0
+ READ INPUT
+ MIX 3 FIXE
+ 0.0 0.0 0.0 0.0 0.0 1.0E+5 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0
+ ;
+*----
+* Geometry MOSTELA : Annular cell with reflective BC
+*----
+MOSTELA := GEO: :: TUBE 3
+ RADIUS 0.0 0.39306 0.45802 0.71206
+ SPLITR 2 1 1
+ MIX 1 2 3
+ R+ REFL ;
+*----
+* Self-Shielding calculation SYBIL
+* Transport calculation SYBIL
+* Flux calculation for keff
+*----
+VOLMAT := SYBILT: MOSTELA ::
+ TITLE 'TCWE10: SYBIL TRACK MOSTELLER BENCHMARK'
+ MAXR 4 QUA1 5 ;
+LIBRARY := SHI: LIBRARY VOLMAT :: EDIT 0 NOLJ ;
+PIJ := ASM: LIBRARY VOLMAT ;
+FLUX := FLU: PIJ LIBRARY VOLMAT ::
+ TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 0.8289741 ;
+OUT := EDI: LIBRARY VOLMAT FLUX ::
+ EDIT 4 MERG MIX 1 2 3 COND 4.0 ;
+FLUX := DELETE: FLUX ;
+FLUX := FLU: PIJ LIBRARY VOLMAT :: EDIT 99
+ TYPE S EXTE 30 UNKT 1.0E-3 ;
+assertV FLUX :: 'FLUX' (*GROUP*) 10 (*REGION*) 4 2.435309E+05 ;
+OUT := EDI: OUT LIBRARY VOLMAT FLUX ::
+ EDIT 4 MERG MIX 1 2 3 COND 4.0 ;
+OUT FLUX PIJ LIBRARY VOLMAT := DELETE:
+ OUT FLUX PIJ LIBRARY VOLMAT ;
+ECHO "test TCWE10 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE11.c2m b/Dragon/data/twimsE_proc/TCWE11.c2m
new file mode 100644
index 0000000..23a73cf
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE11.c2m
@@ -0,0 +1,171 @@
+*----
+* TEST CASE TCWE11
+* CANDU-6 ANNULAR CELL
+* WIMS-E Library
+* TWO GROUP BURNUP
+* POWER (KW) = 615.00000
+* BURN POWER (KW/KG) = 31.97130
+* URANIUM MASS = 19.23600
+* UO2 REAL DENSITY = 10.59300
+* UO2 EFF DENSITY = 10.43750
+* UO2 TEMPERATURE = 941.28998
+* ENRICHMENT = 0.71140
+* COOLANT D2 AT % = 99.222
+* MODERATOR D2 AT % = 99.911
+* NUMBER OF DAYS = 50
+*
+*----
+* Define variables
+* Burnup paremeters
+* a) Power
+* = 31.9713 kw/kg for 0.0 to 300.0 days
+* b) 69 Groups Burnup time interval Delt
+* = 300 day for 0 to 300 day
+* c) 2 Groups Burnup time interval Delt
+* = 1 day for 0 to 1 day
+* = 4 days for 1 to 5 days
+* = 5 days for 5 to 10 days
+* = 10 days for 10 to 50 days
+* = 20 days for 50 to 150 days
+* = 50 days for 150 to 300 days
+* c) Days with burnup interval changes
+* = 1.0, 5.0, 10.0, 50.0, 150.0 and 300.0 days
+* d) Burnup control time variables Timei, Timef
+* Timei = initial time
+* Timef = final time
+*----
+REAL
+ Power Delt Timec Timei Timef :=
+ 31.9713 1.0 300.0 0.0 0.0 ;
+*----
+* Define STRUCTURES and MODULES used
+*----
+LINKED_LIST LIBRARY ;
+LINKED_LIST
+ CANDU6S CANDU6F VOLMATS VOLMATF PIJ FLUX BURNUP EDITION ;
+SEQ_BINARY
+ INTLINS INTLINF ;
+SEQ_ASCII
+ res ;
+MODULE
+ GEO: EXCELT: LIB: SHI: ASM: FLU: EVO: EDI: SPH: DELETE: END: ;
+PROCEDURE assertS ;
+*----
+* Depletion data from file WIMS-E format WIMSE
+* Microscopic cross sections from file WIMS-E format WIMSE
+*----
+PROCEDURE TCWE05Lib ;
+INTEGER iedit := 1 ;
+LIBRARY := TCWE05Lib :: <<iedit>> ;
+*----
+* Geometry CANDU6S : 13 regions annular cluster for self-shielding
+* CANDU6F : 31 regions annular cluster for transport
+*----
+CANDU6S := GEO: :: TUBE 5
+ R+ REFL RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 16.12171
+ MIX 1 2 3 4 5
+ CLUSTER ROD1 ROD2 ROD3 ROD4
+ ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000
+ RADIUS 0.00000 0.6122 0.6540 ;
+ ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ;
+ ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ;
+ ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ;
+ ;
+CANDU6F := GEO: CANDU6S :: SPLITR 6 1 1 1 10
+ ::: ROD1 := GEO: ROD1 SPLITR 2 1 ;
+ ::: ROD2 := GEO: ROD2 SPLITR 2 1 ;
+ ::: ROD3 := GEO: ROD3 SPLITR 2 1 ;
+ ::: ROD4 := GEO: ROD4 SPLITR 2 1 ;
+ ;
+*----
+* Self-Shielding calculation EXCEL
+* Transport calculation EXCEL
+* Flux calculation for keff
+*----
+VOLMATS INTLINS := EXCELT: CANDU6S ::
+ TITLE 'TCWE11: FEW GROUP BURNUP / SELF-SHIELDING TRACKING'
+ EDIT 0 MAXR 13
+ ALLG BATCH 100
+ TRAK TISO 5 10.0 SYMM 12 ;
+LIBRARY := SHI: LIBRARY VOLMATS INTLINS :: EDIT 0 NOLJ ;
+VOLMATF INTLINF := EXCELT: CANDU6F ::
+ TITLE 'TCWE11: FEW GROUP BURNUP / TRANSPORT TRACKING'
+ EDIT 0 MAXR 31
+ ALLG BATCH 100
+ TRAK TISO 5 10.0 SYMM 12 ;
+PIJ := ASM: LIBRARY VOLMATF INTLINF ;
+FLUX := FLU: PIJ LIBRARY VOLMATF ::
+ TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.124889 ;
+EDITION := EDI: LIBRARY VOLMATF FLUX CANDU6F ::
+ MERG REGI 6 6 10 7 7 10 1 1 8 8 10 1
+ 1 9 9 10 1 1 2 3 4 5 5 5
+ 5 5 5 5 5 5 5
+ COND 4.0 MICR ALL SAVE
+ MGEO CANDU6F
+ ;
+EDITION := SPH: EDITION VOLMATF INTLINF ;
+*----
+* 69 group Burnup
+*----
+BURNUP LIBRARY := EVO: LIBRARY FLUX VOLMATF ::
+ EDIT 3 EXPM 200.0 DEPL <<Timei>> <<Timec>> DAY POWR <<Power>> ;
+LIBRARY := SHI: LIBRARY VOLMATS INTLINS :: EDIT 0 NOLJ ;
+PIJ := DELETE: PIJ ;
+PIJ := ASM: LIBRARY VOLMATF INTLINF ;
+FLUX := FLU: FLUX PIJ LIBRARY VOLMATF ::
+ TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 0.9580293 ;
+EDITION := EDI: EDITION LIBRARY VOLMATF FLUX CANDU6F ::
+ MGEO CANDU6F
+ ;
+EDITION := SPH: EDITION VOLMATF INTLINF ;
+BURNUP FLUX PIJ LIBRARY INTLINS VOLMATS CANDU6S := DELETE:
+ BURNUP FLUX PIJ LIBRARY INTLINS VOLMATS CANDU6S ;
+*----
+* 2 group Burnup
+*----
+LIBRARY := EDITION :: STEP UP 'REF-CASE0001' ;
+EDITION := DELETE: EDITION ;
+PIJ := ASM: LIBRARY VOLMATF INTLINF ;
+FLUX := FLU: PIJ LIBRARY VOLMATF ::
+ TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.124906 ;
+EDITION := EDI: LIBRARY VOLMATF FLUX ::
+ EDIT 1 MERGE COMP MICR ALL SAVE ;
+EVALUATE Timec := 1.0 ;
+WHILE Timei Timec < DO
+ EVALUATE Timef := Timei Delt + ;
+ IF Timei 0.0 = THEN
+ BURNUP LIBRARY := EVO: LIBRARY FLUX VOLMATF ::
+ EDIT 3 DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ;
+ ELSE
+ BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX VOLMATF ::
+ EDIT 3 NOEX DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ;
+ ENDIF ;
+ PIJ := DELETE: PIJ ;
+ PIJ := ASM: LIBRARY VOLMATF INTLINF ;
+ FLUX := FLU: FLUX PIJ LIBRARY VOLMATF ::
+ TYPE K ;
+ EDITION := EDI: EDITION LIBRARY VOLMATF FLUX ::
+ SAVE ;
+*----
+* change delta t for burnup and final time if required
+*----
+ IF Timef Timec = THEN
+ IF Timec 5.0 = THEN
+ EVALUATE Delt Timec := 5.0 10.0 ;
+ ENDIF ;
+ IF Timec 1.0 = THEN
+ EVALUATE Delt Timec := 4.0 5.0 ;
+ ENDIF ;
+ ENDIF ;
+ EVALUATE Timei := Timef ;
+ENDWHILE ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.079448 ;
+res := EDITION ;
+EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F := DELETE:
+ EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F ;
+ECHO "test TCWE11 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE12.c2m b/Dragon/data/twimsE_proc/TCWE12.c2m
new file mode 100644
index 0000000..2ea0f32
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE12.c2m
@@ -0,0 +1,107 @@
+*----
+* TEST CASE TCWE12
+* WIMS-E Library
+* GENERATE A LIBRARY USING INFO AND OTHER OPTIONS
+*
+* REF: None
+*
+*----
+* Define variables and initialize
+* Coolant properties
+* a) Input
+* TempCool = Coolant temperature (K)
+* Purity = D2/(D2+H1) Weight % ratio in Coolant
+* b) Output
+* DensCool = Coolant Density (g/cm**3)
+* WH1C = H1 Weight % in Coolant
+* WD2C = D2 Weight % in Coolant
+* WO16C = O16 Weight % in Coolant
+* Fuel properties
+* a) Input
+* TempFuel = Fuel temperature (K)
+* Enrichment = U235/(U235+U238) Weight % ratio in Fuel
+* DensFuel = Fuel Density (g/cm**3)
+* b) Output
+* WU235F = U235 Weight % in Fuel
+* WU238F = U238 Weight % in Fuel
+* WO16F = O16 Weight % in Fuel
+*----
+REAL
+ TempCool Purity TempFuel Enrichment DensFuel :=
+ 560.66 99.95 941.29 0.72 10.437501 ;
+REAL
+ WH1C WD2C WO16C DensCool
+ WU235F WU238F WO16F ;
+*----
+* Define STRUCTURES and MODULES used
+*----
+LINKED_LIST
+ LIBRARY ISOT ;
+MODULE
+ LIB: INFO: END: ;
+PROCEDURE assertS ;
+*----
+* Get Coolant properties
+*----
+ECHO
+ "Input - Coolant temperature (K) " TempCool ;
+ECHO
+ "Input - D2/(D2+H1) Weight % ratio in Coolant" Purity ;
+INFO: ::
+ TMP: <<TempCool>> K
+ PUR: <<Purity>> WGT%
+ CALC DENS WATER >>DensCool<<
+ LIB: WIMSE FIL: WIMS-E
+ ISO: 3 '1001' '1002' '8016'
+ CALC WGT% D2O >>WH1C<< >>WD2C<< >>WO16C<<
+ ;
+ECHO
+ "Output - Coolant Density (g/cm**3)" DensCool ;
+ECHO
+ "Output - H1 Weight % in Coolant " WH1C ;
+ECHO
+ "Output - D2 Weight % in Coolant " WD2C ;
+ECHO
+ "Output - O16 Weight % in Coolant " WO16C ;
+*----
+* Get Fuel properties
+*----
+ECHO
+ "Input - Fuel temperature (K) " TempFuel ;
+ECHO
+ "Input - U235/(U235+U238) Weight % ratio in Fuel" Enrichment ;
+ECHO
+ "Input - Fuel Density (g/cm**3) " DensFuel ;
+INFO: ::
+ ENR: <<Enrichment>> WGT%
+ LIB: WIMSE FIL: WIMS-E
+ ISO: 3 '92235' '92238' '8016'
+ CALC WGT% UO2 >>WU235F<< >>WU238F<< >>WO16F<<
+ ;
+ECHO
+ "Output - U235 Weight % in Fuel" WU235F ;
+ECHO
+ "Output - U238 Weight % in Fuel" WU238F ;
+ECHO
+ "Output - O16 Weight % in Fuel " WO16F ;
+*----
+* Microscopic cross sections from file WIMS-E format WIMSE
+*----
+LIBRARY := LIB: ::
+ NMIX 5 CTRA WIMS
+ MIXS LIB: WIMSE FIL: WIMS-E
+ MIX 1 <<TempCool>> <<DensCool>> O16 = '8016' <<WO16C>>
+ D2D2O = '1002' <<WD2C>> H1H2O = '1001' <<WH1C>>
+ MIX 2 <<TempFuel>> <<DensCool>> O16 = '8016' <<WO16C>>
+ D2D2O = '1002' <<WD2C>> H1H2O = '1001' <<WH1C>>
+ MIX 3 <<TempFuel>> <<DensFuel>> O16 = '8016' <<WO16F>>
+ U235 = '92235' <<WU235F>> 1
+ U238 = '92238' <<WU238F>> 1
+ MIX 4 COMB 2 0.5 0 0.5
+ MIX 5 COMB 3 0.1 3 0.9
+ ;
+ISOT := LIBRARY :: STEP UP ISOTOPESLIST STEP AT 15 ; (*U238 0005*)
+assertS ISOT :: NG 50 4.918399E-01 ;
+ECHO "test TCWE12 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE13.c2m b/Dragon/data/twimsE_proc/TCWE13.c2m
new file mode 100644
index 0000000..21bf25c
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE13.c2m
@@ -0,0 +1,146 @@
+*----
+* TEST CASE TCWE13
+* 17 X 17 PWR ASSEMBLY WITHOUT POISON
+* WIMSE 69 GROUPS LIBRARY FILE WIMS-E
+*
+* REF: none
+*
+*----
+* Define STRUCTURES and MODULES used
+*----
+LINKED_LIST
+ ASSMB DISCR1 DISCR2 LIBRARY CP CALC OUT DATABASE ISOT SPHGEOM
+ MTRACK ;
+SEQ_BINARY FILTRK ;
+SEQ_ASCII
+ res ;
+MODULE
+ LIB: GEO: SYBILT: NXT: BIVACT: SHI: ASM: FLU: EDI: COMPO: SPH:
+ DELETE: END: ;
+PROCEDURE assertS ;
+*----
+* Microscopic cross sections from file WIMS-E format WIMSE
+*----
+LIBRARY := LIB: ::
+ NMIX 8 CTRA WIMS
+ MIXS LIB: WIMSE FIL: WIMS-E
+ MIX 1 579.9
+ H1H2O = '1001' 4.76690E-2 O16H2O = '8016' 2.38345E-2
+ B10 = '5010' 4.73825E-06 B11 = '5011' 1.90721E-05
+ MIX 2 579.9 O16 = '8016' 3.06711E-4
+ Cr52 = '24052' 7.54987E-5 Fe56 = '26056' 1.47624E-4
+ Zr90 = '40090' 2.15381E-02 Zr91 = '40091' 4.69693E-03 2
+ Zr92 = '40092' 7.17935E-03 Zr94 = '40094' 7.27563E-03
+ Zr96 = '40096' 1.17214E-03
+ MIX 3 579.9
+ H1H2O = '1001' 4.65292E-2 O16H2O = '8016' 2.32646E-2
+ Cr52 = '24052' 4.79927E-5
+ Fe56 = '26056' 4.45845E-5 Ni58 = '28058' 1.13521E-4
+ Mo92 = '42092' 5.91461E-07 Mo95 = '42095' 3.70930E-07
+ Mo96 = '42096' 6.40880E-07 Mo97 = '42097' 3.86878E-07
+ Mo98 = '42098' 9.80802E-07 Mo100 = '42100' 3.93419E-07
+ Al27 = '13027' 2.35231E-6
+ Mn55 = '25055' 4.15901E-7
+ B10 = '5010' 4.63194E-06 B11 = '5011' 1.86442E-05
+ Zr90 = '40090' 4.59154E-04 Zr91 = '40091' 1.00130E-04
+ Zr92 = '40092' 1.53051E-04 Zr94 = '40094' 1.55104E-04
+ Zr96 = '40096' 2.49880E-05
+ MIX 4 579.9 O16 = '8016' 2.87335E-4
+ Cr52 = '24052' 7.07291E-5 Fe56 = '26056' 1.38298E-4
+ Zr90 = '40090' 2.01774E-02 Zr91 = '40091' 4.40020E-03 2
+ Zr92 = '40092' 6.72580E-03 Zr94 = '40094' 6.81600E-03
+ Zr96 = '40096' 1.09809E-03
+ MIX 5 579.9
+ H1H2O = '1001' 4.71346E-2 O16H2O = '8016' 2.35673E-2
+ Cr52 = '24052' 2.24991E-5
+ Fe56 = '26056' 2.09013E-5 Ni58 = '28058' 5.32188E-5
+ Mo92 = '42092' 2.77278E-07 Mo95 = '42095' 1.73892E-07
+ Mo96 = '42096' 3.00446E-07 Mo97 = '42097' 1.81369E-07
+ Mo98 = '42098' 4.59801E-07 Mo100 = '42100' 1.84435E-07
+ Al27 = '13027' 1.10277E-6
+ Mn55 = '25055' 1.94976E-7
+ B10 = '5010' 4.68840E-06 B11 = '5011' 1.88714E-05
+ Zr90 = '40090' 2.15252E-04 Zr91 = '40091' 4.69413E-05
+ Zr92 = '40092' 7.17508E-05 Zr94 = '40094' 7.27131E-05
+ Zr96 = '40096' 1.17144E-05
+ MIX 6 579.9
+ H1H2O = '1001' 4.71676E-2 O16H2O = '8016' 2.35838E-2
+ Cr52 = '24052' 2.11122E-5
+ Fe56 = '26056' 1.96130E-5 Ni58 = '28058' 4.99383E-5
+ Mo92 = '42092' 2.60187E-07 Mo95 = '42095' 1.63174E-07
+ Mo96 = '42096' 2.81927E-07 Mo97 = '42097' 1.70190E-07
+ Mo98 = '42098' 4.31460E-07 Mo100 = '42100' 1.73067E-07
+ Al27 = '13027' 1.03479E-6
+ Mn55 = '25055' 1.82957E-7
+ B10 = '5010' 4.69148E-06 B11 = '5011' 1.88838E-05
+ Zr90 = '40090' 2.01984E-04 Zr91 = '40091' 4.40478E-05
+ Zr92 = '40092' 6.73280E-05 Zr94 = '40094' 6.82309E-05
+ Zr96 = '40096' 1.09923E-05
+ MIX 7 579.9
+ H1H2O = '1001' 4.72020E-2 O16H2O = '8016' 2.36010E-2
+ Cr52 = '24052' 1.96591E-5
+ Fe56 = '26056' 1.82630E-5 Ni58 = '28058' 4.65011E-5
+ Mo92 = '42092' 2.42278E-07 Mo95 = '42095' 1.51943E-07
+ Mo96 = '42096' 2.62522E-07 Mo97 = '42097' 1.58476E-07
+ Mo98 = '42098' 4.01763E-07 Mo100 = '42100' 1.61155E-07
+ Al27 = '13027' 9.63569E-7
+ Mn55 = '25055' 1.70365E-7
+ B10 = '5010' 4.69469E-06 B11 = '5011' 1.88967E-05
+ Zr90 = '40090' 1.88082E-04 Zr91 = '40091' 4.10161E-05 2
+ Zr92 = '40092' 6.26939E-05 Zr94 = '40094' 6.35347E-05
+ Zr96 = '40096' 1.02357E-05
+ MIX 8 933.6 O16 = '8016' 4.49355E-2
+ U235 = '92235' 7.39237E-4 1
+ U238 = '92238' 2.17285E-2 1
+ ;
+*----
+* Geometry ASSMB : a 17 X 17 normal PWR assembly
+* contains C1 : cell without fuel
+* C2 : normal fuel cell
+* C3 : peripheral cell
+* C4 : corner cell
+*----
+ASSMB := GEO: :: CAR2D 9 9
+ X- DIAG X+ REFL Y- SYME Y+ DIAG
+ CELL C1 C2 C2 C1 C2 C2 C1 C2 C3
+ C2 C2 C2 C2 C2 C2 C2 C3
+ C2 C2 C2 C2 C2 C2 C3
+ C1 C2 C2 C1 C2 C3
+ C2 C2 C2 C2 C3
+ C1 C2 C2 C3
+ C2 C2 C3
+ C2 C3
+ C4
+ ::: C1 := GEO: CARCEL 2
+ MESHX 0.0 1.26472 MESHY 0.0 1.26472
+ RADIUS 0.0 0.572435 0.613142 MIX 1 2 3 ;
+ ::: C2 := GEO: C1 RADIUS 0.0 0.412660 0.474364 MIX 8 4 5 ;
+ ::: C3 := GEO: C2 MESHX 0.0 1.31472 MIX 8 4 6 ;
+ ::: C4 := GEO: C3 MESHY 0.0 1.31472 MIX 8 4 7 ;
+ ;
+*----
+* Self-Shielding calculation SYBIL
+* Transport calculation SYBIL
+* Flux calculation for B1 homogeneous leakage
+* Editing using SPH model for transport-diffusion
+*----
+DISCR1 := SYBILT: ASSMB ::
+ TITLE 'TCWE13: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON'
+ MAXR 400 QUA2 6 3 ;
+DISCR2 FILTRK := NXT: ASSMB ::
+ TITLE 'TCWE13: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON'
+ ALLG BATCH 100
+ TISO 10 20.0 ;
+LIBRARY := SHI: LIBRARY DISCR1 :: EDIT 0 NOLJ ;
+
+CP := ASM: LIBRARY DISCR2 FILTRK ;
+CALC := FLU: CP LIBRARY DISCR2 ::
+ TYPE B B1 PNL ;
+assertS CALC :: 'K-INFINITY' 1 1.260388 ;
+
+OUT := EDI: LIBRARY DISCR1 CALC ASSMB ::
+ EDIT 3 UPS SAVE MICR RES MERGE CELL COND 4.0
+ ;
+ECHO "test TCWE13 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE14.c2m b/Dragon/data/twimsE_proc/TCWE14.c2m
new file mode 100644
index 0000000..8f626ec
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE14.c2m
@@ -0,0 +1,121 @@
+*----
+* TEST CASE TCWE14
+* CANDU-6 ANNULAR CELL
+* WIMS-E Library
+*----
+*----
+* Define STRUCTURES and MODULES used
+*----
+LINKED_LIST
+ LIBRARY CANDU6S CANDU6F VOLMATS VOLMATF PIJ FLUX EDITION
+ DATABASE ISOT ;
+SEQ_BINARY
+ INTLINS INTLINF ;
+SEQ_ASCII
+ database ;
+MODULE
+ GEO: NXT: SHI: ASM: FLU: EDI: DELETE: END: ;
+PROCEDURE assertS ;
+*----
+* Depletion data from file WIMS-E format WIMSE
+* Microscopic cross sections from file WIMS-E format WIMSE
+*----
+PROCEDURE TCWE05Lib ;
+INTEGER iedit := 1 ;
+LIBRARY := TCWE05Lib :: <<iedit>> ;
+*----
+* Geometry CANDU6S : 13 regions annular cluster for self-shielding
+* CANDU6F : 31 regions annular cluster for transport
+*----
+CANDU6S := GEO: :: CARCEL 5
+ X+ REFL X- REFL
+ Y+ REFL Y- REFL
+ MESHX -14.2875 14.2875
+ MESHY -14.2875 14.2875
+ RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 14.00
+ MIX 1 2 3 4 5 5
+ CLUSTER ROD1 ROD2 ROD3 ROD4
+ ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000
+ RADIUS 0.00000 0.6122 0.6540 ;
+ ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ;
+ ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ;
+ ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ;
+;
+CANDU6F := GEO: :: CAR2D 2 2
+ X+ REFL X- REFL
+ Y+ REFL Y- REFL
+ CELL F1 F2 F3 F4
+ ::: F1 := GEO: CARCEL 5
+ MESHX -14.2875 14.2875
+ MESHY -14.2875 14.2875
+ RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 14.00
+ MIX 1 2 3 4 5 5
+ CLUSTER ROD1 ROD2 ROD3 ROD4
+ ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000
+ RADIUS 0.00000 0.6122 0.6540 ;
+ ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ;
+ ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ;
+ ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ;
+ ;
+ ::: F2 := GEO: CARCEL 5
+ MESHX -14.2875 14.2875
+ MESHY -14.2875 14.2875
+ RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 14.00
+ MIX 1 2 3 4 5 5
+ CLUSTER ROD1 ROD2 ROD3 ROD4
+ ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000
+ RADIUS 0.00000 0.6122 0.6540 ;
+ ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ;
+ ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ;
+ ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ;
+ ;
+ ::: F3 := GEO: CARCEL 5
+ MESHX -14.2875 14.2875
+ MESHY -14.2875 14.2875
+ RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 14.00
+ MIX 1 2 3 4 5 5
+ CLUSTER ROD1 ROD2 ROD3 ROD4
+ ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000
+ RADIUS 0.00000 0.6122 0.6540 ;
+ ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ;
+ ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ;
+ ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ;
+ ;
+ ::: F4 := GEO: CARCEL 5
+ MESHX -14.2875 14.2875
+ MESHY -14.2875 14.2875
+ RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 14.00
+ MIX 1 2 3 4 5 5
+ CLUSTER ROD1 ROD2 ROD3 ROD4
+ ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000
+ RADIUS 0.00000 0.6122 0.6540 ;
+ ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ;
+ ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ;
+ ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ;
+ ;
+;
+*----
+* Self-Shielding calculation NXT
+* Transport calculation NXT
+* Flux calculation for keff
+*----
+VOLMATS INTLINS := NXT: CANDU6S ::
+ TITLE 'TCWE14: CANDU-6 ANNULAR POWER= 31.971 FUEL TEMP= 941.29'
+ ALLG BATCH 100 MERGMIX
+ EDIT 0 TISO 50 10.0 ;
+LIBRARY := SHI: LIBRARY VOLMATS INTLINS ::
+ EDIT 0 ;
+VOLMATF INTLINF := NXT: CANDU6F ::
+ TITLE 'TCWE14: CANDU-6 ANNULAR POWER= 31.971 FUEL TEMP= 941.29'
+ ALLG BATCH 100
+ EDIT 0 TISO 50 10.0 ;
+PIJ := ASM: LIBRARY VOLMATF INTLINF ;
+FLUX := FLU: PIJ LIBRARY VOLMATF ::
+ TYPE K ;
+assertS FLUX :: 'K-INFINITY' 1 1.133684 ;
+EDITION := EDI: LIBRARY VOLMATF FLUX ::
+ COND 4.0 MERGE CELL MICR 1 Xe135 SAVE ON 'fuel' ;
+INTLINF INTLINS := DELETE: INTLINF INTLINS ;
+ECHO "test TCWE14 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE15.c2m b/Dragon/data/twimsE_proc/TCWE15.c2m
new file mode 100644
index 0000000..a2e54f4
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE15.c2m
@@ -0,0 +1,159 @@
+*----
+* TEST CASE TCWE05
+* CANDU-6 ANNULAR CELL
+* WIMS-E Library
+* POWER (KW) = 615.00000
+* BURN POWER (KW/KG) = 31.97130
+* URANIUM MASS = 19.23600
+* UO2 REAL DENSITY = 10.59300
+* UO2 EFF DENSITY = 10.43750
+* UO2 TEMPERATURE = 941.28998
+* ENRICHMENT = 0.71140
+* COOLANT D2 AT % = 99.222
+* MODERATOR D2 AT % = 99.911
+* NUMBER OF DAYS = 50
+*
+*----
+* Define variables and initialize
+* Burnup paremeters
+* a) Power
+* = 31.9713 kw/kg for 0.0 to 300.0 days
+* b) Burnup time interval Delt
+* = 1 day for 0 to 1 day
+* = 4 days for 1 to 5 days
+* = 5 days for 5 to 10 days
+* = 10 days for 10 to 50 days
+* = 20 days for 50 to 150 days
+* = 50 days for 150 to 300 days
+* c) Days with burnup interval changes
+* = 1.0, 5.0, 10.0, 50.0, 150.0 and 300.0 days
+* d) Burnup control time variables Timei, Timef
+* Timei = initial time
+* Timef = final time
+*----
+REAL
+ Power Delt Timec Timei Timef :=
+ 31.9713 1.0 1.0 0.0 0.0 ;
+*----
+* Define STRUCTURES and MODULES used
+*----
+LINKED_LIST
+ LIBRARY CANDU6S CANDU6F VOLMATS VOLMATF PIJ FLUX BURNUP EDITION
+ COMPO1 COMPO2 ;
+SEQ_BINARY
+ INTLINS INTLINF ;
+SEQ_ASCII
+ fuel mode ;
+MODULE
+ GEO: NXT: SHI: ASM: FLU: EVO: EDI: CPO:
+ DELETE: END: ;
+*----
+* Depletion data from file WIMS-E format WIMSE
+* Microscopic cross sections from file WIMS-E format WIMSE
+*----
+PROCEDURE TCWE05Lib ;
+PROCEDURE assertS ;
+INTEGER iedit := 1 ;
+LIBRARY := TCWE05Lib :: <<iedit>> ;
+*----
+* Geometry CANDU6S : 13 regions annular cluster for self-shielding
+* CANDU6F : 31 regions annular cluster for transport
+*----
+CANDU6S := GEO: :: CARCEL 5
+ X+ REFL X- REFL MESHX -14.2875 14.2875
+ Y+ REFL Y- REFL MESHY -14.2875 14.2875
+ RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 14.00
+ MIX 1 2 3 4 5 5
+ CLUSTER ROD1 ROD2 ROD3 ROD4
+ ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000
+ RADIUS 0.00000 0.6122 0.6540 ;
+ ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ;
+ ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ;
+ ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ;
+ ;
+CANDU6F := GEO: CANDU6S :: SPLITR 6 1 1 1 10
+ ::: ROD1 := GEO: ROD1 SPLITR 2 1 ;
+ ::: ROD2 := GEO: ROD2 SPLITR 2 1 ;
+ ::: ROD3 := GEO: ROD3 SPLITR 2 1 ;
+ ::: ROD4 := GEO: ROD4 SPLITR 2 1 ;
+ ;
+*----
+* Self-Shielding calculation EXCEL
+* Transport calculation EXCEL
+* Flux calculation for keff
+*----
+VOLMATS INTLINS := NXT: CANDU6S ::
+ TITLE 'TCWE05: CANDU-6 CARTESIAN POWER= 31.971 FUEL TEMP= 941.29'
+ ALLG BATCH 100 MERGMIX
+ EDIT 0 TRAK TISO 5 10.0 ;
+LIBRARY := SHI: LIBRARY VOLMATS INTLINS ::
+ EDIT 0 ;
+VOLMATF INTLINF := NXT: CANDU6F ::
+ TITLE 'TCWE05: CANDU-6 CARTESIAN POWER= 31.971 FUEL TEMP= 941.29'
+ ALLG BATCH 100
+ EDIT 0 TRAK TISO 5 10.0 ;
+PIJ := ASM: LIBRARY VOLMATF INTLINF :: ;
+FLUX := FLU: PIJ LIBRARY VOLMATF ::
+ TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.124479 ;
+EDITION := EDI: FLUX LIBRARY VOLMATF ::
+ COND 4.0 MERGE MIX 0 0 0 0 1 0 0 0 0 0 SAVE ON 'EDITMOD' ;
+EDITION := EDI: EDITION FLUX LIBRARY VOLMATF ::
+ COND 4.0 MERGE COMP MICR 1 Xe135 SAVE ;
+*----
+* Burnup loop: for first step BURNUP is created
+* while for other steps it is modified
+*----
+WHILE Timei Timec < DO
+ EVALUATE Timef := Timei Delt + ;
+ IF Timei 0.0 = THEN
+ BURNUP LIBRARY := EVO: LIBRARY FLUX VOLMATF ::
+ DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ;
+ ELSE
+ BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX VOLMATF ::
+ NOEX DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ;
+ ENDIF ;
+ LIBRARY := SHI: LIBRARY VOLMATS INTLINS ::
+ EDIT 0 ;
+ PIJ := DELETE: PIJ ;
+ PIJ := ASM: LIBRARY VOLMATF INTLINF :: ;
+ FLUX := FLU: FLUX PIJ LIBRARY VOLMATF ::
+ TYPE K ;
+ EDITION := EDI: EDITION FLUX LIBRARY VOLMATF ::
+ SAVE ;
+*----
+* change delta t for burnup and final time if required
+*----
+ IF Timef Timec = THEN
+ IF Timec 150.0 = THEN
+ EVALUATE Delt Timec := 50.0 300.0 ;
+ ENDIF ;
+ IF Timec 50.0 = THEN
+ EVALUATE Delt Timec := 20.0 150.0 ;
+ ENDIF ;
+ IF Timec 10.0 = THEN
+ EVALUATE Delt Timec := 10.0 50.0 ;
+ ENDIF ;
+ IF Timec 5.0 = THEN
+ EVALUATE Delt Timec := 5.0 10.0 ;
+ ENDIF ;
+ IF Timec 1.0 = THEN
+ EVALUATE Delt Timec := 4.0 5.0 ;
+ ENDIF ;
+ ENDIF ;
+ EVALUATE Timei := Timef ;
+ENDWHILE ;
+assertS FLUX :: 'K-EFFECTIVE' 1 0.9688265 ;
+*----
+* Save calculation results in CPO format file
+*----
+COMPO1 := CPO: BURNUP EDITION ::
+ BURNUP REF-CASE EXTRACT Xe135 Xe135 NAME MIXTRXE ;
+fuel := COMPO1 ;
+COMPO2 := CPO: EDITION ::
+ STEP 'EDITMOD' NAME MIXTMOD ;
+mode := COMPO2 ;
+INTLINF INTLINS := DELETE: INTLINF INTLINS ;
+ECHO "test TCWE15 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE17.c2m b/Dragon/data/twimsE_proc/TCWE17.c2m
new file mode 100644
index 0000000..1d0d955
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE17.c2m
@@ -0,0 +1,337 @@
+*----
+* Exemple of the use of HMIX for cell homogenization
+* 2-D supercell with fuel clusters based on AECL supercell model
+* for G2 SOR and MCA with fuel and reactivity devices parallel
+* References
+* PREPARED BY : G. Marleau on 2013/06/11
+*
+*----
+* modules and data structures
+*----
+SEQ_ASCII MACROLIBF ;
+LINKED_LIST SORINS SORIN TRACK MicLib FLUX EDITION ;
+XSM_FILE ASMPIJ ;
+SEQ_ASCII FigReg.ps FigMix.ps FigHom.ps ;
+SEQ_ASCII HomMix.txt HomHMix.txt ;
+SEQ_BINARY Lines ;
+MODULE GEO: EXCELT: EXCELL: SHI: ASM: LIB: FLU:
+ EDI: UTL: DELETE: FREE: END: INFO: CPO: MAC:
+ NXT: PSP: ;
+PROCEDURE TCWE17Lib ;
+PROCEDURE assertS ;
+INTEGER iedit := 1 ;
+MicLib := TCWE17Lib :: <<iedit>> ;
+*----
+* DEFINE GEOMETRY FOR SUPERCELL CALCULATION
+* SORINS : 2D self-shielding geometry with SHUT-OFF ROD & GT in
+* for annular fuel.
+* SORIN : 2D transport geometry with SHUT-OFF ROD & GT in
+* for annular fuel.
+*----
+SORINS := GEO: :: CAR2D 5 3
+ EDIT 0
+ X- REFL X+ REFL
+ Y- REFL Y+ REFL
+ CELL M MXL MX2 MXR M
+ MY FXYL BXY FXYR MY
+ M MXL MX2 MXR M
+ ::: M := GEO: CAR2D 1 1
+ MESHX 0.0 7.14375
+ MESHY 0.0 7.14375
+ HMIX 0
+ MIX 5 ;
+ ::: MXL := GEO: CAR2D 2 1
+ MESHX -7.14375 0.0 7.14375
+ MESHY 0.0 7.14375
+ HMIX 0 1
+ MIX 5 15 ;
+ ::: MX2 := GEO: CAR2D 1 1
+ MESHX -7.14375 7.14375
+ MESHY 0.0 7.14375
+ HMIX 1
+ MIX 15 ;
+ ::: MXR := GEO: CAR2D 2 1
+ MESHX -7.14375 0.0 7.14375
+ MESHY 0.0 7.14375
+ HMIX 1 0
+ MIX 15 5 ;
+ ::: MY := GEO: CAR2D 1 2
+ MESHY -7.14375 0.0 7.14375
+ MESHX 0.0 7.14375
+ HMIX 0 0
+ MIX 5 5 ;
+ ::: BXY := GEO: CARCEL 2
+ MESHX -7.14375 7.14375
+ MESHY -7.14375 7.14375
+ RADIUS 0.0 6.380 6.530
+ HMIX 1 1 1
+ MIX 15 15 15
+ CLUSTER ROD
+ ::: ROD := GEO: TUBE 4
+ NPIN 1 RPIN 0.0 APIN 0.0
+ RADIUS 0.0 5.4115 5.4877 5.5791 5.6553
+ HMIX 1 1 1 1
+ MIX 15 14 21 14 ; ;
+ ::: FXYL := GEO: CARCEL 5 2 1
+ MESHX -7.14375 0.0 7.14375
+ MESHY -7.14375 7.14375
+ RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 7.00
+ HMIX 0 0 0 0 0 0
+ 1 1 1 1 1 1
+ MIX 1 2 3 4 5 5
+ 11 12 13 14 15 15
+ CLUSTER ROD1 ROD2L ROD2R ROD3L ROD3R ROD4L ROD4R
+ ::: ROD1 := GEO: TUBE 2 1 2 MIX 6 10 16 20 HMIX 0 0 1 1
+ NPIN 1 RPIN 0.0000 APIN 0.0000
+ MESHX -0.6540 0.6540
+ MESHY -0.6540 0.0 0.6540
+ RADIUS 0.00000 0.6122 0.6540 ;
+ ::: ROD2L := GEO: TUBE 2 MIX 7 10 HMIX 0 0
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 3 RPIN 1.4885
+ APIN 2.09439510 3.14159265 4.18879020 ;
+ ::: ROD2R := GEO: TUBE 2 MIX 17 20 HMIX 1 1
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 3 RPIN 1.4885
+ APIN -1.04719755 0.0000 1.04719755 ;
+ ::: ROD3L := GEO: TUBE 2 MIX 8 10 HMIX 0 0
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 6 RPIN 2.8755
+ APIN -2.87979327 -2.35619449 -1.83259571
+ 1.83259571 2.35619449 2.87979327 ;
+ ::: ROD3R := GEO: TUBE 2 MIX 18 20 HMIX 1 1
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 6 RPIN 2.8755
+ APIN -1.30899694 -0.78539816 -0.26179939
+ 0.26179939 0.78539816 1.30899694 ;
+ ::: ROD4L := GEO: TUBE 2 MIX 9 10 HMIX 0 0
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 9 RPIN 4.3305
+ APIN 1.74532925 2.09439510 2.44346095
+ 2.79252680 3.14159265 3.49065850
+ 3.83972435 4.18879020 4.53785606 ;
+ ::: ROD4R := GEO: TUBE 2 MIX 19 20 HMIX 1 1
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 9 RPIN 4.3305
+ APIN -1.39626340 -1.04719755 -0.69813170
+ -0.34906585 0.0 0.34906585
+ 0.69813170 1.04719755 1.39626340 ;
+ ;
+ ::: FXYR := GEO: CARCEL 5 2 1
+ MESHX -7.14375 0.0 7.14375
+ MESHY -7.14375 7.14375
+ RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 7.00
+ HMIX 1 1 1 1 1 1
+ 0 0 0 0 0 0
+ MIX 11 12 13 14 15 15
+ 1 2 3 4 5 5
+ CLUSTER ROD1 ROD2L ROD2R ROD3L ROD3R ROD4L ROD4R
+ ::: ROD1 := GEO: TUBE 2 1 2 MIX 16 20 6 10 HMIX 1 1 0 0
+ NPIN 1 RPIN 0.0000 APIN 0.0000
+ MESHX -0.6540 0.6540
+ MESHY -0.6540 0.0 0.6540
+ RADIUS 0.00000 0.6122 0.6540 ;
+ ::: ROD2L := GEO: TUBE 2 MIX 17 20 HMIX 1 1
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 3 RPIN 1.4885
+ APIN 2.09439510 3.14159265 4.18879020 ;
+ ::: ROD2R := GEO: TUBE 2 MIX 7 10 HMIX 0 0
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 3 RPIN 1.4885
+ APIN -1.04719755 0.0000 1.04719755 ;
+ ::: ROD3L := GEO: TUBE 2 MIX 18 20 HMIX 1 1
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 6 RPIN 2.8755
+ APIN -2.87979327 -2.35619449 -1.83259571
+ 1.83259571 2.35619449 2.87979327 ;
+ ::: ROD3R := GEO: TUBE 2 MIX 8 10 HMIX 0 0
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 6 RPIN 2.8755
+ APIN -1.30899694 -0.78539816 -0.26179939
+ 0.26179939 0.78539816 1.30899694 ;
+ ::: ROD4L := GEO: TUBE 2 MIX 19 20 HMIX 1 1
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 9 RPIN 4.3305
+ APIN 1.74532925 2.09439510 2.44346095
+ 2.79252680 3.14159265 3.49065850
+ 3.83972435 4.18879020 4.53785606 ;
+ ::: ROD4R := GEO: TUBE 2 MIX 9 10 HMIX 0 0
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 9 RPIN 4.3305
+ APIN -1.39626340 -1.04719755 -0.69813170
+ -0.34906585 0.0 0.34906585
+ 0.69813170 1.04719755 1.39626340 ;
+ ;
+;
+SORIN := GEO: :: CAR2D 5 3
+ EDIT 0
+ X- REFL X+ REFL
+ Y- REFL Y+ REFL
+ CELL M MXL MX2 MXR M
+ MY FXYL BXY FXYR MY
+ M MXL MX2 MXR M
+ ::: M := GEO: CAR2D 1 1
+ MESHX 0.0 7.14375 SPLITX 2
+ MESHY 0.0 7.14375 SPLITY 2
+ HMIX 0
+ MIX 5 ;
+ ::: MXL := GEO: CAR2D 2 1
+ MESHX -7.14375 0.0 7.14375 SPLITX 3 3
+ MESHY 0.0 7.14375 SPLITY 3
+ HMIX 0 1
+ MIX 5 15 ;
+ ::: MX2 := GEO: CAR2D 1 1
+ MESHX -7.14375 7.14375 SPLITX 6
+ MESHY 0.0 7.14375 SPLITY 3
+ HMIX 1
+ MIX 15 ;
+ ::: MXR := GEO: CAR2D 2 1
+ MESHX -7.14375 0.0 7.14375 SPLITX 3 3
+ MESHY 0.0 7.14375 SPLITY 3
+ HMIX 1 0
+ MIX 15 5 ;
+ ::: MY := GEO: CAR2D 1 2
+ MESHY -7.14375 0.0 7.14375 SPLITY 3 3
+ MESHX 0.0 7.14375 SPLITX 3
+ HMIX 0 0
+ MIX 5 5 ;
+ ::: BXY := GEO: CARCEL 2
+ MESHX -7.14375 7.14375 SPLITX 6
+ MESHY -7.14375 7.14375 SPLITY 6
+ RADIUS 0.0 6.380 6.530
+ SPLITR 2 2
+ HMIX 1 1 1
+ MIX 15 15 15
+ CLUSTER ROD
+ ::: ROD := GEO: TUBE 4
+ NPIN 1 RPIN 0.0 APIN 0.0
+ RADIUS 0.0 5.4115 5.4877 5.5791 5.6553
+ SPLITR 1 1 1 1
+ HMIX 1 1 1 1
+ MIX 15 14 21 14 ; ;
+ ::: FXYL := GEO: CARCEL 5 2 1
+ MESHX -7.14375 0.0 7.14375 SPLITX 3 3
+ MESHY -7.14375 7.14375 SPLITY 6
+ RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 7.00
+ HMIX 0 0 0 0 0 0
+ 1 1 1 1 1 1
+ MIX 1 2 3 4 5 5
+ 11 12 13 14 15 15
+ CLUSTER ROD1 ROD2L ROD2R ROD3L ROD3R ROD4L ROD4R
+ ::: ROD1 := GEO: TUBE 2 1 2 MIX 6 10 16 20 HMIX 0 0 1 1
+ NPIN 1 RPIN 0.0000 APIN 0.0000
+ MESHX -0.6540 0.6540
+ MESHY -0.6540 0.0 0.6540
+ RADIUS 0.00000 0.6122 0.6540 ;
+ ::: ROD2L := GEO: TUBE 2 MIX 7 10 HMIX 0 0
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 3 RPIN 1.4885
+ APIN 2.09439510 3.14159265 4.18879020 ;
+ ::: ROD2R := GEO: TUBE 2 MIX 17 20 HMIX 1 1
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 3 RPIN 1.4885
+ APIN -1.04719755 0.0000 1.04719755 ;
+ ::: ROD3L := GEO: TUBE 2 MIX 8 10 HMIX 0 0
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 6 RPIN 2.8755
+ APIN -2.87979327 -2.35619449 -1.83259571
+ 1.83259571 2.35619449 2.87979327 ;
+ ::: ROD3R := GEO: TUBE 2 MIX 18 20 HMIX 1 1
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 6 RPIN 2.8755
+ APIN -1.30899694 -0.78539816 -0.26179939
+ 0.26179939 0.78539816 1.30899694 ;
+ ::: ROD4L := GEO: TUBE 2 MIX 9 10 HMIX 0 0
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 9 RPIN 4.3305
+ APIN 1.74532925 2.09439510 2.44346095
+ 2.79252680 3.14159265 3.49065850
+ 3.83972435 4.18879020 4.53785606 ;
+ ::: ROD4R := GEO: TUBE 2 MIX 19 20 HMIX 1 1
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 9 RPIN 4.3305
+ APIN -1.39626340 -1.04719755 -0.69813170
+ -0.34906585 0.0 0.34906585
+ 0.69813170 1.04719755 1.39626340 ;
+ ;
+ ::: FXYR := GEO: CARCEL 5 2 1
+ MESHX -7.14375 0.0 7.14375 SPLITX 3 3
+ MESHY -7.14375 7.14375 SPLITY 6
+ RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 7.00
+ HMIX 1 1 1 1 1 1
+ 0 0 0 0 0 0
+ MIX 11 12 13 14 15 15
+ 1 2 3 4 5 5
+ CLUSTER ROD1 ROD2L ROD2R ROD3L ROD3R ROD4L ROD4R
+ ::: ROD1 := GEO: TUBE 2 1 2 MIX 16 20 6 10 HMIX 1 1 0 0
+ NPIN 1 RPIN 0.0000 APIN 0.0000
+ MESHX -0.6540 0.6540
+ MESHY -0.6540 0.0 0.6540
+ RADIUS 0.00000 0.6122 0.6540 ;
+ ::: ROD2L := GEO: TUBE 2 MIX 17 20 HMIX 1 1
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 3 RPIN 1.4885
+ APIN 2.09439510 3.14159265 4.18879020 ;
+ ::: ROD2R := GEO: TUBE 2 MIX 7 10 HMIX 0 0
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 3 RPIN 1.4885
+ APIN -1.04719755 0.0000 1.04719755 ;
+ ::: ROD3L := GEO: TUBE 2 MIX 18 20 HMIX 1 1
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 6 RPIN 2.8755
+ APIN -2.87979327 -2.35619449 -1.83259571
+ 1.83259571 2.35619449 2.87979327 ;
+ ::: ROD3R := GEO: TUBE 2 MIX 8 10 HMIX 0 0
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 6 RPIN 2.8755
+ APIN -1.30899694 -0.78539816 -0.26179939
+ 0.26179939 0.78539816 1.30899694 ;
+ ::: ROD4L := GEO: TUBE 2 MIX 19 20 HMIX 1 1
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 9 RPIN 4.3305
+ APIN 1.74532925 2.09439510 2.44346095
+ 2.79252680 3.14159265 3.49065850
+ 3.83972435 4.18879020 4.53785606 ;
+ ::: ROD4R := GEO: TUBE 2 MIX 9 10 HMIX 0 0
+ RADIUS 0.00000 0.6122 0.6540
+ NPIN 9 RPIN 4.3305
+ APIN -1.39626340 -1.04719755 -0.69813170
+ -0.34906585 0.0 0.34906585
+ 0.69813170 1.04719755 1.39626340 ;
+ ;
+;
+*----
+* Rod and GT absent
+*----
+TRACK Lines := NXT: SORINS :: EDIT 3
+ ALLG BATCH 100 MERGMIX
+ TISO 40 30.0 ;
+MicLib := SHI: MicLib TRACK Lines ;
+TRACK Lines := DELETE: TRACK Lines ;
+TRACK Lines := NXT: SORIN :: EDIT 3
+ ALLG BATCH 100
+ TISO 40 30.0 ;
+FigReg.ps := PSP: TRACK :: TYPE REGI ;
+FigMix.ps := PSP: TRACK :: TYPE MIXT ;
+FigHom.ps := PSP: TRACK :: TYPE HMIX ;
+ASMPIJ := ASM: MicLib TRACK Lines :: PIJ ;
+FLUX := FLU: ASMPIJ MicLib TRACK :: TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.039269 ;
+EDITION := EDI: FLUX MicLib TRACK ::
+ EDIT 3
+ COND 0.625
+ MERG HMIX
+ SAVE ON SORINHMIX ;
+HomHMix.txt := EDITION ;
+EDITION := DELETE: EDITION ;
+EDITION := EDI: FLUX MicLib TRACK ::
+ EDIT 3
+ COND 0.625
+ MERG MIX 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
+ SAVE ON SORINMIX ;
+HomMix.txt := EDITION ;
+EDITION := DELETE: EDITION ;
+TRACK FLUX ASMPIJ Lines := DELETE: TRACK FLUX ASMPIJ Lines ;
+ECHO "test TCWE17 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE17Lib.c2m b/Dragon/data/twimsE_proc/TCWE17Lib.c2m
new file mode 100644
index 0000000..e085080
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE17Lib.c2m
@@ -0,0 +1,97 @@
+*----
+* Create Library for test CASE TCWE17.x2m
+* Calling :
+* LIBRARY := TCWE17Lib.x2m :: iedit ;
+* with :
+* LIBRARY = Linked list containing the result of LIB: for
+* TCWE17.x2m
+* iprint = print level for LIB: module
+*----
+* Define PARAMETERS,STRUCTURES and MODULES used
+*----
+PARAMETER LIBRARY :: ::: LINKED_LIST LIBRARY ; ;
+MODULE LIB: DELETE: END: ;
+*----
+* Define and read LIB: EDIT option
+INTEGER iedit ;
+:: >>iedit<< ;
+*----
+* Depletion data from file WIMS-E format WIMSE
+* Microscopic cross sections from file WIMS-E format WIMSE
+* All materials are duplicated for left and right cell identification
+* during homogenization
+*----
+LIBRARY := LIB: ::
+ EDIT <<iedit>>
+ NMIX 21 CTRA WIMS
+ DEPL LIB: WIMSE FIL: WIMS-E
+ MIXS LIB: WIMSE FIL: WIMS-E
+ MIX 1 560.66 0.81212 O16 = '8016' 7.99449E-1
+ D2D2O = '1002' 1.99768E-1 H1H2O = '1001' 7.83774E-4
+ MIX 2 560.66 6.57 Nb93 = '41093' 2.50000
+ B10 = '5010' 4.17900E-05 B11 = '5011' 1.68210E-04
+ Zr90 = '40090' 5.01638E+01 Zr91 = '40091' 1.09395E+01 2
+ Zr92 = '40092' 1.67213E+01 Zr94 = '40094' 1.69455E+01
+ Zr96 = '40096' 2.73000E+00
+ MIX 3 345.66 0.0014 He4 = '2004' 1.00000E+2
+ MIX 4 345.66 6.44 Fe56 = '26056' 1.60000E-1
+ Ni58 = '28058' 6.00000E-2 Cr52 = '24052' 1.10000E-1
+ B10 = '5010' 6.16900E-05 B11 = '5011' 2.48310E-04
+ Zr90 = '40090' 5.13008E+01 Zr91 = '40091' 1.11875E+01 2
+ Zr92 = '40092' 1.71003E+01 Zr94 = '40094' 1.73296E+01
+ Zr96 = '40096' 2.79188E+00
+ MIX 5 345.66 1.082885 O16 = '8016' 7.98895E-1
+ D2D2O = '1002' 2.01016E-1 H1H2O = '1001' 8.96000E-5
+ MIX 6 941.29 10.4375010 O16 = '8016' 1.18473E+1
+ Xe135 = '54135' 0.0
+ U235 = '92235' 6.27118E-1 1
+ U238 = '92238' 8.75256E+1 1
+ U236 = '92236' 0.0 1
+ Pu239 = '94239' 0.0 1
+ MIX 7 COMB 6 1.0
+ MIX 8 COMB 6 1.0
+ MIX 9 COMB 6 1.0
+ MIX 10 560.66 6.44 Fe56 = '26056' 1.60000E-1
+ Ni58 = '28058' 6.00000E-2 Cr52 = '24052' 1.10000E-1
+ B10 = '5010' 6.16900E-05 B11 = '5011' 2.48310E-04
+ Zr90 = '40090' 5.13008E+01 Zr91 = '40091' 1.11875E+01 2
+ Zr92 = '40092' 1.71003E+01 Zr94 = '40094' 1.73296E+01
+ Zr96 = '40096' 2.79188E+00
+ MIX 11 560.66 0.81212 O16 = '8016' 7.99449E-1
+ D2D2O = '1002' 1.99768E-1 H1H2O = '1001' 7.83774E-4
+ MIX 12 560.66 6.57 Nb93 = '41093' 2.50000
+ B10 = '5010' 4.17900E-05 B11 = '5011' 1.68210E-04
+ Zr90 = '40090' 5.01638E+01 Zr91 = '40091' 1.09395E+01 2
+ Zr92 = '40092' 1.67213E+01 Zr94 = '40094' 1.69455E+01
+ Zr96 = '40096' 2.73000E+00
+ MIX 13 345.66 0.0014 He4 = '2004' 1.00000E+2
+ MIX 14 345.66 6.44 Fe56 = '26056' 1.60000E-1
+ Ni58 = '28058' 6.00000E-2 Cr52 = '24052' 1.10000E-1
+ B10 = '5010' 6.16900E-05 B11 = '5011' 2.48310E-04
+ Zr90 = '40090' 5.13008E+01 Zr91 = '40091' 1.11875E+01 2
+ Zr92 = '40092' 1.71003E+01 Zr94 = '40094' 1.73296E+01
+ Zr96 = '40096' 2.79188E+00
+ MIX 15 345.66 1.082885 O16 = '8016' 7.98895E-1
+ D2D2O = '1002' 2.01016E-1 H1H2O = '1001' 8.96000E-5
+ MIX 16 941.29 10.4375010 O16 = '8016' 1.18473E+1
+ Xe135 = '54135' 0.0
+ U235 = '92235' 6.27118E-1 1
+ U238 = '92238' 8.75256E+1 1
+ U236 = '92236' 0.0 1
+ Pu239 = '94239' 0.0 1
+ MIX 17 COMB 6 1.0
+ MIX 18 COMB 6 1.0
+ MIX 19 COMB 6 1.0
+ MIX 20 560.66 6.44 Fe56 = '26056' 1.60000E-1
+ Ni58 = '28058' 6.00000E-2 Cr52 = '24052' 1.10000E-1
+ B10 = '5010' 6.16900E-05 B11 = '5011' 2.48310E-04
+ Zr90 = '40090' 5.13008E+01 Zr91 = '40091' 1.11875E+01 2
+ Zr92 = '40092' 1.71003E+01 Zr94 = '40094' 1.73296E+01
+ Zr96 = '40096' 2.79188E+00
+ MIX 21 345.66 Fe56 = '26056' 6.19027E-2
+ Cr52 = '24052' 1.56659E-3 Ni58 = '28058' 6.83337E-3
+ Si29 = '14029' 7.79072E-4 C12 = '6000' 1.46552E-4
+ Mn55 = '25055' 1.25431E-3
+ ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE31.c2m b/Dragon/data/twimsE_proc/TCWE31.c2m
new file mode 100644
index 0000000..084b688
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE31.c2m
@@ -0,0 +1,172 @@
+*----
+* TEST CASE TCWE31
+* CANDU-6 ANNULAR CELL
+* WIMS-E Library
+* MULTI-PARAMETER COMPO ACCESS FOR MACRO-DEPLETION
+* TWO GROUP BURNUP
+* POWER (KW) = 615.00000
+* BURN POWER (KW/KG) = 31.97130
+* URANIUM MASS = 19.23600
+* UO2 REAL DENSITY = 10.59300
+* UO2 EFF DENSITY = 10.43750
+* UO2 TEMPERATURE = 941.28998
+* ENRICHMENT = 0.71140
+* COOLANT D2 AT % = 99.222
+* MODERATOR D2 AT % = 99.911
+* NUMBER OF DAYS = 50
+*
+*----
+* Define variables
+* Burnup paremeters
+* a) Power
+* = 31.9713 kw/kg for 0.0 to 300.0 days
+* b) 69 Groups Burnup time interval Delt
+* = 300 day for 0 to 300 day
+* c) 2 Groups Burnup time interval Delt
+* = 1 day for 0 to 1 day
+* = 4 days for 1 to 5 days
+* = 5 days for 5 to 10 days
+* = 10 days for 10 to 50 days
+* = 20 days for 50 to 150 days
+* = 50 days for 150 to 300 days
+* c) Days with burnup interval changes
+* = 1.0, 5.0, 10.0, 50.0, 150.0 and 300.0 days
+* d) Burnup control time variables Timei, Timef
+* Timei = initial time
+* Timef = final time
+*----
+REAL
+ Power Delt Timec Timei Timef :=
+ 31.9713 1.0 300.0 0.0 0.0 ;
+*----
+* Define STRUCTURES and MODULES used
+*----
+LINKED_LIST LIBRARY ;
+LINKED_LIST
+ CANDU6S CANDU6F VOLMATS VOLMATF PIJ FLUX BURNUP EDITION DATABASE ;
+SEQ_BINARY
+ INTLINS INTLINF ;
+SEQ_ASCII
+ res ;
+MODULE
+ GEO: EXCELT: LIB: SHI: ASM: FLU: EVO: EDI: COMPO: SPH: DELETE: END: ;
+PROCEDURE assertS ;
+*----
+* Depletion data from file WIMS-E format WIMSE
+* Microscopic cross sections from file WIMS-E format WIMSE
+*----
+PROCEDURE TCWE05Lib ;
+INTEGER iedit := 1 ;
+LIBRARY := TCWE05Lib :: <<iedit>> ;
+*----
+* Geometry CANDU6S : 13 regions annular cluster for self-shielding
+* CANDU6F : 31 regions annular cluster for transport
+*----
+CANDU6S := GEO: :: TUBE 5
+ R+ REFL RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 16.12171
+ MIX 1 2 3 4 5
+ CLUSTER ROD1 ROD2 ROD3 ROD4
+ ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000
+ RADIUS 0.00000 0.6122 0.6540 ;
+ ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ;
+ ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ;
+ ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ;
+ ;
+CANDU6F := GEO: CANDU6S :: SPLITR 6 1 1 1 10
+ ::: ROD1 := GEO: ROD1 SPLITR 2 1 ;
+ ::: ROD2 := GEO: ROD2 SPLITR 2 1 ;
+ ::: ROD3 := GEO: ROD3 SPLITR 2 1 ;
+ ::: ROD4 := GEO: ROD4 SPLITR 2 1 ;
+ ;
+*----
+* Self-Shielding calculation EXCEL
+* Transport calculation EXCEL
+* Flux calculation for keff
+*----
+VOLMATS INTLINS := EXCELT: CANDU6S ::
+ TITLE 'TCWE31: FEW GROUP BURNUP / SELF-SHIELDING TRACKING'
+ EDIT 0 MAXR 13 TRAK TISO 5 10.0 SYMM 12 ;
+LIBRARY := SHI: LIBRARY VOLMATS INTLINS :: EDIT 0 NOLJ ;
+VOLMATF INTLINF := EXCELT: CANDU6F ::
+ TITLE 'TCWE31: FEW GROUP BURNUP / TRANSPORT TRACKING'
+ EDIT 0 MAXR 31 TRAK TISO 5 10.0 SYMM 12 ;
+PIJ := ASM: LIBRARY VOLMATF INTLINF ;
+FLUX := FLU: PIJ LIBRARY VOLMATF ::
+ TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.124889 ;
+EDITION := EDI: LIBRARY VOLMATF FLUX CANDU6F ::
+ MERG REGI 6 6 10 7 7 10 1 1 8 8 10 1
+ 1 9 9 10 1 1 2 3 4 5 5 5
+ 5 5 5 5 5 5 5
+ COND 4.0 MICR ALL SAVE
+ MGEO CANDU6F
+ ;
+EDITION := SPH: EDITION VOLMATF INTLINF ;
+FLUX PIJ LIBRARY INTLINS VOLMATS CANDU6S VOLMATF CANDU6F
+INTLINF := DELETE: FLUX PIJ LIBRARY INTLINS VOLMATS CANDU6S
+VOLMATF CANDU6F INTLINF ;
+*----
+* 2 group Burnup (macro-depletion)
+*----
+DATABASE := COMPO: ::
+ EDIT 5
+ COMM 'Multi-parameter reactor database' ENDC
+ INIT
+ ;
+DATABASE := COMPO: DATABASE EDITION ::
+ EDIT 3
+ ;
+EDITION := DELETE: EDITION ;
+*
+* COMPO INTERPOLATION
+LIBRARY := SPH: DATABASE ::
+ EDIT 1
+ STEP UP default STEP AT 1
+ MICRO OFF
+ ;
+CANDU6F := DATABASE :: STEP UP default STEP UP 'GEOMETRIES' STEP AT 1 ;
+VOLMATF INTLINF := EXCELT: CANDU6F ::
+ TITLE 'TCWE31: FEW GROUP BURNUP / TRANSPORT TRACKING'
+ EDIT 0 MAXR 31 TRAK TISO 5 10.0 SYMM 12 ;
+PIJ := ASM: LIBRARY VOLMATF INTLINF ;
+FLUX := FLU: PIJ LIBRARY VOLMATF ::
+ TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.124861 ;
+EDITION := EDI: LIBRARY VOLMATF FLUX ::
+ EDIT 1 MERGE COMP MICR ALL SAVE ;
+EVALUATE Timec := 1.0 ;
+WHILE Timei Timec < DO
+ EVALUATE Timef := Timei Delt + ;
+ IF Timei 0.0 = THEN
+ BURNUP LIBRARY := EVO: LIBRARY FLUX VOLMATF ::
+ EDIT 3 DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ;
+ ELSE
+ BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX VOLMATF ::
+ EDIT 3 NOEX DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ;
+ ENDIF ;
+ PIJ := DELETE: PIJ ;
+ PIJ := ASM: LIBRARY VOLMATF INTLINF ;
+ FLUX := FLU: FLUX PIJ LIBRARY VOLMATF ::
+ TYPE K ;
+ EDITION := EDI: EDITION LIBRARY VOLMATF FLUX ::
+ SAVE ;
+*----
+* change delta t for burnup and final time if required
+*----
+ IF Timef Timec = THEN
+ IF Timec 5.0 = THEN
+ EVALUATE Delt Timec := 5.0 10.0 ;
+ ENDIF ;
+ IF Timec 1.0 = THEN
+ EVALUATE Delt Timec := 4.0 5.0 ;
+ ENDIF ;
+ ENDIF ;
+ EVALUATE Timei := Timef ;
+ENDWHILE ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.079416 ;
+res := EDITION ;
+EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F := DELETE:
+ EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F ;
+ECHO "test TCWE31 completed" ;
+END: ;
+QUIT "LIST" .
diff --git a/Dragon/data/twimsE_proc/TCWE32.c2m b/Dragon/data/twimsE_proc/TCWE32.c2m
new file mode 100644
index 0000000..68096ab
--- /dev/null
+++ b/Dragon/data/twimsE_proc/TCWE32.c2m
@@ -0,0 +1,97 @@
+*----
+* TEST CASE TCWE32
+* VVER Assembly of central water hole with 6 identical fuel cells
+* NXT self-shielding and flux
+* Coarse mesh
+*----
+MODULE
+ LIB: GEO: NXT: SHI: ASM: FLU: DELETE: END: ;
+LINKED_LIST
+ LIBRARY GEOMS GEOMF TRACK PIJMAT FLUX ;
+SEQ_BINARY
+ TRACK_LGT ;
+PROCEDURE assertS ;
+*----
+* Microscopic cross sections from file WNEALIB format WIMSE
+*----
+LIBRARY := LIB: ::
+ EDIT 1
+ NMIX 4
+ CTRA WIMS
+ MIXS LIB: WIMSE FIL: WIMS-E
+ MIX 1 579.9
+ H1H2O = '1001' 4.76690E-2 O16H2O = '8016' 2.38345E-2
+ BNat = '1011' 2.38103E-5
+ MIX 2 609.0 O16 = '8016' 3.06711E-4
+ Cr52 = '24052' 7.54987E-5 Fe56 = '26056' 1.47624E-4
+ Zr90 = '40090' 2.15381E-02 Zr91 = '40091' 4.69693E-03 2
+ Zr92 = '40092' 7.17935E-03 Zr94 = '40094' 7.27563E-03
+ Zr96 = '40096' 1.17214E-03
+ MIX 3 966.0 O16 = '8016' 4.49355E-2
+ U235 = '92235' 7.39237E-4 1
+ U238 = '92238' 2.17285E-2 1
+ MIX 4 579.9
+ H1H2O = '1001' 4.76690E-2 O16H2O = '8016' 2.38345E-2
+ B10 = '5010' 4.73825E-06 B11 = '5011' 1.90721E-05
+ ;
+*----
+* Geometries NXT
+*----
+GEOMS := GEO: :: HEX 7
+ EDIT 2
+ HBC COMPLETE REFL
+ CELL CT CFI CFI CFI CFI CFI CFI
+ ::: CT := GEO: HEXTCEL 2 1
+ SIDE 0.736121593
+ MIX 4 2 1
+ RADIUS 0.0 0.386 0.445
+ ;
+ ::: CFI := GEO: HEXTCEL 2 1
+ SIDE 0.736121593
+ MIX 3 2 1
+ RADIUS 0.0 0.386 0.445
+ ;
+ ;
+GEOMF := GEO: :: HEX 7
+ EDIT 2
+ HBC COMPLETE REFL
+ CELL CT CFI CFI CFI CFI CFI CFI
+ ::: CT := GEO: HEXTCEL 2 1
+ SIDE 0.736121593
+ SPLITR 3 1
+ MIX 1 2 1
+ RADIUS 0.0 0.386 0.445
+ ;
+ ::: CFI := GEO: HEXTCEL 2 1
+ SIDE 0.736121593
+ SPLITR 3 1
+ MIX 3 2 1
+ RADIUS 0.0 0.386 0.445
+ ;
+ ;
+*----
+* Self-Shielding calculation
+* Transport calculation
+* Flux calculation
+*----
+TRACK TRACK_LGT := NXT: GEOMS ::
+ EDIT 3 TITL 'VVER-NXT: Shielding'
+ ALLG BATCH 100
+ TISO 50 20.0 ;
+LIBRARY := SHI: LIBRARY TRACK TRACK_LGT ;
+TRACK TRACK_LGT GEOMS := DELETE: TRACK TRACK_LGT GEOMS ;
+TRACK TRACK_LGT := NXT: GEOMF ::
+ EDIT 3 TITL 'VVER-NXT: Flux'
+ ALLG BATCH 100
+ TISO 50 100.0 ;
+PIJMAT := ASM: LIBRARY TRACK TRACK_LGT ::
+ EDIT 1 ;
+FLUX := FLU: PIJMAT LIBRARY TRACK ::
+ EDIT 1 TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.360705 ;
+TRACK TRACK_LGT GEOMF := DELETE: TRACK TRACK_LGT GEOMF ;
+PIJMAT FLUX := DELETE: PIJMAT FLUX ;
+LIBRARY := DELETE: LIBRARY ;
+ECHO "test TCWE32 completed" ;
+END: ;
+QUIT "LIST" .
generated by cgit v1.2.3 (git 2.46.0) at 2026-02-16 17:29:27 +0000