From 082f41d0f5a806e2f864fb1225fa543ce5d9d565 Mon Sep 17 00:00:00 2001
From: Will Boyd <wbinventor@gmail.com>
Date: Tue, 7 Feb 2017 19:39:18 -0500
Subject: [PATCH] Simplified SMR grid spacers; removed SMR debugging plots

---
 smr/build-xml.py  |   6 +-
 smr/geometry.xml  | 846 ++++++++++++++++++++++++++++++++++++++++++++++
 smr/grids.py      |  36 --
 smr/materials.xml | 203 +++++++++++
 smr/pins.py       |  10 +-
 smr/plots.py      |  50 +--
 smr/plots.xml     |  60 ++++
 smr/settings.xml  |  19 ++
 smr/smr.py        |  59 ----
 smr/surfaces.py   |  98 ++----
 smr/tallies.xml   |  13 +
 11 files changed, 1180 insertions(+), 220 deletions(-)
 create mode 100644 smr/geometry.xml
 delete mode 100644 smr/grids.py
 create mode 100644 smr/materials.xml
 create mode 100644 smr/plots.xml
 create mode 100644 smr/settings.xml
 delete mode 100644 smr/smr.py
 create mode 100644 smr/tallies.xml

diff --git a/smr/build-xml.py b/smr/build-xml.py
index 463b2e3..fcbeb9c 100644
--- a/smr/build-xml.py
+++ b/smr/build-xml.py
@@ -23,8 +23,8 @@ multipole = input('Use multipole cross sections? (y/n): ').lower()
 multipole = True if multipole == 'y' else False
 
 # Construct uniform initial source distribution over fissionable zones
-lower_left = [-8.*lattice_pitch/2., -8.*lattice_pitch/2., bottom_fuel_stack]
-upper_right = [+8.*lattice_pitch/2., +8.*lattice_pitch/2., top_active_core]
+lower_left = [-7.*lattice_pitch/2., -7.*lattice_pitch/2., bottom_fuel_stack]
+upper_right = [+7.*lattice_pitch/2., +7.*lattice_pitch/2., top_active_core]
 source = openmc.source.Source(space=openmc.stats.Box(lower_left, upper_right))
 source.space.only_fissionable = True
 
@@ -52,7 +52,7 @@ plots.export_to_xml()
 # Create a pin-wise mesh
 mesh = openmc.Mesh(name='pin-wise mesh')
 mesh.type = 'regular'
-mesh.dimension = [15*17, 15*17]
+mesh.dimension = [7*17, 7*17]
 mesh.lower_left = lower_left[:2]
 mesh.width = np.array(upper_right[:2]) - np.array(lower_left[:2])
 mesh_filter = openmc.MeshFilter(mesh)
diff --git a/smr/geometry.xml b/smr/geometry.xml
new file mode 100644
index 0000000..7986082
--- /dev/null
+++ b/smr/geometry.xml
@@ -0,0 +1,846 @@
+<?xml version='1.0' encoding='utf-8'?>
+<geometry>
+    <cell id="10000" material="10007" name="water pin 1" universe="10000" />
+    <cell id="10001" material="10007" name="GT empty (0)" region="-10004" universe="10001" />
+    <cell id="10002" material="10005" name="GT empty (1)" region="10004 -10005" universe="10001" />
+    <cell id="10003" material="10007" name="GT empty (last)" region="10005" universe="10001" />
+    <cell id="10004" material="10007" name="GT empty grid (top/bottom) (0)" region="-10004" universe="10002" />
+    <cell id="10005" material="10005" name="GT empty grid (top/bottom) (1)" region="10004 -10005" universe="10002" />
+    <cell id="10006" material="10007" name="GT empty grid (top/bottom) (last)" region="10005 (10019 -10020 10021 -10022)" universe="10002" />
+    <cell id="10007" material="10005" name="GT empty grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10002" />
+    <cell id="10008" material="10007" name="GT empty grid (intermediate) (0)" region="-10004" universe="10003" />
+    <cell id="10009" material="10005" name="GT empty grid (intermediate) (1)" region="10004 -10005" universe="10003" />
+    <cell id="10010" material="10007" name="GT empty grid (intermediate) (last)" region="10005 (10023 -10024 10025 -10026)" universe="10003" />
+    <cell id="10011" material="10005" name="GT empty grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10003" />
+    <cell id="10012" material="10007" name="GT empty nozzle (0)" region="-10004" universe="10004" />
+    <cell id="10013" material="10005" name="GT empty nozzle (1)" region="10004 -10005" universe="10004" />
+    <cell id="10014" material="10007" name="GT empty nozzle (last)" region="10005" universe="10004" />
+    <cell id="10015" material="10007" name="GT empty at dashpot (0)" region="-10006" universe="10005" />
+    <cell id="10016" material="10005" name="GT empty at dashpot (1)" region="10006 -10007" universe="10005" />
+    <cell id="10017" material="10007" name="GT empty at dashpot (last)" region="10007" universe="10005" />
+    <cell id="10018" material="10007" name="GT empty at dashpot grid (top/bottom) (0)" region="-10006" universe="10006" />
+    <cell id="10019" material="10005" name="GT empty at dashpot grid (top/bottom) (1)" region="10006 -10007" universe="10006" />
+    <cell id="10020" material="10007" name="GT empty at dashpot grid (top/bottom) (last)" region="10007 (10019 -10020 10021 -10022)" universe="10006" />
+    <cell id="10021" material="10005" name="GT empty at dashpot grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10006" />
+    <cell id="10022" material="10007" name="GT empty at dashpot grid (intermediate) (0)" region="-10006" universe="10007" />
+    <cell id="10023" material="10005" name="GT empty at dashpot grid (intermediate) (1)" region="10006 -10007" universe="10007" />
+    <cell id="10024" material="10007" name="GT empty at dashpot grid (intermediate) (last)" region="10007 (10023 -10024 10025 -10026)" universe="10007" />
+    <cell id="10025" material="10005" name="GT empty at dashpot grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10007" />
+    <cell id="10026" material="10007" name="GT empty nozzle (0)" region="-10006" universe="10008" />
+    <cell id="10027" material="10005" name="GT empty nozzle (1)" region="10006 -10007" universe="10008" />
+    <cell id="10028" material="10007" name="GT empty nozzle (last)" region="10007" universe="10008" />
+    <cell fill="10000" id="10029" name="GT empty (0)" region="-10036" universe="10009" />
+    <cell fill="10000" id="10030" name="GT empty (1)" region="10036 -10037" universe="10009" />
+    <cell fill="10000" id="10031" name="GT empty (2)" region="10037 -10038" universe="10009" />
+    <cell fill="10005" id="10032" name="GT empty (3)" region="10038 -10039" universe="10009" />
+    <cell fill="10005" id="10033" name="GT empty (4)" region="10039 -10042" universe="10009" />
+    <cell fill="10006" id="10034" name="GT empty (5)" region="10042 -10043" universe="10009" />
+    <cell fill="10005" id="10035" name="GT empty (6)" region="10043 -10044" universe="10009" />
+    <cell fill="10001" id="10036" name="GT empty (7)" region="10044 -10045" universe="10009" />
+    <cell fill="10003" id="10037" name="GT empty (8)" region="10045 -10046" universe="10009" />
+    <cell fill="10001" id="10038" name="GT empty (9)" region="10046 -10047" universe="10009" />
+    <cell fill="10003" id="10039" name="GT empty (10)" region="10047 -10048" universe="10009" />
+    <cell fill="10001" id="10040" name="GT empty (11)" region="10048 -10049" universe="10009" />
+    <cell fill="10001" id="10041" name="GT empty (12)" region="10049 -10050" universe="10009" />
+    <cell fill="10002" id="10042" name="GT empty (13)" region="10050 -10040" universe="10009" />
+    <cell fill="10001" id="10043" name="GT empty (14)" region="10040 -10051" universe="10009" />
+    <cell fill="10001" id="10044" name="GT empty (15)" region="10051 -10052" universe="10009" />
+    <cell fill="10001" id="10045" name="GT empty (16)" region="10052 -10053" universe="10009" />
+    <cell fill="10000" id="10046" name="GT empty (17)" region="10053 -10054" universe="10009" />
+    <cell fill="10000" id="10047" name="GT empty (last)" region="10054" universe="10009" />
+    <cell fill="10000" id="10048" name="GT empty instr (0)" region="-10036" universe="10010" />
+    <cell fill="10000" id="10049" name="GT empty instr (1)" region="10036 -10037" universe="10010" />
+    <cell fill="10000" id="10050" name="GT empty instr (2)" region="10037 -10038" universe="10010" />
+    <cell fill="10009" id="10051" name="GT empty instr (3)" region="10038 -10039" universe="10010" />
+    <cell fill="10009" id="10052" name="GT empty instr (4)" region="10039 -10042" universe="10010" />
+    <cell fill="10002" id="10053" name="GT empty instr (5)" region="10042 -10043" universe="10010" />
+    <cell fill="10009" id="10054" name="GT empty instr (6)" region="10043 -10044" universe="10010" />
+    <cell fill="10009" id="10055" name="GT empty instr (7)" region="10044 -10045" universe="10010" />
+    <cell fill="10003" id="10056" name="GT empty instr (8)" region="10045 -10046" universe="10010" />
+    <cell fill="10009" id="10057" name="GT empty instr (9)" region="10046 -10047" universe="10010" />
+    <cell fill="10003" id="10058" name="GT empty instr (10)" region="10047 -10048" universe="10010" />
+    <cell fill="10009" id="10059" name="GT empty instr (11)" region="10048 -10049" universe="10010" />
+    <cell fill="10009" id="10060" name="GT empty instr (12)" region="10049 -10050" universe="10010" />
+    <cell fill="10002" id="10061" name="GT empty instr (13)" region="10050 -10040" universe="10010" />
+    <cell fill="10009" id="10062" name="GT empty instr (14)" region="10040 -10051" universe="10010" />
+    <cell fill="10009" id="10063" name="GT empty instr (15)" region="10051 -10052" universe="10010" />
+    <cell fill="10009" id="10064" name="GT empty instr (16)" region="10052 -10053" universe="10010" />
+    <cell fill="10000" id="10065" name="GT empty instr (17)" region="10053 -10054" universe="10010" />
+    <cell fill="10000" id="10066" name="GT empty instr (last)" region="10054" universe="10010" />
+    <cell id="10067" material="10001" name="IT (0)" region="-10015" universe="10011" />
+    <cell id="10068" material="10005" name="IT (1)" region="10015 -10016" universe="10011" />
+    <cell id="10069" material="10007" name="IT (2)" region="10016 -10004" universe="10011" />
+    <cell id="10070" material="10005" name="IT (3)" region="10004 -10005" universe="10011" />
+    <cell id="10071" material="10007" name="IT (last)" region="10005" universe="10011" />
+    <cell id="10072" material="10001" name="IT grid (top/bottom) (0)" region="-10015" universe="10012" />
+    <cell id="10073" material="10005" name="IT grid (top/bottom) (1)" region="10015 -10016" universe="10012" />
+    <cell id="10074" material="10007" name="IT grid (top/bottom) (2)" region="10016 -10004" universe="10012" />
+    <cell id="10075" material="10005" name="IT grid (top/bottom) (3)" region="10004 -10005" universe="10012" />
+    <cell id="10076" material="10007" name="IT grid (top/bottom) (last)" region="10005 (10019 -10020 10021 -10022)" universe="10012" />
+    <cell id="10077" material="10005" name="IT grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10012" />
+    <cell id="10078" material="10001" name="IT grid (intermediate) (0)" region="-10015" universe="10013" />
+    <cell id="10079" material="10005" name="IT grid (intermediate) (1)" region="10015 -10016" universe="10013" />
+    <cell id="10080" material="10007" name="IT grid (intermediate) (2)" region="10016 -10004" universe="10013" />
+    <cell id="10081" material="10005" name="IT grid (intermediate) (3)" region="10004 -10005" universe="10013" />
+    <cell id="10082" material="10007" name="IT grid (intermediate) (last)" region="10005 (10023 -10024 10025 -10026)" universe="10013" />
+    <cell id="10083" material="10005" name="IT grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10013" />
+    <cell id="10089" material="10001" name="IT dashpot (0)" region="-10015" universe="10015" />
+    <cell id="10090" material="10005" name="IT dashpot (1)" region="10015 -10016" universe="10015" />
+    <cell id="10091" material="10007" name="IT dashpot (last)" region="10016" universe="10015" />
+    <cell fill="10015" id="10092" name="GT instr (0)" region="-10036" universe="10016" />
+    <cell fill="10015" id="10093" name="GT instr (1)" region="10036 -10037" universe="10016" />
+    <cell fill="10015" id="10094" name="GT instr (2)" region="10037 -10038" universe="10016" />
+    <cell fill="10011" id="10095" name="GT instr (3)" region="10038 -10039" universe="10016" />
+    <cell fill="10011" id="10096" name="GT instr (4)" region="10039 -10042" universe="10016" />
+    <cell fill="10012" id="10097" name="GT instr (5)" region="10042 -10043" universe="10016" />
+    <cell fill="10011" id="10098" name="GT instr (6)" region="10043 -10044" universe="10016" />
+    <cell fill="10011" id="10099" name="GT instr (7)" region="10044 -10045" universe="10016" />
+    <cell fill="10013" id="10100" name="GT instr (8)" region="10045 -10046" universe="10016" />
+    <cell fill="10011" id="10101" name="GT instr (9)" region="10046 -10047" universe="10016" />
+    <cell fill="10013" id="10102" name="GT instr (10)" region="10047 -10048" universe="10016" />
+    <cell fill="10011" id="10103" name="GT instr (11)" region="10048 -10049" universe="10016" />
+    <cell fill="10011" id="10104" name="GT instr (12)" region="10049 -10050" universe="10016" />
+    <cell fill="10012" id="10105" name="GT instr (13)" region="10050 -10040" universe="10016" />
+    <cell fill="10011" id="10106" name="GT instr (14)" region="10040 -10051" universe="10016" />
+    <cell fill="10011" id="10107" name="GT instr (15)" region="10051 -10052" universe="10016" />
+    <cell fill="10011" id="10108" name="GT instr (16)" region="10052 -10053" universe="10016" />
+    <cell fill="10015" id="10109" name="GT instr (17)" region="10053 -10054" universe="10016" />
+    <cell fill="10000" id="10110" name="GT instr (last)" region="10054" universe="10016" />
+    <cell id="10111" material="10006" name="CR (0)" region="-10008" universe="10017" />
+    <cell id="10112" material="10001" name="CR (1)" region="10008 -10009" universe="10017" />
+    <cell id="10113" material="10003" name="CR (2)" region="10009 -10004" universe="10017" />
+    <cell id="10114" material="10007" name="CR (3)" region="10004 -10005" universe="10017" />
+    <cell id="10115" material="10007" name="CR (last)" region="10005" universe="10017" />
+    <cell id="10116" material="10006" name="CR grid (top/bottom) (0)" region="-10008" universe="10018" />
+    <cell id="10117" material="10001" name="CR grid (top/bottom) (1)" region="10008 -10009" universe="10018" />
+    <cell id="10118" material="10003" name="CR grid (top/bottom) (2)" region="10009 -10004" universe="10018" />
+    <cell id="10119" material="10007" name="CR grid (top/bottom) (3)" region="10004 -10005" universe="10018" />
+    <cell id="10120" material="10007" name="CR grid (top/bottom) (last)" region="10005 (10019 -10020 10021 -10022)" universe="10018" />
+    <cell id="10121" material="10005" name="CR grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10018" />
+    <cell id="10122" material="10006" name="CR grid (intermediate) (0)" region="-10008" universe="10019" />
+    <cell id="10123" material="10001" name="CR grid (intermediate) (1)" region="10008 -10009" universe="10019" />
+    <cell id="10124" material="10003" name="CR grid (intermediate) (2)" region="10009 -10004" universe="10019" />
+    <cell id="10125" material="10007" name="CR grid (intermediate) (3)" region="10004 -10005" universe="10019" />
+    <cell id="10126" material="10007" name="CR grid (intermediate) (last)" region="10005 (10023 -10024 10025 -10026)" universe="10019" />
+    <cell id="10127" material="10005" name="CR grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10019" />
+    <cell id="10128" material="10006" name="CR nozzle (0)" region="-10008" universe="10020" />
+    <cell id="10129" material="10001" name="CR nozzle (1)" region="10008 -10009" universe="10020" />
+    <cell id="10130" material="10003" name="CR nozzle (2)" region="10009 -10010" universe="10020" />
+    <cell id="10131" material="10007" name="CR nozzle (last)" region="10010" universe="10020" />
+    <cell id="10132" material="10003" name="CR blank (0)" region="-10008" universe="10021" />
+    <cell id="10133" material="10001" name="CR blank (1)" region="10008 -10009" universe="10021" />
+    <cell id="10134" material="10003" name="CR blank (2)" region="10009 -10010" universe="10021" />
+    <cell id="10135" material="10007" name="CR blank (3)" region="10010 -10004" universe="10021" />
+    <cell id="10136" material="10005" name="CR blank (4)" region="10004 -10005" universe="10021" />
+    <cell id="10137" material="10007" name="CR blank (last)" region="10005" universe="10021" />
+    <cell id="10138" material="10003" name="CR blank grid (top/bottom) (0)" region="-10008" universe="10022" />
+    <cell id="10139" material="10001" name="CR blank grid (top/bottom) (1)" region="10008 -10009" universe="10022" />
+    <cell id="10140" material="10003" name="CR blank grid (top/bottom) (2)" region="10009 -10010" universe="10022" />
+    <cell id="10141" material="10007" name="CR blank grid (top/bottom) (3)" region="10010 -10004" universe="10022" />
+    <cell id="10142" material="10005" name="CR blank grid (top/bottom) (4)" region="10004 -10005" universe="10022" />
+    <cell id="10143" material="10007" name="CR blank grid (top/bottom) (last)" region="10005 (10019 -10020 10021 -10022)" universe="10022" />
+    <cell id="10144" material="10005" name="CR blank grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10022" />
+    <cell id="10145" material="10003" name="CR blank grid (intermediate) (0)" region="-10008" universe="10023" />
+    <cell id="10146" material="10001" name="CR blank grid (intermediate) (1)" region="10008 -10009" universe="10023" />
+    <cell id="10147" material="10003" name="CR blank grid (intermediate) (2)" region="10009 -10010" universe="10023" />
+    <cell id="10148" material="10007" name="CR blank grid (intermediate) (3)" region="10010 -10004" universe="10023" />
+    <cell id="10149" material="10005" name="CR blank grid (intermediate) (4)" region="10004 -10005" universe="10023" />
+    <cell id="10150" material="10007" name="CR blank grid (intermediate) (last)" region="10005 (10023 -10024 10025 -10026)" universe="10023" />
+    <cell id="10151" material="10005" name="CR blank grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10023" />
+    <cell id="10152" material="10003" name="CR blank nozzle (0)" region="-10008" universe="10024" />
+    <cell id="10153" material="10001" name="CR blank nozzle (1)" region="10008 -10009" universe="10024" />
+    <cell id="10154" material="10003" name="CR blank nozzle (2)" region="10009 -10010" universe="10024" />
+    <cell id="10155" material="10007" name="CR blank nozzle (last)" region="10010" universe="10024" />
+    <cell id="10156" material="10003" name="CR blank bare (0)" region="-10008" universe="10025" />
+    <cell id="10157" material="10001" name="CR blank bare (1)" region="10008 -10009" universe="10025" />
+    <cell id="10158" material="10003" name="CR blank bare (2)" region="10009 -10010" universe="10025" />
+    <cell id="10159" material="10007" name="CR blank bare (last)" region="10010" universe="10025" />
+    <cell id="10160" material="10006" name="CR bare (0)" region="-10008" universe="10026" />
+    <cell id="10161" material="10001" name="CR bare (1)" region="10008 -10009" universe="10026" />
+    <cell id="10162" material="10003" name="CR bare (2)" region="10009 -10010" universe="10026" />
+    <cell id="10163" material="10007" name="CR bare (last)" region="10010" universe="10026" />
+    <cell fill="10000" id="10208" name="GT CR bank B dummy (0)" region="-10038" universe="10033" />
+    <cell fill="10005" id="10209" name="GT CR bank B dummy (1)" region="10038 -10044" universe="10033" />
+    <cell fill="10009" id="10210" name="GT CR bank B dummy (2)" region="10044 -10069" universe="10033" />
+    <cell fill="10017" id="10211" name="GT CR bank B dummy (3)" region="10069 -10068" universe="10033" />
+    <cell fill="10021" id="10212" name="GT CR bank B dummy (last)" region="10068" universe="10033" />
+    <cell fill="10000" id="10213" name="GT CR bank B dummy grid (top/bottom) (0)" region="-10038" universe="10034" />
+    <cell fill="10006" id="10214" name="GT CR bank B dummy grid (top/bottom) (1)" region="10038 -10044" universe="10034" />
+    <cell fill="10002" id="10215" name="GT CR bank B dummy grid (top/bottom) (2)" region="10044 -10069" universe="10034" />
+    <cell fill="10018" id="10216" name="GT CR bank B dummy grid (top/bottom) (3)" region="10069 -10068" universe="10034" />
+    <cell fill="10022" id="10217" name="GT CR bank B dummy grid (top/bottom) (last)" region="10068" universe="10034" />
+    <cell fill="10000" id="10218" name="GT CR bank B dummy grid (intermediate) (0)" region="-10038" universe="10035" />
+    <cell fill="10007" id="10219" name="GT CR bank B dummy grid (intermediate) (1)" region="10038 -10044" universe="10035" />
+    <cell fill="10003" id="10220" name="GT CR bank B dummy grid (intermediate) (2)" region="10044 -10069" universe="10035" />
+    <cell fill="10019" id="10221" name="GT CR bank B dummy grid (intermediate) (3)" region="10069 -10068" universe="10035" />
+    <cell fill="10023" id="10222" name="GT CR bank B dummy grid (intermediate) (last)" region="10068" universe="10035" />
+    <cell fill="10000" id="10223" name="GT CR bank B dummy nozzle (0)" region="-10038" universe="10036" />
+    <cell fill="10008" id="10224" name="GT CR bank B dummy nozzle (1)" region="10038 -10044" universe="10036" />
+    <cell fill="10004" id="10225" name="GT CR bank B dummy nozzle (2)" region="10044 -10069" universe="10036" />
+    <cell fill="10020" id="10226" name="GT CR bank B dummy nozzle (3)" region="10069 -10068" universe="10036" />
+    <cell fill="10024" id="10227" name="GT CR bank B dummy nozzle (last)" region="10068" universe="10036" />
+    <cell fill="10000" id="10228" name="GT CR bank B dummy bare (0)" region="-10038" universe="10037" />
+    <cell fill="10008" id="10229" name="GT CR bank B dummy bare (1)" region="10038 -10044" universe="10037" />
+    <cell fill="10004" id="10230" name="GT CR bank B dummy bare (2)" region="10044 -10069" universe="10037" />
+    <cell fill="10026" id="10231" name="GT CR bank B dummy bare (3)" region="10069 -10068" universe="10037" />
+    <cell fill="10025" id="10232" name="GT CR bank B dummy bare (last)" region="10068" universe="10037" />
+    <cell fill="10000" id="10233" name="GT CR bank B (0)" region="-10036" universe="10038" />
+    <cell fill="10036" id="10234" name="GT CR bank B (1)" region="10036 -10037" universe="10038" />
+    <cell fill="10036" id="10235" name="GT CR bank B (2)" region="10037 -10038" universe="10038" />
+    <cell fill="10033" id="10236" name="GT CR bank B (3)" region="10038 -10039" universe="10038" />
+    <cell fill="10033" id="10237" name="GT CR bank B (4)" region="10039 -10042" universe="10038" />
+    <cell fill="10034" id="10238" name="GT CR bank B (5)" region="10042 -10043" universe="10038" />
+    <cell fill="10033" id="10239" name="GT CR bank B (6)" region="10043 -10044" universe="10038" />
+    <cell fill="10033" id="10240" name="GT CR bank B (7)" region="10044 -10045" universe="10038" />
+    <cell fill="10035" id="10241" name="GT CR bank B (8)" region="10045 -10046" universe="10038" />
+    <cell fill="10033" id="10242" name="GT CR bank B (9)" region="10046 -10047" universe="10038" />
+    <cell fill="10035" id="10243" name="GT CR bank B (10)" region="10047 -10048" universe="10038" />
+    <cell fill="10033" id="10244" name="GT CR bank B (11)" region="10048 -10049" universe="10038" />
+    <cell fill="10033" id="10245" name="GT CR bank B (12)" region="10049 -10050" universe="10038" />
+    <cell fill="10034" id="10246" name="GT CR bank B (13)" region="10050 -10040" universe="10038" />
+    <cell fill="10033" id="10247" name="GT CR bank B (14)" region="10040 -10051" universe="10038" />
+    <cell fill="10033" id="10248" name="GT CR bank B (15)" region="10051 -10052" universe="10038" />
+    <cell fill="10033" id="10249" name="GT CR bank B (16)" region="10052 -10053" universe="10038" />
+    <cell fill="10037" id="10250" name="GT CR bank B (17)" region="10053 -10054" universe="10038" />
+    <cell fill="10037" id="10251" name="GT CR bank B (last)" region="10054" universe="10038" />
+    <cell fill="10000" id="10296" name="GT CR bank D dummy (0)" region="-10038" universe="10045" />
+    <cell fill="10005" id="10297" name="GT CR bank D dummy (1)" region="10038 -10044" universe="10045" />
+    <cell fill="10009" id="10298" name="GT CR bank D dummy (2)" region="10044 -10073" universe="10045" />
+    <cell fill="10017" id="10299" name="GT CR bank D dummy (3)" region="10073 -10072" universe="10045" />
+    <cell fill="10021" id="10300" name="GT CR bank D dummy (last)" region="10072" universe="10045" />
+    <cell fill="10000" id="10301" name="GT CR bank D dummy grid (top/bottom) (0)" region="-10038" universe="10046" />
+    <cell fill="10006" id="10302" name="GT CR bank D dummy grid (top/bottom) (1)" region="10038 -10044" universe="10046" />
+    <cell fill="10002" id="10303" name="GT CR bank D dummy grid (top/bottom) (2)" region="10044 -10073" universe="10046" />
+    <cell fill="10018" id="10304" name="GT CR bank D dummy grid (top/bottom) (3)" region="10073 -10072" universe="10046" />
+    <cell fill="10022" id="10305" name="GT CR bank D dummy grid (top/bottom) (last)" region="10072" universe="10046" />
+    <cell fill="10000" id="10306" name="GT CR bank D dummy grid (intermediate) (0)" region="-10038" universe="10047" />
+    <cell fill="10007" id="10307" name="GT CR bank D dummy grid (intermediate) (1)" region="10038 -10044" universe="10047" />
+    <cell fill="10003" id="10308" name="GT CR bank D dummy grid (intermediate) (2)" region="10044 -10073" universe="10047" />
+    <cell fill="10019" id="10309" name="GT CR bank D dummy grid (intermediate) (3)" region="10073 -10072" universe="10047" />
+    <cell fill="10023" id="10310" name="GT CR bank D dummy grid (intermediate) (last)" region="10072" universe="10047" />
+    <cell fill="10000" id="10311" name="GT CR bank D dummy nozzle (0)" region="-10038" universe="10048" />
+    <cell fill="10008" id="10312" name="GT CR bank D dummy nozzle (1)" region="10038 -10044" universe="10048" />
+    <cell fill="10004" id="10313" name="GT CR bank D dummy nozzle (2)" region="10044 -10073" universe="10048" />
+    <cell fill="10020" id="10314" name="GT CR bank D dummy nozzle (3)" region="10073 -10072" universe="10048" />
+    <cell fill="10024" id="10315" name="GT CR bank D dummy nozzle (last)" region="10072" universe="10048" />
+    <cell fill="10000" id="10316" name="GT CR bank D dummy bare (0)" region="-10038" universe="10049" />
+    <cell fill="10008" id="10317" name="GT CR bank D dummy bare (1)" region="10038 -10044" universe="10049" />
+    <cell fill="10004" id="10318" name="GT CR bank D dummy bare (2)" region="10044 -10073" universe="10049" />
+    <cell fill="10026" id="10319" name="GT CR bank D dummy bare (3)" region="10073 -10072" universe="10049" />
+    <cell fill="10025" id="10320" name="GT CR bank D dummy bare (last)" region="10072" universe="10049" />
+    <cell fill="10000" id="10321" name="GT CR bank D (0)" region="-10036" universe="10050" />
+    <cell fill="10048" id="10322" name="GT CR bank D (1)" region="10036 -10037" universe="10050" />
+    <cell fill="10048" id="10323" name="GT CR bank D (2)" region="10037 -10038" universe="10050" />
+    <cell fill="10045" id="10324" name="GT CR bank D (3)" region="10038 -10039" universe="10050" />
+    <cell fill="10045" id="10325" name="GT CR bank D (4)" region="10039 -10042" universe="10050" />
+    <cell fill="10046" id="10326" name="GT CR bank D (5)" region="10042 -10043" universe="10050" />
+    <cell fill="10045" id="10327" name="GT CR bank D (6)" region="10043 -10044" universe="10050" />
+    <cell fill="10045" id="10328" name="GT CR bank D (7)" region="10044 -10045" universe="10050" />
+    <cell fill="10047" id="10329" name="GT CR bank D (8)" region="10045 -10046" universe="10050" />
+    <cell fill="10045" id="10330" name="GT CR bank D (9)" region="10046 -10047" universe="10050" />
+    <cell fill="10047" id="10331" name="GT CR bank D (10)" region="10047 -10048" universe="10050" />
+    <cell fill="10045" id="10332" name="GT CR bank D (11)" region="10048 -10049" universe="10050" />
+    <cell fill="10045" id="10333" name="GT CR bank D (12)" region="10049 -10050" universe="10050" />
+    <cell fill="10046" id="10334" name="GT CR bank D (13)" region="10050 -10040" universe="10050" />
+    <cell fill="10045" id="10335" name="GT CR bank D (14)" region="10040 -10051" universe="10050" />
+    <cell fill="10045" id="10336" name="GT CR bank D (15)" region="10051 -10052" universe="10050" />
+    <cell fill="10045" id="10337" name="GT CR bank D (16)" region="10052 -10053" universe="10050" />
+    <cell fill="10049" id="10338" name="GT CR bank D (17)" region="10053 -10054" universe="10050" />
+    <cell fill="10049" id="10339" name="GT CR bank D (last)" region="10054" universe="10050" />
+    <cell id="10560" material="10001" name="BA (0)" region="-10011" universe="10081" />
+    <cell id="10561" material="10003" name="BA (1)" region="10011 -10012" universe="10081" />
+    <cell id="10562" material="10001" name="BA (2)" region="10012 -10013" universe="10081" />
+    <cell id="10563" material="10008" name="BA (3)" region="10013 -10014" universe="10081" />
+    <cell id="10564" material="10001" name="BA (4)" region="10014 -10015" universe="10081" />
+    <cell id="10565" material="10003" name="BA (5)" region="10015 -10016" universe="10081" />
+    <cell id="10566" material="10007" name="BA (6)" region="10016 -10017" universe="10081" />
+    <cell id="10567" material="10005" name="BA (7)" region="10017 -10018" universe="10081" />
+    <cell id="10568" material="10007" name="BA (last)" region="10018" universe="10081" />
+    <cell id="10579" material="10001" name="BA grid (intermediate) (0)" region="-10011" universe="10083" />
+    <cell id="10580" material="10003" name="BA grid (intermediate) (1)" region="10011 -10012" universe="10083" />
+    <cell id="10581" material="10001" name="BA grid (intermediate) (2)" region="10012 -10013" universe="10083" />
+    <cell id="10582" material="10008" name="BA grid (intermediate) (3)" region="10013 -10014" universe="10083" />
+    <cell id="10583" material="10001" name="BA grid (intermediate) (4)" region="10014 -10015" universe="10083" />
+    <cell id="10584" material="10003" name="BA grid (intermediate) (5)" region="10015 -10016" universe="10083" />
+    <cell id="10585" material="10007" name="BA grid (intermediate) (6)" region="10016 -10017" universe="10083" />
+    <cell id="10586" material="10005" name="BA grid (intermediate) (7)" region="10017 -10018" universe="10083" />
+    <cell id="10587" material="10007" name="BA grid (intermediate) (last)" region="10018 (10023 -10024 10025 -10026)" universe="10083" />
+    <cell id="10588" material="10005" name="BA grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10083" />
+    <cell id="10589" material="10001" name="BA dashpot (0)" region="-10011" universe="10084" />
+    <cell id="10590" material="10003" name="BA dashpot (1)" region="10011 -10012" universe="10084" />
+    <cell id="10591" material="10001" name="BA dashpot (2)" region="10012 -10013" universe="10084" />
+    <cell id="10592" material="10008" name="BA dashpot (3)" region="10013 -10014" universe="10084" />
+    <cell id="10593" material="10001" name="BA dashpot (4)" region="10014 -10015" universe="10084" />
+    <cell id="10594" material="10003" name="BA dashpot (5)" region="10015 -10016" universe="10084" />
+    <cell id="10595" material="10007" name="BA dashpot (6)" region="10016 -10006" universe="10084" />
+    <cell id="10596" material="10005" name="BA dashpot (7)" region="10006 -10007" universe="10084" />
+    <cell id="10597" material="10007" name="BA dashpot (last)" region="10007" universe="10084" />
+    <cell id="10598" material="10001" name="BA dashpot grid (top/bottom) (0)" region="-10011" universe="10085" />
+    <cell id="10599" material="10003" name="BA dashpot grid (top/bottom) (1)" region="10011 -10012" universe="10085" />
+    <cell id="10600" material="10001" name="BA dashpot grid (top/bottom) (2)" region="10012 -10013" universe="10085" />
+    <cell id="10601" material="10008" name="BA dashpot grid (top/bottom) (3)" region="10013 -10014" universe="10085" />
+    <cell id="10602" material="10001" name="BA dashpot grid (top/bottom) (4)" region="10014 -10015" universe="10085" />
+    <cell id="10603" material="10003" name="BA dashpot grid (top/bottom) (5)" region="10015 -10016" universe="10085" />
+    <cell id="10604" material="10007" name="BA dashpot grid (top/bottom) (6)" region="10016 -10006" universe="10085" />
+    <cell id="10605" material="10005" name="BA dashpot grid (top/bottom) (7)" region="10006 -10007" universe="10085" />
+    <cell id="10606" material="10007" name="BA dashpot grid (top/bottom) (last)" region="10007 (10019 -10020 10021 -10022)" universe="10085" />
+    <cell id="10607" material="10005" name="BA dashpot grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10085" />
+    <cell id="10618" material="10003" name="BA blank SS (0)" region="-10016" universe="10087" />
+    <cell id="10619" material="10007" name="BA blank SS (1)" region="10016 -10017" universe="10087" />
+    <cell id="10620" material="10005" name="BA blank SS (2)" region="10017 -10018" universe="10087" />
+    <cell id="10621" material="10007" name="BA blank SS (last)" region="10018" universe="10087" />
+    <cell id="10622" material="10003" name="BA blank SS bare (0)" region="-10016" universe="10088" />
+    <cell id="10623" material="10007" name="BA blank SS bare (last)" region="10016" universe="10088" />
+    <cell fill="10000" id="10624" name="BA stack (0)" region="-10036" universe="10089" />
+    <cell fill="10000" id="10625" name="BA stack (1)" region="10036 -10037" universe="10089" />
+    <cell fill="10000" id="10626" name="BA stack (2)" region="10037 -10038" universe="10089" />
+    <cell fill="10005" id="10627" name="BA stack (3)" region="10038 -10039" universe="10089" />
+    <cell fill="10005" id="10628" name="BA stack (4)" region="10039 -10042" universe="10089" />
+    <cell fill="10006" id="10629" name="BA stack (5)" region="10042 -10041" universe="10089" />
+    <cell fill="10085" id="10630" name="BA stack (6)" region="10041 -10043" universe="10089" />
+    <cell fill="10084" id="10631" name="BA stack (7)" region="10043 -10044" universe="10089" />
+    <cell fill="10081" id="10632" name="BA stack (8)" region="10044 -10045" universe="10089" />
+    <cell fill="10083" id="10633" name="BA stack (9)" region="10045 -10046" universe="10089" />
+    <cell fill="10081" id="10634" name="BA stack (10)" region="10046 -10047" universe="10089" />
+    <cell fill="10083" id="10635" name="BA stack (11)" region="10047 -10048" universe="10089" />
+    <cell fill="10081" id="10636" name="BA stack (12)" region="10048 -10049" universe="10089" />
+    <cell fill="10087" id="10637" name="BA stack (13)" region="10049 -10050" universe="10089" />
+    <cell fill="10087" id="10638" name="BA stack (14)" region="10050 -10040" universe="10089" />
+    <cell fill="10087" id="10639" name="BA stack (15)" region="10040 -10051" universe="10089" />
+    <cell fill="10087" id="10640" name="BA stack (16)" region="10051 -10052" universe="10089" />
+    <cell fill="10087" id="10641" name="BA stack (17)" region="10052 -10053" universe="10089" />
+    <cell fill="10088" id="10642" name="BA stack (18)" region="10053 -10054" universe="10089" />
+    <cell fill="10000" id="10643" name="BA stack (last)" region="10054" universe="10089" />
+    <cell id="10644" material="10003" name="SS pin (0)" region="-10003" universe="10090" />
+    <cell id="10645" material="10007" name="SS pin (last)" region="10003" universe="10090" />
+    <cell id="10646" material="10005" name="end plug (0)" region="-10003" universe="10091" />
+    <cell id="10647" material="10007" name="end plug (last)" region="10003" universe="10091" />
+    <cell id="10648" material="10002" name="pin plenum (0)" region="-10001" universe="10092" />
+    <cell id="10649" material="10000" name="pin plenum (1)" region="10001 -10002" universe="10092" />
+    <cell id="10650" material="10005" name="pin plenum (2)" region="10002 -10003" universe="10092" />
+    <cell id="10651" material="10007" name="pin plenum (last)" region="10003" universe="10092" />
+    <cell id="10652" material="10002" name="pin plenum grid (top/bottom) (0)" region="-10001" universe="10093" />
+    <cell id="10653" material="10000" name="pin plenum grid (top/bottom) (1)" region="10001 -10002" universe="10093" />
+    <cell id="10654" material="10005" name="pin plenum grid (top/bottom) (2)" region="10002 -10003" universe="10093" />
+    <cell id="10655" material="10007" name="pin plenum grid (top/bottom) (last)" region="10003 (10019 -10020 10021 -10022)" universe="10093" />
+    <cell id="10656" material="10005" name="pin plenum grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10093" />
+    <cell id="10657" material="10009" name="Fuel (1.6\%) (0)" region="-10000" universe="10094" />
+    <cell id="10658" material="10000" name="Fuel (1.6\%) (1)" region="10000 -10002" universe="10094" />
+    <cell id="10659" material="10005" name="Fuel (1.6\%) (2)" region="10002 -10003" universe="10094" />
+    <cell id="10660" material="10007" name="Fuel (1.6\%) (last)" region="10003" universe="10094" />
+    <cell id="10661" material="10009" name="Fuel (1.6\%) grid (top/bottom) (0)" region="-10000" universe="10095" />
+    <cell id="10662" material="10000" name="Fuel (1.6\%) grid (top/bottom) (1)" region="10000 -10002" universe="10095" />
+    <cell id="10663" material="10005" name="Fuel (1.6\%) grid (top/bottom) (2)" region="10002 -10003" universe="10095" />
+    <cell id="10664" material="10007" name="Fuel (1.6\%) grid (top/bottom) (last)" region="10003 (10019 -10020 10021 -10022)" universe="10095" />
+    <cell id="10665" material="10005" name="Fuel (1.6\%) grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10095" />
+    <cell id="10666" material="10009" name="Fuel (1.6\%) grid (intermediate) (0)" region="-10000" universe="10096" />
+    <cell id="10667" material="10000" name="Fuel (1.6\%) grid (intermediate) (1)" region="10000 -10002" universe="10096" />
+    <cell id="10668" material="10005" name="Fuel (1.6\%) grid (intermediate) (2)" region="10002 -10003" universe="10096" />
+    <cell id="10669" material="10007" name="Fuel (1.6\%) grid (intermediate) (last)" region="10003 (10023 -10024 10025 -10026)" universe="10096" />
+    <cell id="10670" material="10005" name="Fuel (1.6\%) grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10096" />
+    <cell fill="10000" id="10671" name="Fuel (1.6\%) stack (0)" region="-10036" universe="10097" />
+    <cell fill="10090" id="10672" name="Fuel (1.6\%) stack (1)" region="10036 -10037" universe="10097" />
+    <cell fill="10090" id="10673" name="Fuel (1.6\%) stack (2)" region="10037 -10038" universe="10097" />
+    <cell fill="10091" id="10674" name="Fuel (1.6\%) stack (3)" region="10038 -10039" universe="10097" />
+    <cell fill="10094" id="10675" name="Fuel (1.6\%) stack (4)" region="10039 -10042" universe="10097" />
+    <cell fill="10095" id="10676" name="Fuel (1.6\%) stack (5)" region="10042 -10043" universe="10097" />
+    <cell fill="10094" id="10677" name="Fuel (1.6\%) stack (6)" region="10043 -10044" universe="10097" />
+    <cell fill="10094" id="10678" name="Fuel (1.6\%) stack (7)" region="10044 -10045" universe="10097" />
+    <cell fill="10096" id="10679" name="Fuel (1.6\%) stack (8)" region="10045 -10046" universe="10097" />
+    <cell fill="10094" id="10680" name="Fuel (1.6\%) stack (9)" region="10046 -10047" universe="10097" />
+    <cell fill="10096" id="10681" name="Fuel (1.6\%) stack (10)" region="10047 -10048" universe="10097" />
+    <cell fill="10094" id="10682" name="Fuel (1.6\%) stack (11)" region="10048 -10049" universe="10097" />
+    <cell fill="10092" id="10683" name="Fuel (1.6\%) stack (12)" region="10049 -10050" universe="10097" />
+    <cell fill="10093" id="10684" name="Fuel (1.6\%) stack (13)" region="10050 -10040" universe="10097" />
+    <cell fill="10092" id="10685" name="Fuel (1.6\%) stack (14)" region="10040 -10051" universe="10097" />
+    <cell fill="10091" id="10686" name="Fuel (1.6\%) stack (15)" region="10051 -10052" universe="10097" />
+    <cell fill="10000" id="10687" name="Fuel (1.6\%) stack (16)" region="10052 -10053" universe="10097" />
+    <cell fill="10090" id="10688" name="Fuel (1.6\%) stack (17)" region="10053 -10054" universe="10097" />
+    <cell fill="10000" id="10689" name="Fuel (1.6\%) stack (last)" region="10054" universe="10097" />
+    <cell id="10690" material="10010" name="Fuel (2.4\%) (0)" region="-10000" universe="10098" />
+    <cell id="10691" material="10000" name="Fuel (2.4\%) (1)" region="10000 -10002" universe="10098" />
+    <cell id="10692" material="10005" name="Fuel (2.4\%) (2)" region="10002 -10003" universe="10098" />
+    <cell id="10693" material="10007" name="Fuel (2.4\%) (last)" region="10003" universe="10098" />
+    <cell id="10694" material="10010" name="Fuel (2.4\%) grid (top/bottom) (0)" region="-10000" universe="10099" />
+    <cell id="10695" material="10000" name="Fuel (2.4\%) grid (top/bottom) (1)" region="10000 -10002" universe="10099" />
+    <cell id="10696" material="10005" name="Fuel (2.4\%) grid (top/bottom) (2)" region="10002 -10003" universe="10099" />
+    <cell id="10697" material="10007" name="Fuel (2.4\%) grid (top/bottom) (last)" region="10003 (10019 -10020 10021 -10022)" universe="10099" />
+    <cell id="10698" material="10005" name="Fuel (2.4\%) grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10099" />
+    <cell id="10699" material="10010" name="Fuel (2.4\%) grid (intermediate) (0)" region="-10000" universe="10100" />
+    <cell id="10700" material="10000" name="Fuel (2.4\%) grid (intermediate) (1)" region="10000 -10002" universe="10100" />
+    <cell id="10701" material="10005" name="Fuel (2.4\%) grid (intermediate) (2)" region="10002 -10003" universe="10100" />
+    <cell id="10702" material="10007" name="Fuel (2.4\%) grid (intermediate) (last)" region="10003 (10023 -10024 10025 -10026)" universe="10100" />
+    <cell id="10703" material="10005" name="Fuel (2.4\%) grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10100" />
+    <cell fill="10000" id="10704" name="Fuel (2.4\%) stack (0)" region="-10036" universe="10101" />
+    <cell fill="10090" id="10705" name="Fuel (2.4\%) stack (1)" region="10036 -10037" universe="10101" />
+    <cell fill="10090" id="10706" name="Fuel (2.4\%) stack (2)" region="10037 -10038" universe="10101" />
+    <cell fill="10091" id="10707" name="Fuel (2.4\%) stack (3)" region="10038 -10039" universe="10101" />
+    <cell fill="10098" id="10708" name="Fuel (2.4\%) stack (4)" region="10039 -10042" universe="10101" />
+    <cell fill="10099" id="10709" name="Fuel (2.4\%) stack (5)" region="10042 -10043" universe="10101" />
+    <cell fill="10098" id="10710" name="Fuel (2.4\%) stack (6)" region="10043 -10044" universe="10101" />
+    <cell fill="10098" id="10711" name="Fuel (2.4\%) stack (7)" region="10044 -10045" universe="10101" />
+    <cell fill="10100" id="10712" name="Fuel (2.4\%) stack (8)" region="10045 -10046" universe="10101" />
+    <cell fill="10098" id="10713" name="Fuel (2.4\%) stack (9)" region="10046 -10047" universe="10101" />
+    <cell fill="10100" id="10714" name="Fuel (2.4\%) stack (10)" region="10047 -10048" universe="10101" />
+    <cell fill="10098" id="10715" name="Fuel (2.4\%) stack (11)" region="10048 -10049" universe="10101" />
+    <cell fill="10092" id="10716" name="Fuel (2.4\%) stack (12)" region="10049 -10050" universe="10101" />
+    <cell fill="10093" id="10717" name="Fuel (2.4\%) stack (13)" region="10050 -10040" universe="10101" />
+    <cell fill="10092" id="10718" name="Fuel (2.4\%) stack (14)" region="10040 -10051" universe="10101" />
+    <cell fill="10091" id="10719" name="Fuel (2.4\%) stack (15)" region="10051 -10052" universe="10101" />
+    <cell fill="10000" id="10720" name="Fuel (2.4\%) stack (16)" region="10052 -10053" universe="10101" />
+    <cell fill="10090" id="10721" name="Fuel (2.4\%) stack (17)" region="10053 -10054" universe="10101" />
+    <cell fill="10000" id="10722" name="Fuel (2.4\%) stack (last)" region="10054" universe="10101" />
+    <cell id="10723" material="10011" name="Fuel (3.1\%) (0)" region="-10000" universe="10102" />
+    <cell id="10724" material="10000" name="Fuel (3.1\%) (1)" region="10000 -10002" universe="10102" />
+    <cell id="10725" material="10005" name="Fuel (3.1\%) (2)" region="10002 -10003" universe="10102" />
+    <cell id="10726" material="10007" name="Fuel (3.1\%) (last)" region="10003" universe="10102" />
+    <cell id="10727" material="10011" name="Fuel (3.1\%) grid (top/bottom) (0)" region="-10000" universe="10103" />
+    <cell id="10728" material="10000" name="Fuel (3.1\%) grid (top/bottom) (1)" region="10000 -10002" universe="10103" />
+    <cell id="10729" material="10005" name="Fuel (3.1\%) grid (top/bottom) (2)" region="10002 -10003" universe="10103" />
+    <cell id="10730" material="10007" name="Fuel (3.1\%) grid (top/bottom) (last)" region="10003 (10019 -10020 10021 -10022)" universe="10103" />
+    <cell id="10731" material="10005" name="Fuel (3.1\%) grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10103" />
+    <cell id="10732" material="10011" name="Fuel (3.1\%) grid (intermediate) (0)" region="-10000" universe="10104" />
+    <cell id="10733" material="10000" name="Fuel (3.1\%) grid (intermediate) (1)" region="10000 -10002" universe="10104" />
+    <cell id="10734" material="10005" name="Fuel (3.1\%) grid (intermediate) (2)" region="10002 -10003" universe="10104" />
+    <cell id="10735" material="10007" name="Fuel (3.1\%) grid (intermediate) (last)" region="10003 (10023 -10024 10025 -10026)" universe="10104" />
+    <cell id="10736" material="10005" name="Fuel (3.1\%) grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10104" />
+    <cell fill="10000" id="10737" name="Fuel (3.1\%) stack (0)" region="-10036" universe="10105" />
+    <cell fill="10090" id="10738" name="Fuel (3.1\%) stack (1)" region="10036 -10037" universe="10105" />
+    <cell fill="10090" id="10739" name="Fuel (3.1\%) stack (2)" region="10037 -10038" universe="10105" />
+    <cell fill="10091" id="10740" name="Fuel (3.1\%) stack (3)" region="10038 -10039" universe="10105" />
+    <cell fill="10102" id="10741" name="Fuel (3.1\%) stack (4)" region="10039 -10042" universe="10105" />
+    <cell fill="10103" id="10742" name="Fuel (3.1\%) stack (5)" region="10042 -10043" universe="10105" />
+    <cell fill="10102" id="10743" name="Fuel (3.1\%) stack (6)" region="10043 -10044" universe="10105" />
+    <cell fill="10102" id="10744" name="Fuel (3.1\%) stack (7)" region="10044 -10045" universe="10105" />
+    <cell fill="10104" id="10745" name="Fuel (3.1\%) stack (8)" region="10045 -10046" universe="10105" />
+    <cell fill="10102" id="10746" name="Fuel (3.1\%) stack (9)" region="10046 -10047" universe="10105" />
+    <cell fill="10104" id="10747" name="Fuel (3.1\%) stack (10)" region="10047 -10048" universe="10105" />
+    <cell fill="10102" id="10748" name="Fuel (3.1\%) stack (11)" region="10048 -10049" universe="10105" />
+    <cell fill="10092" id="10749" name="Fuel (3.1\%) stack (12)" region="10049 -10050" universe="10105" />
+    <cell fill="10093" id="10750" name="Fuel (3.1\%) stack (13)" region="10050 -10040" universe="10105" />
+    <cell fill="10092" id="10751" name="Fuel (3.1\%) stack (14)" region="10040 -10051" universe="10105" />
+    <cell fill="10091" id="10752" name="Fuel (3.1\%) stack (15)" region="10051 -10052" universe="10105" />
+    <cell fill="10000" id="10753" name="Fuel (3.1\%) stack (16)" region="10052 -10053" universe="10105" />
+    <cell fill="10090" id="10754" name="Fuel (3.1\%) stack (17)" region="10053 -10054" universe="10105" />
+    <cell fill="10000" id="10755" name="Fuel (3.1\%) stack (last)" region="10054" universe="10105" />
+    <cell fill="10110" id="10778" name="Assembly (1.6\%) CR B lattice" region="10027 -10028 10029 -10030" universe="10111" />
+    <cell fill="10000" id="10779" name="Assembly (1.6\%) CR B lattice outer water" region="~(10031 -10032 10033 -10034)" universe="10111" />
+    <cell id="10780" material="10007" name="Assembly (1.6\%) CR B lattice axial (0)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) -10042" universe="10111" />
+    <cell id="10781" material="10003" name="Assembly (1.6\%) CR B lattice axial (1)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10042 -10043" universe="10111" />
+    <cell id="10782" material="10007" name="Assembly (1.6\%) CR B lattice axial (2)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10043 -10045" universe="10111" />
+    <cell id="10783" material="10005" name="Assembly (1.6\%) CR B lattice axial (3)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10045 -10046" universe="10111" />
+    <cell id="10784" material="10007" name="Assembly (1.6\%) CR B lattice axial (4)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10046 -10047" universe="10111" />
+    <cell id="10785" material="10005" name="Assembly (1.6\%) CR B lattice axial (5)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10047 -10048" universe="10111" />
+    <cell id="10786" material="10007" name="Assembly (1.6\%) CR B lattice axial (6)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10048 -10049" universe="10111" />
+    <cell id="10787" material="10003" name="Assembly (1.6\%) CR B lattice axial (7)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10049 -10050" universe="10111" />
+    <cell id="10788" material="10007" name="Assembly (1.6\%) CR B lattice axial (last)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10050" universe="10111" />
+    <cell fill="10124" id="10855" name="Assembly (1.6\%) no BAs instr lattice" region="10027 -10028 10029 -10030" universe="10125" />
+    <cell fill="10000" id="10856" name="Assembly (1.6\%) no BAs instr lattice outer water" region="~(10031 -10032 10033 -10034)" universe="10125" />
+    <cell id="10857" material="10007" name="Assembly (1.6\%) no BAs instr lattice axial (0)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) -10042" universe="10125" />
+    <cell id="10858" material="10003" name="Assembly (1.6\%) no BAs instr lattice axial (1)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10042 -10043" universe="10125" />
+    <cell id="10859" material="10007" name="Assembly (1.6\%) no BAs instr lattice axial (2)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10043 -10045" universe="10125" />
+    <cell id="10860" material="10005" name="Assembly (1.6\%) no BAs instr lattice axial (3)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10045 -10046" universe="10125" />
+    <cell id="10861" material="10007" name="Assembly (1.6\%) no BAs instr lattice axial (4)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10046 -10047" universe="10125" />
+    <cell id="10862" material="10005" name="Assembly (1.6\%) no BAs instr lattice axial (5)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10047 -10048" universe="10125" />
+    <cell id="10863" material="10007" name="Assembly (1.6\%) no BAs instr lattice axial (6)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10048 -10049" universe="10125" />
+    <cell id="10864" material="10003" name="Assembly (1.6\%) no BAs instr lattice axial (7)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10049 -10050" universe="10125" />
+    <cell id="10865" material="10007" name="Assembly (1.6\%) no BAs instr lattice axial (last)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10050" universe="10125" />
+    <cell fill="10144" id="10965" name="Assembly (2.4\%) CR D lattice" region="10027 -10028 10029 -10030" universe="10145" />
+    <cell fill="10000" id="10966" name="Assembly (2.4\%) CR D lattice outer water" region="~(10031 -10032 10033 -10034)" universe="10145" />
+    <cell id="10967" material="10007" name="Assembly (2.4\%) CR D lattice axial (0)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) -10042" universe="10145" />
+    <cell id="10968" material="10003" name="Assembly (2.4\%) CR D lattice axial (1)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10042 -10043" universe="10145" />
+    <cell id="10969" material="10007" name="Assembly (2.4\%) CR D lattice axial (2)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10043 -10045" universe="10145" />
+    <cell id="10970" material="10005" name="Assembly (2.4\%) CR D lattice axial (3)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10045 -10046" universe="10145" />
+    <cell id="10971" material="10007" name="Assembly (2.4\%) CR D lattice axial (4)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10046 -10047" universe="10145" />
+    <cell id="10972" material="10005" name="Assembly (2.4\%) CR D lattice axial (5)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10047 -10048" universe="10145" />
+    <cell id="10973" material="10007" name="Assembly (2.4\%) CR D lattice axial (6)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10048 -10049" universe="10145" />
+    <cell id="10974" material="10003" name="Assembly (2.4\%) CR D lattice axial (7)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10049 -10050" universe="10145" />
+    <cell id="10975" material="10007" name="Assembly (2.4\%) CR D lattice axial (last)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10050" universe="10145" />
+    <cell fill="10164" id="11075" name="Assembly (3.1\%) 16BA lattice" region="10027 -10028 10029 -10030" universe="10165" />
+    <cell fill="10000" id="11076" name="Assembly (3.1\%) 16BA lattice outer water" region="~(10031 -10032 10033 -10034)" universe="10165" />
+    <cell id="11077" material="10007" name="Assembly (3.1\%) 16BA lattice axial (0)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) -10042" universe="10165" />
+    <cell id="11078" material="10003" name="Assembly (3.1\%) 16BA lattice axial (1)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10042 -10043" universe="10165" />
+    <cell id="11079" material="10007" name="Assembly (3.1\%) 16BA lattice axial (2)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10043 -10045" universe="10165" />
+    <cell id="11080" material="10005" name="Assembly (3.1\%) 16BA lattice axial (3)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10045 -10046" universe="10165" />
+    <cell id="11081" material="10007" name="Assembly (3.1\%) 16BA lattice axial (4)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10046 -10047" universe="10165" />
+    <cell id="11082" material="10005" name="Assembly (3.1\%) 16BA lattice axial (5)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10047 -10048" universe="10165" />
+    <cell id="11083" material="10007" name="Assembly (3.1\%) 16BA lattice axial (6)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10048 -10049" universe="10165" />
+    <cell id="11084" material="10003" name="Assembly (3.1\%) 16BA lattice axial (7)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10049 -10050" universe="10165" />
+    <cell id="11085" material="10007" name="Assembly (3.1\%) 16BA lattice axial (last)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10050" universe="10165" />
+    <cell fill="10182" id="11174" name="Assembly (3.1\%) no BAs instr lattice" region="10027 -10028 10029 -10030" universe="10183" />
+    <cell fill="10000" id="11175" name="Assembly (3.1\%) no BAs instr lattice outer water" region="~(10031 -10032 10033 -10034)" universe="10183" />
+    <cell id="11176" material="10007" name="Assembly (3.1\%) no BAs instr lattice axial (0)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) -10042" universe="10183" />
+    <cell id="11177" material="10003" name="Assembly (3.1\%) no BAs instr lattice axial (1)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10042 -10043" universe="10183" />
+    <cell id="11178" material="10007" name="Assembly (3.1\%) no BAs instr lattice axial (2)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10043 -10045" universe="10183" />
+    <cell id="11179" material="10005" name="Assembly (3.1\%) no BAs instr lattice axial (3)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10045 -10046" universe="10183" />
+    <cell id="11180" material="10007" name="Assembly (3.1\%) no BAs instr lattice axial (4)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10046 -10047" universe="10183" />
+    <cell id="11181" material="10005" name="Assembly (3.1\%) no BAs instr lattice axial (5)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10047 -10048" universe="10183" />
+    <cell id="11182" material="10007" name="Assembly (3.1\%) no BAs instr lattice axial (6)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10048 -10049" universe="10183" />
+    <cell id="11183" material="10003" name="Assembly (3.1\%) no BAs instr lattice axial (7)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10049 -10050" universe="10183" />
+    <cell id="11184" material="10007" name="Assembly (3.1\%) no BAs instr lattice axial (last)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10050" universe="10183" />
+    <cell id="11306" material="10003" name="baffle dummy SS" region="10085" universe="10206" />
+    <cell id="11307" material="10007" name="baffle north water" region="-10085" universe="10206" />
+    <cell fill="10207" id="11308" name="baffle north" universe="10208" />
+    <cell id="11309" material="10007" name="baffle dummy SS" region="10086" universe="10209" />
+    <cell id="11310" material="10003" name="baffle south water" region="-10086" universe="10209" />
+    <cell fill="10210" id="11311" name="baffle south" universe="10211" />
+    <cell id="11312" material="10003" name="baffle dummy SS" region="10087" universe="10212" />
+    <cell id="11313" material="10007" name="baffle east water" region="-10087" universe="10212" />
+    <cell fill="10213" id="11314" name="baffle east" universe="10214" />
+    <cell id="11315" material="10007" name="baffle dummy SS" region="10088" universe="10215" />
+    <cell id="11316" material="10003" name="baffle west water" region="-10088" universe="10215" />
+    <cell fill="10216" id="11317" name="baffle west" universe="10217" />
+    <cell id="11318" material="10007" name="baffle northwest dummy 1" region="10088 -10085" universe="10218" />
+    <cell id="11319" material="10003" name="baffle northwest dummy 2" region="10088 10085" universe="10218" />
+    <cell id="11320" material="10003" name="baffle northwest dummy 3" region="-10088" universe="10218" />
+    <cell fill="10219" id="11321" name="baffle northwest" universe="10220" />
+    <cell id="11322" material="10007" name="baffle northeast dummy 1" region="-10085 -10087" universe="10221" />
+    <cell id="11323" material="10003" name="baffle northeast dummy 2" region="10085 -10087" universe="10221" />
+    <cell id="11324" material="10003" name="baffle northeast dummy 3" region="10087" universe="10221" />
+    <cell fill="10222" id="11325" name="baffle northeast" universe="10223" />
+    <cell id="11326" material="10007" name="baffle southwest dummy 1" region="10086 10088" universe="10224" />
+    <cell id="11327" material="10003" name="baffle southwest dummy 2" region="-10086 10088" universe="10224" />
+    <cell id="11328" material="10003" name="baffle southwest dummy 3" region="-10088" universe="10224" />
+    <cell fill="10225" id="11329" name="baffle southwest" universe="10226" />
+    <cell id="11330" material="10007" name="baffle southeast dummy 1" region="10086 -10087" universe="10227" />
+    <cell id="11331" material="10003" name="baffle southeast dummy 2" region="-10086 -10087" universe="10227" />
+    <cell id="11332" material="10003" name="baffle southeast dummy 3" region="10087" universe="10227" />
+    <cell fill="10228" id="11333" name="baffle southeast" universe="10229" />
+    <cell id="11334" material="10007" name="baffle northwest corner dummy 1" region="-10088 -10085" universe="10230" />
+    <cell id="11335" material="10007" name="baffle northwest corner dummy 2" region="10088" universe="10230" />
+    <cell id="11336" material="10003" name="baffle northwest corner dummy 3" region="10085 -10088" universe="10230" />
+    <cell fill="10231" id="11337" name="baffle northwest corner" universe="10232" />
+    <cell id="11338" material="10007" name="baffle northeast corner dummy 1" region="10087 -10085" universe="10233" />
+    <cell id="11339" material="10007" name="baffle northeast corner dummy 2" region="-10087" universe="10233" />
+    <cell id="11340" material="10003" name="baffle northeast corner dummy 3" region="10085 10087" universe="10233" />
+    <cell fill="10234" id="11341" name="baffle northeast corner" universe="10235" />
+    <cell id="11342" material="10007" name="baffle southeast corner dummy 1" region="10087 10086" universe="10236" />
+    <cell id="11343" material="10007" name="baffle southeast corner dummy 2" region="-10087" universe="10236" />
+    <cell id="11344" material="10003" name="baffle southeast corner dummy 3" region="-10086 10087" universe="10236" />
+    <cell fill="10237" id="11345" name="baffle southeast corner" universe="10238" />
+    <cell id="11346" material="10007" name="baffle southwest corner dummy 1" region="-10088 10086" universe="10239" />
+    <cell id="11347" material="10007" name="baffle southwest corner dummy 2" region="10088" universe="10239" />
+    <cell id="11348" material="10003" name="baffle southwest corner dummy 3" region="-10086 -10088" universe="10239" />
+    <cell fill="10240" id="11349" name="baffle southwest corner" universe="10241" />
+    <cell fill="10242" id="11350" name="Main core" region="-10074 10084 -10083" universe="0" />
+    <cell id="11351" material="10003" name="core barrel" region="10074 -10075 10084 -10083" universe="0" />
+    <cell id="11352" material="10003" name="neutron shield panel NW" region="10075 -10076 10077 -10078 10084 -10083" universe="0" />
+    <cell id="11353" material="10007" name="neutron shield panel N" region="10075 -10076 10078 -10080 10084 -10083" universe="0" />
+    <cell id="11354" material="10003" name="neutron shield panel SE" region="10075 -10076 -10077 10078 10084 -10083" universe="0" />
+    <cell id="11355" material="10007" name="neutron shield panel E" region="10075 -10076 10077 10079 10084 -10083" universe="0" />
+    <cell id="11356" material="10003" name="neutron shield panel NE" region="10075 -10076 10079 -10080 10084 -10083" universe="0" />
+    <cell id="11357" material="10007" name="neutron shield panel S" region="10075 -10076 -10078 10080 10084 -10083" universe="0" />
+    <cell id="11358" material="10003" name="neutron shield panel SW" region="10075 -10076 -10079 10080 10084 -10083" universe="0" />
+    <cell id="11359" material="10007" name="neutron shield panel W" region="10075 -10076 -10077 -10079 10084 -10083" universe="0" />
+    <cell id="11360" material="10007" name="downcomer" region="10076 -10081 10084 -10083" universe="0" />
+    <cell id="11361" material="10004" name="reactor pressure vessel" region="10081 -10082 10084 -10083" universe="0" />
+    <lattice id="10110" name="Assembly (1.6\%) CR B">
+        <pitch>1.25984 1.25984</pitch>
+        <dimension>17 17</dimension>
+        <lower_left>-10.70864 -10.70864</lower_left>
+        <universes>
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10097 10097 10097 10038 10097 10097 10038 10097 10097 10038 10097 10097 10097 10097 10097 
+10097 10097 10097 10038 10097 10097 10097 10097 10097 10097 10097 10097 10097 10038 10097 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10038 10097 10097 10038 10097 10097 10038 10097 10097 10038 10097 10097 10038 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10038 10097 10097 10038 10097 10097 10010 10097 10097 10038 10097 10097 10038 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10038 10097 10097 10038 10097 10097 10038 10097 10097 10038 10097 10097 10038 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10097 10038 10097 10097 10097 10097 10097 10097 10097 10097 10097 10038 10097 10097 10097 
+10097 10097 10097 10097 10097 10038 10097 10097 10038 10097 10097 10038 10097 10097 10097 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 </universes>
+    </lattice>
+    <lattice id="10124" name="Assembly (1.6\%) no BAs instr">
+        <pitch>1.25984 1.25984</pitch>
+        <dimension>17 17</dimension>
+        <lower_left>-10.70864 -10.70864</lower_left>
+        <universes>
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10097 10097 10097 10009 10097 10097 10009 10097 10097 10009 10097 10097 10097 10097 10097 
+10097 10097 10097 10009 10097 10097 10097 10097 10097 10097 10097 10097 10097 10009 10097 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10009 10097 10097 10009 10097 10097 10009 10097 10097 10009 10097 10097 10009 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10009 10097 10097 10009 10097 10097 10016 10097 10097 10009 10097 10097 10009 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10009 10097 10097 10009 10097 10097 10009 10097 10097 10009 10097 10097 10009 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10097 10009 10097 10097 10097 10097 10097 10097 10097 10097 10097 10009 10097 10097 10097 
+10097 10097 10097 10097 10097 10009 10097 10097 10009 10097 10097 10009 10097 10097 10097 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 
+10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 </universes>
+    </lattice>
+    <lattice id="10144" name="Assembly (2.4\%) CR D">
+        <pitch>1.25984 1.25984</pitch>
+        <dimension>17 17</dimension>
+        <lower_left>-10.70864 -10.70864</lower_left>
+        <universes>
+10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 
+10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 
+10101 10101 10101 10101 10101 10050 10101 10101 10050 10101 10101 10050 10101 10101 10101 10101 10101 
+10101 10101 10101 10050 10101 10101 10101 10101 10101 10101 10101 10101 10101 10050 10101 10101 10101 
+10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 
+10101 10101 10050 10101 10101 10050 10101 10101 10050 10101 10101 10050 10101 10101 10050 10101 10101 
+10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 
+10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 
+10101 10101 10050 10101 10101 10050 10101 10101 10010 10101 10101 10050 10101 10101 10050 10101 10101 
+10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 
+10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 
+10101 10101 10050 10101 10101 10050 10101 10101 10050 10101 10101 10050 10101 10101 10050 10101 10101 
+10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 
+10101 10101 10101 10050 10101 10101 10101 10101 10101 10101 10101 10101 10101 10050 10101 10101 10101 
+10101 10101 10101 10101 10101 10050 10101 10101 10050 10101 10101 10050 10101 10101 10101 10101 10101 
+10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 
+10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 </universes>
+    </lattice>
+    <lattice id="10164" name="Assembly (3.1\%) 16BA">
+        <pitch>1.25984 1.25984</pitch>
+        <dimension>17 17</dimension>
+        <lower_left>-10.70864 -10.70864</lower_left>
+        <universes>
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10105 10105 10105 10089 10105 10105 10089 10105 10105 10089 10105 10105 10105 10105 10105 
+10105 10105 10105 10089 10105 10105 10105 10105 10105 10105 10105 10105 10105 10089 10105 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10089 10105 10105 10009 10105 10105 10009 10105 10105 10009 10105 10105 10089 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10089 10105 10105 10009 10105 10105 10010 10105 10105 10009 10105 10105 10089 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10089 10105 10105 10009 10105 10105 10009 10105 10105 10009 10105 10105 10089 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10105 10089 10105 10105 10105 10105 10105 10105 10105 10105 10105 10089 10105 10105 10105 
+10105 10105 10105 10105 10105 10089 10105 10105 10089 10105 10105 10089 10105 10105 10105 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 </universes>
+    </lattice>
+    <lattice id="10182" name="Assembly (3.1\%) no BAs instr">
+        <pitch>1.25984 1.25984</pitch>
+        <dimension>17 17</dimension>
+        <lower_left>-10.70864 -10.70864</lower_left>
+        <universes>
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10105 10105 10105 10009 10105 10105 10009 10105 10105 10009 10105 10105 10105 10105 10105 
+10105 10105 10105 10009 10105 10105 10105 10105 10105 10105 10105 10105 10105 10009 10105 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10009 10105 10105 10009 10105 10105 10009 10105 10105 10009 10105 10105 10009 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10009 10105 10105 10009 10105 10105 10016 10105 10105 10009 10105 10105 10009 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10009 10105 10105 10009 10105 10105 10009 10105 10105 10009 10105 10105 10009 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10105 10009 10105 10105 10105 10105 10105 10105 10105 10105 10105 10009 10105 10105 10105 
+10105 10105 10105 10105 10105 10009 10105 10105 10009 10105 10105 10009 10105 10105 10105 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 
+10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 </universes>
+    </lattice>
+    <lattice id="10207" name="baffle north">
+        <pitch>10.75182 10.75182</pitch>
+        <dimension>2 2</dimension>
+        <lower_left>-10.75182 -10.75182</lower_left>
+        <universes>
+10206 10206 
+10000 10000 </universes>
+    </lattice>
+    <lattice id="10210" name="baffle south">
+        <pitch>10.75182 10.75182</pitch>
+        <dimension>2 2</dimension>
+        <lower_left>-10.75182 -10.75182</lower_left>
+        <universes>
+10000 10000 
+10209 10209 </universes>
+    </lattice>
+    <lattice id="10213" name="baffle east">
+        <pitch>10.75182 10.75182</pitch>
+        <dimension>2 2</dimension>
+        <lower_left>-10.75182 -10.75182</lower_left>
+        <universes>
+10000 10212 
+10000 10212 </universes>
+    </lattice>
+    <lattice id="10216" name="baffle west">
+        <pitch>10.75182 10.75182</pitch>
+        <dimension>2 2</dimension>
+        <lower_left>-10.75182 -10.75182</lower_left>
+        <universes>
+10215 10000 
+10215 10000 </universes>
+    </lattice>
+    <lattice id="10219" name="baffle northwest">
+        <pitch>10.75182 10.75182</pitch>
+        <dimension>2 2</dimension>
+        <lower_left>-10.75182 -10.75182</lower_left>
+        <universes>
+10218 10206 
+10215 10000 </universes>
+    </lattice>
+    <lattice id="10222" name="baffle northeast">
+        <pitch>10.75182 10.75182</pitch>
+        <dimension>2 2</dimension>
+        <lower_left>-10.75182 -10.75182</lower_left>
+        <universes>
+10206 10221 
+10000 10212 </universes>
+    </lattice>
+    <lattice id="10225" name="baffle southwest">
+        <pitch>10.75182 10.75182</pitch>
+        <dimension>2 2</dimension>
+        <lower_left>-10.75182 -10.75182</lower_left>
+        <universes>
+10215 10000 
+10224 10209 </universes>
+    </lattice>
+    <lattice id="10228" name="baffle southeast">
+        <pitch>10.75182 10.75182</pitch>
+        <dimension>2 2</dimension>
+        <lower_left>-10.75182 -10.75182</lower_left>
+        <universes>
+10000 10212 
+10209 10227 </universes>
+    </lattice>
+    <lattice id="10231" name="baffle northwest corner">
+        <pitch>10.75182 10.75182</pitch>
+        <dimension>2 2</dimension>
+        <lower_left>-10.75182 -10.75182</lower_left>
+        <universes>
+10230 10000 
+10000 10000 </universes>
+    </lattice>
+    <lattice id="10234" name="baffle northeast corner">
+        <pitch>10.75182 10.75182</pitch>
+        <dimension>2 2</dimension>
+        <lower_left>-10.75182 -10.75182</lower_left>
+        <universes>
+10000 10233 
+10000 10000 </universes>
+    </lattice>
+    <lattice id="10237" name="baffle southeast corner">
+        <pitch>10.75182 10.75182</pitch>
+        <dimension>2 2</dimension>
+        <lower_left>-10.75182 -10.75182</lower_left>
+        <universes>
+10000 10000 
+10000 10236 </universes>
+    </lattice>
+    <lattice id="10240" name="baffle southwest corner">
+        <pitch>10.75182 10.75182</pitch>
+        <dimension>2 2</dimension>
+        <lower_left>-10.75182 -10.75182</lower_left>
+        <universes>
+10000 10000 
+10239 10000 </universes>
+    </lattice>
+    <lattice id="10242" name="Main core">
+        <pitch>21.50364 21.50364</pitch>
+        <dimension>19 19</dimension>
+        <lower_left>-204.28458 -204.28458</lower_left>
+        <universes>
+10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10000 10000 10238 10211 10211 10211 10241 10000 10000 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10000 10238 10229 10183 10111 10183 10226 10241 10000 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10238 10229 10183 10111 10165 10111 10183 10226 10241 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10214 10183 10111 10165 10145 10165 10111 10183 10217 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10214 10111 10165 10145 10125 10145 10165 10111 10217 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10214 10183 10111 10165 10145 10165 10111 10183 10217 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10235 10223 10183 10111 10165 10111 10183 10220 10232 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10000 10235 10223 10183 10111 10183 10220 10232 10000 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10000 10000 10235 10208 10208 10208 10232 10000 10000 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 
+10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 </universes>
+    </lattice>
+    <surface coeffs="0.0 0.0 0.39218" id="10000" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.06459" id="10001" name="FR Plenum Spring OR" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.40005" id="10002" name="Clad IR" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.4572" id="10003" name="Clad OR" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.56134" id="10004" name="GT IR (above dashpot)" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.60198" id="10005" name="GT OR (above dashpot)" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.50419" id="10006" name="GT IR (at dashpot)" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.5461" id="10007" name="GT OR (at dashpot)" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.4331" id="10008" name="Control Poison OR" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.43688" id="10009" name="CR Clad IR" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.48387" id="10010" name="CR Clad OR" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.214" id="10011" name="BA IR 1" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.23051" id="10012" name="BA IR 2" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.2413" id="10013" name="BA IR 3" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.42672" id="10014" name="BA IR 4" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.43688" id="10015" name="BA IR 5" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.48387" id="10016" name="BA IR 6" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.56134" id="10017" name="BA IR 7" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 0.60198" id="10018" name="BA IR 8" type="z-cylinder" />
+    <surface coeffs="-0.62208" id="10019" name="minimum x" type="x-plane" />
+    <surface coeffs="0.62208" id="10020" name="maximum x" type="x-plane" />
+    <surface coeffs="-0.62208" id="10021" name="minimum y" type="y-plane" />
+    <surface coeffs="0.62208" id="10022" name="maximum y" type="y-plane" />
+    <surface coeffs="-0.60981" id="10023" name="minimum x" type="x-plane" />
+    <surface coeffs="0.60981" id="10024" name="maximum x" type="x-plane" />
+    <surface coeffs="-0.60981" id="10025" name="minimum y" type="y-plane" />
+    <surface coeffs="0.60981" id="10026" name="maximum y" type="y-plane" />
+    <surface coeffs="-10.70864" id="10027" name="minimum x" type="x-plane" />
+    <surface coeffs="10.70864" id="10028" name="maximum x" type="x-plane" />
+    <surface coeffs="-10.70864" id="10029" name="minimum y" type="y-plane" />
+    <surface coeffs="10.70864" id="10030" name="maximum y" type="y-plane" />
+    <surface coeffs="-10.73635" id="10031" name="minimum x" type="x-plane" />
+    <surface coeffs="10.73635" id="10032" name="maximum x" type="x-plane" />
+    <surface coeffs="-10.73635" id="10033" name="minimum y" type="y-plane" />
+    <surface coeffs="10.73635" id="10034" name="maximum y" type="y-plane" />
+    <surface coeffs="20.0" id="10036" name="bot support plate" type="z-plane" />
+    <surface coeffs="25.0" id="10037" name="top support plate" type="z-plane" />
+    <surface coeffs="35.16" id="10038" name="bottom FR" type="z-plane" />
+    <surface coeffs="36.007" id="10039" name="bot active core" type="z-plane" />
+    <surface coeffs="218.887" id="10040" name="top active core" type="z-plane" />
+    <surface coeffs="41.087" id="10041" name="bottom of BA" type="z-plane" />
+    <surface coeffs="37.879" id="10042" name="bottom grid 1" type="z-plane" />
+    <surface coeffs="42.07" id="10043" name="top of grid 1" type="z-plane" />
+    <surface coeffs="45.079" id="10044" name="top dashpot" type="z-plane" />
+    <surface coeffs="99.164" id="10045" name="bottom grid 2" type="z-plane" />
+    <surface coeffs="104.879" id="10046" name="top grid 2" type="z-plane" />
+    <surface coeffs="151.361" id="10047" name="bottom of grid 3" type="z-plane" />
+    <surface coeffs="157.076" id="10048" name="top of grid 3" type="z-plane" />
+    <surface coeffs="203.558" id="10049" name="bottom of grid 4" type="z-plane" />
+    <surface coeffs="209.273" id="10050" name="top grid 4" type="z-plane" />
+    <surface coeffs="221.223" id="10051" name="top pin plenum" type="z-plane" />
+    <surface coeffs="223.272" id="10052" name="top FR" type="z-plane" />
+    <surface coeffs="226.617" id="10053" name="bottom upper nozzle" type="z-plane" />
+    <surface coeffs="235.444" id="10054" name="top upper nozzle" type="z-plane" />
+    <surface coeffs="766.348" id="10068" name="CR bank B top" type="z-plane" />
+    <surface coeffs="405.713" id="10069" name="CR bank B bottom" type="z-plane" />
+    <surface coeffs="766.348" id="10072" name="CR bank D top" type="z-plane" />
+    <surface coeffs="405.713" id="10073" name="CR bank D bottom" type="z-plane" />
+    <surface coeffs="0.0 0.0 85.0" id="10074" name="core barrel IR" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 90.0" id="10075" name="core barrel OR" type="z-cylinder" />
+    <surface coeffs="0.0 0.0 92.0" id="10076" name="neutron shield OR" type="z-cylinder" />
+    <surface coeffs="1.0 1.7320508075688767 0.0 0.0" id="10077" name="neutron shield NWbot SEtop" type="plane" />
+    <surface coeffs="1.0 0.5773502691896257 0.0 0.0" id="10078" name="neutron shield NWtop SEbot" type="plane" />
+    <surface coeffs="1.0 -1.7320508075688767 0.0 0.0" id="10079" name="neutron shield NEbot SWtop" type="plane" />
+    <surface coeffs="1.0 -0.5773502691896257 0.0 0.0" id="10080" name="neutron shield NEtop SWbot" type="plane" />
+    <surface coeffs="0.0 0.0 120.0" id="10081" name="RPV IR" type="z-cylinder" />
+    <surface boundary="vacuum" coeffs="0.0 0.0 135.0" id="10082" name="RPV OR" type="z-cylinder" />
+    <surface boundary="vacuum" coeffs="255.444" id="10083" name="upper problem boundary" type="z-plane" />
+    <surface boundary="vacuum" coeffs="0.0" id="10084" name="lower problem boundary" type="z-plane" />
+    <surface coeffs="3.15341" id="10085" name="baffle north" type="y-plane" />
+    <surface coeffs="-3.15341" id="10086" name="baffle south" type="y-plane" />
+    <surface coeffs="3.15341" id="10087" name="baffle east" type="x-plane" />
+    <surface coeffs="-3.15341" id="10088" name="baffle west" type="x-plane" />
+</geometry>
diff --git a/smr/grids.py b/smr/grids.py
deleted file mode 100644
index d43bdef..0000000
--- a/smr/grids.py
+++ /dev/null
@@ -1,36 +0,0 @@
-from surfaces import surfs
-
-# FIXME: Add docstring explanation
-# FIXME: Put this in surface.py??
-
-# FIXME: Use get_rectangular_prism method instead???
-
-grids = {}
-
-grids['inner grid box (intermediate)'] = (+surfs['rod grid box ybot i'] &
-                                          -surfs['rod grid box ytop i'] &
-                                          +surfs['rod grid box xbot i'] &
-                                          -surfs['rod grid box xtop i'])
-
-grids['outer grid box (intermediate)'] = (-surfs['rod grid box ybot i']) | \
-                                         (+surfs['rod grid box ytop i']) | \
-                                         (+surfs['rod grid box xtop i'] &
-                                          -surfs['rod grid box ytop i'] &
-                                          +surfs['rod grid box ybot i']) | \
-                                         (-surfs['rod grid box xbot i'] &
-                                          -surfs['rod grid box ytop i'] &
-                                          +surfs['rod grid box ybot i'])
-
-grids['inner grid box (top/bottom)'] = (+surfs['rod grid box ybot tb'] &
-                                        -surfs['rod grid box ytop tb'] &
-                                        +surfs['rod grid box xbot tb'] &
-                                        -surfs['rod grid box xtop tb'])
-
-grids['outer grid box (top/bottom)'] = (-surfs['rod grid box ybot tb']) | \
-                                       (+surfs['rod grid box ytop tb']) | \
-                                       (+surfs['rod grid box xtop tb'] &
-                                        -surfs['rod grid box ytop tb'] &
-                                        +surfs['rod grid box ybot tb']) | \
-                                       (-surfs['rod grid box xbot tb'] &
-                                        -surfs['rod grid box ytop tb'] &
-                                        +surfs['rod grid box ybot tb'])
\ No newline at end of file
diff --git a/smr/materials.xml b/smr/materials.xml
new file mode 100644
index 0000000..d392417
--- /dev/null
+++ b/smr/materials.xml
@@ -0,0 +1,203 @@
+<?xml version='1.0' encoding='utf-8'?>
+<materials>
+    <material id="10000" name="Helium">
+        <temperature>300</temperature>
+        <density units="g/cc" value="0.0015981" />
+        <nuclide ao="1.34e-06" name="He3" />
+        <nuclide ao="0.99999866" name="He4" />
+    </material>
+    <material id="10001" name="Air">
+        <temperature>300</temperature>
+        <density units="g/cc" value="0.00616" />
+        <nuclide ao="0.20942038999999998" name="O16" />
+        <nuclide ao="7.961e-05" name="O17" />
+        <nuclide ao="0.7780575240000001" name="N14" />
+        <nuclide ao="0.002842476" name="N15" />
+        <nuclide ao="3.1124879999999996e-05" name="Ar36" />
+        <nuclide ao="5.86857e-06" name="Ar38" />
+        <nuclide ao="0.00929300655" name="Ar40" />
+        <nuclide ao="0.00027" name="C0" />
+    </material>
+    <material id="10002" name="Inconel">
+        <temperature>300</temperature>
+        <density units="g/cc" value="8.2" />
+        <nuclide name="Si28" wo="0.003215327037608538" />
+        <nuclide name="Si29" wo="0.0001691770106155053" />
+        <nuclide name="Si30" wo="0.00011549595177595707" />
+        <nuclide name="Cr50" wo="0.007913309553017726" />
+        <nuclide name="Cr52" wo="0.15869399115925625" />
+        <nuclide name="Cr53" wo="0.018341120727338085" />
+        <nuclide name="Cr54" wo="0.004651578560387908" />
+        <nuclide name="Mn55" wo="0.0087" />
+        <nuclide name="Fe54" wo="0.016163233201258693" />
+        <nuclide name="Fe56" wo="0.26311407687664323" />
+        <nuclide name="Fe57" wo="0.006185135350045743" />
+        <nuclide name="Fe58" wo="0.0008375545720523535" />
+        <nuclide name="Ni58" wo="0.3439855822298418" />
+        <nuclide name="Ni60" wo="0.13706564068516966" />
+        <nuclide name="Ni61" wo="0.006057615566737856" />
+        <nuclide name="Ni62" wo="0.01963099676371214" />
+        <nuclide name="Ni64" wo="0.0051601647545385385" />
+    </material>
+    <material id="10003" name="SS304">
+        <temperature>300</temperature>
+        <density units="g/cc" value="8.03" />
+        <nuclide name="Si28" wo="0.0055119892073289215" />
+        <nuclide name="Si29" wo="0.0002900177324837234" />
+        <nuclide name="Si30" wo="0.00019799306018735498" />
+        <nuclide name="Cr50" wo="0.007930004298910168" />
+        <nuclide name="Cr52" wo="0.15902878860895933" />
+        <nuclide name="Cr53" wo="0.018379815074864116" />
+        <nuclide name="Cr54" wo="0.004661392017266364" />
+        <nuclide name="Mn55" wo="0.02" />
+        <nuclide name="Fe54" wo="0.03861561826636726" />
+        <nuclide name="Fe56" wo="0.6286064568062312" />
+        <nuclide name="Fe57" wo="0.014776921339264018" />
+        <nuclide name="Fe58" wo="0.002001003588137652" />
+        <nuclide name="Ni58" wo="0.06719780860125842" />
+        <nuclide name="Ni60" wo="0.026775862606987628" />
+        <nuclide name="Ni61" wo="0.0011833591652154437" />
+        <nuclide name="Ni62" wo="0.0038349280648001835" />
+        <nuclide name="Ni64" wo="0.0010080415617383354" />
+    </material>
+    <material id="10004" name="Carbon Steel">
+        <temperature>300</temperature>
+        <density units="g/cc" value="7.8" />
+        <nuclide name="C0" wo="0.0027" />
+        <nuclide name="Mn55" wo="0.0075" />
+        <nuclide name="P31" wo="0.00025" />
+        <nuclide name="S32" wo="0.00023678833437985024" />
+        <nuclide name="S33" wo="1.928017942950141e-06" />
+        <nuclide name="S34" wo="1.1255605170437904e-05" />
+        <nuclide name="S36" wo="2.8042506761729612e-08" />
+        <nuclide name="Si28" wo="0.0036746594715526145" />
+        <nuclide name="Si29" wo="0.00019334515498914892" />
+        <nuclide name="Si30" wo="0.00013199537345823665" />
+        <nuclide name="Ni58" wo="0.005039835645094381" />
+        <nuclide name="Ni60" wo="0.002008189695524072" />
+        <nuclide name="Ni61" wo="8.875193739115826e-05" />
+        <nuclide name="Ni62" wo="0.0002876196048600137" />
+        <nuclide name="Ni64" wo="7.560311713037514e-05" />
+        <nuclide name="Cr50" wo="0.00014607902655887153" />
+        <nuclide name="Cr52" wo="0.0029294776849018824" />
+        <nuclide name="Cr53" wo="0.00033857554085276" />
+        <nuclide name="Cr54" wo="8.586774768648565e-05" />
+        <nuclide name="Mo100" wo="0.0006389990196986576" />
+        <nuclide name="Mo92" wo="0.0008697692153084613" />
+        <nuclide name="Mo94" wo="0.0005596299371553791" />
+        <nuclide name="Mo95" wo="0.0009791265884380354" />
+        <nuclide name="Mo96" wo="0.0010412766216951505" />
+        <nuclide name="Mo97" wo="0.0006059164130800636" />
+        <nuclide name="Mo98" wo="0.0015552822046242536" />
+        <nuclide name="V50" wo="1.2256016778573164e-06" />
+        <nuclide name="V51" wo="0.0004987743983221427" />
+        <nuclide name="Nb93" wo="0.0001" />
+        <nuclide name="Cu63" wo="0.001369583906732317" />
+        <nuclide name="Cu65" wo="0.0006304160932676829" />
+        <nuclide name="Ca40" wo="0.00014499268968855714" />
+        <nuclide name="Ca42" wo="1.0160391170196e-06" />
+        <nuclide name="Ca43" wo="2.1705537495761849e-07" />
+        <nuclide name="Ca44" wo="3.4317237524995553e-06" />
+        <nuclide name="Ca46" wo="6.8796341950755285e-09" />
+        <nuclide name="Ca48" wo="3.3561243277098485e-07" />
+        <nuclide name="B10" wo="5.529283001779395e-06" />
+        <nuclide name="B11" wo="2.4470716998220607e-05" />
+        <nuclide name="Ti46" wo="1.1880142852196743e-05" />
+        <nuclide name="Ti47" wo="1.0946673488791549e-05" />
+        <nuclide name="Ti48" wo="0.00011076757837453494" />
+        <nuclide name="Ti49" wo="8.298285799079257e-06" />
+        <nuclide name="Ti50" wo="8.107319485397494e-06" />
+        <nuclide name="Al27" wo="0.00025" />
+        <nuclide name="Fe54" wo="0.054472297655949964" />
+        <nuclide name="Fe56" wo="0.8867302806705092" />
+        <nuclide name="Fe57" wo="0.020844748673414723" />
+        <nuclide name="Fe58" wo="0.0028226730001262813" />
+    </material>
+    <material id="10005" name="Zircaloy-4">
+        <temperature>300</temperature>
+        <density units="g/cc" value="6.55" />
+        <nuclide name="O16" wo="0.001249495189868807" />
+        <nuclide name="O17" wo="5.048101311928357e-07" />
+        <nuclide name="Cr50" wo="4.1736864731106146e-05" />
+        <nuclide name="Cr52" wo="0.0008369936242576807" />
+        <nuclide name="Cr53" wo="9.673586881507429e-05" />
+        <nuclide name="Cr54" wo="2.4533642196138756e-05" />
+        <nuclide name="Fe54" wo="0.00011855672274761877" />
+        <nuclide name="Fe56" wo="0.001929932104229657" />
+        <nuclide name="Fe57" wo="4.536774095388075e-05" />
+        <nuclide name="Fe58" wo="6.143432068843669e-06" />
+        <nuclide name="Zr90" wo="0.49750307249921255" />
+        <nuclide name="Zr91" wo="0.10970127796055709" />
+        <nuclide name="Zr92" wo="0.16952409354767467" />
+        <nuclide name="Zr94" wo="0.17553856942304608" />
+        <nuclide name="Zr96" wo="0.02888298656950975" />
+        <nuclide name="Sn112" wo="0.0001325869644430062" />
+        <nuclide name="Sn114" wo="9.182449637587617e-05" />
+        <nuclide name="Sn115" wo="4.771905922545867e-05" />
+        <nuclide name="Sn116" wo="0.002058423153629443" />
+        <nuclide name="Sn117" wo="0.0010966473429083066" />
+        <nuclide name="Sn118" wo="0.0034879812938438245" />
+        <nuclide name="Sn119" wo="0.001247577110245757" />
+        <nuclide name="Sn120" wo="0.004771539495238715" />
+        <nuclide name="Sn122" wo="0.0006894094798456136" />
+        <nuclide name="Sn124" wo="0.000876291604244001" />
+    </material>
+    <material id="10006" name="aic_rod">
+        <temperature>300</temperature>
+        <density units="g/cc" value="10.16" />
+        <nuclide name="Ag107" wo="0.4110094082785408" />
+        <nuclide name="Ag109" wo="0.3889905917214592" />
+        <nuclide name="In113" wo="0.006327728141638555" />
+        <nuclide name="In115" wo="0.14367227185836143" />
+    </material>
+    <material id="10007" name="Borated Water">
+        <temperature>300</temperature>
+        <density units="g/cc" value="0.7405820675158279" />
+        <nuclide ao="0.00032310954757148997" name="B10" />
+        <nuclide ao="0.001300556520626952" name="B11" />
+        <nuclide ao="1.9965230413067991" name="H1" />
+        <nuclide ao="0.0002296265568043144" name="H2" />
+        <nuclide ao="0.9979969509249076" name="O16" />
+        <nuclide ao="0.00037938300689408465" name="O17" />
+        <sab name="c_H_in_H2O" />
+    </material>
+    <material id="10008" name="Borosilicate Glass">
+        <temperature>300</temperature>
+        <density units="g/cc" value="2.26" />
+        <nuclide ao="0.16923132768520646" name="B10" />
+        <nuclide ao="0.8307686723147936" name="B11" />
+        <nuclide ao="0.6509773086458921" name="O16" />
+        <nuclide ao="0.0002474654141428133" name="O17" />
+        <nuclide ao="0.23638904846495962" name="Si28" />
+        <nuclide ao="0.01200874718951168" name="Si29" />
+        <nuclide ao="0.007925516821765233" name="Si30" />
+        <nuclide ao="0.02423631977691551" name="Al27" />
+    </material>
+    <material id="10009" name="1.6\% Enr. UO2 Fuel">
+        <temperature>300</temperature>
+        <density units="g/cc" value="10.31341" />
+        <nuclide ao="1.9998012531906832" name="O16" />
+        <nuclide ao="0.0007602133572882293" name="O17" />
+        <nuclide ao="1.3101134677091015e-06" name="U234" />
+        <nuclide ao="0.00016306537007819218" name="U235" />
+        <nuclide ao="0.9998356245164544" name="U238" />
+    </material>
+    <material id="10010" name="2.4\% Enr. UO2 Fuel">
+        <temperature>300</temperature>
+        <density units="g/cc" value="10.29748" />
+        <nuclide ao="1.9998012531906837" name="O16" />
+        <nuclide ao="0.0007602133572882294" name="O17" />
+        <nuclide ao="1.9528324431744344e-06" name="U234" />
+        <nuclide ao="0.00024306241626824177" name="U235" />
+        <nuclide ao="0.9997549847512887" name="U238" />
+    </material>
+    <material id="10011" name="3.1\% Enr. UO2 Fuel">
+        <temperature>300</temperature>
+        <density units="g/cc" value="10.30166" />
+        <nuclide ao="1.9998012531906824" name="O16" />
+        <nuclide ao="0.0007602133572882289" name="O17" />
+        <nuclide ao="2.524278138361776e-06" name="U234" />
+        <nuclide ao="0.00031418831952932034" name="U235" />
+        <nuclide ao="0.9996832874023323" name="U238" />
+    </material>
+</materials>
diff --git a/smr/pins.py b/smr/pins.py
index 25d46f7..18b7dfa 100644
--- a/smr/pins.py
+++ b/smr/pins.py
@@ -2,7 +2,6 @@ import openmc
 
 from materials import mats
 from surfaces import surfs
-from grids import grids
 
 # FIXME: Add docstrings, etc.
 # FIXME: Assign names of cells, universes based on dict keys
@@ -35,14 +34,12 @@ def make_pin(name, surfaces, materials, grid=None):
 
     # Add spacer grid cells if specified
     if grid:
-        inner_grid_name = 'inner grid box {}'.format(grid)
-
-        # FIXME: Does this work???
-        cell.region = openmc.Intersection(cell.region, grids[inner_grid_name])
+        grid_name = 'rod grid box {}'.format(grid)
+        cell.region = openmc.Intersection(cell.region, surfs[grid_name])
 
         cell_name = name + ' (grid)'
         cell = openmc.Cell(name=cell_name)
-        cell.region = openmc.Complement(grids[inner_grid_name])
+        cell.region = openmc.Complement(surfs[grid_name])
 
         if grid == 'top/bottom':
             cell.fill = mats['In']
@@ -88,7 +85,6 @@ univs = {}
 
 cell = openmc.Cell(name='water pin 1')
 cell.fill = mats['H2O']
-cell.region = -surfs['dummy outer']
 univs['water pin'] = openmc.Universe(name='Empty water pin cell universe')
 univs['water pin'].add_cell(cell)
 
diff --git a/smr/plots.py b/smr/plots.py
index ccecd83..e895078 100644
--- a/smr/plots.py
+++ b/smr/plots.py
@@ -23,34 +23,12 @@ col_spec = {mats['H2O'].id:     [198, 226, 255],  # light blue
 # Create a collection of plots
 plots = openmc.Plots()
 
-plot = openmc.Plot(name='radial core map 1')
+plot = openmc.Plot(name='radial core map ')
 plot.basis = 'xy'
 plot.color = 'mat'
 plot.origin = [0., 0., (highest_extent-lowest_extent)/2.]
-plot.width = [19*lattice_pitch/2, 19*lattice_pitch/2.]
-plot.filename = 'radial_core_map1'
-plot.col_spec = col_spec
-plot.background = [255, 255, 255]
-plot.pixels = [1000, 1000]
-plots += [plot]
-
-plot = openmc.Plot(name='radial core map 2')
-plot.basis = 'xy'
-plot.color = 'mat'
-plot.origin = [0., 0., 40.]
-plot.width = [19*lattice_pitch/2, 19*lattice_pitch/2.]
-plot.filename = 'radial_core_map2'
-plot.col_spec = col_spec
-plot.background = [255, 255, 255]
-plot.pixels = [1000, 1000]
-plots += [plot]
-
-plot = openmc.Plot(name='radial core map 3')
-plot.basis = 'xy'
-plot.color = 'mat'
-plot.origin = [0., 0., 210.]
-plot.width = [19*lattice_pitch/2, 19*lattice_pitch/2.]
-plot.filename = 'radial_core_map3'
+plot.width = [25*lattice_pitch/2, 25*lattice_pitch/2.]
+plot.filename = 'radial_core_map'
 plot.col_spec = col_spec
 plot.background = [255, 255, 255]
 plot.pixels = [1000, 1000]
@@ -67,17 +45,6 @@ plot.background = [255, 255, 255]
 plot.pixels = [600, 600]
 plots += [plot]
 
-plot = openmc.Plot(name='mats J8 ax bot')
-plot.basis = 'xz'
-plot.color = 'mat'
-plot.origin = [0., lattice_pitch, top_lower_nozzle]
-plot.width = [lattice_pitch, 2.1*(top_lower_nozzle-bottom_support_plate)]
-plot.filename = 'J8_mats_ax_bot'
-plot.col_spec = col_spec
-plot.background = [255, 255, 255]
-plot.pixels = [400, int(400*2.1*(top_lower_nozzle-bottom_support_plate)/lattice_pitch)]
-plots += [plot]
-
 plot = openmc.Plot(name='mats J8 nozzle')
 plot.basis = 'xy'
 plot.color = 'mat'
@@ -87,15 +54,4 @@ plot.filename = 'J8_mats_nozzle'
 plot.col_spec = col_spec
 plot.background = [255, 255, 255]
 plot.pixels = [400, 400]
-plots += [plot]
-
-plot = openmc.Plot(name='mats H8 axial top')
-plot.basis = 'xy'
-plot.color = 'mat'
-plot.origin = [0., 0., top_fuel_rod]
-plot.width = [lattice_pitch, 5*(top_upper_nozzle-top_fuel_rod)]
-plot.filename = 'H8_mats_ax_top'
-plot.col_spec = col_spec
-plot.background = [255, 255, 255]
-plot.pixels = [400, int(400*2.1*(top_upper_nozzle-top_fuel_rod)/lattice_pitch)]
 plots += [plot]
\ No newline at end of file
diff --git a/smr/plots.xml b/smr/plots.xml
new file mode 100644
index 0000000..0041114
--- /dev/null
+++ b/smr/plots.xml
@@ -0,0 +1,60 @@
+<?xml version='1.0' encoding='utf-8'?>
+<plots>
+    <!--radial core map -->
+    <plot basis="xy" color="mat" filename="radial_core_map" id="10000" type="slice">
+        <origin>0.0 0.0 127.722</origin>
+        <width>268.7955 268.7955</width>
+        <pixels>1000 1000</pixels>
+        <background>255 255 255</background>
+        <col_spec id="10000" rgb="255 218 185" />
+        <col_spec id="10001" rgb="255 255 255" />
+        <col_spec id="10002" rgb="101 101 101" />
+        <col_spec id="10003" rgb="0 0 0" />
+        <col_spec id="10004" rgb="0 0 0" />
+        <col_spec id="10005" rgb="201 201 201" />
+        <col_spec id="10006" rgb="255 0 0" />
+        <col_spec id="10007" rgb="198 226 255" />
+        <col_spec id="10008" rgb="0 255 0" />
+        <col_spec id="10009" rgb="142 35 35" />
+        <col_spec id="10010" rgb="255 215 0" />
+        <col_spec id="10011" rgb="0 0 128" />
+    </plot>
+    <!--row 8 axial-->
+    <plot basis="xz" color="mat" filename="row_8_mats_axial" id="10001" type="slice">
+        <origin>0.0 0.0 127.722</origin>
+        <width>127.722 127.722</width>
+        <pixels>600 600</pixels>
+        <background>255 255 255</background>
+        <col_spec id="10000" rgb="255 218 185" />
+        <col_spec id="10001" rgb="255 255 255" />
+        <col_spec id="10002" rgb="101 101 101" />
+        <col_spec id="10003" rgb="0 0 0" />
+        <col_spec id="10004" rgb="0 0 0" />
+        <col_spec id="10005" rgb="201 201 201" />
+        <col_spec id="10006" rgb="255 0 0" />
+        <col_spec id="10007" rgb="198 226 255" />
+        <col_spec id="10008" rgb="0 255 0" />
+        <col_spec id="10009" rgb="142 35 35" />
+        <col_spec id="10010" rgb="255 215 0" />
+        <col_spec id="10011" rgb="0 0 128" />
+    </plot>
+    <!--mats J8 nozzle-->
+    <plot basis="xy" color="mat" filename="J8_mats_nozzle" id="10002" type="slice">
+        <origin>0.0 21.50364 20.0</origin>
+        <width>21.50364 21.50364</width>
+        <pixels>400 400</pixels>
+        <background>255 255 255</background>
+        <col_spec id="10000" rgb="255 218 185" />
+        <col_spec id="10001" rgb="255 255 255" />
+        <col_spec id="10002" rgb="101 101 101" />
+        <col_spec id="10003" rgb="0 0 0" />
+        <col_spec id="10004" rgb="0 0 0" />
+        <col_spec id="10005" rgb="201 201 201" />
+        <col_spec id="10006" rgb="255 0 0" />
+        <col_spec id="10007" rgb="198 226 255" />
+        <col_spec id="10008" rgb="0 255 0" />
+        <col_spec id="10009" rgb="142 35 35" />
+        <col_spec id="10010" rgb="255 215 0" />
+        <col_spec id="10011" rgb="0 0 128" />
+    </plot>
+</plots>
diff --git a/smr/settings.xml b/smr/settings.xml
new file mode 100644
index 0000000..48b2727
--- /dev/null
+++ b/smr/settings.xml
@@ -0,0 +1,19 @@
+<?xml version='1.0' encoding='utf-8'?>
+<settings>
+    <eigenvalue>
+        <particles>10000</particles>
+        <batches>10</batches>
+        <inactive>5</inactive>
+    </eigenvalue>
+    <source strength="1.0">
+        <space type="fission">
+            <parameters>-75.26274000000001 -75.26274000000001 36.007 75.26274000000001 75.26274000000001 218.887</parameters>
+        </space>
+    </source>
+    <output>
+        <tallies>false</tallies>
+    </output>
+    <ptables>true</ptables>
+    <temperature_multipole>True</temperature_multipole>
+    <temperature_tolerance>1000</temperature_tolerance>
+</settings>
diff --git a/smr/smr.py b/smr/smr.py
deleted file mode 100644
index 57da225..0000000
--- a/smr/smr.py
+++ /dev/null
@@ -1,59 +0,0 @@
-  # settings parameters
-  entrX = 15*17
-  entrY = 15*17
-  entrZ = 100
-  xbot = -8*latticePitch/2
-  ybot = -8*latticePitch/2
-  zbot = bottomFuelStack
-  xtop = 8*latticePitch/2
-  ytop = 8*latticePitch/2
-  ztop = topActiveCore
-  if core_D == '2-D':
-    entrZ = 1
-    zbot = twoDlower
-    ztop = twoDhigher 
-  sett = { 
-            'xslib':        '/home/shared/mcnpdata/binary/cross_sections.xml',
-#            'xslib':        '/home/nhorelik/xsdata/cross_sections.xml',
-            'batches':      350,
-            'inactive':     250,
-            'particles':    int(4e4),
-            'verbosity':    7,
-            'entrX':        entrX,
-            'entrY':        entrY,
-            'entrZ':        entrZ,
-            'xbot':   xbot, 'ybot':  ybot, 'zbot': zbot,
-            'xtop':   xtop, 'ytop':  ytop, 'ztop': ztop}
-
-
-
-  tallies = {}
-
-  meshDim = 15*17
-  meshLleft = -15*latticePitch/2
-  tallies['testmesh'] = {'ttype': 'mesh',
-                      'id': 1,
-                      'type': 'rectangular',
-                      'origin': '0.0 0.0',
-                      'width': '{0} {0}'.format(-meshLleft*2/meshDim),
-                      'lleft': '{0} {0}'.format(meshLleft),
-                      'dimension': '{0} {0}'.format(meshDim)}
-  tallies['test'] = { 'ttype': 'tally',
-                      'id': 1,
-                      'mesh':tallies['testmesh']['id'],
-                      'scores':'nu-fission'}
-
-  return mats,surfs,cells,latts,sett,plots,tallies
-
-
-def main():
-
-  mats,surfs,cells,latts,sett,plots,tallies,cmfd = init_data()
-
-  write_materials(mats,"materials.xml")
-  write_geometry(surfs,cells,latts,"geometry.xml")
-  write_settings(sett,"settings.xml")
-  write_plots(plots,"plots.xml")
-
-if __name__ == "__main__":
-  main()
diff --git a/smr/surfaces.py b/smr/surfaces.py
index 91370fc..64db3f7 100644
--- a/smr/surfaces.py
+++ b/smr/surfaces.py
@@ -89,91 +89,53 @@ neutron_shield_NEtop_SWbot = math.tan(-math.pi/6)
 
 surfs = {}
 
-# FIXME: Is this a good idea???
-surfs['dummy outer'] = openmc.Sphere(
-    x0=0., y0=0., R=4000., name='Dummy Outer Surface')
-
-
-surfs['pellet OR'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.39218, name='Pellet OR')
+surfs['pellet OR'] = openmc.ZCylinder(R=0.39218)
 surfs['plenum spring OR'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.06459, name='FR Plenum Spring OR')
+    R=0.06459, name='FR Plenum Spring OR')
 surfs['clad IR'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.40005, name='Clad IR')
+    R=0.40005, name='Clad IR')
 surfs['clad OR'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.45720, name='Clad OR')
+    R=0.45720, name='Clad OR')
 surfs['GT IR'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.56134, name='GT IR (above dashpot)')
+    R=0.56134, name='GT IR (above dashpot)')
 surfs['GT OR'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.60198, name='GT OR (above dashpot)')
+    R=0.60198, name='GT OR (above dashpot)')
 surfs['GT dashpot IR'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.50419, name='GT IR (at dashpot)')
+    R=0.50419, name='GT IR (at dashpot)')
 surfs['GT dashpot OR'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.54610, name='GT OR (at dashpot)')
+    R=0.54610, name='GT OR (at dashpot)')
 surfs['CP OR'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.43310, name='Control Poison OR')
+    R=0.43310, name='Control Poison OR')
 surfs['CR IR'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.43688, name='CR Clad IR')
+    R=0.43688, name='CR Clad IR')
 surfs['CR OR'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.48387, name='CR Clad OR')
+    R=0.48387, name='CR Clad OR')
 surfs['BA IR 1'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.21400, name='BA IR 1')
+    R=0.21400, name='BA IR 1')
 surfs['BA IR 2'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.23051, name='BA IR 2')
+    R=0.23051, name='BA IR 2')
 surfs['BA IR 3'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.24130, name='BA IR 3')
+    R=0.24130, name='BA IR 3')
 surfs['BA IR 4'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.42672, name='BA IR 4')
+    R=0.42672, name='BA IR 4')
 surfs['BA IR 5'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.43688, name='BA IR 5')
+    R=0.43688, name='BA IR 5')
 surfs['BA IR 6'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.48387, name='BA IR 6')
+    R=0.48387, name='BA IR 6')
 surfs['BA IR 7'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.56134, name='BA IR 7')
+    R=0.56134, name='BA IR 7')
 surfs['BA IR 8'] = openmc.ZCylinder(
-    x0=0., y0=0., R=0.60198, name='BA IR 8')
+    R=0.60198, name='BA IR 8')
 surfs['IT IR'] = copy.deepcopy(surfs['BA IR 5'])
 surfs['IT OR'] = copy.deepcopy(surfs['BA IR 6'])
 
-# Planes
-surfs['rod grid box xtop tb'] = openmc.XPlane(
-    x0=rod_grid_side_tb/2., name='X max for grid outside FR in top/bottom spacers')
-surfs['rod grid box xbot tb'] = openmc.XPlane(
-    x0=-rod_grid_side_tb/2., name='X min for grid outside FR in top/bottom spacers')
-surfs['rod grid box ytop tb'] = openmc.YPlane(
-    y0=rod_grid_side_tb/2., name='Y max for grid outside FR in top/bottom spacers')
-surfs['rod grid box ybot tb'] = openmc.YPlane(
-    y0=-rod_grid_side_tb/2., name='Y min for grid outside FR in top/bottom spacers')
+# Rectangular prisms for grid spacers
+surfs['rod grid box (top/bottom)'] = \
+    openmc.get_rectangular_prism(rod_grid_side_tb, rod_grid_side_tb)
+surfs['rod grid box (intermediate)'] = \
+    openmc.get_rectangular_prism(rod_grid_side_i, rod_grid_side_i)
 
-surfs['rod grid box xtop i'] = openmc.XPlane(
-    x0=rod_grid_side_i/2, name='X max for grid outside FR in intermediate spacers')
-surfs['rod grid box xbot i'] = openmc.XPlane(
-    x0=-rod_grid_side_i/2, name='X min for grid outside FR in intermediate spacers')
-surfs['rod grid box ytop i'] = openmc.YPlane(
-    y0=rod_grid_side_i/2, name='Y max for grid outside FR in intermediate spacers')
-surfs['rod grid box ybot i'] = openmc.YPlane(
-    y0=-rod_grid_side_i/2, name='Y min for grid outside FR in intermediate spacers')
-
-# FIXME: Can these be removed??? Use the get_rectangular_prism method instead???
-surfs['lat grid box xtop'] = openmc.XPlane(
-    x0=grid_strap_side/2, name='X max for grid outside FA')
-surfs['lat grid box xbot'] = openmc.XPlane(
-    x0=-grid_strap_side/2, name='Y min for grid outside FA')
-surfs['lat grid box ytop'] = openmc.YPlane(
-    y0=grid_strap_side/2, name='Y max for grid outside FA')
-surfs['lat grid box ybot'] = openmc.YPlane(
-    y0=-grid_strap_side/2, name='Y min for grid outside FA')
-
-surfs['lat box xtop'] = openmc.XPlane(
-    x0=17.*pin_pitch/2., name='lattice X max')
-surfs['lat box xbot'] = openmc.XPlane(
-    x0=-17.*pin_pitch/2., name='lattice X min')
-surfs['lat box ytop'] = openmc.YPlane(
-    y0=17.*pin_pitch/2., name='lattice Y max')
-surfs['lat box ybot'] = openmc.YPlane(
-    y0=-17.*pin_pitch/2., name='lattice Y min')
-
-# FIXME: rectangular prism for lattice grid sleeves
+# Rectangular prisms for lattice grid sleeves
 surfs['lat grid box inner'] = \
     openmc.get_rectangular_prism(17.*pin_pitch, 17.*pin_pitch)
 surfs['lat grid box outer'] = \
@@ -255,11 +217,11 @@ surfs['bankD bot'] = openmc.ZPlane(
 
 # outer radial surfaces
 surfs['core barrel IR'] = openmc.ZCylinder(
-    x0=0., y0=0., R=core_barrel_IR, name='core barrel IR')
+    R=core_barrel_IR, name='core barrel IR')
 surfs['core barrel OR'] = openmc.ZCylinder(
-    x0=0., y0=0., R=core_barrel_OR, name='core barrel OR')
+    R=core_barrel_OR, name='core barrel OR')
 surfs['neutron shield OR'] = openmc.ZCylinder(
-    x0=0., y0=0., R=neutron_shield_OR, name='neutron shield OR')
+    R=neutron_shield_OR, name='neutron shield OR')
 
 # neutron shield planes
 surfs['neutron shield NWbot SEtop'] = openmc.Plane(
@@ -277,9 +239,9 @@ surfs['neutron shield NEtop SWbot'] = openmc.Plane(
 
 # outer radial surfaces
 surfs['RPV IR'] = openmc.ZCylinder(
-    x0=0., y0=0., R=rpv_IR, name='RPV IR')
+    R=rpv_IR, name='RPV IR')
 surfs['RPV OR'] = openmc.ZCylinder(
-    x0=0., y0=0., R=rpv_OR, name='RPV OR', boundary_type='vacuum')
+    R=rpv_OR, name='RPV OR', boundary_type='vacuum')
 
 # outer axial surfaces
 surfs['upper bound'] = openmc.ZPlane(
diff --git a/smr/tallies.xml b/smr/tallies.xml
new file mode 100644
index 0000000..a43c2ac
--- /dev/null
+++ b/smr/tallies.xml
@@ -0,0 +1,13 @@
+<?xml version='1.0' encoding='utf-8'?>
+<tallies>
+    <!--pin-wise mesh-->
+    <mesh id="10000" type="regular">
+        <dimension>119 119</dimension>
+        <lower_left>-75.26274000000001 -75.26274000000001</lower_left>
+        <width>150.52548 150.52548</width>
+    </mesh>
+    <tally id="10000" name="fission rates">
+        <filter bins="10000" type="mesh" />
+        <scores>fission</scores>
+    </tally>
+</tallies>