From 83df2127417dd454376155ebd2e49d725f0e0f4c Mon Sep 17 00:00:00 2001 From: ChasingNeutrons Date: Wed, 22 Jan 2025 17:15:52 +0000 Subject: [PATCH 1/6] Fix to keff for MG scattering --- Tallies/TallyClerks/keffImplicitClerk_class.f90 | 16 ++++++++-------- 1 file changed, 8 insertions(+), 8 deletions(-) diff --git a/Tallies/TallyClerks/keffImplicitClerk_class.f90 b/Tallies/TallyClerks/keffImplicitClerk_class.f90 index 008c50e02..d7f4852dd 100644 --- a/Tallies/TallyClerks/keffImplicitClerk_class.f90 +++ b/Tallies/TallyClerks/keffImplicitClerk_class.f90 @@ -26,7 +26,7 @@ module keffImplicitClerk_class private - !! Locations of diffrent bins wrt memory Address of the clerk + !! Locations of different bins wrt memory Address of the clerk integer(shortInt), parameter :: MEM_SIZE = 5 integer(longInt), parameter :: IMP_PROD = 0 ,& ! Implicit neutron production (from fission) SCATTER_PROD = 1 ,& ! Analog Stattering production (N,XN) @@ -38,7 +38,7 @@ module keffImplicitClerk_class !! and on analog estimators of (N,XN) reactions and leakage !! !! Private Members: - !! targetSTD -> Target Standard Deviation for convergance check + !! targetSTD -> Target Standard Deviation for convergence check !! !! Interface: !! tallyClerk interface @@ -96,10 +96,10 @@ subroutine init(self, dict, name) ! Set name call self % setName(name) - ! Configure convergance trigger + ! Configure convergence trigger call dict % getOrDefault(chr,'trigger','no') - ! Read convergance target + ! Read convergence target if( charCmp(chr,'yes')) then call dict % get(self % targetSTD,'SDtarget') @@ -126,7 +126,7 @@ elemental subroutine kill(self) end subroutine kill !! - !! Returns array of codes that represent diffrent reports + !! Returns array of codes that represent different reports !! !! See tallyClerk_inter for details !! @@ -225,7 +225,7 @@ subroutine reportOutColl(self, p, MT, muL, xsData, mem) score = 2.0_defReal * p % preCollision % wgt case(N_4N) score = 3.0_defReal * p % preCollision % wgt - case(macroAllScatter) ! Catch weight change for MG scattering + case(macroAllScatter, macroIEScatter) ! Catch weight change for MG scattering score = max(p % w - p % preCollision % wgt, ZERO) case default score = ZERO @@ -289,7 +289,7 @@ subroutine reportCycleEnd(self, end, mem) end subroutine reportCycleEnd !! - !! Perform convergance check in the Clerk + !! Perform convergence check in the Clerk !! !! See tallyClerk_inter for details !! @@ -306,7 +306,7 @@ function isConverged(self, mem) result(isIt) end function isConverged !! - !! Display convergance progress on the console + !! Display convergence progress on the console !! !! See tallyClerk_inter for details !! From 0791fe251f443601ed3b11332c7ac187d5ab13c4 Mon Sep 17 00:00:00 2001 From: ChasingNeutrons Date: Thu, 6 Feb 2025 19:15:33 +0000 Subject: [PATCH 2/6] BEAVRS_HZP and geometry fixes Added in the BEAVRS model with the D-bank partially inserted. Also fixed a flaw in the geometry where the outermost cell was defined such that there could be a rare particle lost between it and the geometry boundary. --- InputFiles/Benchmarks/BEAVRS/BEAVRS2D | 3 +- .../Benchmarks/BEAVRS/{BEAVRS => BEAVRS_ARO} | 9 +- InputFiles/Benchmarks/BEAVRS/BEAVRS_HZP | 2410 +++++++++++++++++ 3 files changed, 2416 insertions(+), 6 deletions(-) rename InputFiles/Benchmarks/BEAVRS/{BEAVRS => BEAVRS_ARO} (97%) create mode 100644 InputFiles/Benchmarks/BEAVRS/BEAVRS_HZP diff --git a/InputFiles/Benchmarks/BEAVRS/BEAVRS2D b/InputFiles/Benchmarks/BEAVRS/BEAVRS2D index 0ace8d16c..899a59c71 100644 --- a/InputFiles/Benchmarks/BEAVRS/BEAVRS2D +++ b/InputFiles/Benchmarks/BEAVRS/BEAVRS2D @@ -114,7 +114,8 @@ geometry { thickGrid {type simpleCell; id 55; surfaces (90 -91); filltype mat; material Inconel;} thinGrid {type simpleCell; id 56; surfaces (-92 93); filltype mat; material Inconel;} - pressureVessel { type simpleCell; id 7; surfaces (-1 2); filltype mat; material CarbonSteel;} + ! Does not use surface 1 to define it as that bounds the geometry + pressureVessel { type simpleCell; id 7; surfaces (2); filltype mat; material CarbonSteel;} RPVLiner { type simpleCell; id 8; surfaces (-2 3); filltype mat; material SS304;} outerWater1 {type simpleCell; id 9; surfaces (-3 4 ); filltype mat; material Water;} diff --git a/InputFiles/Benchmarks/BEAVRS/BEAVRS b/InputFiles/Benchmarks/BEAVRS/BEAVRS_ARO similarity index 97% rename from InputFiles/Benchmarks/BEAVRS/BEAVRS rename to InputFiles/Benchmarks/BEAVRS/BEAVRS_ARO index 52e5cb575..dd3ed25c5 100644 --- a/InputFiles/Benchmarks/BEAVRS/BEAVRS +++ b/InputFiles/Benchmarks/BEAVRS/BEAVRS_ARO @@ -3,14 +3,12 @@ !! This is not a fully faithful replica of BEAVRS at HZP: !! all rods are fully withdrawn. !! In reality, some rods are partially inserted. -!! TODO: add rods in the appropriate assemblies at the -!! correct heights! !! type eigenPhysicsPackage; -pop 10000000; +pop 1000000; active 50; -inactive 200; +inactive 150; XSdata ce; dataType ce; @@ -216,7 +214,8 @@ geometry { thickGrid {type simpleCell; id 555; surfaces (90 ); filltype mat; material Inconel;} thinGrid {type simpleCell; id 556; surfaces (92 ); filltype mat; material Zircaloy;} - pressureVessel { type simpleCell; id 7; surfaces (-1 2); filltype mat; material CarbonSteel;} + // Don't need to bound PV by 1 since it is the bounding surface of the geometry. + pressureVessel { type simpleCell; id 7; surfaces ( 2); filltype mat; material CarbonSteel;} RPVLiner { type simpleCell; id 8; surfaces (-2 3); filltype mat; material SS304;} outerWater1 {type simpleCell; id 9; surfaces (-3 4 ); filltype mat; material Water;} diff --git a/InputFiles/Benchmarks/BEAVRS/BEAVRS_HZP b/InputFiles/Benchmarks/BEAVRS/BEAVRS_HZP new file mode 100644 index 000000000..5a740ce28 --- /dev/null +++ b/InputFiles/Benchmarks/BEAVRS/BEAVRS_HZP @@ -0,0 +1,2410 @@ +!! +!! 3D BEAVRS benchmark +!! At HZP, the D bank of rods is partially inserted up +!! to 115 steps inserted / 113 steps withdrawn. +!! A step corresponds to an increment of 1.58193cm +!! +type eigenPhysicsPackage; + +pop 1000000; +active 50; +inactive 250; +XSdata ce; +dataType ce; + +collisionOperator { neutronCE {type neutronCEstd;}} + +transportOperator { + !type transportOperatorDT; + type transportOperatorHT; cache 1; + } + +inactiveTally { + shannon { + type shannonEntropyClerk; + map {type multiMap; + maps (xax yax zax); + xax { type spaceMap; grid lin; min -161.2773; max 161.2773; N 15; axis x;} + yax { type spaceMap; grid lin; min -161.2773; max 161.2773; N 15; axis y;} + zax { type spaceMap; grid lin; min 36.748; max 402.508; N 15; axis z;} + } + cycles 200; + } + +} + +activeTally { + pinFissRadial { type collisionClerk; response (fission); fission { type macroResponse; MT -6;} + map {type multiMap; maps (xax yax); + xax {type spaceMap; axis x; grid lin; N 255; min -161.2773; max 161.2773; } + yax {type spaceMap; axis y; grid lin; N 255; min -161.2773; max 161.2773; } + } + } + assemblyFissRadial { type collisionClerk; response (fission); fission { type macroResponse; MT -6;} + map {type multiMap; maps (xax yax); + xax {type spaceMap; axis x; grid lin; N 15; min -161.2773; max 161.2773; } + yax {type spaceMap; axis y; grid lin; N 15; min -161.2773; max 161.2773; } + } + } + fissionAxial { type collisionClerk; response (fission); fission { type macroResponse; MT -6;} + map {type spaceMap; axis z; grid lin; N 60; min 36.748; max 402.508;} + } + fissionYZ { type collisionClerk; response (fission); fission { type macroResponse; MT -6;} + map {type multiMap; maps (yax zax); + yax {type spaceMap; axis y; grid lin; N 255; min -161.2773; max 161.2773; } + zax {type spaceMap; axis z; grid lin; N 60; min 36.748; max 402.508;} + } + } +} + +geometry { + type geometryStd; + boundary ( 0 0 0 0 0 0); + graph {type shrunk;} + + surfaces { + + // thickness specifications for RPV and RPV liner + outerRPV { id 1; type zTruncCylinder; radius 241.3; origin (0.0 0.0 230.0); halfwidth 230; } + innerRPV { id 2; type zCylinder; radius 219.710; origin (0.0 0.0 0.0); } + innerRPVLiner { id 3; type zCylinder; radius 219.150; origin (0.0 0.0 0.0); } + + // thickness specifications for neutron shield + outerBoundNS { id 4; type zCylinder; radius 201.630; origin (0.0 0.0 0.0); } + innerBoundNS { id 5; type zCylinder; radius 194.84; origin (0.0 0.0 0.0); } + + // thickness specifications for core barrel + outerCoreBarrel { id 6; type zCylinder; radius 193.675; origin (0.0 0.0 0.0); } + innerCoreBarrel { id 7; type zCylinder; radius 187.96; origin (0.0 0.0 0.0); } + + // four planes that intersect to bound the Neutron shield panel + P1 { id 8; type plane; coeffs (-0.48480962025 0.87461970714 0.0 0.0);} + P2 { id 9; type plane; coeffs (-0.87461970714 0.48480962025 0.0 0.0);} + P3 { id 10; type plane; coeffs (-0.87461970714 -0.48480962025 0.0 0.0);} + P4 { id 11; type plane; coeffs (-0.48480962025 -0.87461970714 0.0 0.0);} + + // bounding widths for baffle on various sides + // right & left refers to the side of the reactor that it is on + // close/away refers to its location in relation to the LATTICE it is a part of + // (NOT the reactor itself) + rightClose { id 50; type plane; coeffs (1.0 0.0 0.0 8.36662);} + rightAway { id 51; type plane; coeffs (1.0 0.0 0.0 10.58912);} + leftClose { id 52; type plane; coeffs (-1.0 0.0 0.0 8.36662);} + leftAway { id 53; type plane; coeffs (-1.0 0.0 0.0 10.58912);} + bottomClose { id 54; type plane; coeffs (0.0 -1.0 0.0 8.36662);} + bottomAway { id 55; type plane; coeffs (0.0 -1.0 0.0 10.58912);} + topClose { id 56; type plane; coeffs (0.0 1.0 0.0 8.36662);} + topAway { id 57; type plane; coeffs (0.0 1.0 0.0 10.58912);} + + // thickness specifications for grid with thickness of 0.0198cm (Inconel) + pinThickGridInner { id 90; type zSquareCylinder; origin (0.0 0.0 0.0); halfwidth (0.61015 0.61015 0.0); } + pinThickGridOuter { id 91; type zSquareCylinder; origin (0.0 0.0 0.0); halfwidth (0.62992 0.62992 0.0); } + + // thickness specifications for grid with thickness of 0.0194cm (Zircaloy) + pinThinGridInner { id 92; type zSquareCylinder; origin (0.0 0.0 0.0); halfwidth (0.61049 0.61049 0.0); } + pinThinGridOuter { id 93; type zSquareCylinder; origin (0.0 0.0 0.0); halfwidth (0.62992 0.62992 0.0); } + + // inner and outer surfaces of assembly sleeves (both SS and Zircaloy) + assemblySleeveInner { id 94; type zSquareCylinder; origin (0.0 0.0 0.0); halfwidth (10.70864 10.70864 0.0); } + assemblySleeveOuter { id 95; type zSquareCylinder; origin (0.0 0.0 0.0); halfwidth (10.74798 10.74798 0.0); } + + + // Axial planes across core height + // Names are based on axial heights + plane460 { id 100; type plane; coeffs (0.0 0.0 1.0 460.0 ); } + plane431p876 { id 101; type plane; coeffs (0.0 0.0 1.0 431.876); } + plane423p049 { id 102; type plane; coeffs (0.0 0.0 1.0 423.049); } + plane421p532 { id 103; type plane; coeffs (0.0 0.0 1.0 421.532); } + plane419p704 { id 104; type plane; coeffs (0.0 0.0 1.0 419.704); } + plane417p164 { id 105; type plane; coeffs (0.0 0.0 1.0 417.164); } + plane415p164 { id 106; type plane; coeffs (0.0 0.0 1.0 415.164); } + plane411p806 { id 107; type plane; coeffs (0.0 0.0 1.0 411.806); } + plane403p778 { id 108; type plane; coeffs (0.0 0.0 1.0 403.778); } + plane402p508 { id 109; type plane; coeffs (0.0 0.0 1.0 402.508); } + plane401p238 { id 110; type plane; coeffs (0.0 0.0 1.0 401.238); } + plane364p725 { id 111; type plane; coeffs (0.0 0.0 1.0 364.725); } + plane359p01 { id 112; type plane; coeffs (0.0 0.0 1.0 359.01 ); } + plane312p528 { id 113; type plane; coeffs (0.0 0.0 1.0 312.528); } + plane306p813 { id 114; type plane; coeffs (0.0 0.0 1.0 306.813); } + plane260p331 { id 115; type plane; coeffs (0.0 0.0 1.0 260.331); } + plane254p616 { id 116; type plane; coeffs (0.0 0.0 1.0 254.616); } + plane208p134 { id 117; type plane; coeffs (0.0 0.0 1.0 208.134); } + plane202p419 { id 118; type plane; coeffs (0.0 0.0 1.0 202.419); } + plane155p937 { id 119; type plane; coeffs (0.0 0.0 1.0 155.937); } + plane150p222 { id 120; type plane; coeffs (0.0 0.0 1.0 150.222); } + plane143p428 { id 121; type plane; coeffs (0.0 0.0 1.0 143.428); } + plane103p74 { id 122; type plane; coeffs (0.0 0.0 1.0 103.74 ); } + plane98p025 { id 123; type plane; coeffs (0.0 0.0 1.0 98.025 ); } + plane41p828 { id 124; type plane; coeffs (0.0 0.0 1.0 41.828 ); } + plane40p558 { id 125; type plane; coeffs (0.0 0.0 1.0 40.558 ); } + plane40p52 { id 126; type plane; coeffs (0.0 0.0 1.0 40.52 ); } + plane39p958 { id 127; type plane; coeffs (0.0 0.0 1.0 39.958 ); } + plane38p66 { id 128; type plane; coeffs (0.0 0.0 1.0 38.66 ); } + plane37p1621 { id 129; type plane; coeffs (0.0 0.0 1.0 37.1621); } + plane36p748 { id 130; type plane; coeffs (0.0 0.0 1.0 36.748 ); } + plane35 { id 131; type plane; coeffs (0.0 0.0 1.0 35.0 ); } + plane20 { id 132; type plane; coeffs (0.0 0.0 1.0 20.0 ); } + plane0 { id 133; type plane; coeffs (0.0 0.0 1.0 0.0 ); } + + planeSteelBottom { id 134; type plane; coeffs (0.0 0.0 1.0 400.638); } + planeCRLowerBottom { id 135; type plane; coeffs (0.0 0.0 1.0 402.508); } // same as 109 on withdrawal + planeCRUpperBottom { id 136; type plane; coeffs (0.0 0.0 1.0 504.108); } // out of core on withdrawal + + plane322p1861 {id 137; type plane; coeffs (0 0 1 322.18609);} + plane220p586 {id 138; type plane; coeffs (0 0 1 220.5861);} + plane218p716 {id 139; type plane; coeffs (0 0 1 218.71609);} + + } + + cells { + + // assembly wrappers and surrounding water at various heights + wrapper1 {type simpleCell; id 2001; surfaces (94 -95 129 -126); filltype mat; material SS304;} + wrapper2 {type simpleCell; id 2002; surfaces (94 -95 123 -122); filltype mat; material Zircaloy;} + wrapper3 {type simpleCell; id 2003; surfaces (94 -95 120 -119); filltype mat; material Zircaloy;} + wrapper4 {type simpleCell; id 2004; surfaces (94 -95 118 -117); filltype mat; material Zircaloy;} + wrapper5 {type simpleCell; id 2005; surfaces (94 -95 116 -115); filltype mat; material Zircaloy;} + wrapper6 {type simpleCell; id 2006; surfaces (94 -95 114 -113); filltype mat; material Zircaloy;} + wrapper7 {type simpleCell; id 2007; surfaces (94 -95 112 -111); filltype mat; material Zircaloy;} + wrapper8 {type simpleCell; id 2008; surfaces (94 -95 107 -106); filltype mat; material SS304;} + + assemWater0 {type simpleCell; id 2009; surfaces (94 -129); filltype mat; material Water;} + assemWater1 {type simpleCell; id 2010; surfaces (94 -95 126 -123); filltype mat; material Water;} + assemWater2 {type simpleCell; id 2011; surfaces (94 -95 122 -120); filltype mat; material Water;} + assemWater3 {type simpleCell; id 2012; surfaces (94 -95 119 -118); filltype mat; material Water;} + assemWater4 {type simpleCell; id 2013; surfaces (94 -95 117 -116); filltype mat; material Water;} + assemWater5 {type simpleCell; id 2014; surfaces (94 -95 115 -114); filltype mat; material Water;} + assemWater6 {type simpleCell; id 2015; surfaces (94 -95 113 -112); filltype mat; material Water;} + assemWater7 {type simpleCell; id 2016; surfaces (94 -95 111 -107); filltype mat; material Water;} + assemWater8 {type simpleCell; id 2017; surfaces (94 -95 106); filltype mat; material Water;} + assemWaterEx {type simpleCell; id 2018; surfaces (95); filltype mat; material Water;} + + // assemblies inside the wrappers + assem1424 {type simpleCell; id 2019; surfaces (-94); filltype uni; universe 1424;} + assem1416 {type simpleCell; id 2020; surfaces (-94); filltype uni; universe 1416;} + assem1431 {type simpleCell; id 2021; surfaces (-94); filltype uni; universe 1431;} + assem60316 {type simpleCell; id 2022; surfaces (-94); filltype uni; universe 60316;} + assem603112 {type simpleCell; id 2023; surfaces (-94); filltype uni; universe 603112;} + assem60313 {type simpleCell; id 2024; surfaces (-94); filltype uni; universe 60313;} + assem60319 {type simpleCell; id 2025; surfaces (-94); filltype uni; universe 60319;} + assem15315 {type simpleCell; id 2026; surfaces (-94); filltype uni; universe 15315;} + assem15317 {type simpleCell; id 2027; surfaces (-94); filltype uni; universe 15317;} + assem15311 {type simpleCell; id 2028; surfaces (-94); filltype uni; universe 15311;} + assem153111 {type simpleCell; id 2029; surfaces (-94); filltype uni; universe 153111;} + assem1631 {type simpleCell; id 2030; surfaces (-94); filltype uni; universe 1631;} + assem2031 {type simpleCell; id 2031; surfaces (-94); filltype uni; universe 2031;} + assem1224 {type simpleCell; id 2032; surfaces (-94); filltype uni; universe 1224;} + assem1624 {type simpleCell; id 2033; surfaces (-94); filltype uni; universe 1624;} + + // unrodded assemblies in wrappers + assem2424 {type simpleCell; id 2034; surfaces (-94); filltype uni; universe 2424;} + assem2416 {type simpleCell; id 2035; surfaces (-94); filltype uni; universe 2416;} + assem2431 {type simpleCell; id 2036; surfaces (-94); filltype uni; universe 2431;} + assem70316 {type simpleCell; id 2037; surfaces (-94); filltype uni; universe 70316;} + assem703112 {type simpleCell; id 2038; surfaces (-94); filltype uni; universe 703112;} + assem70313 {type simpleCell; id 2039; surfaces (-94); filltype uni; universe 70313;} + assem70319 {type simpleCell; id 2040; surfaces (-94); filltype uni; universe 70319;} + assem25315 {type simpleCell; id 2041; surfaces (-94); filltype uni; universe 25315;} + assem25317 {type simpleCell; id 2042; surfaces (-94); filltype uni; universe 25317;} + assem25311 {type simpleCell; id 2043; surfaces (-94); filltype uni; universe 25311;} + assem253111 {type simpleCell; id 2044; surfaces (-94); filltype uni; universe 253111;} + assem2631 {type simpleCell; id 2045; surfaces (-94); filltype uni; universe 2631;} + assem3031 {type simpleCell; id 2046; surfaces (-94); filltype uni; universe 3031;} + assem2224 {type simpleCell; id 2047; surfaces (-94); filltype uni; universe 2224;} + assem2624 {type simpleCell; id 2048; surfaces (-94); filltype uni; universe 2624;} + + // rodded assemblies in wrappers + assem1425 {type simpleCell; id 2049; surfaces (-94); filltype uni; universe 1425;} + assem1417 {type simpleCell; id 2050; surfaces (-94); filltype uni; universe 1417;} + + // pin grids - thick at the top and bottom, thin in fuelled region + thickGrid {type simpleCell; id 555; surfaces (90 ); filltype mat; material Inconel;} + thinGrid {type simpleCell; id 556; surfaces (92 ); filltype mat; material Zircaloy;} + + // Don't need to bound PV by 1 since it is the bounding surface of the geometry. + pressureVessel { type simpleCell; id 7; surfaces (2); filltype mat; material CarbonSteel;} + RPVLiner { type simpleCell; id 8; surfaces (-2 3); filltype mat; material SS304;} + outerWater1 {type simpleCell; id 9; surfaces (-3 4 ); filltype mat; material Water;} + + // Neutron shields + NS1 { type simpleCell; id 10; surfaces (-4 5 -8 9); filltype mat; material SS304;} + NS2 { type simpleCell; id 11; surfaces (-4 5 8 -9); filltype mat; material SS304;} + NS3 { type simpleCell; id 12; surfaces (-4 5 -10 11); filltype mat; material SS304;} + NS4 { type simpleCell; id 13; surfaces (-4 5 10 -11); filltype mat; material SS304;} + + // Water in the arc between neutron shields + outerWaterSeg1 {type simpleCell; id 14; surfaces (-4 5 ); filltype mat; material Water;} + outerWater2 {type simpleCell; id 15; surfaces (-5 6 ); filltype mat; material Water;} + + // Outer core + coreBarrel { type simpleCell; id 16; surfaces (-6 7); filltype mat; material SS304;} + core {type simpleCell; id 17; surfaces (-5); filltype uni; universe 9999;} + + + // Gridded pins + + // THICK GRID + // 2.4% in grid + grid24Thick {type simpleCell; id 253; surfaces (-90); filltype uni; universe 24000;} + // guide tube in grid + gridGTThick {type simpleCell; id 254; surfaces (-90); filltype uni; universe 12000;} + // 3.1% in grid + grid31Thick {type simpleCell; id 255; surfaces (-90); filltype uni; universe 31000;} + // 1.6 % in grid + grid16Thick {type simpleCell; id 256; surfaces (-90); filltype uni; universe 16000;} + // instrumentation tube in grid + gridITThick {type simpleCell; id 257; surfaces (-90); filltype uni; universe 14000;} + // empty GT at dashpot in grid + gridGTDPThick {type simpleCell; id 258; surfaces (-90); filltype uni; universe 1010;} + // pin upper fuel plenum in grid + gridFPPThick {type simpleCell; id 259; surfaces (-90); filltype uni; universe 1008;} + // stainless steel in guide tube in grid + gridSSGThick {type simpleCell; id 260; surfaces (-90); filltype uni; universe 1023;} + // stainless steel in dash pot in grid + gridSSDPThick {type simpleCell; id 261; surfaces (-90); filltype uni; universe 1024;} + // BP plenum in grid + gridBPPThick {type simpleCell; id 262; surfaces (-90); filltype uni; universe 1012;} + // Lower rodded GT in grid + gridLRGTThick {type simpleCell; id 263; surfaces (-90); filltype uni; universe 1014;} + // Upper rodded GT in grid + gridURGTThick {type simpleCell; id 463; surfaces (-90); filltype uni; universe 1013;} + + // THIN GRID + // 2.4% in grid + grid24Thin {type simpleCell; id 264; surfaces (-92); filltype uni; universe 24000;} + // guide tube in grid + gridGTThin {type simpleCell; id 265; surfaces (-92); filltype uni; universe 12000;} + // burnable poison in grid + gridBPThin {type simpleCell; id 266; surfaces (-92); filltype uni; universe 1000;} + // 3.1% in grid + grid31Thin {type simpleCell; id 267; surfaces (-92); filltype uni; universe 31000;} + // 1.6 % in grid + grid16Thin {type simpleCell; id 268; surfaces (-92); filltype uni; universe 16000;} + // instrumentation tube in grid + gridITThin {type simpleCell; id 269; surfaces (-92); filltype uni; universe 14000;} + // empty GT at dashpot in grid + gridGTDPThin {type simpleCell; id 270; surfaces (-92); filltype uni; universe 1010;} + // pin upper fuel plenum in grid + gridFPPThin {type simpleCell; id 271; surfaces (-92); filltype uni; universe 1008;} + // stainless steel in guide tube in grid + gridSSGTThin {type simpleCell; id 272; surfaces (-92); filltype uni; universe 1023;} + // stainless steel in dash pot in grid + gridSSDPThin {type simpleCell; id 273; surfaces (-92); filltype uni; universe 1024;} + // BP plenum in grid + gridBPPThin {type simpleCell; id 274; surfaces (-92); filltype uni; universe 1012;} + // Lower rodded GT in grid (not used when rods fully withdrawn) + gridLRGTThin {type simpleCell; id 275; surfaces (-92); filltype uni; universe 1014;} + // Upper rodded GT in grid (not used when rods fully withdrawn) + gridURGTThin {type simpleCell; id 475; surfaces (-92); filltype uni; universe 1013;} + + + // 3.1% enriched pins, axial layering + 31FP460 {type simpleCell; id 100; surfaces ( 101); filltype uni; universe 1001;} + 31FP431 {type simpleCell; id 101; surfaces (-101 102); filltype uni; universe 1003;} + 31FP423 {type simpleCell; id 102; surfaces (-102 104); filltype uni; universe 1001;} + 31FP419 {type simpleCell; id 103; surfaces (-104 105); filltype uni; universe 1006;} + 31FP417 {type simpleCell; id 104; surfaces (-105 106); filltype uni; universe 1008;} + 31FP415 {type simpleCell; id 105; surfaces (-106 107); filltype uni; universe 1017;} + 31FP411 {type simpleCell; id 106; surfaces (-107 109); filltype uni; universe 1008;} + 31FP402 {type simpleCell; id 107; surfaces (-109 111); filltype uni; universe 31000;} + 31FP364 {type simpleCell; id 108; surfaces (-111 112); filltype uni; universe 9231;} + 31FP359 {type simpleCell; id 109; surfaces (-112 113); filltype uni; universe 31000;} + 31FP312 {type simpleCell; id 110; surfaces (-113 114); filltype uni; universe 9231;} + 31FP306 {type simpleCell; id 111; surfaces (-114 115); filltype uni; universe 31000;} + 31FP260 {type simpleCell; id 112; surfaces (-115 116); filltype uni; universe 9231;} + 31FP254 {type simpleCell; id 113; surfaces (-116 117); filltype uni; universe 31000;} + 31FP208 {type simpleCell; id 114; surfaces (-117 118); filltype uni; universe 9231;} + 31FP202 {type simpleCell; id 115; surfaces (-118 119); filltype uni; universe 31000;} + 31FP155 {type simpleCell; id 116; surfaces (-119 120); filltype uni; universe 9231;} + 31FP150 {type simpleCell; id 117; surfaces (-120 122); filltype uni; universe 31000;} + 31FP103 {type simpleCell; id 118; surfaces (-122 123); filltype uni; universe 9231;} + 31FP98 {type simpleCell; id 119; surfaces (-123 126); filltype uni; universe 31000;} + 31FP4052 {type simpleCell; id 120; surfaces (-126 129); filltype uni; universe 1131;} + 31FP37 {type simpleCell; id 121; surfaces (-129 130); filltype uni; universe 31000;} + 31FP36 {type simpleCell; id 122; surfaces (-130 131); filltype uni; universe 1006;} + 31FP35 {type simpleCell; id 123; surfaces (-131 132); filltype uni; universe 1003;} + 31FP20 {type simpleCell; id 124; surfaces (-132 ); filltype uni; universe 1001;} + + + //2.4% enriched pins, axial layering + 24FP460 {type simpleCell; id 125; surfaces ( 101); filltype uni; universe 1001;} + 24FP431 {type simpleCell; id 126; surfaces (-101 102); filltype uni; universe 1003;} + 24FP423 {type simpleCell; id 127; surfaces (-102 104); filltype uni; universe 1001;} + 24FP419 {type simpleCell; id 128; surfaces (-104 105); filltype uni; universe 1006;} + 24FP417 {type simpleCell; id 129; surfaces (-105 106); filltype uni; universe 1008;} + 24FP415 {type simpleCell; id 130; surfaces (-106 107); filltype uni; universe 1017;} + 24FP411 {type simpleCell; id 131; surfaces (-107 109); filltype uni; universe 1008;} + 24FP402 {type simpleCell; id 132; surfaces (-109 111); filltype uni; universe 24000;} + 24FP364 {type simpleCell; id 133; surfaces (-111 112); filltype uni; universe 9224;} + 24FP359 {type simpleCell; id 134; surfaces (-112 113); filltype uni; universe 24000;} + 24FP312 {type simpleCell; id 135; surfaces (-113 114); filltype uni; universe 9224;} + 24FP306 {type simpleCell; id 136; surfaces (-114 115); filltype uni; universe 24000;} + 24FP260 {type simpleCell; id 137; surfaces (-115 116); filltype uni; universe 9224;} + 24FP254 {type simpleCell; id 138; surfaces (-116 117); filltype uni; universe 24000;} + 24FP208 {type simpleCell; id 139; surfaces (-117 118); filltype uni; universe 9224;} + 24FP202 {type simpleCell; id 140; surfaces (-118 119); filltype uni; universe 24000;} + 24FP155 {type simpleCell; id 141; surfaces (-119 120); filltype uni; universe 9224;} + 24FP150 {type simpleCell; id 142; surfaces (-120 122); filltype uni; universe 24000;} + 24FP103 {type simpleCell; id 143; surfaces (-122 123); filltype uni; universe 9224;} + 24FP98 {type simpleCell; id 144; surfaces (-123 126); filltype uni; universe 24000;} + 24FP4052 {type simpleCell; id 145; surfaces (-126 129); filltype uni; universe 1124;} + 24FP37 {type simpleCell; id 146; surfaces (-129 130); filltype uni; universe 24000;} + 24FP36 {type simpleCell; id 147; surfaces (-130 131); filltype uni; universe 1006;} + 24FP35 {type simpleCell; id 148; surfaces (-131 132); filltype uni; universe 1003;} + 24FP20 {type simpleCell; id 149; surfaces (-132); filltype uni; universe 1001;} + + //1.6% enriched pins, axial layering + 16FP460 {type simpleCell; id 150; surfaces ( 101); filltype uni; universe 1001;} + 16FP431 {type simpleCell; id 151; surfaces (-101 102); filltype uni; universe 1003;} + 16FP423 {type simpleCell; id 152; surfaces (-102 104); filltype uni; universe 1001;} + 16FP419 {type simpleCell; id 153; surfaces (-104 105); filltype uni; universe 1006;} + 16FP417 {type simpleCell; id 154; surfaces (-105 106); filltype uni; universe 1008;} + 16FP415 {type simpleCell; id 155; surfaces (-106 107); filltype uni; universe 1017;} + 16FP411 {type simpleCell; id 156; surfaces (-107 109); filltype uni; universe 1008;} + 16FP402 {type simpleCell; id 157; surfaces (-109 111); filltype uni; universe 16000;} + 16FP364 {type simpleCell; id 158; surfaces (-111 112); filltype uni; universe 9216;} + 16FP359 {type simpleCell; id 159; surfaces (-112 113); filltype uni; universe 16000;} + 16FP312 {type simpleCell; id 160; surfaces (-113 114); filltype uni; universe 9216;} + 16FP306 {type simpleCell; id 161; surfaces (-114 115); filltype uni; universe 16000;} + 16FP260 {type simpleCell; id 162; surfaces (-115 116); filltype uni; universe 9216;} + 16FP254 {type simpleCell; id 163; surfaces (-116 117); filltype uni; universe 16000;} + 16FP208 {type simpleCell; id 164; surfaces (-117 118); filltype uni; universe 9216;} + 16FP202 {type simpleCell; id 165; surfaces (-118 119); filltype uni; universe 16000;} + 16FP155 {type simpleCell; id 166; surfaces (-119 120); filltype uni; universe 9216;} + 16FP150 {type simpleCell; id 167; surfaces (-120 122); filltype uni; universe 16000;} + 16FP103 {type simpleCell; id 168; surfaces (-122 123); filltype uni; universe 9216;} + 16FP98 {type simpleCell; id 169; surfaces (-123 126); filltype uni; universe 16000;} + 16FP4052 {type simpleCell; id 170; surfaces (-126 129); filltype uni; universe 1116;} + 16FP37 {type simpleCell; id 171; surfaces (-129 130); filltype uni; universe 16000;} + 16FP36 {type simpleCell; id 172; surfaces (-130 131); filltype uni; universe 1006;} + 16FP35 {type simpleCell; id 173; surfaces (-131 132); filltype uni; universe 1003;} + 16FP20 {type simpleCell; id 174; surfaces (-132); filltype uni; universe 1001;} + + + //guide tube, with CR, axial layering + GC460 {type simpleCell; id 175; surfaces ( 101); filltype uni; universe 1014;} + GC431 {type simpleCell; id 176; surfaces (-101 106); filltype uni; universe 1014;} // Nozzle/support plate BW? Replaced just with rod + GC415 {type simpleCell; id 177; surfaces (-106 107); filltype uni; universe 1018;} // Rodded w/ grid + GC411 {type simpleCell; id 178; surfaces (-107 109); filltype uni; universe 1014;} + GC402 {type simpleCell; id 179; surfaces (-109 134); filltype uni; universe 1015;} + GC400 {type simpleCell; id 180; surfaces (-134 111); filltype uni; universe 12000;} + GC364 {type simpleCell; id 181; surfaces (-111 112); filltype uni; universe 9112;} + GC359 {type simpleCell; id 182; surfaces (-112 113); filltype uni; universe 12000;} + GC312 {type simpleCell; id 183; surfaces (-113 114); filltype uni; universe 9112;} + GC306 {type simpleCell; id 184; surfaces (-114 115); filltype uni; universe 12000;} + GC260 {type simpleCell; id 185; surfaces (-115 116); filltype uni; universe 9112;} + GC254 {type simpleCell; id 186; surfaces (-116 117); filltype uni; universe 12000;} + GC208 {type simpleCell; id 187; surfaces (-117 118); filltype uni; universe 9112;} + GC202 {type simpleCell; id 188; surfaces (-118 119); filltype uni; universe 12000;} + GC155 {type simpleCell; id 189; surfaces (-119 120); filltype uni; universe 9112;} + GC150 {type simpleCell; id 190; surfaces (-120 122); filltype uni; universe 12000;} + GC103 {type simpleCell; id 191; surfaces (-122 123); filltype uni; universe 9112;} + GC98 {type simpleCell; id 192; surfaces (-123 126); filltype uni; universe 1010;} + GC4052 {type simpleCell; id 193; surfaces (-126 127); filltype uni; universe 12000;} + GC39 {type simpleCell; id 194; surfaces (-127 129); filltype uni; universe 1019;} + GC37 {type simpleCell; id 195; surfaces (-129 131); filltype uni; universe 1010;} + GC35 {type simpleCell; id 196; surfaces (-131 132); filltype uni; universe 1005;} + GC20 {type simpleCell; id 197; surfaces (-132); filltype uni; universe 1001;} + + + // instrumentation tube, axial layering + IT460 {type simpleCell; id 198; surfaces ( 102); filltype uni; universe 1001;} + IT423 {type simpleCell; id 199; surfaces (-102 106); filltype uni; universe 14000;} + IT415 {type simpleCell; id 200; surfaces (-106 107); filltype uni; universe 1114;} + IT411 {type simpleCell; id 201; surfaces (-107 111); filltype uni; universe 14000;} + IT364 {type simpleCell; id 202; surfaces (-111 112); filltype uni; universe 9114;} + IT359 {type simpleCell; id 203; surfaces (-112 113); filltype uni; universe 14000;} + IT312 {type simpleCell; id 204; surfaces (-113 114); filltype uni; universe 9114;} + IT306 {type simpleCell; id 205; surfaces (-114 115); filltype uni; universe 14000;} + IT260 {type simpleCell; id 206; surfaces (-115 116); filltype uni; universe 9114;} + IT254 {type simpleCell; id 207; surfaces (-116 117); filltype uni; universe 14000;} + IT208 {type simpleCell; id 208; surfaces (-117 118); filltype uni; universe 9114;} + IT202 {type simpleCell; id 209; surfaces (-118 119); filltype uni; universe 14000;} + IT155 {type simpleCell; id 210; surfaces (-119 120); filltype uni; universe 9114;} + IT150 {type simpleCell; id 211; surfaces (-120 122); filltype uni; universe 14000;} + IT103 {type simpleCell; id 212; surfaces (-122 123); filltype uni; universe 9114;} + IT98 {type simpleCell; id 213; surfaces (-123 126); filltype uni; universe 14000;} + IT4052 {type simpleCell; id 214; surfaces (-126 129); filltype uni; universe 1114;} + IT37 {type simpleCell; id 215; surfaces (-129 131); filltype uni; universe 14000;} + IT35 {type simpleCell; id 216; surfaces (-131 132); filltype uni; universe 1005;} + IT20 {type simpleCell; id 217; surfaces (-132 ); filltype uni; universe 1011;} + + + // burnable absorber, axial layering + BA460 {type simpleCell; id 218; surfaces ( 101); filltype uni; universe 1001;} + BA431 {type simpleCell; id 219; surfaces (-101 102); filltype uni; universe 1002;} + BA423 {type simpleCell; id 220; surfaces (-102 103); filltype uni; universe 1023;} + BA421 {type simpleCell; id 230; surfaces (-103 106); filltype uni; universe 1012;} + BA415 {type simpleCell; id 231; surfaces (-106 107); filltype uni; universe 1027;} + BA411 {type simpleCell; id 232; surfaces (-107 110); filltype uni; universe 1012;} + BA401 {type simpleCell; id 233; surfaces (-110 111); filltype uni; universe 1000;} + BA364 {type simpleCell; id 234; surfaces (-111 112); filltype uni; universe 1110;} + BA359 {type simpleCell; id 235; surfaces (-112 113); filltype uni; universe 1000;} + BA312 {type simpleCell; id 236; surfaces (-113 114); filltype uni; universe 1110;} + BA306 {type simpleCell; id 237; surfaces (-114 115); filltype uni; universe 1000;} + BA260 {type simpleCell; id 238; surfaces (-115 116); filltype uni; universe 1110;} + BA254 {type simpleCell; id 239; surfaces (-116 117); filltype uni; universe 1000;} + BA208 {type simpleCell; id 240; surfaces (-117 118); filltype uni; universe 1110;} + BA202 {type simpleCell; id 241; surfaces (-118 119); filltype uni; universe 1000;} + BA155 {type simpleCell; id 242; surfaces (-119 120); filltype uni; universe 1110;} + BA150 {type simpleCell; id 243; surfaces (-120 122); filltype uni; universe 1000;} + BA103 {type simpleCell; id 244; surfaces (-122 123); filltype uni; universe 1110;} + BA98 {type simpleCell; id 245; surfaces (-123 125); filltype uni; universe 1000;} + BA4055 {type simpleCell; id 246; surfaces (-125 126); filltype uni; universe 1023;} + BA4052 {type simpleCell; id 247; surfaces (-126 127); filltype uni; universe 1021;} + BA39 {type simpleCell; id 248; surfaces (-127 128); filltype uni; universe 1025;} + BA38 {type simpleCell; id 249; surfaces (-128 129); filltype uni; universe 1019;} + BA37 {type simpleCell; id 250; surfaces (-129 131); filltype uni; universe 1010;} + BA35 {type simpleCell; id 251; surfaces (-131 132); filltype uni; universe 1005;} + BA20 {type simpleCell; id 252; surfaces (-132 ); filltype uni; universe 1001;} + + //guide tube, no CR, axial layering + GT460 {type simpleCell; id 353; surfaces ( 101); filltype uni; universe 1001;} + GT431 {type simpleCell; id 354; surfaces (-101 102); filltype uni; universe 1005;} + GT423 {type simpleCell; id 355; surfaces (-102 106); filltype uni; universe 12000;} + GT415 {type simpleCell; id 356; surfaces (-106 107); filltype uni; universe 1112;} + GT411 {type simpleCell; id 357; surfaces (-107 111); filltype uni; universe 12000;} + GT364 {type simpleCell; id 358; surfaces (-111 112); filltype uni; universe 9112;} + GT359 {type simpleCell; id 359; surfaces (-112 113); filltype uni; universe 12000;} + GT312 {type simpleCell; id 360; surfaces (-113 114); filltype uni; universe 9112;} + GT306 {type simpleCell; id 361; surfaces (-114 115); filltype uni; universe 12000;} + GT260 {type simpleCell; id 362; surfaces (-115 116); filltype uni; universe 9112;} + GT254 {type simpleCell; id 363; surfaces (-116 117); filltype uni; universe 12000;} + GT208 {type simpleCell; id 364; surfaces (-117 118); filltype uni; universe 9112;} + GT202 {type simpleCell; id 365; surfaces (-118 119); filltype uni; universe 12000;} + GT155 {type simpleCell; id 366; surfaces (-119 120); filltype uni; universe 9112;} + GT150 {type simpleCell; id 367; surfaces (-120 122); filltype uni; universe 12000;} + GT103 {type simpleCell; id 368; surfaces (-122 123); filltype uni; universe 9112;} + GT98 {type simpleCell; id 369; surfaces (-123 126); filltype uni; universe 1010;} + GT4052 {type simpleCell; id 370; surfaces (-126 127); filltype uni; universe 12000;} + GT39 {type simpleCell; id 371; surfaces (-127 129); filltype uni; universe 1019;} + GT37 {type simpleCell; id 372; surfaces (-129 131); filltype uni; universe 1010;} + GT35 {type simpleCell; id 373; surfaces (-131 132); filltype uni; universe 1005;} + GT20 {type simpleCell; id 374; surfaces (-132); filltype uni; universe 1001;} + + //guide tube, partially inserted CR, axial layering + // Note: universe is basically a combo of fully and partially inserted rods + GP460 {type simpleCell; id 375; surfaces ( 101); filltype uni; universe 1013;} + GP431 {type simpleCell; id 376; surfaces (-101 102); filltype uni; universe 1013;} + GP423 {type simpleCell; id 377; surfaces (-102 106); filltype uni; universe 1013;} + GP415 {type simpleCell; id 378; surfaces (-106 107); filltype uni; universe 1133;} + GP411 {type simpleCell; id 379; surfaces (-107 111); filltype uni; universe 1013;} + GP364 {type simpleCell; id 380; surfaces (-111 112); filltype uni; universe 9233;} + GP359 {type simpleCell; id 381; surfaces (-112 137); filltype uni; universe 1013;} + GP322 {type simpleCell; id 382; surfaces (-137 113); filltype uni; universe 1014;} + GP312 {type simpleCell; id 383; surfaces (-113 114); filltype uni; universe 9232;} + GP306 {type simpleCell; id 384; surfaces (-114 115); filltype uni; universe 1014;} + GP260 {type simpleCell; id 385; surfaces (-115 116); filltype uni; universe 9232;} + GP254 {type simpleCell; id 386; surfaces (-116 138); filltype uni; universe 1014;} + GP220 {type simpleCell; id 387; surfaces (-138 139); filltype uni; universe 1023;} + GP219 {type simpleCell; id 388; surfaces (-139 117); filltype uni; universe 12000;} + GP208 {type simpleCell; id 389; surfaces (-117 118); filltype uni; universe 9112;} + GP202 {type simpleCell; id 390; surfaces (-118 119); filltype uni; universe 12000;} + GP155 {type simpleCell; id 391; surfaces (-119 120); filltype uni; universe 9112;} + GP150 {type simpleCell; id 392; surfaces (-120 122); filltype uni; universe 12000;} + GP103 {type simpleCell; id 393; surfaces (-122 123); filltype uni; universe 9112;} + GP98 {type simpleCell; id 394; surfaces (-123 126); filltype uni; universe 1010;} + GP4052 {type simpleCell; id 395; surfaces (-126 127); filltype uni; universe 12000;} + GP39 {type simpleCell; id 396; surfaces (-127 129); filltype uni; universe 1019;} + GP37 {type simpleCell; id 397; surfaces (-129 131); filltype uni; universe 1010;} + GP35 {type simpleCell; id 398; surfaces (-131 132); filltype uni; universe 1005;} + GP20 {type simpleCell; id 399; surfaces (-132); filltype uni; universe 1001;} + + // control rod, axial layering + // Used (probably with some modification) only when fully inserted + //CR460 {type simpleCell; id 448; surfaces (-120 121); filltype uni; universe 1002;} + //CR415 {type simpleCell; id 449; surfaces (-126 403); filltype uni; universe 1013;} + //CR403 {type simpleCell; id 450; surfaces (-403 402); filltype uni; universe 1015;} + //CR402 {type simpleCell; id 451; surfaces (-402 1430); filltype uni; universe 1013;} + //CR143 {type simpleCell; id 452; surfaces (-1430 41); filltype uni; universe 1014;} + //CR41 {type simpleCell; id 453; surfaces (-41 143); filltype uni; universe 1002;} + //CR39 {type simpleCell; id 454; surfaces (-143 147); filltype uni; universe 1001;} + //CR35 {type simpleCell; id 455; surfaces (-147 148); filltype uni; universe 1005;} + //CR20 {type simpleCell; id 456; surfaces (-148 149); filltype uni; universe 1001;} + + + + outsideLeftBaffle { type simpleCell; id 52; surfaces (-50); filltype mat; material Water;} + leftBaffle { type simpleCell; id 53; surfaces (50 -51); filltype mat; material SS304;} + insideLeftBaffle { type simpleCell; id 54; surfaces (51); filltype mat; material Water;} + + outsideRightBaffle { type simpleCell; id 55; surfaces (-52); filltype mat; material Water;} + RightBaffle { type simpleCell; id 56; surfaces (52 -53); filltype mat; material SS304;} + insideRightBaffle { type simpleCell; id 57; surfaces (53); filltype mat; material Water;} + + outsideTopBaffle { type simpleCell; id 58; surfaces (-54); filltype mat; material Water;} + TopBaffle { type simpleCell; id 59; surfaces (54 -55); filltype mat; material SS304;} + insideTopBaffle { type simpleCell; id 60; surfaces (55); filltype mat; material Water;} + + outsideBottomBaffle { type simpleCell; id 61; surfaces (-56); filltype mat; material Water;} + BottomBaffle { type simpleCell; id 62; surfaces (56 -57); filltype mat; material SS304;} + insideBottomBaffle { type simpleCell; id 63; surfaces (57); filltype mat; material Water;} + + topLeftCornerBaffle1 { type simpleCell; id 64; surfaces (52 -53 -57); filltype mat; material SS304;} + topLeftCornerBaffle2 { type simpleCell; id 65; surfaces (56 -57 -52); filltype mat; material SS304;} + topLeftCornerGap1 { type simpleCell; id 66; surfaces (57); filltype mat; material Water;} + topLeftCornerGap2 { type simpleCell; id 67; surfaces (53); filltype mat; material Water;} + topLeftMajorGap { type simpleCell; id 68; surfaces (-56 -52); filltype mat; material Water;} + + topRightCornerBaffle1 { type simpleCell; id 69; surfaces (-57 50 -51); filltype mat; material SS304;} + topRightCornerBaffle2 { type simpleCell; id 70; surfaces (-50 56 -57); filltype mat; material SS304;} + topRightCornerGap1 { type simpleCell; id 71; surfaces (57); filltype mat; material Water;} + topRightCornerGap2 { type simpleCell; id 72; surfaces (51); filltype mat; material Water;} + topRightMajorGap { type simpleCell; id 73; surfaces (-56 -50); filltype mat; material Water;} + + bottomLeftCornerBaffle1 { type simpleCell; id 74; surfaces (-55 52 -53); filltype mat; material SS304;} + bottomLeftCornerBaffle2 { type simpleCell; id 75; surfaces (-55 54 -52); filltype mat; material SS304;} + bottomLeftCornerGap1 { type simpleCell; id 76; surfaces (55); filltype mat; material Water;} + bottomLeftCornerGap2 { type simpleCell; id 77; surfaces (53); filltype mat; material Water;} + bottomLeftMajorGap { type simpleCell; id 78; surfaces (-54 -52); filltype mat; material Water;} + + bottomRightCornerBaffle1 { type simpleCell; id 79; surfaces (-51 50 -55); filltype mat; material SS304;} + bottomRightCornerBaffle2 { type simpleCell; id 80; surfaces (-55 54 -50); filltype mat; material SS304;} + bottomRightCornerGap1 { type simpleCell; id 81; surfaces (51); filltype mat; material Water;} + bottomRightCornerGap2 { type simpleCell; id 82; surfaces (55); filltype mat; material Water;} + bottomRightMajorGap { type simpleCell; id 83; surfaces (-50 -54); filltype mat; material Water;} + + + TLSG1 { type simpleCell; id 84; surfaces (-56 -52); filltype mat; material Water;} + TLSG2 { type simpleCell; id 85; surfaces (-56 52); filltype mat; material Water;} + TLSG3 { type simpleCell; id 86; surfaces (56 -52); filltype mat; material Water;} + topLeftSquare { type simpleCell; id 87; surfaces (56 52); filltype mat; material SS304;} + + TRSG1 { type simpleCell; id 88; surfaces (-56 50); filltype mat; material Water;} + TRSG2 { type simpleCell; id 89; surfaces (-56 -50); filltype mat; material Water;} + TRSG3 { type simpleCell; id 90; surfaces (56 -50); filltype mat; material Water;} + topRightSquare { type simpleCell; id 91; surfaces (56 50); filltype mat; material SS304;} + + BLSG1 { type simpleCell; id 92; surfaces (54 -52); filltype mat; material Water;} + BLSG2 { type simpleCell; id 93; surfaces (-54 52); filltype mat; material Water;} + BLSG3 { type simpleCell; id 94; surfaces (-54 -52); filltype mat; material Water;} + bottomLeftSquare { type simpleCell; id 95; surfaces (54 52); filltype mat; material SS304;} + + BRSG1 { type simpleCell; id 96; surfaces (-54 50); filltype mat; material Water;} + BRSG2 { type simpleCell; id 97; surfaces (54 -50); filltype mat; material Water;} + BRSG3 { type simpleCell; id 98; surfaces (-54 -50); filltype mat; material Water;} + bottomRightSquare { type simpleCell; id 99; surfaces (54 50); filltype mat; material SS304;} + } + + universes { + root { id 1; type rootUniverse; border 1; fill u<8888>; } + + // Pin universes + + //Burnable poison + pinBPaboveDP { id 1000; type pinUniverse; radii (0.21400 0.23051 0.24130 0.42672 0.43688 0.48387 0.56134 0.60198 0.0); + fills (Air SS304 Helium BorosilicateGlass Helium SS304 Water Zircaloy Water);} + pinBPPlenumGeometry { id 1012; type pinUniverse; radii ( 0.21400 0.23051 0.43688 0.48387 0.50419 0.54610 0.0); + fills (Air SS304 Helium SS304 Water Zircaloy Water);} + + //guide tubes + pinGTaboveDP { id 12000; type pinUniverse; radii (0.56134 0.60198 0.0 ); fills (Water Zircaloy Water);} + pinGTatDP { id 1010; type pinUniverse; radii (0.50419 0.54610 0.0); fills (Water Zircaloy Water);} + + //INST Tube + pinIT { id 14000; type pinUniverse; radii (0.43688 0.48387 0.56134 0.60198 0.0 ); + fills (Air Zircaloy Water Zircaloy Water);} + pinBareInstrumentThimble { id 1011; type pinUniverse; radii (0.43688 0.48387 0.0); fills (Air Zircaloy Water);} + + // Fuel pins + pin16 { id 16000; type pinUniverse; radii (0.39218 0.40005 0.45720 0.0); + fills (UO2-16 Helium Zircaloy Water);} + pin24 { id 24000; type pinUniverse; radii (0.39218 0.40005 0.45720 0.0); + fills (UO2-24 Helium Zircaloy Water);} + pin31 { id 31000; type pinUniverse; radii (0.39218 0.40005 0.45720 0.0); + fills (UO2-31 Helium Zircaloy Water);} + // Higher enrichments not used + //pin32 { id 32000; type pinUniverse; radii (0.39218 0.40005 0.45720 0.0); + // fills (UO2-32 Helium Zircaloy Water);} + //pin34 { id 34000; type pinUniverse; radii (0.39218 0.40005 0.45720 0.0); + // fills (UO2-34 Helium Zircaloy Water);} + + pinWater { id 1001; type pinUniverse; radii ( 0.0); fills (Water);} + + + // Solid pins, assumed radius to be that of a fuel pin (0.45720) + pinNozzle_SupportSteel { id 1003; type pinUniverse; radii ( 0.45720 0.0); fills (SupportPlateSS Water);} + pinSupportPlateBW { id 1005; type pinUniverse; radii ( 0.45720 0.0); fills (SupportPlateBW Water);} + pinZircaloy { id 1006; type pinUniverse; radii ( 0.45720 0.0); fills (Zircaloy Water);} + + + SSinDashPot { id 1024; type pinUniverse; radii (0.50419 0.54610 0.0); fills (SS304 Zircaloy Water);} + SSinGuideTube { id 1023; type pinUniverse; radii ( 0.56134 0.60198 0.0); fills (SS304 Zircaloy Water);} + SSnoGuideTube { id 1002; type pinUniverse; radii ( 0.56134 0.0); fills (SS304 Water);} + + + pinUpperFuelPlenum { id 1008; type pinUniverse; radii ( 0.06459 0.40005 0.45720 0.0); + fills (Inconel Helium Zircaloy Water);} + + // Control rod pins + pinControlRodUpper { id 1013; type pinUniverse; radii ( 0.37338 0.38608 0.48387 0.56134 0.60198 0.0); + fills (B4C Helium SS304 Water Zircaloy Water);} + pinControlRodLower { id 1014; type pinUniverse; radii ( 0.38227 0.38608 0.48387 0.56134 0.60198 0.0); + fills (Ag-In-Cd Helium SS304 Water Zircaloy Water);} + pinControlRodSpacer { id 1015; type pinUniverse; radii ( 0.37845 0.38608 0.48387 0.56134 0.60198 0.0); + fills (SS304 Helium SS304 Water Zircaloy Water);} + pinControlRodPlenum { id 1016; type pinUniverse; radii ( 0.06459 0.38608 0.48387 0.56134 0.60198 0.0); + fills (Inconel Helium SS304 Water Zircaloy Water);} + + // pins that have grids + fuelRodPlenumWithGridThick { + id 1017; + type cellUniverse; + cells ( 259 555);} + + GTRodThick { + id 1018; + type cellUniverse; + cells (263 555);} + + dashPotGuideTubeGridThick { + id 1019; + type cellUniverse; + cells ( 258 555);} + + dashPotGuideTubeGridThin { + id 1020; + type cellUniverse; + cells ( 270 556);} + + SSinGuideTubeThick { + id 1021; + type cellUniverse; + cells ( 260 555);} + + SSinGuideTubeThin { + id 1022; + type cellUniverse; + cells ( 272 556);} + + SSinDashPotThick { + id 1025; + type cellUniverse; + cells ( 261 555);} + + SSinDashPotThin { + id 1026; + type cellUniverse; + cells ( 273 556);} + + BPPlenumThick { + id 1027; + type cellUniverse; + cells ( 262 555);} + + BPPlenumThin { + id 1028; + type cellUniverse; + cells ( 274 556);} + + BPaboveDPThin { + id 1110; + type cellUniverse; + cells (266 556);} + + GTThick { + id 1112; + type cellUniverse; + cells (254 555);} + + ITThick { + id 1114; + type cellUniverse; + cells (257 555);} + + pin16Thick { + id 1116; + type cellUniverse; + cells (256 555);} + + pin24Thick { + id 1124; + type cellUniverse; + cells (253 555);} + + pin31Thick { + id 1131; + type cellUniverse; + cells (255 555);} + + LowerRodGTThick { + id 1132; + type cellUniverse; + cells (263 555);} + + UpperRodGTThick { + id 1133; + type cellUniverse; + cells (463 555);} + + BPThin { // Is this necessary given 1110, BPaboveDPThin??? + id 9110; + type cellUniverse; + cells (266 556);} + + GTThin { + id 9112; + type cellUniverse; + cells (265 556);} + + ITThin { + id 9114; + type cellUniverse; + cells (269 556);} + + pin16Thin { + id 9216; + type cellUniverse; + cells (268 556);} + + pin24Thin { + id 9224; + type cellUniverse; + cells (264 556);} + + pin31Thin { + id 9231; + type cellUniverse; + cells (267 556);} + + LowerRodGTThin { + id 9232; + type cellUniverse; + cells (275 556);} + + UpperRodGTThin { + id 9233; + type cellUniverse; + cells (475 556);} + + // Axial stacks of universes to make up full pins + + // 3.1 % + fuelPin31 { + id 31; + type cellUniverse; + cells ( 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124);} + + // 2.4 % + fuelPin24 { + id 24; + type cellUniverse; + cells ( 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149);} + + //1.6 % + fuelPin16 { + id 16; + type cellUniverse; + cells ( 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174);} + + //burnable absorber + BP { + id 10; + type cellUniverse; + cells (218 219 220 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252);} + + // guide tube, with CR + GuideTubeRodded { + id 12; + type cellUniverse; + cells (175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197);} + + //control rod, fully retracted + GuideTubeEmpty { + id 13; + type cellUniverse; + cells (353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374);} + + //instr. tube + instrumentTube { + id 14; + type cellUniverse; + cells (198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217);} + + + // control rod, partially inserted + GuideTubePartial { + id 15; + type cellUniverse; + cells (375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399); + + } + + + + // Lattices w/o grid + // Names represent AE + A0E24 { + id 1424; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 12 24 24 12 24 24 12 24 24 24 24 24 + 24 24 24 12 24 24 24 24 24 24 24 24 24 12 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 12 24 24 12 24 24 12 24 24 12 24 24 12 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 12 24 24 12 24 24 14 24 24 12 24 24 12 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 12 24 24 12 24 24 12 24 24 12 24 24 12 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 12 24 24 24 24 24 24 24 24 24 12 24 24 24 + 24 24 24 24 24 12 24 24 12 24 24 12 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); } + + // assembly with sleeves at different heights + A0E24Sleeve { + id 14240; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2019 + );} + + + A0E16 { + id 1416; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 16 16 12 16 16 12 16 16 12 16 16 16 16 16 + 16 16 16 12 16 16 16 16 16 16 16 16 16 12 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 12 16 16 12 16 16 12 16 16 12 16 16 12 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 12 16 16 12 16 16 14 16 16 12 16 16 12 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 12 16 16 12 16 16 12 16 16 12 16 16 12 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 12 16 16 16 16 16 16 16 16 16 12 16 16 16 + 16 16 16 16 16 12 16 16 12 16 16 12 16 16 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16); } + + A0E16Sleeve { + id 14160; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2020 + ); + } + + A0E31 { + id 1431; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 12 31 31 12 31 31 12 31 31 31 31 31 + 31 31 31 12 31 31 31 31 31 31 31 31 31 12 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 12 31 31 12 31 31 12 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 12 31 31 14 31 31 12 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 12 31 31 12 31 31 12 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 12 31 31 31 31 31 31 31 31 31 12 31 31 31 + 31 31 31 31 31 12 31 31 12 31 31 12 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A0E31Sleeve { + id 14310; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2021 + ); + } + + + A6BE31B { + id 60316; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 10 31 31 12 31 31 10 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 12 31 31 12 31 31 12 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 12 31 31 14 31 31 12 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 12 31 31 12 31 31 12 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 12 31 31 31 31 31 31 31 31 31 12 31 31 31 + 31 31 31 31 31 12 31 31 12 31 31 12 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + + A6BE31BSleeve { + id 603160; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2022 + ); + } + + A6BE31T { + id 603112; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 12 31 31 12 31 31 12 31 31 31 31 31 + 31 31 31 12 31 31 31 31 31 31 31 31 31 12 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 12 31 31 12 31 31 12 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 12 31 31 14 31 31 12 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 12 31 31 12 31 31 12 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 10 31 31 12 31 31 10 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A6BE31TSleeve { + id 6031120; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2023 + ); + } + + A6BE31R { + id 60313; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 10 31 31 12 31 31 12 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 12 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 12 31 31 12 31 31 12 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 12 31 31 14 31 31 12 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 12 31 31 12 31 31 12 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 12 31 31 31 + 31 31 31 31 31 10 31 31 12 31 31 12 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A6BE31RSleeve { + id 603130; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2024 + ); + } + + A6BE31L { + id 60319; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 12 31 31 12 31 31 10 31 31 31 31 31 + 31 31 31 12 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 12 31 31 12 31 31 12 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 12 31 31 14 31 31 12 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 12 31 31 12 31 31 12 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 12 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 12 31 31 12 31 31 10 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A6BE31LSleeve { + id 603190; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2025 + ); + } + + A15BE31BR { + id 15315; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 12 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 10 31 31 10 31 31 10 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 10 31 31 14 31 31 10 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 10 31 31 10 31 31 10 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 12 31 31 31 31 31 31 31 31 31 12 31 31 31 + 31 31 31 31 31 12 31 31 12 31 31 12 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A15BE31BRSleeve { + id 153150; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2026 + ); + } + + A15BE31BL { + id 15317; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 + 31 31 31 12 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 10 31 31 10 31 31 10 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 10 31 31 14 31 31 10 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 10 31 31 10 31 31 10 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 12 31 31 31 31 31 31 31 31 31 12 31 31 31 + 31 31 31 31 31 12 31 31 12 31 31 12 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A15BE31BLSleeve { + id 153170; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2027 + ); + } + + A15BE31TR { + id 15311; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 12 31 31 12 31 31 12 31 31 31 31 31 + 31 31 31 12 31 31 31 31 31 31 31 31 31 12 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 10 31 31 10 31 31 10 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 10 31 31 14 31 31 10 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 10 31 31 10 31 31 10 31 31 12 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 12 31 31 31 + 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A15BE31TRSleeve { + id 153110; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2028 + ); + } + + A15BE31TL { + id 153111; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 12 31 31 12 31 31 12 31 31 31 31 31 + 31 31 31 12 31 31 31 31 31 31 31 31 31 12 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 10 31 31 10 31 31 10 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 10 31 31 14 31 31 10 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 12 31 31 10 31 31 10 31 31 10 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 12 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A15BE31TLSleeve { + id 1531110; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2029 + ); + } + + A16BE31 { + id 1631; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 12 31 31 12 31 31 12 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 12 31 31 14 31 31 12 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 12 31 31 12 31 31 12 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A16BE31Sleeve { + id 16310; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2030 + ); + } + + A20BE31 { + id 2031; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 10 31 31 12 31 31 10 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 12 31 31 14 31 31 12 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 10 31 31 12 31 31 10 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A20BE31Sleeve { + id 20310; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2031 + ); + } + + A12BE24 { + id 1224; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 10 24 24 12 24 24 10 24 24 24 24 24 + 24 24 24 10 24 24 24 24 24 24 24 24 24 10 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 10 24 24 12 24 24 12 24 24 12 24 24 10 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 12 24 24 12 24 24 14 24 24 12 24 24 12 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 10 24 24 12 24 24 12 24 24 12 24 24 10 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 10 24 24 24 24 24 24 24 24 24 10 24 24 24 + 24 24 24 24 24 10 24 24 12 24 24 10 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); } + + A12BE24Sleeve { + id 12240; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2032 + ); + } + + A16BE24 { + id 1624; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 10 24 24 10 24 24 10 24 24 24 24 24 + 24 24 24 10 24 24 24 24 24 24 24 24 24 10 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 10 24 24 12 24 24 12 24 24 12 24 24 10 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 10 24 24 12 24 24 14 24 24 12 24 24 10 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 10 24 24 12 24 24 12 24 24 12 24 24 10 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 10 24 24 24 24 24 24 24 24 24 10 24 24 24 + 24 24 24 24 24 10 24 24 10 24 24 10 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); } + + A16BE24Sleeve { + id 16240; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2033 + ); + } + + // Unrodded assemblies + A0E24U { + id 2424; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 13 24 24 13 24 24 13 24 24 24 24 24 + 24 24 24 13 24 24 24 24 24 24 24 24 24 13 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 13 24 24 13 24 24 13 24 24 13 24 24 13 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 13 24 24 13 24 24 14 24 24 13 24 24 13 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 13 24 24 13 24 24 13 24 24 13 24 24 13 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 13 24 24 24 24 24 24 24 24 24 13 24 24 24 + 24 24 24 24 24 13 24 24 13 24 24 13 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); } + + A0E24USleeve { + id 24240; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2034 + );} + + A0E16U { + id 2416; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 16 16 13 16 16 13 16 16 13 16 16 16 16 16 + 16 16 16 13 16 16 16 16 16 16 16 16 16 13 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 13 16 16 13 16 16 13 16 16 13 16 16 13 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 13 16 16 13 16 16 14 16 16 13 16 16 13 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 13 16 16 13 16 16 13 16 16 13 16 16 13 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 13 16 16 16 16 16 16 16 16 16 13 16 16 16 + 16 16 16 16 16 13 16 16 13 16 16 13 16 16 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16); } + + // sleeved + A0E16USleeve { + id 24160; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2035 + ); + } + + A0E31U { + id 2431; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 13 31 31 13 31 31 13 31 31 31 31 31 + 31 31 31 13 31 31 31 31 31 31 31 31 31 13 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 13 31 31 13 31 31 13 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 13 31 31 14 31 31 13 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 13 31 31 13 31 31 13 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 13 31 31 31 31 31 31 31 31 31 13 31 31 31 + 31 31 31 31 31 13 31 31 13 31 31 13 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A0E31USleeve { + id 24310; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2036 + ); + } + + + A6BE31BU { + id 70316; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 10 31 31 13 31 31 10 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 13 31 31 13 31 31 13 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 13 31 31 14 31 31 13 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 13 31 31 13 31 31 13 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 13 31 31 31 31 31 31 31 31 31 13 31 31 31 + 31 31 31 31 31 13 31 31 13 31 31 13 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + + A6BE31BUSleeve { + id 703160; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2037 + ); + } + + A6BE31TU { + id 703112; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 13 31 31 13 31 31 13 31 31 31 31 31 + 31 31 31 13 31 31 31 31 31 31 31 31 31 13 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 13 31 31 13 31 31 13 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 13 31 31 14 31 31 13 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 13 31 31 13 31 31 13 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 10 31 31 13 31 31 10 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A6BE31TUSleeve { + id 7031120; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2038 + ); + } + + A6BE31RU { + id 70313; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 10 31 31 13 31 31 13 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 13 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 13 31 31 13 31 31 13 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 13 31 31 14 31 31 13 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 13 31 31 13 31 31 13 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 13 31 31 31 + 31 31 31 31 31 10 31 31 13 31 31 13 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A6BE31RUSleeve { + id 703130; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2039 + ); + } + + A6BE31LU { + id 70319; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 13 31 31 13 31 31 10 31 31 31 31 31 + 31 31 31 13 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 13 31 31 13 31 31 13 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 13 31 31 14 31 31 13 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 13 31 31 13 31 31 13 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 13 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 13 31 31 13 31 31 10 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A6BE31LUSleeve { + id 703190; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2040 + ); + } + + A15BE31BRU { + id 25315; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 13 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 10 31 31 10 31 31 10 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 10 31 31 14 31 31 10 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 10 31 31 10 31 31 10 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 13 31 31 31 31 31 31 31 31 31 13 31 31 31 + 31 31 31 31 31 13 31 31 13 31 31 13 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A15BE31BRUSleeve { + id 253150; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2041 + ); + } + + A15BE31BLU { + id 25317; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 + 31 31 31 13 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 10 31 31 10 31 31 10 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 10 31 31 14 31 31 10 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 10 31 31 10 31 31 10 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 13 31 31 31 31 31 31 31 31 31 13 31 31 31 + 31 31 31 31 31 13 31 31 13 31 31 13 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A15BE31BLUSleeve { + id 253170; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2042 + ); + } + + A15BE31TRU { + id 25311; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 13 31 31 13 31 31 13 31 31 31 31 31 + 31 31 31 13 31 31 31 31 31 31 31 31 31 13 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 10 31 31 10 31 31 10 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 10 31 31 14 31 31 10 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 10 31 31 10 31 31 10 31 31 13 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 13 31 31 31 + 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A15BE31TRUSleeve { + id 253110; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2043 + ); + } + + A15BE31TLU { + id 253111; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 13 31 31 13 31 31 13 31 31 31 31 31 + 31 31 31 13 31 31 31 31 31 31 31 31 31 13 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 10 31 31 10 31 31 10 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 10 31 31 14 31 31 10 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 13 31 31 10 31 31 10 31 31 10 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 13 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A15BE31TLUSleeve { + id 2531110; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2044 + ); + } + + A16BE31U { + id 2631; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 13 31 31 13 31 31 13 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 13 31 31 14 31 31 13 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 13 31 31 13 31 31 13 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A16BE31USleeve { + id 26310; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2045 + ); + } + + A20BE31U { + id 3031; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 10 31 31 13 31 31 10 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 13 31 31 14 31 31 13 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 10 31 31 10 31 31 13 31 31 10 31 31 10 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 + 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + + A20BE31USleeve { + id 30310; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2046 + ); + } + + A12BE24U { + id 2224; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 10 24 24 13 24 24 10 24 24 24 24 24 + 24 24 24 10 24 24 24 24 24 24 24 24 24 10 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 10 24 24 13 24 24 13 24 24 13 24 24 10 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 13 24 24 13 24 24 14 24 24 13 24 24 13 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 10 24 24 13 24 24 13 24 24 13 24 24 10 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 10 24 24 24 24 24 24 24 24 24 10 24 24 24 + 24 24 24 24 24 10 24 24 13 24 24 10 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); } + + A12BE24USleeve { + id 22240; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2047 + ); + } + + A16BE24U { + id 2624; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 10 24 24 10 24 24 10 24 24 24 24 24 + 24 24 24 10 24 24 24 24 24 24 24 24 24 10 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 10 24 24 13 24 24 13 24 24 13 24 24 10 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 10 24 24 13 24 24 14 24 24 13 24 24 10 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 10 24 24 13 24 24 13 24 24 13 24 24 10 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 10 24 24 24 24 24 24 24 24 24 10 24 24 24 + 24 24 24 24 24 10 24 24 10 24 24 10 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); } + + A16BE24USleeve { + id 26240; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2048 + ); + } + + A0E24DBank { + id 1425; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 15 24 24 15 24 24 15 24 24 24 24 24 + 24 24 24 15 24 24 24 24 24 24 24 24 24 15 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 15 24 24 15 24 24 15 24 24 15 24 24 15 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 15 24 24 15 24 24 14 24 24 15 24 24 15 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 15 24 24 15 24 24 15 24 24 15 24 24 15 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 15 24 24 24 24 24 24 24 24 24 15 24 24 24 + 24 24 24 24 24 15 24 24 15 24 24 15 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); } + + // assembly with sleeves at different heights + A0E24SleeveDBank { + id 14250; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2049 + );} + + A0E16DBank { + id 1417; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.26 1.26 0.0); + shape (17 17 0); + padMat Water; + map ( + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 16 16 15 16 16 15 16 16 15 16 16 16 16 16 + 16 16 16 15 16 16 16 16 16 16 16 16 16 15 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 15 16 16 15 16 16 15 16 16 15 16 16 15 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 15 16 16 15 16 16 14 16 16 15 16 16 15 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 15 16 16 15 16 16 15 16 16 15 16 16 15 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 15 16 16 16 16 16 16 16 16 16 15 16 16 16 + 16 16 16 16 16 15 16 16 15 16 16 15 16 16 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16); } + + A0E16SleeveDBank { + id 14170; + type cellUniverse; + cells ( + 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 + 2016 2017 2018 2050 + ); + } + + + leftBaffleUni { + id 5222; + type cellUniverse; + cells (52 53 54);} + + + rightBaffleUni { + id 5223; + type cellUniverse; + cells (55 56 57);} + + topBaffleUni { + id 5224; + type cellUniverse; + cells (58 59 60);} + + bottomBaffleUni { + id 5225; + type cellUniverse; + cells (61 62 63);} + + + topLeft { + id 5226; + type cellUniverse; + cells (64 65 66 67 68);} + + topRight { + id 5227; + type cellUniverse; + cells ( 69 70 71 72 73);} + + BottomLeft { + id 5228; + type cellUniverse; + cells ( 74 75 76 77 78);} + + BottomRight { + id 5229; + type cellUniverse; + cells ( 79 80 81 82 83);} + + + SQTL { + id 1500; + type cellUniverse; + cells (84 85 86 87);} + + + SQTR { + id 1600; + type cellUniverse; + cells (88 89 90 91);} + + SQBL { + id 1700; + type cellUniverse; + cells (92 93 94 95);} + + SQBR { + id 1800; + type cellUniverse; + cells (96 97 98 99);} + + + latCore { + id 9999; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (21.50364 21.50364 0.0); + shape (17 17 0); + padMat Water; + map ( + 1001 1001 1001 1001 1800 5224 5224 5224 5224 5224 5224 5224 1700 1001 1001 1001 1001 + 1001 1001 1800 5224 5229 24310 7031120 24310 7031120 24310 7031120 24310 5228 5224 1700 1001 1001 + 1001 1800 5229 24310 14310 26310 14160 30310 14160 30310 14160 26310 14310 24310 5228 1700 1001 + 1001 5222 24310 2531110 26240 14160 26240 14160 26240 14160 26240 14160 26240 25311 24310 5223 1001 + 1800 5229 14310 26240 14250 26240 24160 22240 14160 22240 24160 26240 14250 26240 14310 5228 1700 + 5222 24310 26310 14160 26240 24160 22240 24160 22240 24160 22240 24160 26240 14160 26310 24310 5223 + 5222 703190 14160 26240 24160 22240 14160 22240 14160 22240 14160 22240 24160 26240 14160 703130 5223 + 5222 24310 30310 14160 22240 24160 22240 24160 26240 24160 22240 24160 22240 14160 30310 24310 5223 + 5222 703190 14160 26240 14160 22240 14160 26240 14170 26240 14160 22240 14160 26240 14160 703130 5223 + 5222 24310 30310 14160 22240 24160 22240 24160 26240 24160 22240 24160 22240 14160 30310 24310 5223 + 5222 703190 14160 26240 24160 22240 14160 22240 14160 22240 14160 22240 24160 26240 14160 703130 5223 + 5222 24310 26310 14160 26240 24160 22240 24160 22240 24160 22240 24160 26240 14160 26310 24310 5223 + 1600 5227 14310 26240 14250 26240 24160 22240 14160 22240 24160 26240 14250 26240 14310 5226 1500 + 1001 5222 24310 253170 26240 14160 26240 14160 26240 14160 26240 14160 26240 253150 24310 5223 1001 + 1001 1600 5227 24310 14310 26310 14160 30310 14160 30310 14160 26310 14310 24310 5226 1500 1001 + 1001 1001 1600 5225 5227 24310 703160 24310 703160 24310 703160 24310 5226 5225 1500 1001 1001 + 1001 1001 1001 1001 1600 5225 5225 5225 5225 5225 5225 5225 1500 1001 1001 1001 1001 ); } + +! Note partial rodded assemblies end with a 1, i.e., 14170 and 14250 + + coreAndStructures { + id 8888; + type cellUniverse; + cells (7 8 9 10 11 12 13 14 15 16 17);} + + + + } +} + + +viz { + bmpZ { + type bmp; + output imgXY; + what material; + centre (-17.13 240.69 167.74); + width (50 50); + axis z; + res (1000 1000); + } + bmpYZ { + type bmp; + output imgYZ; + what material; + centre (0.0 0.0 232.0); + width (100.0 200.0); + axis x; + res (1000 2000); + } +} + + +nuclearData { + handles { + ce {type aceNeutronDatabase; aceLibrary $SCONE_ACE; ures 0; } + } + + materials { + + // Note that commented nuclide densities are included in the specification + // but are not available in the JEFF-3.11 data library + + Air { + temp 566; + composition { + 18036.06 7.8730E-09; + 18038.06 1.4844E-09; + 18040.06 2.3506E-06; + 6012.06 6.7539E-08; + //6013.06 7.5658E-10; + 7014.06 1.9680E-04; + 7015.06 7.2354E-07; + 8016.06 5.2866E-05; + 8017.06 2.0084E-08; + //8018.06 1.0601E-07; + } + } + + SS304 { + temp 566; + composition { + 24050.06 7.6778E-04; + 24052.06 1.4806E-02; + 24053.06 1.6789E-03; + 24054.06 4.1791E-04; + 26054.06 3.4620E-03; + 26056.06 5.4345E-02; + 26057.06 1.2551E-03; + 26058.06 1.6703E-04; + 25055.06 1.7604E-03; + 28058.06 5.6089E-03; + 28060.06 2.1605E-03; + 28061.06 9.3917E-05; + 28062.06 2.9945E-04; + 28064.06 7.6261E-05; + 14028.06 9.5281E-04; + 14029.06 4.8381E-05; + 14030.06 3.1893E-05; } + } + + Helium { + temp 566; + composition { + 2003.06 4.8089E-10; + 2004.06 2.4044E-04; } + } + + BorosilicateGlass { + temp 566; + composition { + 13027.06 1.7352E-03; + 5010.06 9.6506E-04; + 5011.06 3.9189E-03; + 8016.06 4.6514E-02; + 8017.06 1.7671E-05; + //8018.06 9.3268E-05; + 14028.06 1.6926E-02; + 14029.06 8.5944E-04; + 14030.06 5.6654E-04; } + } + + Water { + temp 566; + moder {1001.06 (lwj3.11 lwj3.09); } + composition { + 5010.06 7.9714E-06; + 5011.06 3.2247E-05; + 1001.06 4.9456E-02; + 1002.06 7.7035E-06; + 8016.06 2.4673E-02; + 8017.06 9.3734E-06; + //8018.06 4.9474E-05; + } + } + + Zircaloy { + temp 566; + composition { + 24050.06 3.2962E-06; + 24052.06 6.3564E-05; + 24053.06 7.2076E-06; + 24054.06 1.7941E-06; + 26054.06 8.6698E-06; + 26056.06 1.3610E-04; + 26057.06 3.1431E-06; + 26058.06 4.1829E-07; + 8016.06 3.0744E-04; + 8017.06 1.1680E-07; + //8018.03 6.1648E-07; + 50112.06 4.6735E-06; + 50114.06 3.1799E-06; + 50115.06 1.6381E-06; + 50116.06 7.0055E-05; + 50117.06 3.7003E-05; + 50118.06 1.1669E-04; + 50119.06 4.1387E-05; + 50120.06 1.5697E-04; + 50122.06 2.2308E-05; + 50124.06 2.7897E-05; + 40090.06 2.1828E-02; + 40091.06 4.7601E-03; + 40092.06 7.2759E-03; + 40094.06 7.3734E-03; + 40096.06 1.1879E-03; } + } + + Inconel{ + temp 566; + composition { + 24050.06 7.8239E-04; + 24052.06 1.5088E-02; + 24053.06 1.7108E-03; + 24054.06 4.2586E-04; + 26054.06 1.4797E-03; + 26056.06 2.3229E-02; + 26057.06 5.3645E-04; + 26058.06 7.1392E-05; + 25055.06 7.8201E-04; + 28058.06 2.9320E-02; + 28060.06 1.1294E-02; + 28061.06 4.9094E-04; + 28062.06 1.5653E-03; + 28064.06 3.9864E-04; + 14028.06 5.6757E-04; + 14029.06 2.8820E-05; + 14030.06 1.8998E-05; } + } + + B4C{ + temp 566; + composition { + 5010.06 1.5206E-02; + 5011.06 6.1514E-02; + 6012.06 1.8972E-02; + //6013.06 2.1252E-04; + } + } + + Ag-In-Cd{ + temp 566; + composition { + 47107.06 2.3523E-02; + 47109.06 2.1854E-02; + 48106.06 3.3882E-05; + 48108.06 2.4166E-05; + 48110.06 3.3936E-04; + 48111.06 3.4821E-04; + 48112.06 6.5611E-04; + 48113.06 3.3275E-04; + 48114.06 7.8252E-04; + 48116.06 2.0443E-04; + 49113.06 3.4219E-04; + 49115.06 7.6511E-03; } + } + + UO2-16 { + temp 566; + tms 1; + composition { + 8016.03 4.5897E-02; + 8017.03 1.7436E-05; + //8018.03 9.2032E-05; + 92234.03 3.0131E-06; + 92235.03 3.7503E-04; + 92238.03 2.2625E-02;} + } + + UO2-24 { + temp 566; + tms 1; + composition { + 8016.03 4.5830E-02; + 8017.03 1.7411E-05; + //8018.03 9.1898E-05; + 92234.03 4.4842E-06; + 92235.03 5.5814E-04; + 92238.03 2.2407E-02;} + } + + UO2-31 { + temp 566; + tms 1; + composition { + 8016.03 4.5853E-02; + 8017.03 1.7420E-05; + //8018.03 9.1942E-05; + 92234.03 5.7987E-06; + 92235.03 7.2175E-04; + 92238.03 2.2253E-02;} + } + + UO2-32 { + temp 566; + tms 1; + composition { + 8016.03 4.6029E-02; + 8017.03 1.7487E-05; + //8018.03 9.2296E-05; + 92234.03 5.9959E-06; + 92235.03 7.4630E-04; + 92238.03 2.2317E-02; + } + } + + UO2-34 { + temp 566; + tms 1; + composition { + 8016.03 4.6110E-02; + 8017.03 1.7517E-05; + //8018.03 9.2459E-05; + 92234.03 6.4018E-06; + 92235.03 7.9681E-04; + 92238.03 2.2307E-02;} + } + + // vanadium51 was stated twice in carbonsteel below + // in the beavrs pdf - typo? + CarbonSteel { + temp 566; + composition { + 13027.06 4.3523E-05; + 5010.06 2.5833E-06; + 5011.06 1.0450E-05; + 6012.06 1.0442E-03; + //6013.06 1.1697E-05 ; + 20040.06 1.7043E-05; + 20042.06 1.1375E-07; + 20043.06 2.3734E-08; + 20044.06 3.6673E-07; + 20046.06 7.0322E-10; + 20048.06 3.2875E-08; + 24050.06 1.3738E-05; + 24052.06 2.6493E-04; + 24053.06 3.0041E-05; + 24054.06 7.4778E-06; + 29063.06 1.0223E-04; + 29065.06 4.5608E-05; + 26054.06 4.7437E-03; + 26056.06 7.4465E-02; + 26057.06 1.7197E-03; + 26058.06 2.2886E-04; + 25055.06 6.4126E-04; + 42100.06 2.9814E-05; + 42092.06 4.4822E-05; + 42094.06 2.8110E-05; + 42095.06 4.8567E-05; + 42096.06 5.1015E-05; + 42097.06 2.9319E-05; + 42098.06 7.4327E-05; + 41093.06 5.0559E-06; + 28058.06 4.0862E-04; + 28060.06 1.5740E-04; + 28061.06 6.8420E-06; + 28062.06 2.1815E-05; + 28064.06 5.5557E-06; + 15031.06 3.7913E-05; + 16032.06 3.4808E-05; + 16033.06 2.7420E-07; + 16034.06 1.5368E-06; + 16036.06 5.3398E-09; + 14028.06 6.1702E-04; + 14029.06 3.1330E-05; + 14030.06 2.0653E-05; + 22046.06 1.2144E-06; + 22047.06 1.0952E-06; + 22048.06 1.0851E-05; + 22049.06 7.9634E-07; + 22050.06 7.6249E-07; + //23050.06 1.1526E-07; + 23051.06 4.5989E-05; + } + } + + SupportPlateSS { + temp 566; + composition { + 24050.06 3.5223E-04; + 24052.06 6.7924E-03; + 24053.06 7.7020E-04; + 24054.06 1.9172E-04; + 26054.06 1.5882E-03; + 26056.06 2.4931E-02; + 26057.06 5.7578E-04; + 26058.06 7.6625E-05; + 25055.06 8.0762E-04; + 28058.06 2.5731E-03; + 28060.06 9.9117E-04; + 28061.06 4.3085E-05; + 28062.06 1.3738E-04; + 28064.06 3.4985E-05; + 14028.06 4.3711E-04; + 14029.06 2.2195E-05; + 14030.06 1.4631E-05;} + } + + SupportPlateBW { + temp 566; + moder {1001.06 (lwj3.11 lwj3.09); } + composition { + 5010.06 1.0559E-05; + 5011.06 4.2716E-05; + 1001.06 6.5512E-02; + 1002.06 1.0204E-05; + 8016.06 3.2683E-02; + 8017.06 1.2416E-05; + //8018.06 6.5535E-05; + } + } + + +} +} From f20df29f3122c22dc243643225515851ea020656 Mon Sep 17 00:00:00 2001 From: ChasingNeutrons Date: Mon, 23 Jun 2025 01:17:39 +0100 Subject: [PATCH 3/6] Allow offset to be disabled in nested universes It is desirable for input to allow some universes not to be offset as might be done in a lattice. ALternatively, it may be preferable to evaluate their locations in the global reference frame. The former has been done by allowing lattice universes to have an 'offset map', the latter has been done by allowing a chosen universe to be evaluated in global coordinates. This allowed a fix to the BEAVRS inputs to correctly position the instrumentation tubes. There seems to be a dictionary error with the offsetMap. Occasionally it is not recognised as an integer list. This has prevented the creation of a new test in latUniverse. This has been commented out until the dictionary has been fixed. --- DataStructures/dictionary_class.f90 | 144 +++- Geometry/Surfaces/box_class.f90 | 4 +- Geometry/Surfaces/squareCylinder_class.f90 | 4 +- Geometry/Universes/Tests/latUniverse_test.f90 | 17 +- Geometry/Universes/latUniverse_class.f90 | 113 +++- Geometry/Universes/universe_inter.f90 | 45 +- Geometry/coord_class.f90 | 5 +- Geometry/geometryStd_class.f90 | 14 +- InputFiles/Benchmarks/BEAVRS/BEAVRS_ARO | 596 +++++++++-------- InputFiles/Benchmarks/BEAVRS/BEAVRS_HZP | 615 ++++++++++-------- SharedModules/genericProcedures.f90 | 42 +- SharedModules/universalVariables.f90 | 2 +- 12 files changed, 1019 insertions(+), 582 deletions(-) diff --git a/DataStructures/dictionary_class.f90 b/DataStructures/dictionary_class.f90 index 4434dda6a..e07fa7cf6 100644 --- a/DataStructures/dictionary_class.f90 +++ b/DataStructures/dictionary_class.f90 @@ -32,6 +32,8 @@ module dictionary_class ! -> 1D array of shortInt - integer(shortInt) :: i(:) ! -> scalar character of length "charLen" defined below - character(charLen) :: c ! -> array of characters of length "charLen" defined below - character(charLen) :: c(:) + ! -> boolean logical - logical(defBool) :: b + ! -> array of boolean logicals - logical(defBool) :: b(:) ! -> another dictionary - type(dictionary) :: dict ! ! To add additional structure to store it is necessary to : @@ -100,7 +102,7 @@ module dictionary_class ! *** Note that dictionary is defined as pointer not allocatable ! *** This is because gfortran < 7.0 does not supports circular derived types with ! *** allocatable keyword. This line may change in the future - type(dictionary), pointer :: dict0_alloc => null() + type(dictionary), pointer :: dict0_alloc => null() ! dictContent type ID integer(shortInt) :: type = empty @@ -149,7 +151,9 @@ module dictionary_class getCharArray_alloc_new,& getCharArray_ptr_new,& getDict_new,& - getBool_new + getBool_new, & + getBoolArray_alloc_new,& + getBoolArray_ptr_new procedure,private :: getReal_new procedure,private :: getRealArray_alloc_new @@ -162,6 +166,8 @@ module dictionary_class procedure,private :: getCharArray_ptr_new procedure,private :: getDict_new procedure,private :: getBool_new + procedure,private :: getBoolArray_alloc_new + procedure,private :: getBoolArray_ptr_new generic :: getOrDefault => getOrDefault_real ,& getOrDefault_realArray_alloc ,& @@ -172,7 +178,9 @@ module dictionary_class getOrDefault_char ,& getOrDefault_charArray_alloc ,& getOrDefault_charArray_ptr, & - getOrDefault_bool + getOrDefault_bool, & + getOrDefault_boolArray_alloc, & + getOrDefault_boolArray_ptr procedure,private :: getOrDefault_real procedure,private :: getOrDefault_realArray_alloc @@ -184,6 +192,8 @@ module dictionary_class procedure,private :: getOrDefault_charArray_alloc procedure,private :: getOrDefault_charArray_ptr procedure,private :: getOrDefault_bool + procedure,private :: getOrDefault_boolArray_alloc + procedure,private :: getOrDefault_boolArray_ptr ! Keys inquiry procedures procedure :: keys @@ -779,6 +789,83 @@ subroutine getBool_new(self,value,keyword) call fatalError(Here,'Entry under keyword ' // keyword // ' is not an integer') end select end subroutine getBool_new + + !! + !! Loads a boolean rank 1 from dictionary into provided variable + !! If keyword is associated with an integer which is not 1 or 0, it returns an error. + !! Variable needs to be allocatable. It will be deallocated before assignment + !! + subroutine getBoolArray_alloc_new(self,value,keyword) + class(dictionary), intent(in) :: self + logical(defBool),dimension(:),allocatable,intent(inout) :: value + character(*),intent(in) :: keyword + integer(shortInt) :: idx + integer(shortInt),dimension(:),allocatable :: tempInt + character(100),parameter :: Here='getBoolArray_alloc_new (dictionary_class.f90)' + + idx = self % search(keyword, Here, fatal =.true.) + + if(allocated(value)) deallocate(value) + + select case (self % entries(idx) % getType()) + case(arrInt) + tempInt = self % entries(idx) % int1_alloc + allocate(value(size(tempInt))) + + if (any(tempInt < 0) .or. any(tempInt > 1)) then + call fatalError(Here,'Entry under keyword ' // keyword // ' has non-logical values (not 0 or 1).') + else + do idx = 1, size(tempInt) + value(idx) = tempInt(idx) == 1 + end do + end if + + case default + call fatalError(Here,'Entry under keyword ' // keyword // ' is not an int array.') + + end select + + end subroutine getBoolArray_alloc_new + + !! + !! Loads a boolean rank 1 from dictionary into provided variable. + !! If keyword is associated with an integer which is not 1 or 0, it returns an error. + !! Variable needs to be pointer. It will be deallocated before assignment + !! + subroutine getBoolArray_ptr_new(self,value,keyword) + class(dictionary), intent(in) :: self + logical(defBool),dimension(:),pointer,intent(inout) :: value + character(*),intent(in) :: keyword + integer(shortInt) :: idx, N + integer(shortInt),dimension(:),allocatable :: tempInt + character(100),parameter :: Here='getBoolArray_ptr_new (dictionary_class.f90)' + + idx = self % search(keyword, Here, fatal =.true.) + + if(associated(value)) deallocate(value) + + select case (self % entries(idx) % getType()) + case(arrInt) + + N = size(self % entries(idx) % int1_alloc) + allocate(value(N)) + allocate(tempInt(N)) + tempInt = self % entries(idx) % int1_alloc + + if (any(tempInt < 0) .or. any(tempInt > 1)) then + call fatalError(Here,'Entry under keyword ' // keyword // ' has non-logical values (not 0 or 1).') + else + do idx = 1, N + value(idx) = tempInt(idx) == 1 + end do + end if + + case default + call fatalError(Here,'Entry under keyword ' // keyword // ' is not an int array.') + + end select + + end subroutine getBoolArray_ptr_new !! !! Loads a real rank 0 from a dictionary into provided variable @@ -1044,6 +1131,57 @@ subroutine getOrDefault_bool(self,value,keyword,default) end if end subroutine getOrDefault_bool + !! + !! Loads a boolean rank 1 from a dictionary. + !! Variable needs to be allocatable. It will be deallocated before assignment + !! + subroutine getOrDefault_boolArray_alloc(self, value, keyword, default) + class(dictionary), intent(in) :: self + logical(defBool),dimension(:),allocatable,intent(inout) :: value + character(*),intent(in) :: keyword + logical(defBool),dimension(:),intent(in) :: default + integer(shortInt) :: idx + character(100),parameter :: Here='getOrDefault_boolArray_alloc (dictionary_class.f90)' + + idx = self % search(keyword, Here, fatal =.false.) + + if(allocated(value)) deallocate(value) + + if (idx == targetNotFound) then + value = default + else + call self % get(value, keyword) + end if + end subroutine getOrDefault_boolArray_alloc + + + !! + !! Loads a boolean rank 1 from a dictionary. + !! Variable needs to be pointer. It will be deallocated before assignment + !! + !! For further details refer to doc of getOrDefault_real + !! + subroutine getOrDefault_boolArray_ptr(self, value, keyword, default) + class(dictionary), intent(in) :: self + logical(defBool),dimension(:),pointer,intent(inout) :: value + character(*),intent(in) :: keyword + logical(defBool),dimension(:),intent(in) :: default + integer(shortInt) :: idx + character(100),parameter :: Here='getOrDefault_boolArray_ptr (dictionary_class.f90)' + + idx = self % search(keyword, Here, fatal =.false.) + + if(associated(value)) deallocate(value) + + if (idx == targetNotFound) then + allocate(value(size(default))) + value = default + else + call self % get(value, keyword) + end if + + end subroutine getOrDefault_boolArray_ptr + !! !! Returns an array of all keywords !! diff --git a/Geometry/Surfaces/box_class.f90 b/Geometry/Surfaces/box_class.f90 index 221903c9b..21e7f0d6a 100644 --- a/Geometry/Surfaces/box_class.f90 +++ b/Geometry/Surfaces/box_class.f90 @@ -347,7 +347,7 @@ subroutine explicitBC(self, r, u) r0 = r(ax) - self % origin(ax) ! Skip if particle is well inside the domain - if (abs(r0) <= self % halfwidth(ax) - self % surfTol()) cycle axis + if (abs(r0) <= self % halfwidth(ax) * (ONE - self % surfTol())) cycle axis ! Choose correct BC if (r0 < ZERO) then @@ -396,7 +396,7 @@ subroutine transformBC(self, r, u) ! Calculate halfwidth reduced by the surface_tolerance ! Necessary to capture particles at the boundary - a_bar = self % halfwidth - self % surfTol() + a_bar = self % halfwidth * (ONE - self % surfTol()) ! Calculate distance (in # of transformations) in each direction Ri = ceiling(abs(r - self % origin) / a_bar) / 2 diff --git a/Geometry/Surfaces/squareCylinder_class.f90 b/Geometry/Surfaces/squareCylinder_class.f90 index 5e41d95ec..196c17027 100644 --- a/Geometry/Surfaces/squareCylinder_class.f90 +++ b/Geometry/Surfaces/squareCylinder_class.f90 @@ -402,7 +402,7 @@ subroutine explicitBC(self, r, u) r0 = r(ax) - self % origin(i) ! Skip if particle is well inside the domain - if (abs(r0) <= self % halfwidth(i) - self % surfTol()) cycle axis + if (abs(r0) <= self % halfwidth(i) * (ONE- self % surfTol())) cycle axis ! Choose correct BC if (r0 < ZERO) then @@ -452,7 +452,7 @@ subroutine transformBC(self, r, u) ! Calculate halfwidth reduced by the surface_tolerance ! Necessary to capture particles at the boundary - a_bar = self % halfwidth - self % surfTol() + a_bar = self % halfwidth * (ONE - self % surfTol()) ! Calculate distance (in # of transformations) in each direction Ri = ceiling(abs(r(self % plane) - self % origin) / a_bar) / 2 diff --git a/Geometry/Universes/Tests/latUniverse_test.f90 b/Geometry/Universes/Tests/latUniverse_test.f90 index 6384dac7f..99b005133 100644 --- a/Geometry/Universes/Tests/latUniverse_test.f90 +++ b/Geometry/Universes/Tests/latUniverse_test.f90 @@ -18,11 +18,15 @@ module latUniverse_test character(*), parameter :: UNI1_DEF = & "id 1; type latUniverse; origin (0.0 0.0 0.0); rotation (0.0 0.0 0.0); & &pitch (1.0 2.0 3.0); shape (3 2 2); padMat void; & - &map ( 3 4 5 & + & map ( 3 4 5 & & 7 4 8 & & & & 1 2 3 & & 4 5 6); " + !& offsetMap ( 1 1 1 & + !& 0 1 1 & + !& 1 1 1 & + !& 1 1 1 ); & character(*), parameter :: UNI2_DEF = & "id 2; type latUniverse; pitch (1.0 2.0 0.0); shape (2 1 0); padMat u<1>; & @@ -389,6 +393,17 @@ subroutine test_cellOffset() ref = [0.0_defReal, 1.0_defReal, 1.5_defReal] @assertEqual(ref, uni1 % cellOffset(pos), TOL) + + !! Inside but at the position without an offset + !pos % r = [-1.5_defReal, 1.0_defReal, -0.5_defReal] + !pos % dir = [-ONE, ONE, -ONE] + !pos % dir = pos % dir / norm2(pos % dir) + !pos % uniIdx = 8 + !pos % cellIdx = 0 + !pos % localId = 4 + + !ref = [0.0_defReal, 0.0_defReal, 0.0_defReal] + !@assertEqual(ref, uni1 % cellOffset(pos), TOL) ! Outside pos % r = [-7.0_defReal, 0.0_defReal, 0.5_defReal] diff --git a/Geometry/Universes/latUniverse_class.f90 b/Geometry/Universes/latUniverse_class.f90 index ca950e9cd..5dea8b0a1 100644 --- a/Geometry/Universes/latUniverse_class.f90 +++ b/Geometry/Universes/latUniverse_class.f90 @@ -20,6 +20,9 @@ module latUniverse_class integer(shortInt), parameter :: X_MIN = -1, X_MAX = -2, Y_MIN = -3, Y_MAX = -4, Z_MIN = -5, & Z_MAX = -6, OUTLINE_SURF = -7 + ! Options for the offset map + integer(shortInt), parameter :: local = 1, noOffset = 0 + !! !! 2D or 3D Cartesian lattice with constant pitch !! @@ -32,7 +35,9 @@ module latUniverse_class !! Background cell can have any filling given by keyword (material or universe) !! !! Every lattice cell has an offset to its centre (so the centre of the nested universe - !! is in the center of the lattice cell). + !! is in the center of the lattice cell). Optionally an offset map can be provided, determining + !! whether to apply an offset in a given cell position. This can disables the local universe + !! offset. Alternatively a single offset flag can be provided, disabling offset in all cells. !! !! Minimum lattice pitch is set to 10 * SURF_TOL !! @@ -47,7 +52,12 @@ module latUniverse_class !! map ( 1 2 3 // Top layer !! 4 5 6 // Lower Y row !! 7 8 9 // Bottom layer - !! 10 11 12 ) + !! 10 11 12 ); + !! #offsetMap ( 1 1 1 + !! # 1 0 1 + !! # 1 0 1 + !! # 1 1 1 ); + !! # offset 1; !! } !! !! Sample Input Dictionary (2D): @@ -70,18 +80,22 @@ module latUniverse_class !! a_bar -> Halfwidth of lattice cell reduced by surface tolerance !! outline -> Box type surface that is a boundary between lattice & background !! outLocalID -> LocalID of the background cell + !! offset -> Flag to disable all offsets + !! offsetMap -> Map determining which cells have a lattice offset or not !! !! Interface: !! universe interface !! type, public, extends(universe) :: latUniverse private - real(defReal), dimension(3) :: pitch = ZERO - integer(shortInt), dimension(3) :: sizeN = 0 - real(defReal), dimension(3) :: corner = ZERO - real(defReal), dimension(3) :: a_bar = ZERO - type(box) :: outline - integer(shortInt) :: outLocalID = 0 + real(defReal), dimension(3) :: pitch = ZERO + integer(shortInt), dimension(3) :: sizeN = 0 + real(defReal), dimension(3) :: corner = ZERO + real(defReal), dimension(3) :: a_bar = ZERO + type(box) :: outline + integer(shortInt) :: outLocalID = 0 + logical(defBool) :: offset = .true. + integer(shortInt), dimension(:), allocatable :: offsetMap contains ! Superclass procedures procedure :: init @@ -111,7 +125,7 @@ subroutine init(self, fill, dict, cells, surfs, mats) type(charMap), intent(in) :: mats real(defReal), dimension(:), allocatable :: temp integer(shortInt), dimension(:), allocatable :: tempI - integer(shortInt) :: N, i, j, outFill + integer(shortInt) :: N, i, j, outFill, val type(dictionary) :: tempDict integer(shortInt), dimension(:,:), allocatable :: tempMap character(nameLen) :: name @@ -120,6 +134,9 @@ subroutine init(self, fill, dict, cells, surfs, mats) ! Setup the base class ! With: id, origin rotations... call self % setupBase(dict) + + ! Perform offsets on every cell? + call dict % getOrDefault(self % offset, 'offset', .true.) ! Load pitch call dict % get(temp, 'pitch') @@ -152,11 +169,11 @@ subroutine init(self, fill, dict, cells, surfs, mats) ! Check for invalid pitch if (any(self % pitch < 10 * SURF_TOL)) then - call fatalError(Here, 'Pitch size must be larger than: '//numToChar( 10 * SURF_TOL)) - end if + call fatalError(Here, 'Pitch size must be larger than: '//numToChar( 10 * SURF_TOL)) + end if ! Calculate halfwidth and corner - self % a_bar = self % pitch * HALF - SURF_TOL + self % a_bar = self % pitch * HALF - self % pitch * SURF_TOL self % corner = -(self % sizeN * HALF * self % pitch) ! Calculate local ID of the background @@ -173,6 +190,11 @@ subroutine init(self, fill, dict, cells, surfs, mats) ! Construct fill array call dict % get(tempI, 'map') + ! Ensure size matches sizeN + if (size(tempI) /= product(self % sizeN)) call fatalError(Here, & + 'Lattice map size not equal to size implied by shape. Respectively: '//& + numToChar(size(tempI))//' '//numToChar(product(self % sizeN))) + ! Flip array up-down for more natural input ! Reshape into rank 2 array tempMap = reshape(tempI, [self % sizeN(1), self % sizeN(2) * self % sizeN(3)]) @@ -195,6 +217,46 @@ subroutine init(self, fill, dict, cells, surfs, mats) end do end do fill(self % outLocalID) = outFill + deallocate(tempI, tempMap) + + ! Check whether there is an offset map + if (dict % isPresent('offsetMap')) then + + if (.not. self % offset) call fatalError(Here, 'Cannot have both an offset map '//& + 'and no offset.') + + call dict % get(tempI, 'offsetMap') + + ! Ensure size matches sizeN + if (size(tempI) /= product(self % sizeN)) call fatalError(Here, & + 'Offset map size not equal to size implied by shape. Respectively: '//& + numToChar(size(tempI))//' '//numToChar(product(self % sizeN))) + + ! Flip array up-down for more natural input + ! Reshape into rank 2 array + tempMap = reshape(tempI, [self % sizeN(1), self % sizeN(2) * self % sizeN(3)]) + N = size(tempMap, 2) + do i = 1, N/2 + call swap(tempMap(:,i), tempMap(:,N - i + 1)) + end do + + allocate(self % offsetMap(product(self % sizeN) + 1)) + N = size(tempMap, 1) + do j = 1, size(tempMap, 2) + do i = 1, N + val = tempMap(i,j) + self % offsetMap(i + (j-1) * N) = val + + ! Check that the entries are valid + if (val /= local .and. val /= noOffset) call fatalError(Here,& + 'Invalid entry to the offset map. Must be one of '//numToChar(local)//& + ' '//numToChar(noOffset)//'. Contains: '//numToChar(val)) + end do + end do + ! Add an entry for the padMat + self % offsetMap(self % outLocalID) = noOffset + + end if end subroutine init @@ -254,12 +316,12 @@ end subroutine findCell !! subroutine distance(self, d, surfIdx, coords) class(latUniverse), intent(inout) :: self - real(defReal), intent(out) :: d - integer(shortInt), intent(out) :: surfIdx - type(coord), intent(in) :: coords - real(defReal), dimension(3) :: r_bar, u, bounds - real(defReal) :: test_d - integer(shortInt) :: i, ax + real(defReal), intent(out) :: d + integer(shortInt), intent(out) :: surfIdx + type(coord), intent(in) :: coords + real(defReal), dimension(3) :: r_bar, u, bounds + real(defReal) :: test_d + integer(shortInt) :: i, ax ! Catch case if particle is outside the lattice if (coords % localID == self % outLocalID) then @@ -331,13 +393,20 @@ function cellOffset(self, coords) result (offset) class(latUniverse), intent(in) :: self type(coord), intent(in) :: coords real(defReal), dimension(3) :: offset + logical(defBool) :: doOffset - if (coords % localID == self % outLocalID) then - offset = ZERO - + if (allocated(self % offsetMap)) then + doOffset = self % offsetMap(coords % localID) == noOffset else + doOffset = self % offset + end if + + if (doOffset .and. coords % localID /= self % outLocalID) then offset = (get_ijk(coords % localID, self % sizeN) - HALF) * self % pitch + self % corner + else + offset = ZERO + end if end function cellOffset @@ -358,6 +427,8 @@ elemental subroutine kill(self) self % a_bar = ZERO call self % outline % kill() self % outLocalID = 0 + if (allocated(self % offsetMap)) deallocate(self % offsetMap) + self % offset = .true. end subroutine kill diff --git a/Geometry/Universes/universe_inter.f90 b/Geometry/Universes/universe_inter.f90 index c903a50bd..47afcca57 100644 --- a/Geometry/Universes/universe_inter.f90 +++ b/Geometry/Universes/universe_inter.f90 @@ -38,11 +38,12 @@ module universe_inter !! } !! !! Private Members: - !! uniId -> Id of the universe - !! uniIdx -> Index of the universe - !! origin -> Location of the origin of the universe co-ordinates in the frame of higher universe - !! rotMat -> Rotation matrix for rotation with respect to the higher universe - !! rot -> rotation flag. True is universe is rotated + !! uniId -> Id of the universe + !! uniIdx -> Index of the universe + !! origin -> Location of the origin of the universe co-ordinates in the frame of higher universe + !! rotMat -> Rotation matrix for rotation with respect to the higher universe + !! rot -> rotation flag. True if universe is rotated + !! globalTrans -> global transformation flag. True if universe coordinates should be evaluated in global frame !! !! Interface: !! id -> Get Id of the universe @@ -60,11 +61,12 @@ module universe_inter !! type, public, abstract :: universe private - integer(shortInt) :: uniId = 0 + integer(shortInt), public :: uniId = 0 integer(shortInt) :: uniIdx = 0 real(defReal), dimension(3) :: origin = ZERO real(defReal), dimension(3,3) :: rotMat = ZERO logical(defBool) :: rot = .false. + logical(defBool) :: globalTrans = .false. contains ! Build procedures procedure, non_overridable :: id @@ -81,6 +83,7 @@ module universe_inter procedure(distance), deferred :: distance procedure(cross), deferred :: cross procedure(cellOffset), deferred :: cellOffset + procedure :: transformToGlobal end type universe abstract interface @@ -316,6 +319,11 @@ subroutine setupBase(self, dict) call self % setTransform(rotation=temp) end if + ! Load global translation + if (dict % isPresent('global')) then + call dict % get(self % globalTrans, 'global') + end if + end subroutine setupBase !! @@ -401,13 +409,28 @@ end subroutine enter elemental subroutine kill(self) class(universe), intent(inout) :: self - self % uniIdx = 0 - self % uniId = 0 - self % origin = ZERO - self % rotMat = ZERO - self % rot = .false. + self % uniIdx = 0 + self % uniId = 0 + self % origin = ZERO + self % rotMat = ZERO + self % rot = .false. + self % globalTrans = .false. end subroutine kill + + !! + !! Should co-ordinates be transformed to the global frame? + !! + !! Args: + !! doTrans -> Bool deciding whether to transform the co-ordinates + !! + function transformToGlobal(self) result(doTrans) + class(universe), intent(in) :: self + logical(defBool) :: doTrans + + doTrans = self % globalTrans + + end function transformToGlobal !!<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><> !! Utility Functions diff --git a/Geometry/coord_class.f90 b/Geometry/coord_class.f90 index d7900d9a7..5bc9b2daa 100644 --- a/Geometry/coord_class.f90 +++ b/Geometry/coord_class.f90 @@ -18,6 +18,7 @@ module coord_class !! r -> Position !! dir -> Direction !! isRotated -> Is rotated wrt previous (higher by 1) level + !! isGlobal -> Is the coordinate reverted to level 1? !! rotMat -> Rotation matrix wrt previous level !! uniIdx -> Index of the occupied universe !! uniRootID -> Location of the occupied universe in geometry graph @@ -33,6 +34,7 @@ module coord_class real(defReal), dimension(3) :: r = ZERO real(defReal), dimension(3) :: dir = ZERO logical(defBool) :: isRotated = .false. + logical(defBool) :: isGlobal = .false. real(defReal), dimension(3,3) :: rotMat = ZERO integer(shortInt) :: uniIdx = 0 integer(shortInt) :: uniRootID = 0 @@ -45,7 +47,7 @@ module coord_class end type coord !! - !! List of co-ordinates at diffrent level of a geometry + !! List of co-ordinates at different level of a geometry !! !! Specifies the position of a particle in space !! @@ -160,6 +162,7 @@ elemental subroutine kill_coord(self) self % r = ZERO self % dir = ZERO self % isRotated = .false. + self % isGlobal = .false. self % rotMat = ZERO self % uniIdx = 0 self % uniRootID = 0 diff --git a/Geometry/geometryStd_class.f90 b/Geometry/geometryStd_class.f90 index d53e33020..8a10311b7 100644 --- a/Geometry/geometryStd_class.f90 +++ b/Geometry/geometryStd_class.f90 @@ -438,7 +438,7 @@ subroutine diveToMat(self, coords, start) integer(shortInt), intent(in) :: start integer(shortInt) :: rootID, localID, fill, id, i class(universe), pointer :: uni - real(defReal), dimension(3) :: offset + real(defReal), dimension(3) :: offset, r character(100), parameter :: Here = 'diveToMat (geometryStd_class.f90)' do i = start, HARDCODED_MAX_NEST @@ -462,13 +462,21 @@ subroutine diveToMat(self, coords, start) ! Get cell offset offset = uni % cellOffset(coords % lvl(i)) - + ! Get nested universe uni => self % geom % unis % getPtr_fast(fill) + + ! Does this level revert to the global frame? + ! Unrotate direction as well? + if (uni % transformToGlobal()) then + r = coords % lvl(1) % r + else + r = coords % lvl(i) % r - offset + end if ! Enter nested univers call coords % addLevel() - call uni % enter(coords % lvl(i+1), coords % lvl(i) % r - offset, coords % lvl(i) % dir) + call uni % enter(coords % lvl(i+1), r, coords % lvl(i) % dir) coords % lvl(i+1) % uniRootID = id ! Must be after enter where coord has intent out end if diff --git a/InputFiles/Benchmarks/BEAVRS/BEAVRS_ARO b/InputFiles/Benchmarks/BEAVRS/BEAVRS_ARO index dd3ed25c5..89434cfe9 100644 --- a/InputFiles/Benchmarks/BEAVRS/BEAVRS_ARO +++ b/InputFiles/Benchmarks/BEAVRS/BEAVRS_ARO @@ -1,14 +1,13 @@ !! !! 3D BEAVRS benchmark -!! This is not a fully faithful replica of BEAVRS at HZP: !! all rods are fully withdrawn. -!! In reality, some rods are partially inserted. +!! This could definitely be better written !! type eigenPhysicsPackage; -pop 1000000; +pop 10000000; active 50; -inactive 150; +inactive 250; XSdata ce; dataType ce; @@ -471,6 +470,9 @@ geometry { GT35 {type simpleCell; id 373; surfaces (-131 132); filltype uni; universe 1005;} GT20 {type simpleCell; id 374; surfaces (-132); filltype uni; universe 1001;} + // Extra GT cells without dashpot for empty instrumentation tubes + CGT98 {type simpleCell; id 379; surfaces (-123 126); filltype uni; universe 9112;} + // control rod, axial layering // Used (probably with some modification) only when fully inserted //CR460 {type simpleCell; id 448; surfaces (-120 121); filltype uni; universe 1002;} @@ -752,9 +754,15 @@ geometry { //instr. tube instrumentTube { - id 14; + id 15; type cellUniverse; cells (198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217);} + + //empty instr. tube + emptyInstrumentTube { + id 11; + type cellUniverse; + cells (353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 379 370 371 372 373 374);} //burnable absorber BP { @@ -763,6 +771,33 @@ geometry { cells (218 219 220 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252);} + // Lattice of instrumentation tubes + instLattice { + type latUniverse; + id 14; + origin (0 0 0); + pitch (21.50364 21.50364 0.0); + shape (15 15 0); + global 1; + padMat Water; + map ( + 11 11 11 11 11 11 15 11 11 15 11 11 11 11 11 + 11 11 15 11 11 15 11 15 11 11 11 11 11 11 11 + 11 11 11 11 11 11 11 15 11 15 11 15 11 15 11 + 11 15 15 11 11 11 11 15 11 11 11 11 11 11 11 + 11 11 11 11 15 11 11 11 15 11 15 11 15 11 11 + 15 11 15 11 11 15 11 15 11 11 11 11 11 15 11 + 11 11 11 15 11 11 15 11 11 15 11 11 15 11 11 + 15 11 15 11 15 11 15 11 11 15 11 15 15 15 11 + 11 15 11 11 11 11 11 11 15 11 15 11 11 11 15 + 11 11 11 11 15 11 15 11 11 11 11 15 11 11 11 + 15 11 11 11 15 11 11 15 11 11 15 11 11 11 15 + 11 11 11 11 11 15 11 11 15 11 11 15 11 11 11 + 11 11 15 11 15 11 11 15 11 11 11 11 11 15 11 + 11 11 15 11 11 11 15 11 11 15 11 15 11 11 11 + 11 11 11 11 15 11 11 15 11 11 11 11 11 11 11 + ); + } // Lattices w/o grid @@ -791,7 +826,8 @@ geometry { 24 24 24 12 24 24 24 24 24 24 24 24 24 12 24 24 24 24 24 24 24 24 12 24 24 12 24 24 12 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 - 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); } + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); + } // assembly with sleeves at different heights A0E24Sleeve { @@ -826,7 +862,8 @@ geometry { 16 16 16 12 16 16 16 16 16 16 16 16 16 12 16 16 16 16 16 16 16 16 12 16 16 12 16 16 12 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 - 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16); } + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16); + } A0E16Sleeve { id 14160; @@ -861,7 +898,9 @@ geometry { 31 31 31 12 31 31 31 31 31 31 31 31 31 12 31 31 31 31 31 31 31 31 12 31 31 12 31 31 12 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A0E31Sleeve { id 14310; @@ -897,7 +936,9 @@ geometry { 31 31 31 12 31 31 31 31 31 31 31 31 31 12 31 31 31 31 31 31 31 31 12 31 31 12 31 31 12 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A6BE31BSleeve { @@ -933,7 +974,9 @@ geometry { 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 31 31 31 31 31 10 31 31 12 31 31 10 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A6BE31TSleeve { id 6031120; @@ -968,7 +1011,9 @@ geometry { 31 31 31 10 31 31 31 31 31 31 31 31 31 12 31 31 31 31 31 31 31 31 10 31 31 12 31 31 12 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A6BE31RSleeve { id 603130; @@ -1003,7 +1048,9 @@ geometry { 31 31 31 12 31 31 31 31 31 31 31 31 31 10 31 31 31 31 31 31 31 31 12 31 31 12 31 31 10 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A6BE31LSleeve { id 603190; @@ -1038,7 +1085,9 @@ geometry { 31 31 31 12 31 31 31 31 31 31 31 31 31 12 31 31 31 31 31 31 31 31 12 31 31 12 31 31 12 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A15BE31BRSleeve { id 153150; @@ -1073,7 +1122,9 @@ geometry { 31 31 31 12 31 31 31 31 31 31 31 31 31 12 31 31 31 31 31 31 31 31 12 31 31 12 31 31 12 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A15BE31BLSleeve { id 153170; @@ -1108,7 +1159,9 @@ geometry { 31 31 31 10 31 31 31 31 31 31 31 31 31 12 31 31 31 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A15BE31TRSleeve { id 153110; @@ -1143,7 +1196,9 @@ geometry { 31 31 31 12 31 31 31 31 31 31 31 31 31 10 31 31 31 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A15BE31TLSleeve { id 1531110; @@ -1178,7 +1233,9 @@ geometry { 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A16BE31Sleeve { id 16310; @@ -1213,7 +1270,9 @@ geometry { 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A20BE31Sleeve { id 20310; @@ -1248,7 +1307,9 @@ geometry { 24 24 24 10 24 24 24 24 24 24 24 24 24 10 24 24 24 24 24 24 24 24 10 24 24 12 24 24 10 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 - 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); } + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); + + } A12BE24Sleeve { id 12240; @@ -1283,7 +1344,9 @@ geometry { 24 24 24 10 24 24 24 24 24 24 24 24 24 10 24 24 24 24 24 24 24 24 10 24 24 10 24 24 10 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 - 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); } + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); + + } A16BE24Sleeve { id 16240; @@ -1319,7 +1382,9 @@ geometry { 24 24 24 13 24 24 24 24 24 24 24 24 24 13 24 24 24 24 24 24 24 24 13 24 24 13 24 24 13 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 - 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); } + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); + + } A0E24USleeve { id 24240; @@ -1353,7 +1418,8 @@ geometry { 16 16 16 13 16 16 16 16 16 16 16 16 16 13 16 16 16 16 16 16 16 16 13 16 16 13 16 16 13 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 - 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16); } + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16); + } // sleeved A0E16USleeve { @@ -1389,7 +1455,9 @@ geometry { 31 31 31 13 31 31 31 31 31 31 31 31 31 13 31 31 31 31 31 31 31 31 13 31 31 13 31 31 13 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A0E31USleeve { id 24310; @@ -1425,7 +1493,9 @@ geometry { 31 31 31 13 31 31 31 31 31 31 31 31 31 13 31 31 31 31 31 31 31 31 13 31 31 13 31 31 13 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A6BE31BUSleeve { @@ -1461,7 +1531,9 @@ geometry { 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 31 31 31 31 31 10 31 31 13 31 31 10 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A6BE31TUSleeve { id 7031120; @@ -1496,7 +1568,9 @@ geometry { 31 31 31 10 31 31 31 31 31 31 31 31 31 13 31 31 31 31 31 31 31 31 10 31 31 13 31 31 13 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A6BE31RUSleeve { id 703130; @@ -1531,7 +1605,9 @@ geometry { 31 31 31 13 31 31 31 31 31 31 31 31 31 10 31 31 31 31 31 31 31 31 13 31 31 13 31 31 10 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A6BE31LUSleeve { id 703190; @@ -1566,7 +1642,9 @@ geometry { 31 31 31 13 31 31 31 31 31 31 31 31 31 13 31 31 31 31 31 31 31 31 13 31 31 13 31 31 13 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A15BE31BRUSleeve { id 253150; @@ -1601,7 +1679,9 @@ geometry { 31 31 31 13 31 31 31 31 31 31 31 31 31 13 31 31 31 31 31 31 31 31 13 31 31 13 31 31 13 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A15BE31BLUSleeve { id 253170; @@ -1636,7 +1716,10 @@ geometry { 31 31 31 10 31 31 31 31 31 31 31 31 31 13 31 31 31 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + + } A15BE31TRUSleeve { id 253110; @@ -1671,7 +1754,9 @@ geometry { 31 31 31 13 31 31 31 31 31 31 31 31 31 10 31 31 31 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A15BE31TLUSleeve { id 2531110; @@ -1706,7 +1791,9 @@ geometry { 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A16BE31USleeve { id 26310; @@ -1741,7 +1828,9 @@ geometry { 31 31 31 10 31 31 31 31 31 31 31 31 31 10 31 31 31 31 31 31 31 31 10 31 31 10 31 31 10 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 - 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); } + 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 ); + + } A20BE31USleeve { id 30310; @@ -1776,7 +1865,9 @@ geometry { 24 24 24 10 24 24 24 24 24 24 24 24 24 10 24 24 24 24 24 24 24 24 10 24 24 13 24 24 10 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 - 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); } + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); + + } A12BE24USleeve { id 22240; @@ -1811,7 +1902,8 @@ geometry { 24 24 24 10 24 24 24 24 24 24 24 24 24 10 24 24 24 24 24 24 24 24 10 24 24 10 24 24 10 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 - 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); } + 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); + } A16BE24USleeve { id 26240; @@ -1932,7 +2024,7 @@ viz { output imgZ; what material; centre (0.0 0.0 232.0); - //width (400.0 400.0); + width (100.0 100.0); axis z; res (4000 4000); } @@ -1945,222 +2037,219 @@ nuclearData { } materials { - - // Note that commented nuclide densities are included in the specification - // but are not available in the JEFF-3.11 data library - - Air { + + Air { temp 566; composition { - 18036.06 7.8730E-09; - 18038.06 1.4844E-09; - 18040.06 2.3506E-06; - 6012.06 6.7539E-08; - //6013.06 7.5658E-10; - 7014.06 1.9680E-04; - 7015.06 7.2354E-07; - 8016.06 5.2866E-05; - 8017.06 2.0084E-08; - //8018.06 1.0601E-07; - } + 18036.01 7.8730E-09; + 18038.01 1.4844E-09; + 18040.01 2.3506E-06; + 6012.01 6.7539E-08; + 6013.01 7.5658E-10; + 7014.01 1.9680E-04; + 7015.01 7.2354E-07; + 8016.01 5.2866E-05; + 8017.01 2.0084E-08; + 8018.01 1.0601E-07; } + } SS304 { temp 566; composition { - 24050.06 7.6778E-04; - 24052.06 1.4806E-02; - 24053.06 1.6789E-03; - 24054.06 4.1791E-04; - 26054.06 3.4620E-03; - 26056.06 5.4345E-02; - 26057.06 1.2551E-03; - 26058.06 1.6703E-04; - 25055.06 1.7604E-03; - 28058.06 5.6089E-03; - 28060.06 2.1605E-03; - 28061.06 9.3917E-05; - 28062.06 2.9945E-04; - 28064.06 7.6261E-05; - 14028.06 9.5281E-04; - 14029.06 4.8381E-05; - 14030.06 3.1893E-05; } + 24050.01 7.6778E-04; + 24052.01 1.4806E-02; + 24053.01 1.6789E-03; + 24054.01 4.1791E-04; + 26054.01 3.4620E-03; + 26056.01 5.4345E-02; + 26057.01 1.2551E-03; + 26058.01 1.6703E-04; + 25055.01 1.7604E-03; + 28058.01 5.6089E-03; + 28060.01 2.1605E-03; + 28061.01 9.3917E-05; + 28062.01 2.9945E-04; + 28064.01 7.6261E-05; + 14028.01 9.5281E-04; + 14029.01 4.8381E-05; + 14030.01 3.1893E-05; } } Helium { temp 566; composition { - 2003.06 4.8089E-10; - 2004.06 2.4044E-04; } + 2003.01 4.8089E-10; + 2004.01 2.4044E-04; } } BorosilicateGlass { temp 566; composition { - 13027.06 1.7352E-03; - 5010.06 9.6506E-04; - 5011.06 3.9189E-03; - 8016.06 4.6514E-02; - 8017.06 1.7671E-05; - //8018.06 9.3268E-05; - 14028.06 1.6926E-02; - 14029.06 8.5944E-04; - 14030.06 5.6654E-04; } + 13027.01 1.7352E-03; + 5010.01 9.6506E-04; + 5011.01 3.9189E-03; + 8016.01 4.6514E-02; + 8017.01 1.7671E-05; + 8018.01 9.3268E-05; + 14028.01 1.6926E-02; + 14029.01 8.5944E-04; + 14030.01 5.6654E-04; } } Water { temp 566; - moder {1001.06 (lwj3.11 lwj3.09); } + moder {1001.01 (h-h2o.52 h-h2o.53); } composition { - 5010.06 7.9714E-06; - 5011.06 3.2247E-05; - 1001.06 4.9456E-02; - 1002.06 7.7035E-06; - 8016.06 2.4673E-02; - 8017.06 9.3734E-06; - //8018.06 4.9474E-05; + 5010.01 7.9714E-06; + 5011.01 3.2247E-05; + 1001.01 4.9456E-02; + 1002.01 7.7035E-06; + 8016.01 2.4673E-02; + 8017.01 9.3734E-06; + 8018.01 4.9474E-05; } } Zircaloy { temp 566; composition { - 24050.06 3.2962E-06; - 24052.06 6.3564E-05; - 24053.06 7.2076E-06; - 24054.06 1.7941E-06; - 26054.06 8.6698E-06; - 26056.06 1.3610E-04; - 26057.06 3.1431E-06; - 26058.06 4.1829E-07; - 8016.06 3.0744E-04; - 8017.06 1.1680E-07; - //8018.03 6.1648E-07; - 50112.06 4.6735E-06; - 50114.06 3.1799E-06; - 50115.06 1.6381E-06; - 50116.06 7.0055E-05; - 50117.06 3.7003E-05; - 50118.06 1.1669E-04; - 50119.06 4.1387E-05; - 50120.06 1.5697E-04; - 50122.06 2.2308E-05; - 50124.06 2.7897E-05; - 40090.06 2.1828E-02; - 40091.06 4.7601E-03; - 40092.06 7.2759E-03; - 40094.06 7.3734E-03; - 40096.06 1.1879E-03; } + 24050.01 3.2962E-06; + 24052.01 6.3564E-05; + 24053.01 7.2076E-06; + 24054.01 1.7941E-06; + 26054.01 8.6698E-06; + 26056.01 1.3610E-04; + 26057.01 3.1431E-06; + 26058.01 4.1829E-07; + 8016.01 3.0744E-04; + 8017.01 1.1680E-07; + 8018.01 6.1648E-07; + 50112.01 4.6735E-06; + 50114.01 3.1799E-06; + 50115.01 1.6381E-06; + 50116.01 7.0055E-05; + 50117.01 3.7003E-05; + 50118.01 1.1669E-04; + 50119.01 4.1387E-05; + 50120.01 1.5697E-04; + 50122.01 2.2308E-05; + 50124.01 2.7897E-05; + 40090.01 2.1828E-02; + 40091.01 4.7601E-03; + 40092.01 7.2759E-03; + 40094.01 7.3734E-03; + 40096.01 1.1879E-03; } } Inconel{ temp 566; composition { - 24050.06 7.8239E-04; - 24052.06 1.5088E-02; - 24053.06 1.7108E-03; - 24054.06 4.2586E-04; - 26054.06 1.4797E-03; - 26056.06 2.3229E-02; - 26057.06 5.3645E-04; - 26058.06 7.1392E-05; - 25055.06 7.8201E-04; - 28058.06 2.9320E-02; - 28060.06 1.1294E-02; - 28061.06 4.9094E-04; - 28062.06 1.5653E-03; - 28064.06 3.9864E-04; - 14028.06 5.6757E-04; - 14029.06 2.8820E-05; - 14030.06 1.8998E-05; } + 24050.01 7.8239E-04; + 24052.01 1.5088E-02; + 24053.01 1.7108E-03; + 24054.01 4.2586E-04; + 26054.01 1.4797E-03; + 26056.01 2.3229E-02; + 26057.01 5.3645E-04; + 26058.01 7.1392E-05; + 25055.01 7.8201E-04; + 28058.01 2.9320E-02; + 28060.01 1.1294E-02; + 28061.01 4.9094E-04; + 28062.01 1.5653E-03; + 28064.01 3.9864E-04; + 14028.01 5.6757E-04; + 14029.01 2.8820E-05; + 14030.01 1.8998E-05; } } B4C{ temp 566; composition { - 5010.06 1.5206E-02; - 5011.06 6.1514E-02; - 6012.06 1.8972E-02; - //6013.06 2.1252E-04; + 5010.01 1.5206E-02; + 5011.01 6.1514E-02; + 6012.01 1.8972E-02; + 6013.01 2.1252E-04; } } Ag-In-Cd{ temp 566; composition { - 47107.06 2.3523E-02; - 47109.06 2.1854E-02; - 48106.06 3.3882E-05; - 48108.06 2.4166E-05; - 48110.06 3.3936E-04; - 48111.06 3.4821E-04; - 48112.06 6.5611E-04; - 48113.06 3.3275E-04; - 48114.06 7.8252E-04; - 48116.06 2.0443E-04; - 49113.06 3.4219E-04; - 49115.06 7.6511E-03; } + 47107.01 2.3523E-02; + 47109.01 2.1854E-02; + 48106.01 3.3882E-05; + 48108.01 2.4166E-05; + 48110.01 3.3936E-04; + 48111.01 3.4821E-04; + 48112.01 6.5611E-04; + 48113.01 3.3275E-04; + 48114.01 7.8252E-04; + 48116.01 2.0443E-04; + 49113.01 3.4219E-04; + 49115.01 7.6511E-03; } } UO2-16 { temp 566; tms 1; composition { - 8016.03 4.5897E-02; - 8017.03 1.7436E-05; - //8018.03 9.2032E-05; - 92234.03 3.0131E-06; - 92235.03 3.7503E-04; - 92238.03 2.2625E-02;} + 8016.00 4.5897E-02; + 8017.00 1.7436E-05; + 8018.00 9.2032E-05; + 92234.00 3.0131E-06; + 92235.00 3.7503E-04; + 92238.00 2.2625E-02;} } UO2-24 { temp 566; tms 1; composition { - 8016.03 4.5830E-02; - 8017.03 1.7411E-05; - //8018.03 9.1898E-05; - 92234.03 4.4842E-06; - 92235.03 5.5814E-04; - 92238.03 2.2407E-02;} + 8016.00 4.5830E-02; + 8017.00 1.7411E-05; + 8018.00 9.1898E-05; + 92234.00 4.4842E-06; + 92235.00 5.5814E-04; + 92238.00 2.2407E-02;} } UO2-31 { temp 566; tms 1; composition { - 8016.03 4.5853E-02; - 8017.03 1.7420E-05; - //8018.03 9.1942E-05; - 92234.03 5.7987E-06; - 92235.03 7.2175E-04; - 92238.03 2.2253E-02;} + 8016.00 4.5853E-02; + 8017.00 1.7420E-05; + 8018.00 9.1942E-05; + 92234.00 5.7987E-06; + 92235.00 7.2175E-04; + 92238.00 2.2253E-02;} } UO2-32 { temp 566; tms 1; composition { - 8016.03 4.6029E-02; - 8017.03 1.7487E-05; - //8018.03 9.2296E-05; - 92234.03 5.9959E-06; - 92235.03 7.4630E-04; - 92238.03 2.2317E-02; + 8016.00 4.6029E-02; + 8017.00 1.7487E-05; + 8018.00 9.2296E-05; + 92234.00 5.9959E-06; + 92235.00 7.4630E-04; + 92238.00 2.2317E-02; } } - + UO2-34 { temp 566; tms 1; composition { - 8016.03 4.6110E-02; - 8017.03 1.7517E-05; - //8018.03 9.2459E-05; - 92234.03 6.4018E-06; - 92235.03 7.9681E-04; - 92238.03 2.2307E-02;} + 8016.00 4.6110E-02; + 8017.00 1.7517E-05; + 8018.00 9.2459E-05; + 92234.00 6.4018E-06; + 92235.00 7.9681E-04; + 92238.00 2.2307E-02;} } // vanadium51 was stated twice in carbonsteel below @@ -2168,95 +2257,94 @@ nuclearData { CarbonSteel { temp 566; composition { - 13027.06 4.3523E-05; - 5010.06 2.5833E-06; - 5011.06 1.0450E-05; - 6012.06 1.0442E-03; - //6013.06 1.1697E-05 ; - 20040.06 1.7043E-05; - 20042.06 1.1375E-07; - 20043.06 2.3734E-08; - 20044.06 3.6673E-07; - 20046.06 7.0322E-10; - 20048.06 3.2875E-08; - 24050.06 1.3738E-05; - 24052.06 2.6493E-04; - 24053.06 3.0041E-05; - 24054.06 7.4778E-06; - 29063.06 1.0223E-04; - 29065.06 4.5608E-05; - 26054.06 4.7437E-03; - 26056.06 7.4465E-02; - 26057.06 1.7197E-03; - 26058.06 2.2886E-04; - 25055.06 6.4126E-04; - 42100.06 2.9814E-05; - 42092.06 4.4822E-05; - 42094.06 2.8110E-05; - 42095.06 4.8567E-05; - 42096.06 5.1015E-05; - 42097.06 2.9319E-05; - 42098.06 7.4327E-05; - 41093.06 5.0559E-06; - 28058.06 4.0862E-04; - 28060.06 1.5740E-04; - 28061.06 6.8420E-06; - 28062.06 2.1815E-05; - 28064.06 5.5557E-06; - 15031.06 3.7913E-05; - 16032.06 3.4808E-05; - 16033.06 2.7420E-07; - 16034.06 1.5368E-06; - 16036.06 5.3398E-09; - 14028.06 6.1702E-04; - 14029.06 3.1330E-05; - 14030.06 2.0653E-05; - 22046.06 1.2144E-06; - 22047.06 1.0952E-06; - 22048.06 1.0851E-05; - 22049.06 7.9634E-07; - 22050.06 7.6249E-07; - //23050.06 1.1526E-07; - 23051.06 4.5989E-05; + 13027.01 4.3523E-05; + 5010.01 2.5833E-06; + 5011.01 1.0450E-05; + 6012.01 1.0442E-03; + 6013.01 1.1697E-05 ; + 20040.01 1.7043E-05; + 20042.01 1.1375E-07; + 20043.01 2.3734E-08; + 20044.01 3.6673E-07; + 20046.01 7.0322E-10; + 20048.01 3.2875E-08; + 24050.01 1.3738E-05; + 24052.01 2.6493E-04; + 24053.01 3.0041E-05; + 24054.01 7.4778E-06; + 29063.01 1.0223E-04; + 29065.01 4.5608E-05; + 26054.01 4.7437E-03; + 26056.01 7.4465E-02; + 26057.01 1.7197E-03; + 26058.01 2.2886E-04; + 25055.01 6.4126E-04; + 42100.01 2.9814E-05; + 42092.01 4.4822E-05; + 42094.01 2.8110E-05; + 42095.01 4.8567E-05; + 42096.01 5.1015E-05; + 42097.01 2.9319E-05; + 42098.01 7.4327E-05; + 41093.01 5.0559E-06; + 28058.01 4.0862E-04; + 28060.01 1.5740E-04; + 28061.01 6.8420E-06; + 28062.01 2.1815E-05; + 28064.01 5.5557E-06; + 15031.01 3.7913E-05; + 16032.01 3.4808E-05; + 16033.01 2.7420E-07; + 16034.01 1.5368E-06; + 16036.01 5.3398E-09; + 14028.01 6.1702E-04; + 14029.01 3.1330E-05; + 14030.01 2.0653E-05; + 22046.01 1.2144E-06; + 22047.01 1.0952E-06; + 22048.01 1.0851E-05; + 22049.01 7.9634E-07; + 22050.01 7.6249E-07; + 23050.01 1.1526E-07; + 23051.01 4.5989E-05; } } SupportPlateSS { temp 566; composition { - 24050.06 3.5223E-04; - 24052.06 6.7924E-03; - 24053.06 7.7020E-04; - 24054.06 1.9172E-04; - 26054.06 1.5882E-03; - 26056.06 2.4931E-02; - 26057.06 5.7578E-04; - 26058.06 7.6625E-05; - 25055.06 8.0762E-04; - 28058.06 2.5731E-03; - 28060.06 9.9117E-04; - 28061.06 4.3085E-05; - 28062.06 1.3738E-04; - 28064.06 3.4985E-05; - 14028.06 4.3711E-04; - 14029.06 2.2195E-05; - 14030.06 1.4631E-05;} + 24050.01 3.5223E-04; + 24052.01 6.7924E-03; + 24053.01 7.7020E-04; + 24054.01 1.9172E-04; + 26054.01 1.5882E-03; + 26056.01 2.4931E-02; + 26057.01 5.7578E-04; + 26058.01 7.6625E-05; + 25055.01 8.0762E-04; + 28058.01 2.5731E-03; + 28060.01 9.9117E-04; + 28061.01 4.3085E-05; + 28062.01 1.3738E-04; + 28064.01 3.4985E-05; + 14028.01 4.3711E-04; + 14029.01 2.2195E-05; + 14030.01 1.4631E-05;} } SupportPlateBW { temp 566; - moder {1001.06 (lwj3.11 lwj3.09); } + moder {1001.01 (h-h2o.52 h-h2o.53); } composition { - 5010.06 1.0559E-05; - 5011.06 4.2716E-05; - 1001.06 6.5512E-02; - 1002.06 1.0204E-05; - 8016.06 3.2683E-02; - 8017.06 1.2416E-05; - //8018.06 6.5535E-05; + 5010.01 1.0559E-05; + 5011.01 4.2716E-05; + 1001.01 6.5512E-02; + 1002.01 1.0204E-05; + 8016.01 3.2683E-02; + 8017.01 1.2416E-05; + 8018.01 6.5535E-05; } } - } } diff --git a/InputFiles/Benchmarks/BEAVRS/BEAVRS_HZP b/InputFiles/Benchmarks/BEAVRS/BEAVRS_HZP index 5a740ce28..23c0e5dab 100644 --- a/InputFiles/Benchmarks/BEAVRS/BEAVRS_HZP +++ b/InputFiles/Benchmarks/BEAVRS/BEAVRS_HZP @@ -1,13 +1,18 @@ !! !! 3D BEAVRS benchmark !! At HZP, the D bank of rods is partially inserted up -!! to 115 steps inserted / 113 steps withdrawn. -!! A step corresponds to an increment of 1.58193cm +!! to 15 steps inserted / 213 steps withdrawn. +!! A step corresponds to an increment of 1.58193cm. +!! Hence the D bank is inserted by 23.7255cm. +!! +!! This could probably be written in a much more efficient and +!! understandable way! For example, by being more clever with +!! transformations. !! type eigenPhysicsPackage; -pop 1000000; -active 50; +pop 10000000; +active 100; inactive 250; XSdata ce; dataType ce; @@ -16,7 +21,7 @@ collisionOperator { neutronCE {type neutronCEstd;}} transportOperator { !type transportOperatorDT; - type transportOperatorHT; cache 1; + type transportOperatorST; cache 1; } inactiveTally { @@ -28,7 +33,7 @@ inactiveTally { yax { type spaceMap; grid lin; min -161.2773; max 161.2773; N 15; axis y;} zax { type spaceMap; grid lin; min 36.748; max 402.508; N 15; axis z;} } - cycles 200; + cycles 250; } } @@ -49,12 +54,12 @@ activeTally { fissionAxial { type collisionClerk; response (fission); fission { type macroResponse; MT -6;} map {type spaceMap; axis z; grid lin; N 60; min 36.748; max 402.508;} } - fissionYZ { type collisionClerk; response (fission); fission { type macroResponse; MT -6;} - map {type multiMap; maps (yax zax); - yax {type spaceMap; axis y; grid lin; N 255; min -161.2773; max 161.2773; } - zax {type spaceMap; axis z; grid lin; N 60; min 36.748; max 402.508;} - } - } + !fissionYZ { type collisionClerk; response (fission); fission { type macroResponse; MT -6;} + ! map {type multiMap; maps (yax zax); + ! yax {type spaceMap; axis y; grid lin; N 255; min -161.2773; max 161.2773; } + ! zax {type spaceMap; axis z; grid lin; N 60; min 36.748; max 402.508;} + ! } + ! } } geometry { @@ -150,9 +155,12 @@ geometry { planeCRLowerBottom { id 135; type plane; coeffs (0.0 0.0 1.0 402.508); } // same as 109 on withdrawal planeCRUpperBottom { id 136; type plane; coeffs (0.0 0.0 1.0 504.108); } // out of core on withdrawal - plane322p1861 {id 137; type plane; coeffs (0 0 1 322.18609);} - plane220p586 {id 138; type plane; coeffs (0 0 1 220.5861);} - plane218p716 {id 139; type plane; coeffs (0 0 1 218.71609);} + !plane322p1861 {id 137; type plane; coeffs (0 0 1 322.18609);} + plane478p512 {id 137; type plane; coeffs (0 0 1 478.51249);} + !plane220p586 {id 138; type plane; coeffs (0 0 1 220.5861);} + plane378p783 {id 138; type plane; coeffs (0 0 1 378.783);} + !plane218p716 {id 139; type plane; coeffs (0 0 1 218.71609);} + planep376p913 {id 139; type plane; coeffs (0 0 1 376.9125);} } @@ -481,34 +489,36 @@ geometry { GT37 {type simpleCell; id 372; surfaces (-129 131); filltype uni; universe 1010;} GT35 {type simpleCell; id 373; surfaces (-131 132); filltype uni; universe 1005;} GT20 {type simpleCell; id 374; surfaces (-132); filltype uni; universe 1001;} + + // Extra GT cells without dashpot for empty instrumentation tubes + CGT98 {type simpleCell; id 469; surfaces (-123 126); filltype uni; universe 9112;} //guide tube, partially inserted CR, axial layering // Note: universe is basically a combo of fully and partially inserted rods - GP460 {type simpleCell; id 375; surfaces ( 101); filltype uni; universe 1013;} - GP431 {type simpleCell; id 376; surfaces (-101 102); filltype uni; universe 1013;} - GP423 {type simpleCell; id 377; surfaces (-102 106); filltype uni; universe 1013;} - GP415 {type simpleCell; id 378; surfaces (-106 107); filltype uni; universe 1133;} - GP411 {type simpleCell; id 379; surfaces (-107 111); filltype uni; universe 1013;} - GP364 {type simpleCell; id 380; surfaces (-111 112); filltype uni; universe 9233;} - GP359 {type simpleCell; id 381; surfaces (-112 137); filltype uni; universe 1013;} - GP322 {type simpleCell; id 382; surfaces (-137 113); filltype uni; universe 1014;} - GP312 {type simpleCell; id 383; surfaces (-113 114); filltype uni; universe 9232;} - GP306 {type simpleCell; id 384; surfaces (-114 115); filltype uni; universe 1014;} - GP260 {type simpleCell; id 385; surfaces (-115 116); filltype uni; universe 9232;} - GP254 {type simpleCell; id 386; surfaces (-116 138); filltype uni; universe 1014;} - GP220 {type simpleCell; id 387; surfaces (-138 139); filltype uni; universe 1023;} - GP219 {type simpleCell; id 388; surfaces (-139 117); filltype uni; universe 12000;} - GP208 {type simpleCell; id 389; surfaces (-117 118); filltype uni; universe 9112;} - GP202 {type simpleCell; id 390; surfaces (-118 119); filltype uni; universe 12000;} - GP155 {type simpleCell; id 391; surfaces (-119 120); filltype uni; universe 9112;} - GP150 {type simpleCell; id 392; surfaces (-120 122); filltype uni; universe 12000;} - GP103 {type simpleCell; id 393; surfaces (-122 123); filltype uni; universe 9112;} - GP98 {type simpleCell; id 394; surfaces (-123 126); filltype uni; universe 1010;} - GP4052 {type simpleCell; id 395; surfaces (-126 127); filltype uni; universe 12000;} - GP39 {type simpleCell; id 396; surfaces (-127 129); filltype uni; universe 1019;} - GP37 {type simpleCell; id 397; surfaces (-129 131); filltype uni; universe 1010;} - GP35 {type simpleCell; id 398; surfaces (-131 132); filltype uni; universe 1005;} - GP20 {type simpleCell; id 399; surfaces (-132); filltype uni; universe 1001;} + GP460 {type simpleCell; id 375; surfaces ( 101); filltype uni; universe 1014;} + GP431 {type simpleCell; id 376; surfaces (-101 102); filltype uni; universe 1014;} + GP423 {type simpleCell; id 377; surfaces (-102 106); filltype uni; universe 1014;} + GP415 {type simpleCell; id 378; surfaces (-106 107); filltype uni; universe 1132;} + GP411 {type simpleCell; id 379; surfaces (-107 138); filltype uni; universe 1014;} + GP378 {type simpleCell; id 380; surfaces (-138 139); filltype uni; universe 1023;} + GP377 {type simpleCell; id 381; surfaces (-139 111); filltype uni; universe 12000;} + GP364 {type simpleCell; id 382; surfaces (-111 112); filltype uni; universe 9112;} + GP359 {type simpleCell; id 383; surfaces (-112 113); filltype uni; universe 12000;} + GP312 {type simpleCell; id 384; surfaces (-113 114); filltype uni; universe 9112;} + GP306 {type simpleCell; id 385; surfaces (-114 115); filltype uni; universe 12000;} + GP260 {type simpleCell; id 386; surfaces (-115 116); filltype uni; universe 9112;} + GP254 {type simpleCell; id 387; surfaces (-116 117); filltype uni; universe 12000;} + GP208 {type simpleCell; id 388; surfaces (-117 118); filltype uni; universe 9112;} + GP202 {type simpleCell; id 389; surfaces (-118 119); filltype uni; universe 12000;} + GP155 {type simpleCell; id 390; surfaces (-119 120); filltype uni; universe 9112;} + GP150 {type simpleCell; id 391; surfaces (-120 122); filltype uni; universe 12000;} + GP103 {type simpleCell; id 392; surfaces (-122 123); filltype uni; universe 9112;} + GP98 {type simpleCell; id 393; surfaces (-123 126); filltype uni; universe 1010;} + GP4052 {type simpleCell; id 394; surfaces (-126 127); filltype uni; universe 12000;} + GP39 {type simpleCell; id 395; surfaces (-127 129); filltype uni; universe 1019;} + GP37 {type simpleCell; id 396; surfaces (-129 131); filltype uni; universe 1010;} + GP35 {type simpleCell; id 397; surfaces (-131 132); filltype uni; universe 1005;} + GP20 {type simpleCell; id 398; surfaces (-132); filltype uni; universe 1001;} // control rod, axial layering // Used (probably with some modification) only when fully inserted @@ -817,19 +827,53 @@ geometry { //instr. tube instrumentTube { - id 14; + id 15; type cellUniverse; cells (198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217);} + + //empty instr. tube + emptyInstrumentTube { + id 11; + type cellUniverse; + cells (353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 469 370 371 372 373 374);} // control rod, partially inserted GuideTubePartial { - id 15; + id 17; type cellUniverse; - cells (375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399); + cells (375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398); } + // Lattice of instrumentation tubes + instLattice { + type latUniverse; + id 14; + origin (0 0 0); + pitch (21.50364 21.50364 0.0); + shape (15 15 0); + padMat Water; + global 1; + map ( + 11 11 11 11 11 11 15 11 11 15 11 11 11 11 11 + 11 11 15 11 11 15 11 15 11 11 11 11 11 11 11 + 11 11 11 11 11 11 11 15 11 15 11 15 11 15 11 + 11 15 15 11 11 11 11 15 11 11 11 11 11 11 11 + 11 11 11 11 15 11 11 11 15 11 15 11 15 11 11 + 15 11 15 11 11 15 11 15 11 11 11 11 11 15 11 + 11 11 11 15 11 11 15 11 11 15 11 11 15 11 11 + 15 11 15 11 15 11 15 11 11 15 11 15 15 15 11 + 11 15 11 11 11 11 11 11 15 11 15 11 11 11 15 + 11 11 11 11 15 11 15 11 11 11 11 15 11 11 11 + 15 11 11 11 15 11 11 15 11 11 15 11 11 11 15 + 11 11 11 11 11 15 11 11 15 11 11 15 11 11 11 + 11 11 15 11 15 11 11 15 11 11 11 11 11 15 11 + 11 11 15 11 11 11 15 11 11 15 11 15 11 11 11 + 11 11 11 11 15 11 11 15 11 11 11 11 11 11 11 + ); + } + // Lattices w/o grid @@ -1900,19 +1944,19 @@ geometry { map ( 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 - 24 24 24 24 24 15 24 24 15 24 24 15 24 24 24 24 24 - 24 24 24 15 24 24 24 24 24 24 24 24 24 15 24 24 24 + 24 24 24 24 24 17 24 24 17 24 24 17 24 24 24 24 24 + 24 24 24 17 24 24 24 24 24 24 24 24 24 17 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 - 24 24 15 24 24 15 24 24 15 24 24 15 24 24 15 24 24 + 24 24 17 24 24 17 24 24 17 24 24 17 24 24 17 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 - 24 24 15 24 24 15 24 24 14 24 24 15 24 24 15 24 24 + 24 24 17 24 24 17 24 24 14 24 24 17 24 24 17 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 - 24 24 15 24 24 15 24 24 15 24 24 15 24 24 15 24 24 + 24 24 17 24 24 17 24 24 17 24 24 17 24 24 17 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 - 24 24 24 15 24 24 24 24 24 24 24 24 24 15 24 24 24 - 24 24 24 24 24 15 24 24 15 24 24 15 24 24 24 24 24 + 24 24 24 17 24 24 24 24 24 24 24 24 24 17 24 24 24 + 24 24 24 24 24 17 24 24 17 24 24 17 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 ); } @@ -1935,21 +1979,23 @@ geometry { map ( 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 - 16 16 16 16 16 15 16 16 15 16 16 15 16 16 16 16 16 - 16 16 16 15 16 16 16 16 16 16 16 16 16 15 16 16 16 + 16 16 16 16 16 17 16 16 17 16 16 17 16 16 16 16 16 + 16 16 16 17 16 16 16 16 16 16 16 16 16 17 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 - 16 16 15 16 16 15 16 16 15 16 16 15 16 16 15 16 16 + 16 16 17 16 16 17 16 16 17 16 16 17 16 16 17 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 - 16 16 15 16 16 15 16 16 14 16 16 15 16 16 15 16 16 + 16 16 17 16 16 17 16 16 14 16 16 17 16 16 17 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 - 16 16 15 16 16 15 16 16 15 16 16 15 16 16 15 16 16 + 16 16 17 16 16 17 16 16 17 16 16 17 16 16 17 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 - 16 16 16 15 16 16 16 16 16 16 16 16 16 15 16 16 16 - 16 16 16 16 16 15 16 16 15 16 16 15 16 16 16 16 16 + 16 16 16 17 16 16 16 16 16 16 16 16 16 17 16 16 16 + 16 16 16 16 16 17 16 16 17 16 16 17 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 - 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16); } + 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16); + + } A0E16SleeveDBank { id 14170; @@ -2037,7 +2083,7 @@ geometry { 1001 1001 1001 1001 1800 5224 5224 5224 5224 5224 5224 5224 1700 1001 1001 1001 1001 1001 1001 1800 5224 5229 24310 7031120 24310 7031120 24310 7031120 24310 5228 5224 1700 1001 1001 1001 1800 5229 24310 14310 26310 14160 30310 14160 30310 14160 26310 14310 24310 5228 1700 1001 - 1001 5222 24310 2531110 26240 14160 26240 14160 26240 14160 26240 14160 26240 25311 24310 5223 1001 + 1001 5222 24310 2531110 26240 14160 26240 14160 26240 14160 26240 14160 26240 253110 24310 5223 1001 1800 5229 14310 26240 14250 26240 24160 22240 14160 22240 24160 26240 14250 26240 14310 5228 1700 5222 24310 26310 14160 26240 24160 22240 24160 22240 24160 22240 24160 26240 14160 26310 24310 5223 5222 703190 14160 26240 24160 22240 14160 22240 14160 22240 14160 22240 24160 26240 14160 703130 5223 @@ -2052,8 +2098,6 @@ geometry { 1001 1001 1600 5225 5227 24310 703160 24310 703160 24310 703160 24310 5226 5225 1500 1001 1001 1001 1001 1001 1001 1600 5225 5225 5225 5225 5225 5225 5225 1500 1001 1001 1001 1001 ); } -! Note partial rodded assemblies end with a 1, i.e., 14170 and 14250 - coreAndStructures { id 8888; type cellUniverse; @@ -2070,19 +2114,28 @@ viz { type bmp; output imgXY; what material; - centre (-17.13 240.69 167.74); - width (50 50); + centre (0 0 167.74); + width (100 100); axis z; - res (1000 1000); + res (2000 2000); } bmpYZ { type bmp; output imgYZ; what material; - centre (0.0 0.0 232.0); - width (100.0 200.0); + centre (0.0 0.0 230.0); + !width (400.0 400.0); + axis x; + res (2000 2000); + } + bmpXZ { + type bmp; + output imgXZ; + what material; + centre (0.0 0.0 230.0); + !width (400.0 400.0); axis x; - res (1000 2000); + res (2000 2000); } } @@ -2094,208 +2147,205 @@ nuclearData { materials { - // Note that commented nuclide densities are included in the specification - // but are not available in the JEFF-3.11 data library - - Air { + Air { temp 566; composition { - 18036.06 7.8730E-09; - 18038.06 1.4844E-09; - 18040.06 2.3506E-06; - 6012.06 6.7539E-08; - //6013.06 7.5658E-10; - 7014.06 1.9680E-04; - 7015.06 7.2354E-07; - 8016.06 5.2866E-05; - 8017.06 2.0084E-08; - //8018.06 1.0601E-07; + 18036.01 7.8730E-09; + 18038.01 1.4844E-09; + 18040.01 2.3506E-06; + 6012.01 6.7539E-08; + 6013.01 7.5658E-10; + 7014.01 1.9680E-04; + 7015.01 7.2354E-07; + 8016.01 5.2866E-05; + 8017.01 2.0084E-08; + //8018.01 1.0601E-07; } } SS304 { temp 566; composition { - 24050.06 7.6778E-04; - 24052.06 1.4806E-02; - 24053.06 1.6789E-03; - 24054.06 4.1791E-04; - 26054.06 3.4620E-03; - 26056.06 5.4345E-02; - 26057.06 1.2551E-03; - 26058.06 1.6703E-04; - 25055.06 1.7604E-03; - 28058.06 5.6089E-03; - 28060.06 2.1605E-03; - 28061.06 9.3917E-05; - 28062.06 2.9945E-04; - 28064.06 7.6261E-05; - 14028.06 9.5281E-04; - 14029.06 4.8381E-05; - 14030.06 3.1893E-05; } + 24050.01 7.6778E-04; + 24052.01 1.4806E-02; + 24053.01 1.6789E-03; + 24054.01 4.1791E-04; + 26054.01 3.4620E-03; + 26056.01 5.4345E-02; + 26057.01 1.2551E-03; + 26058.01 1.6703E-04; + 25055.01 1.7604E-03; + 28058.01 5.6089E-03; + 28060.01 2.1605E-03; + 28061.01 9.3917E-05; + 28062.01 2.9945E-04; + 28064.01 7.6261E-05; + 14028.01 9.5281E-04; + 14029.01 4.8381E-05; + 14030.01 3.1893E-05; } } Helium { temp 566; composition { - 2003.06 4.8089E-10; - 2004.06 2.4044E-04; } + 2003.01 4.8089E-10; + 2004.01 2.4044E-04; } } BorosilicateGlass { temp 566; composition { - 13027.06 1.7352E-03; - 5010.06 9.6506E-04; - 5011.06 3.9189E-03; - 8016.06 4.6514E-02; - 8017.06 1.7671E-05; - //8018.06 9.3268E-05; - 14028.06 1.6926E-02; - 14029.06 8.5944E-04; - 14030.06 5.6654E-04; } + 13027.01 1.7352E-03; + 5010.01 9.6506E-04; + 5011.01 3.9189E-03; + 8016.01 4.6514E-02; + 8017.01 1.7671E-05; + 8018.01 9.3268E-05; + 14028.01 1.6926E-02; + 14029.01 8.5944E-04; + 14030.01 5.6654E-04; } } Water { temp 566; - moder {1001.06 (lwj3.11 lwj3.09); } + moder {1001.01 (h-h2o.52 h-h2o.53); } composition { - 5010.06 7.9714E-06; - 5011.06 3.2247E-05; - 1001.06 4.9456E-02; - 1002.06 7.7035E-06; - 8016.06 2.4673E-02; - 8017.06 9.3734E-06; - //8018.06 4.9474E-05; + 5010.01 7.9714E-06; + 5011.01 3.2247E-05; + 1001.01 4.9456E-02; + 1002.01 7.7035E-06; + 8016.01 2.4673E-02; + 8017.01 9.3734E-06; + 8018.01 4.9474E-05; } } Zircaloy { temp 566; composition { - 24050.06 3.2962E-06; - 24052.06 6.3564E-05; - 24053.06 7.2076E-06; - 24054.06 1.7941E-06; - 26054.06 8.6698E-06; - 26056.06 1.3610E-04; - 26057.06 3.1431E-06; - 26058.06 4.1829E-07; - 8016.06 3.0744E-04; - 8017.06 1.1680E-07; - //8018.03 6.1648E-07; - 50112.06 4.6735E-06; - 50114.06 3.1799E-06; - 50115.06 1.6381E-06; - 50116.06 7.0055E-05; - 50117.06 3.7003E-05; - 50118.06 1.1669E-04; - 50119.06 4.1387E-05; - 50120.06 1.5697E-04; - 50122.06 2.2308E-05; - 50124.06 2.7897E-05; - 40090.06 2.1828E-02; - 40091.06 4.7601E-03; - 40092.06 7.2759E-03; - 40094.06 7.3734E-03; - 40096.06 1.1879E-03; } + 24050.01 3.2962E-06; + 24052.01 6.3564E-05; + 24053.01 7.2076E-06; + 24054.01 1.7941E-06; + 26054.01 8.6698E-06; + 26056.01 1.3610E-04; + 26057.01 3.1431E-06; + 26058.01 4.1829E-07; + 8016.01 3.0744E-04; + 8017.01 1.1680E-07; + 8018.01 6.1648E-07; + 50112.01 4.6735E-06; + 50114.01 3.1799E-06; + 50115.01 1.6381E-06; + 50116.01 7.0055E-05; + 50117.01 3.7003E-05; + 50118.01 1.1669E-04; + 50119.01 4.1387E-05; + 50120.01 1.5697E-04; + 50122.01 2.2308E-05; + 50124.01 2.7897E-05; + 40090.01 2.1828E-02; + 40091.01 4.7601E-03; + 40092.01 7.2759E-03; + 40094.01 7.3734E-03; + 40096.01 1.1879E-03; } } Inconel{ temp 566; composition { - 24050.06 7.8239E-04; - 24052.06 1.5088E-02; - 24053.06 1.7108E-03; - 24054.06 4.2586E-04; - 26054.06 1.4797E-03; - 26056.06 2.3229E-02; - 26057.06 5.3645E-04; - 26058.06 7.1392E-05; - 25055.06 7.8201E-04; - 28058.06 2.9320E-02; - 28060.06 1.1294E-02; - 28061.06 4.9094E-04; - 28062.06 1.5653E-03; - 28064.06 3.9864E-04; - 14028.06 5.6757E-04; - 14029.06 2.8820E-05; - 14030.06 1.8998E-05; } + 24050.01 7.8239E-04; + 24052.01 1.5088E-02; + 24053.01 1.7108E-03; + 24054.01 4.2586E-04; + 26054.01 1.4797E-03; + 26056.01 2.3229E-02; + 26057.01 5.3645E-04; + 26058.01 7.1392E-05; + 25055.01 7.8201E-04; + 28058.01 2.9320E-02; + 28060.01 1.1294E-02; + 28061.01 4.9094E-04; + 28062.01 1.5653E-03; + 28064.01 3.9864E-04; + 14028.01 5.6757E-04; + 14029.01 2.8820E-05; + 14030.01 1.8998E-05; } } B4C{ temp 566; composition { - 5010.06 1.5206E-02; - 5011.06 6.1514E-02; - 6012.06 1.8972E-02; - //6013.06 2.1252E-04; + 5010.01 1.5206E-02; + 5011.01 6.1514E-02; + 6012.01 1.8972E-02; + 6013.01 2.1252E-04; } } Ag-In-Cd{ temp 566; composition { - 47107.06 2.3523E-02; - 47109.06 2.1854E-02; - 48106.06 3.3882E-05; - 48108.06 2.4166E-05; - 48110.06 3.3936E-04; - 48111.06 3.4821E-04; - 48112.06 6.5611E-04; - 48113.06 3.3275E-04; - 48114.06 7.8252E-04; - 48116.06 2.0443E-04; - 49113.06 3.4219E-04; - 49115.06 7.6511E-03; } + 47107.01 2.3523E-02; + 47109.01 2.1854E-02; + 48106.01 3.3882E-05; + 48108.01 2.4166E-05; + 48110.01 3.3936E-04; + 48111.01 3.4821E-04; + 48112.01 6.5611E-04; + 48113.01 3.3275E-04; + 48114.01 7.8252E-04; + 48116.01 2.0443E-04; + 49113.01 3.4219E-04; + 49115.01 7.6511E-03; } } UO2-16 { temp 566; tms 1; composition { - 8016.03 4.5897E-02; - 8017.03 1.7436E-05; - //8018.03 9.2032E-05; - 92234.03 3.0131E-06; - 92235.03 3.7503E-04; - 92238.03 2.2625E-02;} + 8016.00 4.5897E-02; + 8017.00 1.7436E-05; + 8018.00 9.2032E-05; + 92234.00 3.0131E-06; + 92235.00 3.7503E-04; + 92238.00 2.2625E-02;} } UO2-24 { temp 566; tms 1; composition { - 8016.03 4.5830E-02; - 8017.03 1.7411E-05; - //8018.03 9.1898E-05; - 92234.03 4.4842E-06; - 92235.03 5.5814E-04; - 92238.03 2.2407E-02;} + 8016.00 4.5830E-02; + 8017.00 1.7411E-05; + 8018.00 9.1898E-05; + 92234.00 4.4842E-06; + 92235.00 5.5814E-04; + 92238.00 2.2407E-02;} } UO2-31 { temp 566; tms 1; composition { - 8016.03 4.5853E-02; - 8017.03 1.7420E-05; - //8018.03 9.1942E-05; - 92234.03 5.7987E-06; - 92235.03 7.2175E-04; - 92238.03 2.2253E-02;} + 8016.00 4.5853E-02; + 8017.00 1.7420E-05; + 8018.00 9.1942E-05; + 92234.00 5.7987E-06; + 92235.00 7.2175E-04; + 92238.00 2.2253E-02;} } UO2-32 { temp 566; tms 1; composition { - 8016.03 4.6029E-02; - 8017.03 1.7487E-05; - //8018.03 9.2296E-05; - 92234.03 5.9959E-06; - 92235.03 7.4630E-04; - 92238.03 2.2317E-02; + 8016.00 4.6029E-02; + 8017.00 1.7487E-05; + 8018.00 9.2296E-05; + 92234.00 5.9959E-06; + 92235.00 7.4630E-04; + 92238.00 2.2317E-02; } } @@ -2303,12 +2353,12 @@ nuclearData { temp 566; tms 1; composition { - 8016.03 4.6110E-02; - 8017.03 1.7517E-05; - //8018.03 9.2459E-05; - 92234.03 6.4018E-06; - 92235.03 7.9681E-04; - 92238.03 2.2307E-02;} + 8016.00 4.6110E-02; + 8017.00 1.7517E-05; + 8018.00 9.2459E-05; + 92234.00 6.4018E-06; + 92235.00 7.9681E-04; + 92238.00 2.2307E-02;} } // vanadium51 was stated twice in carbonsteel below @@ -2316,95 +2366,96 @@ nuclearData { CarbonSteel { temp 566; composition { - 13027.06 4.3523E-05; - 5010.06 2.5833E-06; - 5011.06 1.0450E-05; - 6012.06 1.0442E-03; - //6013.06 1.1697E-05 ; - 20040.06 1.7043E-05; - 20042.06 1.1375E-07; - 20043.06 2.3734E-08; - 20044.06 3.6673E-07; - 20046.06 7.0322E-10; - 20048.06 3.2875E-08; - 24050.06 1.3738E-05; - 24052.06 2.6493E-04; - 24053.06 3.0041E-05; - 24054.06 7.4778E-06; - 29063.06 1.0223E-04; - 29065.06 4.5608E-05; - 26054.06 4.7437E-03; - 26056.06 7.4465E-02; - 26057.06 1.7197E-03; - 26058.06 2.2886E-04; - 25055.06 6.4126E-04; - 42100.06 2.9814E-05; - 42092.06 4.4822E-05; - 42094.06 2.8110E-05; - 42095.06 4.8567E-05; - 42096.06 5.1015E-05; - 42097.06 2.9319E-05; - 42098.06 7.4327E-05; - 41093.06 5.0559E-06; - 28058.06 4.0862E-04; - 28060.06 1.5740E-04; - 28061.06 6.8420E-06; - 28062.06 2.1815E-05; - 28064.06 5.5557E-06; - 15031.06 3.7913E-05; - 16032.06 3.4808E-05; - 16033.06 2.7420E-07; - 16034.06 1.5368E-06; - 16036.06 5.3398E-09; - 14028.06 6.1702E-04; - 14029.06 3.1330E-05; - 14030.06 2.0653E-05; - 22046.06 1.2144E-06; - 22047.06 1.0952E-06; - 22048.06 1.0851E-05; - 22049.06 7.9634E-07; - 22050.06 7.6249E-07; - //23050.06 1.1526E-07; - 23051.06 4.5989E-05; + 13027.01 4.3523E-05; + 5010.01 2.5833E-06; + 5011.01 1.0450E-05; + 6012.01 1.0442E-03; + 6013.01 1.1697E-05 ; + 20040.01 1.7043E-05; + 20042.01 1.1375E-07; + 20043.01 2.3734E-08; + 20044.01 3.6673E-07; + 20046.01 7.0322E-10; + 20048.01 3.2875E-08; + 24050.01 1.3738E-05; + 24052.01 2.6493E-04; + 24053.01 3.0041E-05; + 24054.01 7.4778E-06; + 29063.01 1.0223E-04; + 29065.01 4.5608E-05; + 26054.01 4.7437E-03; + 26056.01 7.4465E-02; + 26057.01 1.7197E-03; + 26058.01 2.2886E-04; + 25055.01 6.4126E-04; + 42100.01 2.9814E-05; + 42092.01 4.4822E-05; + 42094.01 2.8110E-05; + 42095.01 4.8567E-05; + 42096.01 5.1015E-05; + 42097.01 2.9319E-05; + 42098.01 7.4327E-05; + 41093.01 5.0559E-06; + 28058.01 4.0862E-04; + 28060.01 1.5740E-04; + 28061.01 6.8420E-06; + 28062.01 2.1815E-05; + 28064.01 5.5557E-06; + 15031.01 3.7913E-05; + 16032.01 3.4808E-05; + 16033.01 2.7420E-07; + 16034.01 1.5368E-06; + 16036.01 5.3398E-09; + 14028.01 6.1702E-04; + 14029.01 3.1330E-05; + 14030.01 2.0653E-05; + 22046.01 1.2144E-06; + 22047.01 1.0952E-06; + 22048.01 1.0851E-05; + 22049.01 7.9634E-07; + 22050.01 7.6249E-07; + 23050.01 1.1526E-07; + 23051.01 4.5989E-05; } } SupportPlateSS { temp 566; composition { - 24050.06 3.5223E-04; - 24052.06 6.7924E-03; - 24053.06 7.7020E-04; - 24054.06 1.9172E-04; - 26054.06 1.5882E-03; - 26056.06 2.4931E-02; - 26057.06 5.7578E-04; - 26058.06 7.6625E-05; - 25055.06 8.0762E-04; - 28058.06 2.5731E-03; - 28060.06 9.9117E-04; - 28061.06 4.3085E-05; - 28062.06 1.3738E-04; - 28064.06 3.4985E-05; - 14028.06 4.3711E-04; - 14029.06 2.2195E-05; - 14030.06 1.4631E-05;} + 24050.01 3.5223E-04; + 24052.01 6.7924E-03; + 24053.01 7.7020E-04; + 24054.01 1.9172E-04; + 26054.01 1.5882E-03; + 26056.01 2.4931E-02; + 26057.01 5.7578E-04; + 26058.01 7.6625E-05; + 25055.01 8.0762E-04; + 28058.01 2.5731E-03; + 28060.01 9.9117E-04; + 28061.01 4.3085E-05; + 28062.01 1.3738E-04; + 28064.01 3.4985E-05; + 14028.01 4.3711E-04; + 14029.01 2.2195E-05; + 14030.01 1.4631E-05;} } SupportPlateBW { temp 566; - moder {1001.06 (lwj3.11 lwj3.09); } + moder {1001.01 (h-h2o.52 h-h2o.53); } composition { - 5010.06 1.0559E-05; - 5011.06 4.2716E-05; - 1001.06 6.5512E-02; - 1002.06 1.0204E-05; - 8016.06 3.2683E-02; - 8017.06 1.2416E-05; - //8018.06 6.5535E-05; + 5010.01 1.0559E-05; + 5011.01 4.2716E-05; + 1001.01 6.5512E-02; + 1002.01 1.0204E-05; + 8016.01 3.2683E-02; + 8017.01 1.2416E-05; + 8018.01 6.5535E-05; } } - + + } } diff --git a/SharedModules/genericProcedures.f90 b/SharedModules/genericProcedures.f90 index ddcf9f286..086c6ed69 100644 --- a/SharedModules/genericProcedures.f90 +++ b/SharedModules/genericProcedures.f90 @@ -13,8 +13,10 @@ module genericProcedures interface swap module procedure swap_shortInt module procedure swap_defReal + module procedure swap_defBool module procedure swap_char_nameLen module procedure swap_defReal_defReal + module procedure swap_defBool_defBool end interface interface quickSort @@ -1348,7 +1350,7 @@ elemental subroutine swap_shortInt(i1,i2) end subroutine swap_shortInt !! - !! Swap to reals + !! Swap two reals !! elemental subroutine swap_defReal(r1,r2) real(defReal), intent(inout) :: r1 @@ -1384,6 +1386,44 @@ elemental subroutine swap_defReal_defReal(r1_1, r1_2, r2_1, r2_2) r2_2 = temp2 end subroutine swap_defReal_defReal + + !! + !! Swap two bools + !! + elemental subroutine swap_defBool(r1,r2) + logical(defBool), intent(inout) :: r1 + logical(defBool), intent(inout) :: r2 + logical(defBool) :: temp + + temp = r1 + r1 = r2 + r2 = temp + + end subroutine swap_defBool + + !! + !! Swap two pair of bools + !! + elemental subroutine swap_defBool_defBool(r1_1, r1_2, r2_1, r2_2) + logical(defBool), intent(inout) :: r1_1 + logical(defBool), intent(inout) :: r1_2 + logical(defBool), intent(inout) :: r2_1 + logical(defBool), intent(inout) :: r2_2 + logical(defBool) :: temp1, temp2 + + ! Load first pair into temps + temp1 = r1_1 + temp2 = r1_2 + + ! Assign values of 2nd pair to 1st pair + r1_1 = r2_1 + r1_2 = r2_2 + + ! Assign values of 1st pair to 2nd pair + r2_1 = temp1 + r2_2 = temp2 + + end subroutine swap_defBool_defBool !! !! Swap character of length nameLen diff --git a/SharedModules/universalVariables.f90 b/SharedModules/universalVariables.f90 index fdb8a2d57..a75393dca 100644 --- a/SharedModules/universalVariables.f90 +++ b/SharedModules/universalVariables.f90 @@ -7,7 +7,7 @@ module universalVariables ! *** DON't CHANGE THIS. HARDCODED IS FINE ! CHANGE THIS: NUMBER MUST BE CALCULATED DURING INITIAL GEOMETRY PROCESSING ! Problematic for separating modules! - integer(shortInt), parameter, public :: HARDCODED_MAX_NEST = 8 + integer(shortInt), parameter, public :: HARDCODED_MAX_NEST = 9 integer(shortInt), parameter, public :: MAX_OUTGOING_PARTICLES = 5 ! CHANGE THIS: NUMBER WILL DEPEND ON SYSTEM ARCHITECTURE From fe0f8452f9ed5519a181c7c646f22a5bcb9d5919 Mon Sep 17 00:00:00 2001 From: ChasingNeutrons Date: Tue, 24 Jun 2025 18:23:28 +0100 Subject: [PATCH 4/6] Add docs, small fix, and added example input Added description of offsets in latUniverse and global in universe to docs. Also added an example of the offsetMap in the form of a 2D CROCUS input. --- Geometry/Universes/latUniverse_class.f90 | 2 +- InputFiles/CROCUS_2D | 201 +++++++++++++++++++++++ docs/User Manual.rst | 11 ++ 3 files changed, 213 insertions(+), 1 deletion(-) create mode 100644 InputFiles/CROCUS_2D diff --git a/Geometry/Universes/latUniverse_class.f90 b/Geometry/Universes/latUniverse_class.f90 index 5dea8b0a1..1060e36c4 100644 --- a/Geometry/Universes/latUniverse_class.f90 +++ b/Geometry/Universes/latUniverse_class.f90 @@ -396,7 +396,7 @@ function cellOffset(self, coords) result (offset) logical(defBool) :: doOffset if (allocated(self % offsetMap)) then - doOffset = self % offsetMap(coords % localID) == noOffset + doOffset = self % offsetMap(coords % localID) == local else doOffset = self % offset end if diff --git a/InputFiles/CROCUS_2D b/InputFiles/CROCUS_2D new file mode 100644 index 000000000..f4a3b69a7 --- /dev/null +++ b/InputFiles/CROCUS_2D @@ -0,0 +1,201 @@ +!! +!! 2D model of the CROCUS reactor +!! Composition and geometry taken from IRPhE +!! +type eigenPhysicsPackage; + +pop 100000; +active 500; +inactive 100; +XSdata ce; +dataType ce; +outputFile CROCUS_2D; + +// Specify output format default asciiMATLAB +//outputFormat asciiJSON; + +collisionOperator { neutronCE {type neutronCEstd;} } + +transportOperator {// type transportOperatorST; + cache 1; + type transportOperatorHT; + } + +inactiveTally {} + +activeTally { + fluxRadial { type collisionClerk; response (flx); flx { type fluxResponse;} + map {type multiMap; maps (xax yax emap); + xax {type spaceMap; axis x; grid lin; N 240; min -60.0; max 60.0; } + yax {type spaceMap; axis y; grid lin; N 240; min -60.0; max 60.0; } + emap {type energyMap; grid unstruct; bins (1.0E-11 0.625E-6 20); } + } + } +} + +geometry { + type geometryStd; + boundary ( 0 0 0 0 1 1); + graph {type shrunk;} + + surfaces { + boundary { id 1; type zTruncCylinder; origin (0.0 0.0 0.0); radius 60; halfwidth 60;} + } + + cells { + // in {id 1; surfaces ( -1); filltype uni; universe 401;} + // out {id 2; surfaces (1); filltype outside; } + } + + universes { + root { id 1; type rootUniverse; border 1; fill u<10>; } + + // Pin universes + pinUO2 { id 2; type pinUniverse; radii (0.526 0.545 0.63 0.0 ); fills (uo2 void clad water);} + pinUMet {id 3; type pinUniverse; radii (0.85 0.8675 0.965 0.0 ); fills (umet void clad water);} + water {id 4; type pinUniverse; radii (0.0 ); fills (water);} + + // Lattices + inner { + id 5; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (1.837 1.837 0.0); + shape (22 22 0); + padMat water; + map ( + + 4 4 4 4 4 4 4 4 2 2 2 2 2 2 4 4 4 4 4 4 4 4 + 4 4 4 4 4 4 4 4 2 2 2 2 2 2 4 4 4 4 4 4 4 4 + 4 4 4 4 4 2 2 2 2 2 2 2 2 2 2 2 2 4 4 4 4 4 + 4 4 4 4 4 2 2 2 2 2 2 2 2 2 2 2 2 4 4 4 4 4 + 4 4 4 4 4 2 2 2 2 2 2 2 2 2 2 2 2 4 4 4 4 4 + 4 4 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 4 4 + 4 4 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 4 4 + 4 4 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 4 4 + 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 + 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 + 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 + 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 + 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 + 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 + 4 4 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 4 4 + 4 4 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 4 4 + 4 4 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 4 4 + 4 4 4 4 4 2 2 2 2 2 2 2 2 2 2 2 2 4 4 4 4 4 + 4 4 4 4 4 2 2 2 2 2 2 2 2 2 2 2 2 4 4 4 4 4 + 4 4 4 4 4 2 2 2 2 2 2 2 2 2 2 2 2 4 4 4 4 4 + 4 4 4 4 4 4 4 4 2 2 2 2 2 2 4 4 4 4 4 4 4 4 + 4 4 4 4 4 4 4 4 2 2 2 2 2 2 4 4 4 4 4 4 4 4 + ); + } + + outer { + id 10; + type latUniverse; + origin (0.0 0.0 0.0); + pitch (2.917 2.917 0.0); + shape (20 20 0); + padMat water; + map ( + 4 4 4 4 4 4 4 3 3 3 3 3 3 4 4 4 4 4 4 4 + 4 4 4 4 3 3 3 3 3 3 3 3 3 3 3 4 4 4 4 4 + 4 4 4 3 3 3 3 3 3 3 3 3 3 3 3 3 3 4 4 4 + 4 4 3 3 3 3 3 3 5 5 5 5 3 3 3 3 3 3 4 4 + 4 3 3 3 3 3 5 5 5 5 5 5 5 5 3 3 3 3 4 4 + 4 3 3 3 3 3 5 5 5 5 5 5 5 5 3 3 3 3 3 4 + 3 3 3 3 5 5 5 5 5 5 5 5 5 5 5 5 3 3 3 4 + 3 3 3 3 5 5 5 5 5 5 5 5 5 5 5 5 3 3 3 3 + 3 3 3 5 5 5 5 5 5 5 5 5 5 5 5 5 5 3 3 3 + 3 3 3 5 5 5 5 5 5 5 5 5 5 5 5 5 5 3 3 3 + 3 3 3 5 5 5 5 5 5 5 5 5 5 5 5 5 5 3 3 3 + 3 3 3 5 5 5 5 5 5 5 5 5 5 5 5 5 5 3 3 3 + 3 3 3 3 5 5 5 5 5 5 5 5 5 5 5 5 3 3 3 3 + 4 3 3 3 5 5 5 5 5 5 5 5 5 5 5 5 3 3 3 3 + 4 3 3 3 3 3 5 5 5 5 5 5 5 5 3 3 3 3 3 4 + 4 4 3 3 3 3 5 5 5 5 5 5 5 5 3 3 3 3 3 4 + 4 4 3 3 3 3 3 3 5 5 5 5 3 3 3 3 3 3 4 4 + 4 4 4 3 3 3 3 3 3 3 3 3 3 3 3 3 3 4 4 4 + 4 4 4 4 4 3 3 3 3 3 3 3 3 3 3 3 4 4 4 4 + 4 4 4 4 4 4 4 3 3 3 3 3 3 4 4 4 4 4 4 4 + ); + offsetMap ( + 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 0 0 0 0 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 0 0 0 0 0 0 0 0 1 1 1 1 1 1 + 1 1 1 1 1 1 0 0 0 0 0 0 0 0 1 1 1 1 1 1 + 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 + 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 + 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 + 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 + 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 + 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 + 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 + 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 + 1 1 1 1 1 1 0 0 0 0 0 0 0 0 1 1 1 1 1 1 + 1 1 1 1 1 1 0 0 0 0 0 0 0 0 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 0 0 0 0 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 + ); + } + } + +} + + +viz { + bmp { + type bmp; + output img; + what material; + centre (0.0 0.0 0.0); + !width (80.0 80.0); + axis z; + offset -17; + res (1500 1500); + } +} + + +nuclearData { + handles { + ce {type aceNeutronDatabase; aceLibrary $SCONE_ACE;} + } + materials { + + water { + temp 300; + rgb (0 0 139); + moder { 1001.03 (lwj3.00); } + composition { + 1001.03 6.67578E-2; + 8016.03 3.33789E-2; } + } + + clad { + temp 300; + composition { + 13027.03 6.02611E-2;} + } + + uo2 { + temp 300; + composition { + 92235.03 4.30565E-4; + 92238.03 2.31145E-2; + 8016.03 4.70902E-2; } + } + + umet { + temp 300; + composition { + 92235.03 4.5316E-4; + 92238.03 4.68003E-2;} + } + +} +} diff --git a/docs/User Manual.rst b/docs/User Manual.rst index 56cc5eb9f..bccd0c926 100644 --- a/docs/User Manual.rst +++ b/docs/User Manual.rst @@ -613,6 +613,10 @@ Similarly to the surfaces and cells, the **universes** in the geometry can be de { id ; type ; *keywords* } } +One can disable the translation of a given universe, should it be nested in others, by the use +of the keyword ``global 1;``. This evaluates a particle's position in the universe using the +global frame of reference. + Several ``universeTypes`` are possible: * cellUniverse, composed of the union of different cells. Note that overlaps are @@ -667,6 +671,13 @@ Example: :: origin of the universe. [cm] - rotation (*optional*, default = (0.0 0.0 0.0)): (x y z) array with the rotation angles in degrees applied to the universe. [°] + - offsetMap (*optional*, default = all elements offset): map that specifies which elements + of the lattice are offset with respect to the lattice origin. Elements with 1 are offset, + while elements with a 0 are not. Must have the same size as the map. Allows creating, e.g., + BWR assemblies with water rods covering multiple lattice elements. + - offset (*optional*, default = true): enables/disables the offset of all entries in + the latUniverse. Has relatively specialised use cases, e.g., imposing a discretisation + by placing another universe inside the lattice. Example: :: From 3db89076bcf0bb382155cb16c9dd8a16b21ac97c Mon Sep 17 00:00:00 2001 From: ChasingNeutrons Date: Tue, 24 Jun 2025 19:44:38 +0100 Subject: [PATCH 5/6] Added flux tallies to BEAVRS HZP --- InputFiles/Benchmarks/BEAVRS/BEAVRS_HZP | 9 ++++++++- 1 file changed, 8 insertions(+), 1 deletion(-) diff --git a/InputFiles/Benchmarks/BEAVRS/BEAVRS_HZP b/InputFiles/Benchmarks/BEAVRS/BEAVRS_HZP index 23c0e5dab..09102661b 100644 --- a/InputFiles/Benchmarks/BEAVRS/BEAVRS_HZP +++ b/InputFiles/Benchmarks/BEAVRS/BEAVRS_HZP @@ -21,7 +21,7 @@ collisionOperator { neutronCE {type neutronCEstd;}} transportOperator { !type transportOperatorDT; - type transportOperatorST; cache 1; + type transportOperatorHT; cache 1; } inactiveTally { @@ -45,6 +45,13 @@ activeTally { yax {type spaceMap; axis y; grid lin; N 255; min -161.2773; max 161.2773; } } } + pinFlxRadial { type collisionClerk; response (flx); flx { type fluxResponse;} + map {type multiMap; maps (xax yax emap); + xax {type spaceMap; axis x; grid lin; N 255; min -161.2773; max 161.2773; } + yax {type spaceMap; axis y; grid lin; N 255; min -161.2773; max 161.2773; } + emap {type energyMap; grid unstruct; bins (1.0E-11 0.625E-6 20); } + } + } assemblyFissRadial { type collisionClerk; response (fission); fission { type macroResponse; MT -6;} map {type multiMap; maps (xax yax); xax {type spaceMap; axis x; grid lin; N 15; min -161.2773; max 161.2773; } From 387e6ce8cb1a5d2f910c4929089eb5e3f08f7c2d Mon Sep 17 00:00:00 2001 From: Paul Cosgrove Date: Tue, 12 Aug 2025 12:41:35 +0100 Subject: [PATCH 6/6] Update Geometry/Universes/latUniverse_class.f90 Co-authored-by: valeriaRaffuzzi <108435337+valeriaRaffuzzi@users.noreply.github.com> --- Geometry/Universes/latUniverse_class.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/Geometry/Universes/latUniverse_class.f90 b/Geometry/Universes/latUniverse_class.f90 index 1060e36c4..7f890bbc4 100644 --- a/Geometry/Universes/latUniverse_class.f90 +++ b/Geometry/Universes/latUniverse_class.f90 @@ -36,7 +36,7 @@ module latUniverse_class !! !! Every lattice cell has an offset to its centre (so the centre of the nested universe !! is in the center of the lattice cell). Optionally an offset map can be provided, determining - !! whether to apply an offset in a given cell position. This can disables the local universe + !! whether to apply an offset in a given cell position. This can disable the local universe !! offset. Alternatively a single offset flag can be provided, disabling offset in all cells. !! !! Minimum lattice pitch is set to 10 * SURF_TOL