From e9ec9c0f9b1777f39b18d9e5e91652698939b9da Mon Sep 17 00:00:00 2001
From: Adam Nelson <nelsonag@umich.edu>
Date: Sun, 18 May 2014 06:49:15 -0400
Subject: [PATCH] Adding tests to harmonics branch; need to update results
 after obtaining NNDC data

---
 tests/test_score_flux_yn/geometry.xml         | 185 ++++++++++++
 tests/test_score_flux_yn/materials.xml        | 272 ++++++++++++++++++
 tests/test_score_flux_yn/results.py           |  44 +++
 tests/test_score_flux_yn/settings.xml         |  19 ++
 tests/test_score_flux_yn/tallies.xml          |  14 +
 .../test_score_flux_yn/test_score_flux_yn.py  |  62 ++++
 tests/test_score_nuscatter_n/geometry.xml     | 185 ++++++++++++
 tests/test_score_nuscatter_n/materials.xml    | 272 ++++++++++++++++++
 tests/test_score_nuscatter_n/results.py       |  37 +++
 tests/test_score_nuscatter_n/settings.xml     |  19 ++
 tests/test_score_nuscatter_n/tallies.xml      |   9 +
 .../test_score_nuscatter_n.py                 |  62 ++++
 tests/test_score_nuscatter_pn/geometry.xml    | 185 ++++++++++++
 tests/test_score_nuscatter_pn/materials.xml   | 272 ++++++++++++++++++
 tests/test_score_nuscatter_pn/results.py      |  44 +++
 tests/test_score_nuscatter_pn/settings.xml    |  19 ++
 tests/test_score_nuscatter_pn/tallies.xml     |  14 +
 .../test_score_nuscatter_pn.py                |  62 ++++
 tests/test_score_nuscatter_yn/geometry.xml    | 185 ++++++++++++
 tests/test_score_nuscatter_yn/materials.xml   | 272 ++++++++++++++++++
 tests/test_score_nuscatter_yn/results.py      |  44 +++
 tests/test_score_nuscatter_yn/settings.xml    |  19 ++
 tests/test_score_nuscatter_yn/tallies.xml     |  14 +
 .../test_score_nuscatter_yn.py                |  62 ++++
 tests/test_score_scatter_yn/geometry.xml      | 185 ++++++++++++
 tests/test_score_scatter_yn/materials.xml     | 272 ++++++++++++++++++
 tests/test_score_scatter_yn/results.py        |  44 +++
 tests/test_score_scatter_yn/settings.xml      |  19 ++
 tests/test_score_scatter_yn/tallies.xml       |  14 +
 .../test_score_scatter_yn.py                  |  62 ++++
 tests/test_score_total_yn/geometry.xml        | 185 ++++++++++++
 tests/test_score_total_yn/materials.xml       | 272 ++++++++++++++++++
 tests/test_score_total_yn/results.py          |  44 +++
 tests/test_score_total_yn/settings.xml        |  19 ++
 tests/test_score_total_yn/tallies.xml         |  14 +
 .../test_score_total_yn.py                    |  62 ++++
 36 files changed, 3564 insertions(+)
 create mode 100644 tests/test_score_flux_yn/geometry.xml
 create mode 100644 tests/test_score_flux_yn/materials.xml
 create mode 100644 tests/test_score_flux_yn/results.py
 create mode 100644 tests/test_score_flux_yn/settings.xml
 create mode 100644 tests/test_score_flux_yn/tallies.xml
 create mode 100644 tests/test_score_flux_yn/test_score_flux_yn.py
 create mode 100644 tests/test_score_nuscatter_n/geometry.xml
 create mode 100644 tests/test_score_nuscatter_n/materials.xml
 create mode 100644 tests/test_score_nuscatter_n/results.py
 create mode 100644 tests/test_score_nuscatter_n/settings.xml
 create mode 100644 tests/test_score_nuscatter_n/tallies.xml
 create mode 100644 tests/test_score_nuscatter_n/test_score_nuscatter_n.py
 create mode 100644 tests/test_score_nuscatter_pn/geometry.xml
 create mode 100644 tests/test_score_nuscatter_pn/materials.xml
 create mode 100644 tests/test_score_nuscatter_pn/results.py
 create mode 100644 tests/test_score_nuscatter_pn/settings.xml
 create mode 100644 tests/test_score_nuscatter_pn/tallies.xml
 create mode 100644 tests/test_score_nuscatter_pn/test_score_nuscatter_pn.py
 create mode 100644 tests/test_score_nuscatter_yn/geometry.xml
 create mode 100644 tests/test_score_nuscatter_yn/materials.xml
 create mode 100644 tests/test_score_nuscatter_yn/results.py
 create mode 100644 tests/test_score_nuscatter_yn/settings.xml
 create mode 100644 tests/test_score_nuscatter_yn/tallies.xml
 create mode 100644 tests/test_score_nuscatter_yn/test_score_nuscatter_yn.py
 create mode 100644 tests/test_score_scatter_yn/geometry.xml
 create mode 100644 tests/test_score_scatter_yn/materials.xml
 create mode 100644 tests/test_score_scatter_yn/results.py
 create mode 100644 tests/test_score_scatter_yn/settings.xml
 create mode 100644 tests/test_score_scatter_yn/tallies.xml
 create mode 100644 tests/test_score_scatter_yn/test_score_scatter_yn.py
 create mode 100644 tests/test_score_total_yn/geometry.xml
 create mode 100644 tests/test_score_total_yn/materials.xml
 create mode 100644 tests/test_score_total_yn/results.py
 create mode 100644 tests/test_score_total_yn/settings.xml
 create mode 100644 tests/test_score_total_yn/tallies.xml
 create mode 100644 tests/test_score_total_yn/test_score_total_yn.py

diff --git a/tests/test_score_flux_yn/geometry.xml b/tests/test_score_flux_yn/geometry.xml
new file mode 100644
index 0000000000..2407bf74b3
--- /dev/null
+++ b/tests/test_score_flux_yn/geometry.xml
@@ -0,0 +1,185 @@
+<?xml version="1.0"?>
+<geometry>
+
+  <surface id="1" type="z-cylinder" coeffs="0. 0. 0.41" />
+  <surface id="2" type="z-cylinder" coeffs="0. 0. 0.475" />
+  <surface id="3" type="z-cylinder" coeffs="0. 0. 0.56" />
+  <surface id="4" type="z-cylinder" coeffs="0. 0. 0.62" />
+  <surface id="5" type="z-cylinder" coeffs="0. 0. 187.6" />
+  <surface id="6" type="z-cylinder" coeffs="0. 0. 209.0" />
+  <surface id="7" type="z-cylinder" coeffs="0. 0. 229.0" />
+  <surface id="8" type="z-cylinder" coeffs="0. 0. 249.0" boundary="vacuum" />
+
+  <surface id="31" type="z-plane" coeffs="-229.0" boundary="vacuum" />
+  <surface id="32" type="z-plane" coeffs="-199.0" />
+  <surface id="33" type="z-plane" coeffs="-193.0" />
+  <surface id="34" type="z-plane" coeffs="-183.0" />
+  <surface id="35" type="z-plane" coeffs="0.0" />
+  <surface id="36" type="z-plane" coeffs="183.0" />
+  <surface id="37" type="z-plane" coeffs="203.0" />
+  <surface id="38" type="z-plane" coeffs="215.0" />
+  <surface id="39" type="z-plane" coeffs="223.0" boundary="vacuum" />
+
+  <!-- All geometry on base universe -->
+  <cell id="1"  fill="200"    surfaces="  -6 34 -35" /> <!-- Lower core -->
+  <cell id="2"  fill="201"    surfaces="  -6 35 -36" /> <!-- Upper core -->
+  <cell id="3"  material="8"  surfaces="  -7 31 -32" /> <!-- Lower core plate region -->
+  <cell id="4"  material="9"  surfaces="  -5 32 -33" /> <!-- Bottom nozzle region -->
+  <cell id="5"  material="12" surfaces="  -5 33 -34" /> <!-- Bottom FA region -->
+  <cell id="6"  material="11" surfaces="  -5 36 -37" /> <!-- Top FA region -->
+  <cell id="7"  material="10" surfaces="  -5 37 -38" /> <!-- Top nozzle region -->
+  <cell id="8"  material="7"  surfaces="  -7 38 -39" /> <!-- Upper plate region -->
+  <cell id="9"  material="4"  surfaces="6 -7 32 -38" /> <!-- Downcomer -->
+  <cell id="10" material="5"  surfaces="7 -8 31 -39" /> <!-- RPV -->
+  <cell id="11" material="6"  surfaces="5 -6 32 -34" /> <!-- Bottom of radial reflector -->
+  <cell id="12" material="7"  surfaces="5 -6 36 -38" /> <!-- Top of radial reflector -->
+
+  <!-- Fuel pin, cladding, cold water -->
+  <cell id="21" universe="1" material="1" surfaces="-1"   />
+  <cell id="22" universe="1" material="2" surfaces="1 -2" />
+  <cell id="23" universe="1" material="3" surfaces="2"    />
+
+  <!-- Instrumentation guide tube -->
+  <cell id="24" universe="2" material="3" surfaces="-3"   />
+  <cell id="25" universe="2" material="2" surfaces="3 -4" />
+  <cell id="26" universe="2" material="3" surfaces="4"    />
+
+  <!-- Fuel pin, cladding, hot water -->
+  <cell id="27" universe="3" material="1" surfaces="-1"   />
+  <cell id="28" universe="3" material="2" surfaces="1 -2" />
+  <cell id="29" universe="3" material="4" surfaces="2"    />
+
+  <!-- Instrumentation guide tube -->
+  <cell id="30" universe="4" material="4" surfaces="-3"   />
+  <cell id="31" universe="4" material="2" surfaces="3 -4" />
+  <cell id="32" universe="4" material="4" surfaces="4"    />
+
+  <!-- cell for water assembly (cold) -->
+  <cell id="50" universe="5" material="4" surfaces="34 -35" />
+
+  <!-- containing cell for fuel assembly -->
+  <cell id="60" universe="6" fill="100"   surfaces="34 -35" />
+
+  <!-- cell for water assembly (hot) -->
+  <cell id="70" universe="7" material="3" surfaces="35 -36" />
+
+  <!-- containing cell for fuel assembly -->
+  <cell id="80" universe="8" fill="101"   surfaces="35 -36" />
+
+  <!-- Fuel Assembly (Lower Half) -->
+  <lattice id="100">
+    <type>rectangular</type>
+    <dimension>17 17</dimension>
+    <lower_left>-10.71 -10.71</lower_left>
+    <width>1.26 1.26</width>
+    <universes>
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1
+      1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
+      1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+    </universes>
+  </lattice>
+
+  <!-- Fuel Assembly (Upper Half) -->
+  <lattice id="101">
+    <type>rectangular</type>
+    <dimension>17 17</dimension>
+    <lower_left>-10.71 -10.71</lower_left>
+    <width>1.26 1.26</width>
+    <universes>
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 4 3 3 4 3 3 4 3 3 3 3 3
+      3 3 3 4 3 3 3 3 3 3 3 3 3 4 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 4 3 3 3 3 3 3 3 3 3 4 3 3 3
+      3 3 3 3 3 4 3 3 4 3 3 4 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+    </universes>
+  </lattice>
+
+  <!-- Core Lattice (Lower Half) -->
+  <lattice id="200">
+    <type>rectangular</type>
+    <dimension>21 21</dimension>
+    <lower_left>-224.91 -224.91</lower_left>
+    <width>21.42 21.42</width>
+    <universes>
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 6 6 6 6 6 6 6 5 5 5 5 5 5 5
+      5 5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5
+      5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5
+      5 5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5
+      5 5 5 5 5 5 5 6 6 6 6 6 6 6 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+    </universes>
+  </lattice>
+
+  <!-- Core Lattice (Upper Half) -->
+  <lattice id="201">
+    <type>rectangular</type>
+    <dimension>21 21</dimension>
+    <lower_left>-224.91 -224.91</lower_left>
+    <width>21.42 21.42</width>
+    <universes>
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 8 8 8 8 8 8 8 7 7 7 7 7 7 7
+      7 7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7 7
+      7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7
+      7 7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7 7
+      7 7 7 7 7 7 7 8 8 8 8 8 8 8 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+    </universes>
+  </lattice>
+
+</geometry>
\ No newline at end of file
diff --git a/tests/test_score_flux_yn/materials.xml b/tests/test_score_flux_yn/materials.xml
new file mode 100644
index 0000000000..9c0b74f3f1
--- /dev/null
+++ b/tests/test_score_flux_yn/materials.xml
@@ -0,0 +1,272 @@
+<?xml version="1.0"?>
+<materials>
+
+  <default_xs>71c</default_xs>
+
+  <!-- Fuel composition -->
+  <material id="1">
+    <density value="10.062" units="g/cm3" />
+    <nuclide name="U-234"   ao="4.9476e-6" />
+    <nuclide name="U-235"   ao="4.8218e-4" />
+    <nuclide name="U-236"   ao="9.0402e-5" />
+    <nuclide name="U-238"   ao="2.1504e-2" />
+    <nuclide name="Np-237"  ao="7.3733e-6" />
+    <nuclide name="Pu-238"  ao="1.5148e-6" />
+    <nuclide name="Pu-239"  ao="1.3955e-4" />
+    <nuclide name="Pu-240"  ao="3.4405e-5" />
+    <nuclide name="Pu-241"  ao="2.1439e-5" />
+    <nuclide name="Pu-242"  ao="3.7422e-6" />
+    <nuclide name="Am-241"  ao="4.5041e-7" />
+    <nuclide name="Am-242m" ao="9.2301e-9" />
+    <nuclide name="Am-243"  ao="4.7878e-7" />
+    <nuclide name="Cm-242"  ao="1.0485e-7" />
+    <nuclide name="Cm-243"  ao="1.4268e-9" />
+    <nuclide name="Cm-244"  ao="8.8756e-8" />
+    <nuclide name="Cm-245"  ao="3.5285e-9" />
+    <nuclide name="Mo-95"   ao="2.6497e-5" />
+    <nuclide name="Tc-99"   ao="3.2772e-5" />
+    <nuclide name="Ru-101"  ao="3.0742e-5" />
+    <nuclide name="Ru-103"  ao="2.3505e-6" />
+    <nuclide name="Ag-109"  ao="2.0009e-6" />
+    <nuclide name="Xe-135"  ao="1.0801e-8" />
+    <nuclide name="Cs-133"  ao="3.4612e-5" />
+    <nuclide name="Nd-143"  ao="2.6078e-5" />
+    <nuclide name="Nd-145"  ao="1.9898e-5" />
+    <nuclide name="Sm-147"  ao="1.6128e-6" />
+    <nuclide name="Sm-149"  ao="1.1627e-7" />
+    <nuclide name="Sm-150"  ao="7.1727e-6" />
+    <nuclide name="Sm-151"  ao="5.4947e-7" />
+    <nuclide name="Sm-152"  ao="3.0221e-6" />
+    <nuclide name="Eu-153"  ao="2.6209e-6" />
+    <nuclide name="Gd-155"  ao="1.5369e-9" />
+    <nuclide name="O-16"    ao="4.5737e-2" />
+  </material>
+
+  <!-- Cladding composition -->
+  <material id="2">
+    <density value="5.77" units="g/cm3" />
+    <nuclide name="Zr-90" ao="0.5145" />
+    <nuclide name="Zr-91" ao="0.1122" />
+    <nuclide name="Zr-92" ao="0.1715" />
+    <nuclide name="Zr-94" ao="0.1738" />
+    <nuclide name="Zr-96" ao="0.0280" />
+  </material>
+
+  <!-- Cold borated water -->
+  <material id="3">
+    <density value="0.07416" units="atom/b-cm" />
+    <nuclide name="H-1" ao="2.0" />
+    <nuclide name="O-16" ao="1.0" />
+    <nuclide name="B-10" ao="6.490e-4" />
+    <nuclide name="B-11" ao="2.689e-3" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Hot borated water -->
+  <material id="4">
+    <density value="0.06614" units="atom/b-cm" />
+    <nuclide name="H-1"  ao="2.0" />
+    <nuclide name="O-16" ao="1.0" />
+    <nuclide name="B-10" ao="6.490e-4" />
+    <nuclide name="B-11" ao="2.689e-3" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- RPV Composition -->
+  <material id="5">
+    <density value="7.9" units="g/cm3" />
+    <nuclide name="Fe-54"  wo="0.05437098" />
+    <nuclide name="Fe-56"  wo="0.88500663" />
+    <nuclide name="Fe-57"  wo="0.0208008" />
+    <nuclide name="Fe-58"  wo="0.00282159" />
+    <nuclide name="Ni-58"  wo="0.0067198" />
+    <nuclide name="Ni-60"  wo="0.0026776" />
+    <nuclide name="Ni-61"  wo="0.0001183" />
+    <nuclide name="Ni-62"  wo="0.0003835" />
+    <nuclide name="Ni-64"  wo="0.0001008" />
+    <nuclide name="Mn-55"  wo="0.01" />
+    <nuclide name="Mo-92"  wo="0.000849" />
+    <nuclide name="Mo-94"  wo="0.0005418" />
+    <nuclide name="Mo-95"  wo="0.0009438" />
+    <nuclide name="Mo-96"  wo="0.0010002" />
+    <nuclide name="Mo-97"  wo="0.0005796" />
+    <nuclide name="Mo-98"  wo="0.0014814" />
+    <nuclide name="Mo-100" wo="0.0006042" />
+    <nuclide name="Si-28"  wo="0.00367464" />
+    <nuclide name="Si-29"  wo="0.00019336" />
+    <nuclide name="Si-30"  wo="0.000132" />
+    <nuclide name="Cr-50"  wo="0.00010435" />
+    <nuclide name="Cr-52"  wo="0.002092475" />
+    <nuclide name="Cr-53"  wo="0.00024185" />
+    <nuclide name="Cr-54"  wo="6.1325e-05" />
+    <nuclide name="C-Nat"  wo="0.0025" />
+    <nuclide name="Cu-63"  wo="0.0013696" />
+    <nuclide name="Cu-65"  wo="0.0006304" />
+  </material>
+
+  <!-- Lower radial reflector -->
+  <material id="6">
+    <density value="4.32" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0095661" />
+    <nuclide name="O-16"  wo="0.0759107" />
+    <nuclide name="B-10"  wo="3.08409e-5" />
+    <nuclide name="B-11"  wo="1.40499e-4" />
+    <nuclide name="Fe-54" wo="0.035620772088" />
+    <nuclide name="Fe-56" wo="0.579805982228" />
+    <nuclide name="Fe-57" wo="0.01362750048" />
+    <nuclide name="Fe-58" wo="0.001848545204" />
+    <nuclide name="Ni-58" wo="0.055298376566" />
+    <nuclide name="Ni-60" wo="0.022034425592" />
+    <nuclide name="Ni-61" wo="0.000973510811" />
+    <nuclide name="Ni-62" wo="0.003155886695" />
+    <nuclide name="Ni-64" wo="0.000829500336" />
+    <nuclide name="Mn-55" wo="0.0182870" />
+    <nuclide name="Si-28" wo="0.00839976771" />
+    <nuclide name="Si-29" wo="0.00044199679" />
+    <nuclide name="Si-30" wo="0.0003017355" />
+    <nuclide name="Cr-50" wo="0.007251360806" />
+    <nuclide name="Cr-52" wo="0.145407678031" />
+    <nuclide name="Cr-53" wo="0.016806340306" />
+    <nuclide name="Cr-54" wo="0.004261520857" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Upper radial reflector / Top plate region -->
+  <material id="7">
+    <density value="4.28" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0086117" />
+    <nuclide name="O-16"  wo="0.0683369" />
+    <nuclide name="B-10"  wo="2.77638e-5" />
+    <nuclide name="B-11"  wo="1.26481e-4" />
+    <nuclide name="Fe-54" wo="0.035953677186" />
+    <nuclide name="Fe-56" wo="0.585224740891" />
+    <nuclide name="Fe-57" wo="0.01375486056" />
+    <nuclide name="Fe-58" wo="0.001865821363" />
+    <nuclide name="Ni-58" wo="0.055815129186" />
+    <nuclide name="Ni-60" wo="0.022240333032" />
+    <nuclide name="Ni-61" wo="0.000982608081" />
+    <nuclide name="Ni-62" wo="0.003185377845" />
+    <nuclide name="Ni-64" wo="0.000837251856" />
+    <nuclide name="Mn-55" wo="0.0184579" />
+    <nuclide name="Si-28" wo="0.00847831314" />
+    <nuclide name="Si-29" wo="0.00044612986" />
+    <nuclide name="Si-30" wo="0.000304557" />
+    <nuclide name="Cr-50" wo="0.00731912987" />
+    <nuclide name="Cr-52" wo="0.146766614995" />
+    <nuclide name="Cr-53" wo="0.01696340737" />
+    <nuclide name="Cr-54" wo="0.004301347765" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom plate region -->
+  <material id="8">
+    <density value="7.184" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0011505" />
+    <nuclide name="O-16"  wo="0.0091296" />
+    <nuclide name="B-10"  wo="3.70915e-6" />
+    <nuclide name="B-11"  wo="1.68974e-5" />
+    <nuclide name="Fe-54" wo="0.03855611055" />
+    <nuclide name="Fe-56" wo="0.627585036425" />
+    <nuclide name="Fe-57" wo="0.014750478" />
+    <nuclide name="Fe-58" wo="0.002000875025" />
+    <nuclide name="Ni-58" wo="0.059855207342" />
+    <nuclide name="Ni-60" wo="0.023850159704" />
+    <nuclide name="Ni-61" wo="0.001053732407" />
+    <nuclide name="Ni-62" wo="0.003415945715" />
+    <nuclide name="Ni-64" wo="0.000897854832" />
+    <nuclide name="Mn-55" wo="0.0197940" />
+    <nuclide name="Si-28" wo="0.00909197802" />
+    <nuclide name="Si-29" wo="0.00047842098" />
+    <nuclide name="Si-30" wo="0.000326601" />
+    <nuclide name="Cr-50" wo="0.007848910646" />
+    <nuclide name="Cr-52" wo="0.157390026871" />
+    <nuclide name="Cr-53" wo="0.018191270146" />
+    <nuclide name="Cr-54" wo="0.004612692337" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom nozzle region -->
+  <material id="9">
+    <density value="2.53" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0245014" />
+    <nuclide name="O-16"  wo="0.1944274" />
+    <nuclide name="B-10"  wo="7.89917e-5" />
+    <nuclide name="B-11"  wo="3.59854e-4" />
+    <nuclide name="Fe-54" wo="0.030411411144" />
+    <nuclide name="Fe-56" wo="0.495012237964" />
+    <nuclide name="Fe-57" wo="0.01163454624" />
+    <nuclide name="Fe-58" wo="0.001578204652" />
+    <nuclide name="Ni-58" wo="0.047211231662" />
+    <nuclide name="Ni-60" wo="0.018811987544" />
+    <nuclide name="Ni-61" wo="0.000831139127" />
+    <nuclide name="Ni-62" wo="0.002694352115" />
+    <nuclide name="Ni-64" wo="0.000708189552" />
+    <nuclide name="Mn-55" wo="0.0156126" />
+    <nuclide name="Si-28" wo="0.007171335558" />
+    <nuclide name="Si-29" wo="0.000377356542" />
+    <nuclide name="Si-30" wo="0.0002576079" />
+    <nuclide name="Cr-50" wo="0.006190885148" />
+    <nuclide name="Cr-52" wo="0.124142524198" />
+    <nuclide name="Cr-53" wo="0.014348496148" />
+    <nuclide name="Cr-54" wo="0.003638294506" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Top nozzle region -->
+  <material id="10">
+    <density value="1.746" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0358870" />
+    <nuclide name="O-16"  wo="0.2847761" />
+    <nuclide name="B-10"  wo="1.15699e-4" />
+    <nuclide name="B-11"  wo="5.27075e-4" />
+    <nuclide name="Fe-54" wo="0.02644016154" />
+    <nuclide name="Fe-56" wo="0.43037146399" />
+    <nuclide name="Fe-57" wo="0.0101152584" />
+    <nuclide name="Fe-58" wo="0.00137211607" />
+    <nuclide name="Ni-58" wo="0.04104621835" />
+    <nuclide name="Ni-60" wo="0.0163554502" />
+    <nuclide name="Ni-61" wo="0.000722605975" />
+    <nuclide name="Ni-62" wo="0.002342513875" />
+    <nuclide name="Ni-64" wo="0.0006157116" />
+    <nuclide name="Mn-55" wo="0.0135739" />
+    <nuclide name="Si-28" wo="0.006234853554" />
+    <nuclide name="Si-29" wo="0.000328078746" />
+    <nuclide name="Si-30" wo="0.0002239677" />
+    <nuclide name="Cr-50" wo="0.005382452306" />
+    <nuclide name="Cr-52" wo="0.107931450781" />
+    <nuclide name="Cr-53" wo="0.012474806806" />
+    <nuclide name="Cr-54" wo="0.003163190107" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Top of Fuel Assemblies -->
+  <material id="11">
+    <density value="3.044" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0162913" />
+    <nuclide name="O-16"  wo="0.1292776" />
+    <nuclide name="B-10"  wo="5.25228e-5" />
+    <nuclide name="B-11"  wo="2.39272e-4" />
+    <nuclide name="Zr-90" wo="0.43313403903" />
+    <nuclide name="Zr-91" wo="0.09549277374" />
+    <nuclide name="Zr-92" wo="0.14759527104" />
+    <nuclide name="Zr-94" wo="0.15280552077" />
+    <nuclide name="Zr-96" wo="0.02511169542" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom of Fuel Assemblies -->
+  <material id="12">
+    <density value="1.762" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0292856" />
+    <nuclide name="O-16"  wo="0.2323919" />
+    <nuclide name="B-10"  wo="9.44159e-5" />
+    <nuclide name="B-11"  wo="4.30120e-4" />
+    <nuclide name="Zr-90" wo="0.3741373658" />
+    <nuclide name="Zr-91" wo="0.0824858164" />
+    <nuclide name="Zr-92" wo="0.1274914944" />
+    <nuclide name="Zr-94" wo="0.1319920622" />
+    <nuclide name="Zr-96" wo="0.0216912612" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+</materials>
diff --git a/tests/test_score_flux_yn/results.py b/tests/test_score_flux_yn/results.py
new file mode 100644
index 0000000000..f4c32a893f
--- /dev/null
+++ b/tests/test_score_flux_yn/results.py
@@ -0,0 +1,44 @@
+#!/usr/bin/env python
+
+import sys
+import numpy as np
+
+# import statepoint
+sys.path.append('../../src/utils')
+import statepoint
+
+# read in statepoint file
+if len(sys.argv) > 1:
+    sp = statepoint.StatePoint(sys.argv[1])
+else:
+    sp = statepoint.StatePoint('statepoint.10.binary')
+sp.read_results()
+
+# extract tally results and convert to vector
+results1 = sp.tallies[0].results
+shape1 = results1.shape
+size1 = (np.product(shape1))
+results1 = np.reshape(results1, size1)
+results2 = sp.tallies[1].results
+shape2 = results2.shape
+size2 = (np.product(shape2))
+results2 = np.reshape(results2, size2)
+
+# set up output string
+outstr = ''
+ 
+# write out k-combined
+outstr += 'k-combined:\n'
+outstr += "{0:12.6E} {1:12.6E}\n".format(sp.k_combined[0], sp.k_combined[1])
+
+# write out tally results
+outstr += 'tally 1:\n'
+for item in results1:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'tally 2:\n'
+for item in results2:
+  outstr += "{0:12.6E}\n".format(item)
+
+# write results to file
+with open('results_test.dat','w') as fh:
+    fh.write(outstr)
diff --git a/tests/test_score_flux_yn/settings.xml b/tests/test_score_flux_yn/settings.xml
new file mode 100644
index 0000000000..517637a59f
--- /dev/null
+++ b/tests/test_score_flux_yn/settings.xml
@@ -0,0 +1,19 @@
+<?xml version="1.0"?>
+<settings>
+
+  <eigenvalue>
+    <batches>10</batches>
+    <inactive>5</inactive>
+    <particles>100</particles>
+  </eigenvalue>
+
+  <source>
+    <space type="box">
+      <parameters>
+        -160 -160 -183
+         160  160  183
+      </parameters>
+    </space>
+  </source>
+
+</settings>
diff --git a/tests/test_score_flux_yn/tallies.xml b/tests/test_score_flux_yn/tallies.xml
new file mode 100644
index 0000000000..71301d0e7f
--- /dev/null
+++ b/tests/test_score_flux_yn/tallies.xml
@@ -0,0 +1,14 @@
+<?xml version="1.0"?>
+<tallies>
+
+  <tally id="1">
+    <filter type="cell" bins="21 22 23 27 28 29" />
+    <scores>flux</scores>
+  </tally>
+  
+  <tally id="2">
+    <filter type="cell" bins="21 22 23 27 28 29" />
+    <scores>flux-y5</scores>
+  </tally>
+
+</tallies>
diff --git a/tests/test_score_flux_yn/test_score_flux_yn.py b/tests/test_score_flux_yn/test_score_flux_yn.py
new file mode 100644
index 0000000000..6e54389fc1
--- /dev/null
+++ b/tests/test_score_flux_yn/test_score_flux_yn.py
@@ -0,0 +1,62 @@
+#!/usr/bin/env python
+
+import os
+from subprocess import Popen, STDOUT, PIPE, call
+import filecmp
+import glob
+from optparse import OptionParser
+
+parser = OptionParser()
+parser.add_option('--mpi_exec', dest='mpi_exec', default='')
+parser.add_option('--mpi_np', dest='mpi_np', default='3')
+parser.add_option('--exe', dest='exe')
+(opts, args) = parser.parse_args()
+cwd = os.getcwd()
+
+def test_run():
+    if opts.mpi_exec != '':
+        proc = Popen([opts.mpi_exec, '-np', opts.mpi_np, opts.exe, cwd],
+               stderr=STDOUT, stdout=PIPE)
+    else:
+        proc = Popen([opts.exe, cwd], stderr=STDOUT, stdout=PIPE)
+    print(proc.communicate()[0])
+    returncode = proc.returncode
+    assert returncode == 0, 'OpenMC did not exit successfully.'
+
+def test_created_statepoint():
+    statepoint = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    assert len(statepoint) == 1, 'Either multiple or no statepoint files exist.'
+    assert statepoint[0].endswith('binary') or statepoint[0].endswith('h5'),\
+        'Statepoint file is not a binary or hdf5 file.'
+
+def test_output_exists():
+    assert os.path.exists(os.path.join(cwd, 'tallies.out')), 'Tally output file does not exist.'
+
+def test_results():
+    statepoint = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    call(['python', 'results.py', statepoint[0]])
+    compare = filecmp.cmp('results_test.dat', 'results_true.dat')
+    if not compare:
+      os.rename('results_test.dat', 'results_error.dat')
+    assert compare, 'Results do not agree.'
+
+def teardown():
+    output = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    output.append(os.path.join(cwd, 'tallies.out'))
+    output.append(os.path.join(cwd, 'results_test.dat'))
+    for f in output:
+        if os.path.exists(f):
+            os.remove(f)
+
+if __name__ == '__main__':
+
+    # test for openmc executable
+    if opts.exe is None:
+        raise Exception('Must specify OpenMC executable from command line with --exe.')
+
+    # run tests
+    test_run()
+    test_created_statepoint()
+    test_output_exists()
+    test_results()
+    teardown()
diff --git a/tests/test_score_nuscatter_n/geometry.xml b/tests/test_score_nuscatter_n/geometry.xml
new file mode 100644
index 0000000000..2407bf74b3
--- /dev/null
+++ b/tests/test_score_nuscatter_n/geometry.xml
@@ -0,0 +1,185 @@
+<?xml version="1.0"?>
+<geometry>
+
+  <surface id="1" type="z-cylinder" coeffs="0. 0. 0.41" />
+  <surface id="2" type="z-cylinder" coeffs="0. 0. 0.475" />
+  <surface id="3" type="z-cylinder" coeffs="0. 0. 0.56" />
+  <surface id="4" type="z-cylinder" coeffs="0. 0. 0.62" />
+  <surface id="5" type="z-cylinder" coeffs="0. 0. 187.6" />
+  <surface id="6" type="z-cylinder" coeffs="0. 0. 209.0" />
+  <surface id="7" type="z-cylinder" coeffs="0. 0. 229.0" />
+  <surface id="8" type="z-cylinder" coeffs="0. 0. 249.0" boundary="vacuum" />
+
+  <surface id="31" type="z-plane" coeffs="-229.0" boundary="vacuum" />
+  <surface id="32" type="z-plane" coeffs="-199.0" />
+  <surface id="33" type="z-plane" coeffs="-193.0" />
+  <surface id="34" type="z-plane" coeffs="-183.0" />
+  <surface id="35" type="z-plane" coeffs="0.0" />
+  <surface id="36" type="z-plane" coeffs="183.0" />
+  <surface id="37" type="z-plane" coeffs="203.0" />
+  <surface id="38" type="z-plane" coeffs="215.0" />
+  <surface id="39" type="z-plane" coeffs="223.0" boundary="vacuum" />
+
+  <!-- All geometry on base universe -->
+  <cell id="1"  fill="200"    surfaces="  -6 34 -35" /> <!-- Lower core -->
+  <cell id="2"  fill="201"    surfaces="  -6 35 -36" /> <!-- Upper core -->
+  <cell id="3"  material="8"  surfaces="  -7 31 -32" /> <!-- Lower core plate region -->
+  <cell id="4"  material="9"  surfaces="  -5 32 -33" /> <!-- Bottom nozzle region -->
+  <cell id="5"  material="12" surfaces="  -5 33 -34" /> <!-- Bottom FA region -->
+  <cell id="6"  material="11" surfaces="  -5 36 -37" /> <!-- Top FA region -->
+  <cell id="7"  material="10" surfaces="  -5 37 -38" /> <!-- Top nozzle region -->
+  <cell id="8"  material="7"  surfaces="  -7 38 -39" /> <!-- Upper plate region -->
+  <cell id="9"  material="4"  surfaces="6 -7 32 -38" /> <!-- Downcomer -->
+  <cell id="10" material="5"  surfaces="7 -8 31 -39" /> <!-- RPV -->
+  <cell id="11" material="6"  surfaces="5 -6 32 -34" /> <!-- Bottom of radial reflector -->
+  <cell id="12" material="7"  surfaces="5 -6 36 -38" /> <!-- Top of radial reflector -->
+
+  <!-- Fuel pin, cladding, cold water -->
+  <cell id="21" universe="1" material="1" surfaces="-1"   />
+  <cell id="22" universe="1" material="2" surfaces="1 -2" />
+  <cell id="23" universe="1" material="3" surfaces="2"    />
+
+  <!-- Instrumentation guide tube -->
+  <cell id="24" universe="2" material="3" surfaces="-3"   />
+  <cell id="25" universe="2" material="2" surfaces="3 -4" />
+  <cell id="26" universe="2" material="3" surfaces="4"    />
+
+  <!-- Fuel pin, cladding, hot water -->
+  <cell id="27" universe="3" material="1" surfaces="-1"   />
+  <cell id="28" universe="3" material="2" surfaces="1 -2" />
+  <cell id="29" universe="3" material="4" surfaces="2"    />
+
+  <!-- Instrumentation guide tube -->
+  <cell id="30" universe="4" material="4" surfaces="-3"   />
+  <cell id="31" universe="4" material="2" surfaces="3 -4" />
+  <cell id="32" universe="4" material="4" surfaces="4"    />
+
+  <!-- cell for water assembly (cold) -->
+  <cell id="50" universe="5" material="4" surfaces="34 -35" />
+
+  <!-- containing cell for fuel assembly -->
+  <cell id="60" universe="6" fill="100"   surfaces="34 -35" />
+
+  <!-- cell for water assembly (hot) -->
+  <cell id="70" universe="7" material="3" surfaces="35 -36" />
+
+  <!-- containing cell for fuel assembly -->
+  <cell id="80" universe="8" fill="101"   surfaces="35 -36" />
+
+  <!-- Fuel Assembly (Lower Half) -->
+  <lattice id="100">
+    <type>rectangular</type>
+    <dimension>17 17</dimension>
+    <lower_left>-10.71 -10.71</lower_left>
+    <width>1.26 1.26</width>
+    <universes>
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1
+      1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
+      1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+    </universes>
+  </lattice>
+
+  <!-- Fuel Assembly (Upper Half) -->
+  <lattice id="101">
+    <type>rectangular</type>
+    <dimension>17 17</dimension>
+    <lower_left>-10.71 -10.71</lower_left>
+    <width>1.26 1.26</width>
+    <universes>
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 4 3 3 4 3 3 4 3 3 3 3 3
+      3 3 3 4 3 3 3 3 3 3 3 3 3 4 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 4 3 3 3 3 3 3 3 3 3 4 3 3 3
+      3 3 3 3 3 4 3 3 4 3 3 4 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+    </universes>
+  </lattice>
+
+  <!-- Core Lattice (Lower Half) -->
+  <lattice id="200">
+    <type>rectangular</type>
+    <dimension>21 21</dimension>
+    <lower_left>-224.91 -224.91</lower_left>
+    <width>21.42 21.42</width>
+    <universes>
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 6 6 6 6 6 6 6 5 5 5 5 5 5 5
+      5 5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5
+      5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5
+      5 5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5
+      5 5 5 5 5 5 5 6 6 6 6 6 6 6 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+    </universes>
+  </lattice>
+
+  <!-- Core Lattice (Upper Half) -->
+  <lattice id="201">
+    <type>rectangular</type>
+    <dimension>21 21</dimension>
+    <lower_left>-224.91 -224.91</lower_left>
+    <width>21.42 21.42</width>
+    <universes>
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 8 8 8 8 8 8 8 7 7 7 7 7 7 7
+      7 7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7 7
+      7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7
+      7 7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7 7
+      7 7 7 7 7 7 7 8 8 8 8 8 8 8 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+    </universes>
+  </lattice>
+
+</geometry>
\ No newline at end of file
diff --git a/tests/test_score_nuscatter_n/materials.xml b/tests/test_score_nuscatter_n/materials.xml
new file mode 100644
index 0000000000..9c0b74f3f1
--- /dev/null
+++ b/tests/test_score_nuscatter_n/materials.xml
@@ -0,0 +1,272 @@
+<?xml version="1.0"?>
+<materials>
+
+  <default_xs>71c</default_xs>
+
+  <!-- Fuel composition -->
+  <material id="1">
+    <density value="10.062" units="g/cm3" />
+    <nuclide name="U-234"   ao="4.9476e-6" />
+    <nuclide name="U-235"   ao="4.8218e-4" />
+    <nuclide name="U-236"   ao="9.0402e-5" />
+    <nuclide name="U-238"   ao="2.1504e-2" />
+    <nuclide name="Np-237"  ao="7.3733e-6" />
+    <nuclide name="Pu-238"  ao="1.5148e-6" />
+    <nuclide name="Pu-239"  ao="1.3955e-4" />
+    <nuclide name="Pu-240"  ao="3.4405e-5" />
+    <nuclide name="Pu-241"  ao="2.1439e-5" />
+    <nuclide name="Pu-242"  ao="3.7422e-6" />
+    <nuclide name="Am-241"  ao="4.5041e-7" />
+    <nuclide name="Am-242m" ao="9.2301e-9" />
+    <nuclide name="Am-243"  ao="4.7878e-7" />
+    <nuclide name="Cm-242"  ao="1.0485e-7" />
+    <nuclide name="Cm-243"  ao="1.4268e-9" />
+    <nuclide name="Cm-244"  ao="8.8756e-8" />
+    <nuclide name="Cm-245"  ao="3.5285e-9" />
+    <nuclide name="Mo-95"   ao="2.6497e-5" />
+    <nuclide name="Tc-99"   ao="3.2772e-5" />
+    <nuclide name="Ru-101"  ao="3.0742e-5" />
+    <nuclide name="Ru-103"  ao="2.3505e-6" />
+    <nuclide name="Ag-109"  ao="2.0009e-6" />
+    <nuclide name="Xe-135"  ao="1.0801e-8" />
+    <nuclide name="Cs-133"  ao="3.4612e-5" />
+    <nuclide name="Nd-143"  ao="2.6078e-5" />
+    <nuclide name="Nd-145"  ao="1.9898e-5" />
+    <nuclide name="Sm-147"  ao="1.6128e-6" />
+    <nuclide name="Sm-149"  ao="1.1627e-7" />
+    <nuclide name="Sm-150"  ao="7.1727e-6" />
+    <nuclide name="Sm-151"  ao="5.4947e-7" />
+    <nuclide name="Sm-152"  ao="3.0221e-6" />
+    <nuclide name="Eu-153"  ao="2.6209e-6" />
+    <nuclide name="Gd-155"  ao="1.5369e-9" />
+    <nuclide name="O-16"    ao="4.5737e-2" />
+  </material>
+
+  <!-- Cladding composition -->
+  <material id="2">
+    <density value="5.77" units="g/cm3" />
+    <nuclide name="Zr-90" ao="0.5145" />
+    <nuclide name="Zr-91" ao="0.1122" />
+    <nuclide name="Zr-92" ao="0.1715" />
+    <nuclide name="Zr-94" ao="0.1738" />
+    <nuclide name="Zr-96" ao="0.0280" />
+  </material>
+
+  <!-- Cold borated water -->
+  <material id="3">
+    <density value="0.07416" units="atom/b-cm" />
+    <nuclide name="H-1" ao="2.0" />
+    <nuclide name="O-16" ao="1.0" />
+    <nuclide name="B-10" ao="6.490e-4" />
+    <nuclide name="B-11" ao="2.689e-3" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Hot borated water -->
+  <material id="4">
+    <density value="0.06614" units="atom/b-cm" />
+    <nuclide name="H-1"  ao="2.0" />
+    <nuclide name="O-16" ao="1.0" />
+    <nuclide name="B-10" ao="6.490e-4" />
+    <nuclide name="B-11" ao="2.689e-3" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- RPV Composition -->
+  <material id="5">
+    <density value="7.9" units="g/cm3" />
+    <nuclide name="Fe-54"  wo="0.05437098" />
+    <nuclide name="Fe-56"  wo="0.88500663" />
+    <nuclide name="Fe-57"  wo="0.0208008" />
+    <nuclide name="Fe-58"  wo="0.00282159" />
+    <nuclide name="Ni-58"  wo="0.0067198" />
+    <nuclide name="Ni-60"  wo="0.0026776" />
+    <nuclide name="Ni-61"  wo="0.0001183" />
+    <nuclide name="Ni-62"  wo="0.0003835" />
+    <nuclide name="Ni-64"  wo="0.0001008" />
+    <nuclide name="Mn-55"  wo="0.01" />
+    <nuclide name="Mo-92"  wo="0.000849" />
+    <nuclide name="Mo-94"  wo="0.0005418" />
+    <nuclide name="Mo-95"  wo="0.0009438" />
+    <nuclide name="Mo-96"  wo="0.0010002" />
+    <nuclide name="Mo-97"  wo="0.0005796" />
+    <nuclide name="Mo-98"  wo="0.0014814" />
+    <nuclide name="Mo-100" wo="0.0006042" />
+    <nuclide name="Si-28"  wo="0.00367464" />
+    <nuclide name="Si-29"  wo="0.00019336" />
+    <nuclide name="Si-30"  wo="0.000132" />
+    <nuclide name="Cr-50"  wo="0.00010435" />
+    <nuclide name="Cr-52"  wo="0.002092475" />
+    <nuclide name="Cr-53"  wo="0.00024185" />
+    <nuclide name="Cr-54"  wo="6.1325e-05" />
+    <nuclide name="C-Nat"  wo="0.0025" />
+    <nuclide name="Cu-63"  wo="0.0013696" />
+    <nuclide name="Cu-65"  wo="0.0006304" />
+  </material>
+
+  <!-- Lower radial reflector -->
+  <material id="6">
+    <density value="4.32" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0095661" />
+    <nuclide name="O-16"  wo="0.0759107" />
+    <nuclide name="B-10"  wo="3.08409e-5" />
+    <nuclide name="B-11"  wo="1.40499e-4" />
+    <nuclide name="Fe-54" wo="0.035620772088" />
+    <nuclide name="Fe-56" wo="0.579805982228" />
+    <nuclide name="Fe-57" wo="0.01362750048" />
+    <nuclide name="Fe-58" wo="0.001848545204" />
+    <nuclide name="Ni-58" wo="0.055298376566" />
+    <nuclide name="Ni-60" wo="0.022034425592" />
+    <nuclide name="Ni-61" wo="0.000973510811" />
+    <nuclide name="Ni-62" wo="0.003155886695" />
+    <nuclide name="Ni-64" wo="0.000829500336" />
+    <nuclide name="Mn-55" wo="0.0182870" />
+    <nuclide name="Si-28" wo="0.00839976771" />
+    <nuclide name="Si-29" wo="0.00044199679" />
+    <nuclide name="Si-30" wo="0.0003017355" />
+    <nuclide name="Cr-50" wo="0.007251360806" />
+    <nuclide name="Cr-52" wo="0.145407678031" />
+    <nuclide name="Cr-53" wo="0.016806340306" />
+    <nuclide name="Cr-54" wo="0.004261520857" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Upper radial reflector / Top plate region -->
+  <material id="7">
+    <density value="4.28" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0086117" />
+    <nuclide name="O-16"  wo="0.0683369" />
+    <nuclide name="B-10"  wo="2.77638e-5" />
+    <nuclide name="B-11"  wo="1.26481e-4" />
+    <nuclide name="Fe-54" wo="0.035953677186" />
+    <nuclide name="Fe-56" wo="0.585224740891" />
+    <nuclide name="Fe-57" wo="0.01375486056" />
+    <nuclide name="Fe-58" wo="0.001865821363" />
+    <nuclide name="Ni-58" wo="0.055815129186" />
+    <nuclide name="Ni-60" wo="0.022240333032" />
+    <nuclide name="Ni-61" wo="0.000982608081" />
+    <nuclide name="Ni-62" wo="0.003185377845" />
+    <nuclide name="Ni-64" wo="0.000837251856" />
+    <nuclide name="Mn-55" wo="0.0184579" />
+    <nuclide name="Si-28" wo="0.00847831314" />
+    <nuclide name="Si-29" wo="0.00044612986" />
+    <nuclide name="Si-30" wo="0.000304557" />
+    <nuclide name="Cr-50" wo="0.00731912987" />
+    <nuclide name="Cr-52" wo="0.146766614995" />
+    <nuclide name="Cr-53" wo="0.01696340737" />
+    <nuclide name="Cr-54" wo="0.004301347765" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom plate region -->
+  <material id="8">
+    <density value="7.184" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0011505" />
+    <nuclide name="O-16"  wo="0.0091296" />
+    <nuclide name="B-10"  wo="3.70915e-6" />
+    <nuclide name="B-11"  wo="1.68974e-5" />
+    <nuclide name="Fe-54" wo="0.03855611055" />
+    <nuclide name="Fe-56" wo="0.627585036425" />
+    <nuclide name="Fe-57" wo="0.014750478" />
+    <nuclide name="Fe-58" wo="0.002000875025" />
+    <nuclide name="Ni-58" wo="0.059855207342" />
+    <nuclide name="Ni-60" wo="0.023850159704" />
+    <nuclide name="Ni-61" wo="0.001053732407" />
+    <nuclide name="Ni-62" wo="0.003415945715" />
+    <nuclide name="Ni-64" wo="0.000897854832" />
+    <nuclide name="Mn-55" wo="0.0197940" />
+    <nuclide name="Si-28" wo="0.00909197802" />
+    <nuclide name="Si-29" wo="0.00047842098" />
+    <nuclide name="Si-30" wo="0.000326601" />
+    <nuclide name="Cr-50" wo="0.007848910646" />
+    <nuclide name="Cr-52" wo="0.157390026871" />
+    <nuclide name="Cr-53" wo="0.018191270146" />
+    <nuclide name="Cr-54" wo="0.004612692337" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom nozzle region -->
+  <material id="9">
+    <density value="2.53" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0245014" />
+    <nuclide name="O-16"  wo="0.1944274" />
+    <nuclide name="B-10"  wo="7.89917e-5" />
+    <nuclide name="B-11"  wo="3.59854e-4" />
+    <nuclide name="Fe-54" wo="0.030411411144" />
+    <nuclide name="Fe-56" wo="0.495012237964" />
+    <nuclide name="Fe-57" wo="0.01163454624" />
+    <nuclide name="Fe-58" wo="0.001578204652" />
+    <nuclide name="Ni-58" wo="0.047211231662" />
+    <nuclide name="Ni-60" wo="0.018811987544" />
+    <nuclide name="Ni-61" wo="0.000831139127" />
+    <nuclide name="Ni-62" wo="0.002694352115" />
+    <nuclide name="Ni-64" wo="0.000708189552" />
+    <nuclide name="Mn-55" wo="0.0156126" />
+    <nuclide name="Si-28" wo="0.007171335558" />
+    <nuclide name="Si-29" wo="0.000377356542" />
+    <nuclide name="Si-30" wo="0.0002576079" />
+    <nuclide name="Cr-50" wo="0.006190885148" />
+    <nuclide name="Cr-52" wo="0.124142524198" />
+    <nuclide name="Cr-53" wo="0.014348496148" />
+    <nuclide name="Cr-54" wo="0.003638294506" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Top nozzle region -->
+  <material id="10">
+    <density value="1.746" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0358870" />
+    <nuclide name="O-16"  wo="0.2847761" />
+    <nuclide name="B-10"  wo="1.15699e-4" />
+    <nuclide name="B-11"  wo="5.27075e-4" />
+    <nuclide name="Fe-54" wo="0.02644016154" />
+    <nuclide name="Fe-56" wo="0.43037146399" />
+    <nuclide name="Fe-57" wo="0.0101152584" />
+    <nuclide name="Fe-58" wo="0.00137211607" />
+    <nuclide name="Ni-58" wo="0.04104621835" />
+    <nuclide name="Ni-60" wo="0.0163554502" />
+    <nuclide name="Ni-61" wo="0.000722605975" />
+    <nuclide name="Ni-62" wo="0.002342513875" />
+    <nuclide name="Ni-64" wo="0.0006157116" />
+    <nuclide name="Mn-55" wo="0.0135739" />
+    <nuclide name="Si-28" wo="0.006234853554" />
+    <nuclide name="Si-29" wo="0.000328078746" />
+    <nuclide name="Si-30" wo="0.0002239677" />
+    <nuclide name="Cr-50" wo="0.005382452306" />
+    <nuclide name="Cr-52" wo="0.107931450781" />
+    <nuclide name="Cr-53" wo="0.012474806806" />
+    <nuclide name="Cr-54" wo="0.003163190107" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Top of Fuel Assemblies -->
+  <material id="11">
+    <density value="3.044" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0162913" />
+    <nuclide name="O-16"  wo="0.1292776" />
+    <nuclide name="B-10"  wo="5.25228e-5" />
+    <nuclide name="B-11"  wo="2.39272e-4" />
+    <nuclide name="Zr-90" wo="0.43313403903" />
+    <nuclide name="Zr-91" wo="0.09549277374" />
+    <nuclide name="Zr-92" wo="0.14759527104" />
+    <nuclide name="Zr-94" wo="0.15280552077" />
+    <nuclide name="Zr-96" wo="0.02511169542" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom of Fuel Assemblies -->
+  <material id="12">
+    <density value="1.762" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0292856" />
+    <nuclide name="O-16"  wo="0.2323919" />
+    <nuclide name="B-10"  wo="9.44159e-5" />
+    <nuclide name="B-11"  wo="4.30120e-4" />
+    <nuclide name="Zr-90" wo="0.3741373658" />
+    <nuclide name="Zr-91" wo="0.0824858164" />
+    <nuclide name="Zr-92" wo="0.1274914944" />
+    <nuclide name="Zr-94" wo="0.1319920622" />
+    <nuclide name="Zr-96" wo="0.0216912612" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+</materials>
diff --git a/tests/test_score_nuscatter_n/results.py b/tests/test_score_nuscatter_n/results.py
new file mode 100644
index 0000000000..93371960ae
--- /dev/null
+++ b/tests/test_score_nuscatter_n/results.py
@@ -0,0 +1,37 @@
+#!/usr/bin/env python
+
+import sys
+import numpy as np
+
+# import statepoint
+sys.path.append('../../src/utils')
+import statepoint
+
+# read in statepoint file
+if len(sys.argv) > 1:
+    sp = statepoint.StatePoint(sys.argv[1])
+else:
+    sp = statepoint.StatePoint('statepoint.10.binary')
+sp.read_results()
+
+# extract tally results and convert to vector
+results = sp.tallies[0].results
+shape = results.shape
+size = (np.product(shape))
+results = np.reshape(results, size)
+
+# set up output string
+outstr = ''
+ 
+# write out k-combined
+outstr += 'k-combined:\n'
+outstr += "{0:12.6E} {1:12.6E}\n".format(sp.k_combined[0], sp.k_combined[1])
+
+# write out tally results
+outstr += 'tallies:\n'
+for item in results:
+  outstr += "{0:12.6E}\n".format(item)
+
+# write results to file
+with open('results_test.dat','w') as fh:
+    fh.write(outstr)
diff --git a/tests/test_score_nuscatter_n/settings.xml b/tests/test_score_nuscatter_n/settings.xml
new file mode 100644
index 0000000000..517637a59f
--- /dev/null
+++ b/tests/test_score_nuscatter_n/settings.xml
@@ -0,0 +1,19 @@
+<?xml version="1.0"?>
+<settings>
+
+  <eigenvalue>
+    <batches>10</batches>
+    <inactive>5</inactive>
+    <particles>100</particles>
+  </eigenvalue>
+
+  <source>
+    <space type="box">
+      <parameters>
+        -160 -160 -183
+         160  160  183
+      </parameters>
+    </space>
+  </source>
+
+</settings>
diff --git a/tests/test_score_nuscatter_n/tallies.xml b/tests/test_score_nuscatter_n/tallies.xml
new file mode 100644
index 0000000000..b962a81cc1
--- /dev/null
+++ b/tests/test_score_nuscatter_n/tallies.xml
@@ -0,0 +1,9 @@
+<?xml version="1.0"?>
+<tallies>
+
+  <tally id="1">
+    <filter type="cell" bins="21 22 23" />
+    <scores>nuscatter nuscatter-1 nuscatter-2 nuscatter-3 nuscatter-4</scores>
+  </tally>
+
+</tallies>
diff --git a/tests/test_score_nuscatter_n/test_score_nuscatter_n.py b/tests/test_score_nuscatter_n/test_score_nuscatter_n.py
new file mode 100644
index 0000000000..6e54389fc1
--- /dev/null
+++ b/tests/test_score_nuscatter_n/test_score_nuscatter_n.py
@@ -0,0 +1,62 @@
+#!/usr/bin/env python
+
+import os
+from subprocess import Popen, STDOUT, PIPE, call
+import filecmp
+import glob
+from optparse import OptionParser
+
+parser = OptionParser()
+parser.add_option('--mpi_exec', dest='mpi_exec', default='')
+parser.add_option('--mpi_np', dest='mpi_np', default='3')
+parser.add_option('--exe', dest='exe')
+(opts, args) = parser.parse_args()
+cwd = os.getcwd()
+
+def test_run():
+    if opts.mpi_exec != '':
+        proc = Popen([opts.mpi_exec, '-np', opts.mpi_np, opts.exe, cwd],
+               stderr=STDOUT, stdout=PIPE)
+    else:
+        proc = Popen([opts.exe, cwd], stderr=STDOUT, stdout=PIPE)
+    print(proc.communicate()[0])
+    returncode = proc.returncode
+    assert returncode == 0, 'OpenMC did not exit successfully.'
+
+def test_created_statepoint():
+    statepoint = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    assert len(statepoint) == 1, 'Either multiple or no statepoint files exist.'
+    assert statepoint[0].endswith('binary') or statepoint[0].endswith('h5'),\
+        'Statepoint file is not a binary or hdf5 file.'
+
+def test_output_exists():
+    assert os.path.exists(os.path.join(cwd, 'tallies.out')), 'Tally output file does not exist.'
+
+def test_results():
+    statepoint = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    call(['python', 'results.py', statepoint[0]])
+    compare = filecmp.cmp('results_test.dat', 'results_true.dat')
+    if not compare:
+      os.rename('results_test.dat', 'results_error.dat')
+    assert compare, 'Results do not agree.'
+
+def teardown():
+    output = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    output.append(os.path.join(cwd, 'tallies.out'))
+    output.append(os.path.join(cwd, 'results_test.dat'))
+    for f in output:
+        if os.path.exists(f):
+            os.remove(f)
+
+if __name__ == '__main__':
+
+    # test for openmc executable
+    if opts.exe is None:
+        raise Exception('Must specify OpenMC executable from command line with --exe.')
+
+    # run tests
+    test_run()
+    test_created_statepoint()
+    test_output_exists()
+    test_results()
+    teardown()
diff --git a/tests/test_score_nuscatter_pn/geometry.xml b/tests/test_score_nuscatter_pn/geometry.xml
new file mode 100644
index 0000000000..2407bf74b3
--- /dev/null
+++ b/tests/test_score_nuscatter_pn/geometry.xml
@@ -0,0 +1,185 @@
+<?xml version="1.0"?>
+<geometry>
+
+  <surface id="1" type="z-cylinder" coeffs="0. 0. 0.41" />
+  <surface id="2" type="z-cylinder" coeffs="0. 0. 0.475" />
+  <surface id="3" type="z-cylinder" coeffs="0. 0. 0.56" />
+  <surface id="4" type="z-cylinder" coeffs="0. 0. 0.62" />
+  <surface id="5" type="z-cylinder" coeffs="0. 0. 187.6" />
+  <surface id="6" type="z-cylinder" coeffs="0. 0. 209.0" />
+  <surface id="7" type="z-cylinder" coeffs="0. 0. 229.0" />
+  <surface id="8" type="z-cylinder" coeffs="0. 0. 249.0" boundary="vacuum" />
+
+  <surface id="31" type="z-plane" coeffs="-229.0" boundary="vacuum" />
+  <surface id="32" type="z-plane" coeffs="-199.0" />
+  <surface id="33" type="z-plane" coeffs="-193.0" />
+  <surface id="34" type="z-plane" coeffs="-183.0" />
+  <surface id="35" type="z-plane" coeffs="0.0" />
+  <surface id="36" type="z-plane" coeffs="183.0" />
+  <surface id="37" type="z-plane" coeffs="203.0" />
+  <surface id="38" type="z-plane" coeffs="215.0" />
+  <surface id="39" type="z-plane" coeffs="223.0" boundary="vacuum" />
+
+  <!-- All geometry on base universe -->
+  <cell id="1"  fill="200"    surfaces="  -6 34 -35" /> <!-- Lower core -->
+  <cell id="2"  fill="201"    surfaces="  -6 35 -36" /> <!-- Upper core -->
+  <cell id="3"  material="8"  surfaces="  -7 31 -32" /> <!-- Lower core plate region -->
+  <cell id="4"  material="9"  surfaces="  -5 32 -33" /> <!-- Bottom nozzle region -->
+  <cell id="5"  material="12" surfaces="  -5 33 -34" /> <!-- Bottom FA region -->
+  <cell id="6"  material="11" surfaces="  -5 36 -37" /> <!-- Top FA region -->
+  <cell id="7"  material="10" surfaces="  -5 37 -38" /> <!-- Top nozzle region -->
+  <cell id="8"  material="7"  surfaces="  -7 38 -39" /> <!-- Upper plate region -->
+  <cell id="9"  material="4"  surfaces="6 -7 32 -38" /> <!-- Downcomer -->
+  <cell id="10" material="5"  surfaces="7 -8 31 -39" /> <!-- RPV -->
+  <cell id="11" material="6"  surfaces="5 -6 32 -34" /> <!-- Bottom of radial reflector -->
+  <cell id="12" material="7"  surfaces="5 -6 36 -38" /> <!-- Top of radial reflector -->
+
+  <!-- Fuel pin, cladding, cold water -->
+  <cell id="21" universe="1" material="1" surfaces="-1"   />
+  <cell id="22" universe="1" material="2" surfaces="1 -2" />
+  <cell id="23" universe="1" material="3" surfaces="2"    />
+
+  <!-- Instrumentation guide tube -->
+  <cell id="24" universe="2" material="3" surfaces="-3"   />
+  <cell id="25" universe="2" material="2" surfaces="3 -4" />
+  <cell id="26" universe="2" material="3" surfaces="4"    />
+
+  <!-- Fuel pin, cladding, hot water -->
+  <cell id="27" universe="3" material="1" surfaces="-1"   />
+  <cell id="28" universe="3" material="2" surfaces="1 -2" />
+  <cell id="29" universe="3" material="4" surfaces="2"    />
+
+  <!-- Instrumentation guide tube -->
+  <cell id="30" universe="4" material="4" surfaces="-3"   />
+  <cell id="31" universe="4" material="2" surfaces="3 -4" />
+  <cell id="32" universe="4" material="4" surfaces="4"    />
+
+  <!-- cell for water assembly (cold) -->
+  <cell id="50" universe="5" material="4" surfaces="34 -35" />
+
+  <!-- containing cell for fuel assembly -->
+  <cell id="60" universe="6" fill="100"   surfaces="34 -35" />
+
+  <!-- cell for water assembly (hot) -->
+  <cell id="70" universe="7" material="3" surfaces="35 -36" />
+
+  <!-- containing cell for fuel assembly -->
+  <cell id="80" universe="8" fill="101"   surfaces="35 -36" />
+
+  <!-- Fuel Assembly (Lower Half) -->
+  <lattice id="100">
+    <type>rectangular</type>
+    <dimension>17 17</dimension>
+    <lower_left>-10.71 -10.71</lower_left>
+    <width>1.26 1.26</width>
+    <universes>
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1
+      1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
+      1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+    </universes>
+  </lattice>
+
+  <!-- Fuel Assembly (Upper Half) -->
+  <lattice id="101">
+    <type>rectangular</type>
+    <dimension>17 17</dimension>
+    <lower_left>-10.71 -10.71</lower_left>
+    <width>1.26 1.26</width>
+    <universes>
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 4 3 3 4 3 3 4 3 3 3 3 3
+      3 3 3 4 3 3 3 3 3 3 3 3 3 4 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 4 3 3 3 3 3 3 3 3 3 4 3 3 3
+      3 3 3 3 3 4 3 3 4 3 3 4 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+    </universes>
+  </lattice>
+
+  <!-- Core Lattice (Lower Half) -->
+  <lattice id="200">
+    <type>rectangular</type>
+    <dimension>21 21</dimension>
+    <lower_left>-224.91 -224.91</lower_left>
+    <width>21.42 21.42</width>
+    <universes>
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 6 6 6 6 6 6 6 5 5 5 5 5 5 5
+      5 5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5
+      5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5
+      5 5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5
+      5 5 5 5 5 5 5 6 6 6 6 6 6 6 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+    </universes>
+  </lattice>
+
+  <!-- Core Lattice (Upper Half) -->
+  <lattice id="201">
+    <type>rectangular</type>
+    <dimension>21 21</dimension>
+    <lower_left>-224.91 -224.91</lower_left>
+    <width>21.42 21.42</width>
+    <universes>
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 8 8 8 8 8 8 8 7 7 7 7 7 7 7
+      7 7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7 7
+      7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7
+      7 7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7 7
+      7 7 7 7 7 7 7 8 8 8 8 8 8 8 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+    </universes>
+  </lattice>
+
+</geometry>
\ No newline at end of file
diff --git a/tests/test_score_nuscatter_pn/materials.xml b/tests/test_score_nuscatter_pn/materials.xml
new file mode 100644
index 0000000000..9c0b74f3f1
--- /dev/null
+++ b/tests/test_score_nuscatter_pn/materials.xml
@@ -0,0 +1,272 @@
+<?xml version="1.0"?>
+<materials>
+
+  <default_xs>71c</default_xs>
+
+  <!-- Fuel composition -->
+  <material id="1">
+    <density value="10.062" units="g/cm3" />
+    <nuclide name="U-234"   ao="4.9476e-6" />
+    <nuclide name="U-235"   ao="4.8218e-4" />
+    <nuclide name="U-236"   ao="9.0402e-5" />
+    <nuclide name="U-238"   ao="2.1504e-2" />
+    <nuclide name="Np-237"  ao="7.3733e-6" />
+    <nuclide name="Pu-238"  ao="1.5148e-6" />
+    <nuclide name="Pu-239"  ao="1.3955e-4" />
+    <nuclide name="Pu-240"  ao="3.4405e-5" />
+    <nuclide name="Pu-241"  ao="2.1439e-5" />
+    <nuclide name="Pu-242"  ao="3.7422e-6" />
+    <nuclide name="Am-241"  ao="4.5041e-7" />
+    <nuclide name="Am-242m" ao="9.2301e-9" />
+    <nuclide name="Am-243"  ao="4.7878e-7" />
+    <nuclide name="Cm-242"  ao="1.0485e-7" />
+    <nuclide name="Cm-243"  ao="1.4268e-9" />
+    <nuclide name="Cm-244"  ao="8.8756e-8" />
+    <nuclide name="Cm-245"  ao="3.5285e-9" />
+    <nuclide name="Mo-95"   ao="2.6497e-5" />
+    <nuclide name="Tc-99"   ao="3.2772e-5" />
+    <nuclide name="Ru-101"  ao="3.0742e-5" />
+    <nuclide name="Ru-103"  ao="2.3505e-6" />
+    <nuclide name="Ag-109"  ao="2.0009e-6" />
+    <nuclide name="Xe-135"  ao="1.0801e-8" />
+    <nuclide name="Cs-133"  ao="3.4612e-5" />
+    <nuclide name="Nd-143"  ao="2.6078e-5" />
+    <nuclide name="Nd-145"  ao="1.9898e-5" />
+    <nuclide name="Sm-147"  ao="1.6128e-6" />
+    <nuclide name="Sm-149"  ao="1.1627e-7" />
+    <nuclide name="Sm-150"  ao="7.1727e-6" />
+    <nuclide name="Sm-151"  ao="5.4947e-7" />
+    <nuclide name="Sm-152"  ao="3.0221e-6" />
+    <nuclide name="Eu-153"  ao="2.6209e-6" />
+    <nuclide name="Gd-155"  ao="1.5369e-9" />
+    <nuclide name="O-16"    ao="4.5737e-2" />
+  </material>
+
+  <!-- Cladding composition -->
+  <material id="2">
+    <density value="5.77" units="g/cm3" />
+    <nuclide name="Zr-90" ao="0.5145" />
+    <nuclide name="Zr-91" ao="0.1122" />
+    <nuclide name="Zr-92" ao="0.1715" />
+    <nuclide name="Zr-94" ao="0.1738" />
+    <nuclide name="Zr-96" ao="0.0280" />
+  </material>
+
+  <!-- Cold borated water -->
+  <material id="3">
+    <density value="0.07416" units="atom/b-cm" />
+    <nuclide name="H-1" ao="2.0" />
+    <nuclide name="O-16" ao="1.0" />
+    <nuclide name="B-10" ao="6.490e-4" />
+    <nuclide name="B-11" ao="2.689e-3" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Hot borated water -->
+  <material id="4">
+    <density value="0.06614" units="atom/b-cm" />
+    <nuclide name="H-1"  ao="2.0" />
+    <nuclide name="O-16" ao="1.0" />
+    <nuclide name="B-10" ao="6.490e-4" />
+    <nuclide name="B-11" ao="2.689e-3" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- RPV Composition -->
+  <material id="5">
+    <density value="7.9" units="g/cm3" />
+    <nuclide name="Fe-54"  wo="0.05437098" />
+    <nuclide name="Fe-56"  wo="0.88500663" />
+    <nuclide name="Fe-57"  wo="0.0208008" />
+    <nuclide name="Fe-58"  wo="0.00282159" />
+    <nuclide name="Ni-58"  wo="0.0067198" />
+    <nuclide name="Ni-60"  wo="0.0026776" />
+    <nuclide name="Ni-61"  wo="0.0001183" />
+    <nuclide name="Ni-62"  wo="0.0003835" />
+    <nuclide name="Ni-64"  wo="0.0001008" />
+    <nuclide name="Mn-55"  wo="0.01" />
+    <nuclide name="Mo-92"  wo="0.000849" />
+    <nuclide name="Mo-94"  wo="0.0005418" />
+    <nuclide name="Mo-95"  wo="0.0009438" />
+    <nuclide name="Mo-96"  wo="0.0010002" />
+    <nuclide name="Mo-97"  wo="0.0005796" />
+    <nuclide name="Mo-98"  wo="0.0014814" />
+    <nuclide name="Mo-100" wo="0.0006042" />
+    <nuclide name="Si-28"  wo="0.00367464" />
+    <nuclide name="Si-29"  wo="0.00019336" />
+    <nuclide name="Si-30"  wo="0.000132" />
+    <nuclide name="Cr-50"  wo="0.00010435" />
+    <nuclide name="Cr-52"  wo="0.002092475" />
+    <nuclide name="Cr-53"  wo="0.00024185" />
+    <nuclide name="Cr-54"  wo="6.1325e-05" />
+    <nuclide name="C-Nat"  wo="0.0025" />
+    <nuclide name="Cu-63"  wo="0.0013696" />
+    <nuclide name="Cu-65"  wo="0.0006304" />
+  </material>
+
+  <!-- Lower radial reflector -->
+  <material id="6">
+    <density value="4.32" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0095661" />
+    <nuclide name="O-16"  wo="0.0759107" />
+    <nuclide name="B-10"  wo="3.08409e-5" />
+    <nuclide name="B-11"  wo="1.40499e-4" />
+    <nuclide name="Fe-54" wo="0.035620772088" />
+    <nuclide name="Fe-56" wo="0.579805982228" />
+    <nuclide name="Fe-57" wo="0.01362750048" />
+    <nuclide name="Fe-58" wo="0.001848545204" />
+    <nuclide name="Ni-58" wo="0.055298376566" />
+    <nuclide name="Ni-60" wo="0.022034425592" />
+    <nuclide name="Ni-61" wo="0.000973510811" />
+    <nuclide name="Ni-62" wo="0.003155886695" />
+    <nuclide name="Ni-64" wo="0.000829500336" />
+    <nuclide name="Mn-55" wo="0.0182870" />
+    <nuclide name="Si-28" wo="0.00839976771" />
+    <nuclide name="Si-29" wo="0.00044199679" />
+    <nuclide name="Si-30" wo="0.0003017355" />
+    <nuclide name="Cr-50" wo="0.007251360806" />
+    <nuclide name="Cr-52" wo="0.145407678031" />
+    <nuclide name="Cr-53" wo="0.016806340306" />
+    <nuclide name="Cr-54" wo="0.004261520857" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Upper radial reflector / Top plate region -->
+  <material id="7">
+    <density value="4.28" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0086117" />
+    <nuclide name="O-16"  wo="0.0683369" />
+    <nuclide name="B-10"  wo="2.77638e-5" />
+    <nuclide name="B-11"  wo="1.26481e-4" />
+    <nuclide name="Fe-54" wo="0.035953677186" />
+    <nuclide name="Fe-56" wo="0.585224740891" />
+    <nuclide name="Fe-57" wo="0.01375486056" />
+    <nuclide name="Fe-58" wo="0.001865821363" />
+    <nuclide name="Ni-58" wo="0.055815129186" />
+    <nuclide name="Ni-60" wo="0.022240333032" />
+    <nuclide name="Ni-61" wo="0.000982608081" />
+    <nuclide name="Ni-62" wo="0.003185377845" />
+    <nuclide name="Ni-64" wo="0.000837251856" />
+    <nuclide name="Mn-55" wo="0.0184579" />
+    <nuclide name="Si-28" wo="0.00847831314" />
+    <nuclide name="Si-29" wo="0.00044612986" />
+    <nuclide name="Si-30" wo="0.000304557" />
+    <nuclide name="Cr-50" wo="0.00731912987" />
+    <nuclide name="Cr-52" wo="0.146766614995" />
+    <nuclide name="Cr-53" wo="0.01696340737" />
+    <nuclide name="Cr-54" wo="0.004301347765" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom plate region -->
+  <material id="8">
+    <density value="7.184" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0011505" />
+    <nuclide name="O-16"  wo="0.0091296" />
+    <nuclide name="B-10"  wo="3.70915e-6" />
+    <nuclide name="B-11"  wo="1.68974e-5" />
+    <nuclide name="Fe-54" wo="0.03855611055" />
+    <nuclide name="Fe-56" wo="0.627585036425" />
+    <nuclide name="Fe-57" wo="0.014750478" />
+    <nuclide name="Fe-58" wo="0.002000875025" />
+    <nuclide name="Ni-58" wo="0.059855207342" />
+    <nuclide name="Ni-60" wo="0.023850159704" />
+    <nuclide name="Ni-61" wo="0.001053732407" />
+    <nuclide name="Ni-62" wo="0.003415945715" />
+    <nuclide name="Ni-64" wo="0.000897854832" />
+    <nuclide name="Mn-55" wo="0.0197940" />
+    <nuclide name="Si-28" wo="0.00909197802" />
+    <nuclide name="Si-29" wo="0.00047842098" />
+    <nuclide name="Si-30" wo="0.000326601" />
+    <nuclide name="Cr-50" wo="0.007848910646" />
+    <nuclide name="Cr-52" wo="0.157390026871" />
+    <nuclide name="Cr-53" wo="0.018191270146" />
+    <nuclide name="Cr-54" wo="0.004612692337" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom nozzle region -->
+  <material id="9">
+    <density value="2.53" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0245014" />
+    <nuclide name="O-16"  wo="0.1944274" />
+    <nuclide name="B-10"  wo="7.89917e-5" />
+    <nuclide name="B-11"  wo="3.59854e-4" />
+    <nuclide name="Fe-54" wo="0.030411411144" />
+    <nuclide name="Fe-56" wo="0.495012237964" />
+    <nuclide name="Fe-57" wo="0.01163454624" />
+    <nuclide name="Fe-58" wo="0.001578204652" />
+    <nuclide name="Ni-58" wo="0.047211231662" />
+    <nuclide name="Ni-60" wo="0.018811987544" />
+    <nuclide name="Ni-61" wo="0.000831139127" />
+    <nuclide name="Ni-62" wo="0.002694352115" />
+    <nuclide name="Ni-64" wo="0.000708189552" />
+    <nuclide name="Mn-55" wo="0.0156126" />
+    <nuclide name="Si-28" wo="0.007171335558" />
+    <nuclide name="Si-29" wo="0.000377356542" />
+    <nuclide name="Si-30" wo="0.0002576079" />
+    <nuclide name="Cr-50" wo="0.006190885148" />
+    <nuclide name="Cr-52" wo="0.124142524198" />
+    <nuclide name="Cr-53" wo="0.014348496148" />
+    <nuclide name="Cr-54" wo="0.003638294506" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Top nozzle region -->
+  <material id="10">
+    <density value="1.746" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0358870" />
+    <nuclide name="O-16"  wo="0.2847761" />
+    <nuclide name="B-10"  wo="1.15699e-4" />
+    <nuclide name="B-11"  wo="5.27075e-4" />
+    <nuclide name="Fe-54" wo="0.02644016154" />
+    <nuclide name="Fe-56" wo="0.43037146399" />
+    <nuclide name="Fe-57" wo="0.0101152584" />
+    <nuclide name="Fe-58" wo="0.00137211607" />
+    <nuclide name="Ni-58" wo="0.04104621835" />
+    <nuclide name="Ni-60" wo="0.0163554502" />
+    <nuclide name="Ni-61" wo="0.000722605975" />
+    <nuclide name="Ni-62" wo="0.002342513875" />
+    <nuclide name="Ni-64" wo="0.0006157116" />
+    <nuclide name="Mn-55" wo="0.0135739" />
+    <nuclide name="Si-28" wo="0.006234853554" />
+    <nuclide name="Si-29" wo="0.000328078746" />
+    <nuclide name="Si-30" wo="0.0002239677" />
+    <nuclide name="Cr-50" wo="0.005382452306" />
+    <nuclide name="Cr-52" wo="0.107931450781" />
+    <nuclide name="Cr-53" wo="0.012474806806" />
+    <nuclide name="Cr-54" wo="0.003163190107" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Top of Fuel Assemblies -->
+  <material id="11">
+    <density value="3.044" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0162913" />
+    <nuclide name="O-16"  wo="0.1292776" />
+    <nuclide name="B-10"  wo="5.25228e-5" />
+    <nuclide name="B-11"  wo="2.39272e-4" />
+    <nuclide name="Zr-90" wo="0.43313403903" />
+    <nuclide name="Zr-91" wo="0.09549277374" />
+    <nuclide name="Zr-92" wo="0.14759527104" />
+    <nuclide name="Zr-94" wo="0.15280552077" />
+    <nuclide name="Zr-96" wo="0.02511169542" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom of Fuel Assemblies -->
+  <material id="12">
+    <density value="1.762" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0292856" />
+    <nuclide name="O-16"  wo="0.2323919" />
+    <nuclide name="B-10"  wo="9.44159e-5" />
+    <nuclide name="B-11"  wo="4.30120e-4" />
+    <nuclide name="Zr-90" wo="0.3741373658" />
+    <nuclide name="Zr-91" wo="0.0824858164" />
+    <nuclide name="Zr-92" wo="0.1274914944" />
+    <nuclide name="Zr-94" wo="0.1319920622" />
+    <nuclide name="Zr-96" wo="0.0216912612" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+</materials>
diff --git a/tests/test_score_nuscatter_pn/results.py b/tests/test_score_nuscatter_pn/results.py
new file mode 100644
index 0000000000..f4c32a893f
--- /dev/null
+++ b/tests/test_score_nuscatter_pn/results.py
@@ -0,0 +1,44 @@
+#!/usr/bin/env python
+
+import sys
+import numpy as np
+
+# import statepoint
+sys.path.append('../../src/utils')
+import statepoint
+
+# read in statepoint file
+if len(sys.argv) > 1:
+    sp = statepoint.StatePoint(sys.argv[1])
+else:
+    sp = statepoint.StatePoint('statepoint.10.binary')
+sp.read_results()
+
+# extract tally results and convert to vector
+results1 = sp.tallies[0].results
+shape1 = results1.shape
+size1 = (np.product(shape1))
+results1 = np.reshape(results1, size1)
+results2 = sp.tallies[1].results
+shape2 = results2.shape
+size2 = (np.product(shape2))
+results2 = np.reshape(results2, size2)
+
+# set up output string
+outstr = ''
+ 
+# write out k-combined
+outstr += 'k-combined:\n'
+outstr += "{0:12.6E} {1:12.6E}\n".format(sp.k_combined[0], sp.k_combined[1])
+
+# write out tally results
+outstr += 'tally 1:\n'
+for item in results1:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'tally 2:\n'
+for item in results2:
+  outstr += "{0:12.6E}\n".format(item)
+
+# write results to file
+with open('results_test.dat','w') as fh:
+    fh.write(outstr)
diff --git a/tests/test_score_nuscatter_pn/settings.xml b/tests/test_score_nuscatter_pn/settings.xml
new file mode 100644
index 0000000000..517637a59f
--- /dev/null
+++ b/tests/test_score_nuscatter_pn/settings.xml
@@ -0,0 +1,19 @@
+<?xml version="1.0"?>
+<settings>
+
+  <eigenvalue>
+    <batches>10</batches>
+    <inactive>5</inactive>
+    <particles>100</particles>
+  </eigenvalue>
+
+  <source>
+    <space type="box">
+      <parameters>
+        -160 -160 -183
+         160  160  183
+      </parameters>
+    </space>
+  </source>
+
+</settings>
diff --git a/tests/test_score_nuscatter_pn/tallies.xml b/tests/test_score_nuscatter_pn/tallies.xml
new file mode 100644
index 0000000000..9914c7aa4b
--- /dev/null
+++ b/tests/test_score_nuscatter_pn/tallies.xml
@@ -0,0 +1,14 @@
+<?xml version="1.0"?>
+<tallies>
+
+  <tally id="1">
+    <filter type="cell" bins="21" />
+    <scores>nuscatter-0 nuscatter-1 nuscatter-2 nuscatter-3 nuscatter-4</scores>
+  </tally>
+
+  <tally id="2">
+    <filter type="cell" bins="21" />
+    <scores>nuscatter-p4</scores>
+  </tally>
+
+</tallies>
diff --git a/tests/test_score_nuscatter_pn/test_score_nuscatter_pn.py b/tests/test_score_nuscatter_pn/test_score_nuscatter_pn.py
new file mode 100644
index 0000000000..6e54389fc1
--- /dev/null
+++ b/tests/test_score_nuscatter_pn/test_score_nuscatter_pn.py
@@ -0,0 +1,62 @@
+#!/usr/bin/env python
+
+import os
+from subprocess import Popen, STDOUT, PIPE, call
+import filecmp
+import glob
+from optparse import OptionParser
+
+parser = OptionParser()
+parser.add_option('--mpi_exec', dest='mpi_exec', default='')
+parser.add_option('--mpi_np', dest='mpi_np', default='3')
+parser.add_option('--exe', dest='exe')
+(opts, args) = parser.parse_args()
+cwd = os.getcwd()
+
+def test_run():
+    if opts.mpi_exec != '':
+        proc = Popen([opts.mpi_exec, '-np', opts.mpi_np, opts.exe, cwd],
+               stderr=STDOUT, stdout=PIPE)
+    else:
+        proc = Popen([opts.exe, cwd], stderr=STDOUT, stdout=PIPE)
+    print(proc.communicate()[0])
+    returncode = proc.returncode
+    assert returncode == 0, 'OpenMC did not exit successfully.'
+
+def test_created_statepoint():
+    statepoint = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    assert len(statepoint) == 1, 'Either multiple or no statepoint files exist.'
+    assert statepoint[0].endswith('binary') or statepoint[0].endswith('h5'),\
+        'Statepoint file is not a binary or hdf5 file.'
+
+def test_output_exists():
+    assert os.path.exists(os.path.join(cwd, 'tallies.out')), 'Tally output file does not exist.'
+
+def test_results():
+    statepoint = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    call(['python', 'results.py', statepoint[0]])
+    compare = filecmp.cmp('results_test.dat', 'results_true.dat')
+    if not compare:
+      os.rename('results_test.dat', 'results_error.dat')
+    assert compare, 'Results do not agree.'
+
+def teardown():
+    output = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    output.append(os.path.join(cwd, 'tallies.out'))
+    output.append(os.path.join(cwd, 'results_test.dat'))
+    for f in output:
+        if os.path.exists(f):
+            os.remove(f)
+
+if __name__ == '__main__':
+
+    # test for openmc executable
+    if opts.exe is None:
+        raise Exception('Must specify OpenMC executable from command line with --exe.')
+
+    # run tests
+    test_run()
+    test_created_statepoint()
+    test_output_exists()
+    test_results()
+    teardown()
diff --git a/tests/test_score_nuscatter_yn/geometry.xml b/tests/test_score_nuscatter_yn/geometry.xml
new file mode 100644
index 0000000000..2407bf74b3
--- /dev/null
+++ b/tests/test_score_nuscatter_yn/geometry.xml
@@ -0,0 +1,185 @@
+<?xml version="1.0"?>
+<geometry>
+
+  <surface id="1" type="z-cylinder" coeffs="0. 0. 0.41" />
+  <surface id="2" type="z-cylinder" coeffs="0. 0. 0.475" />
+  <surface id="3" type="z-cylinder" coeffs="0. 0. 0.56" />
+  <surface id="4" type="z-cylinder" coeffs="0. 0. 0.62" />
+  <surface id="5" type="z-cylinder" coeffs="0. 0. 187.6" />
+  <surface id="6" type="z-cylinder" coeffs="0. 0. 209.0" />
+  <surface id="7" type="z-cylinder" coeffs="0. 0. 229.0" />
+  <surface id="8" type="z-cylinder" coeffs="0. 0. 249.0" boundary="vacuum" />
+
+  <surface id="31" type="z-plane" coeffs="-229.0" boundary="vacuum" />
+  <surface id="32" type="z-plane" coeffs="-199.0" />
+  <surface id="33" type="z-plane" coeffs="-193.0" />
+  <surface id="34" type="z-plane" coeffs="-183.0" />
+  <surface id="35" type="z-plane" coeffs="0.0" />
+  <surface id="36" type="z-plane" coeffs="183.0" />
+  <surface id="37" type="z-plane" coeffs="203.0" />
+  <surface id="38" type="z-plane" coeffs="215.0" />
+  <surface id="39" type="z-plane" coeffs="223.0" boundary="vacuum" />
+
+  <!-- All geometry on base universe -->
+  <cell id="1"  fill="200"    surfaces="  -6 34 -35" /> <!-- Lower core -->
+  <cell id="2"  fill="201"    surfaces="  -6 35 -36" /> <!-- Upper core -->
+  <cell id="3"  material="8"  surfaces="  -7 31 -32" /> <!-- Lower core plate region -->
+  <cell id="4"  material="9"  surfaces="  -5 32 -33" /> <!-- Bottom nozzle region -->
+  <cell id="5"  material="12" surfaces="  -5 33 -34" /> <!-- Bottom FA region -->
+  <cell id="6"  material="11" surfaces="  -5 36 -37" /> <!-- Top FA region -->
+  <cell id="7"  material="10" surfaces="  -5 37 -38" /> <!-- Top nozzle region -->
+  <cell id="8"  material="7"  surfaces="  -7 38 -39" /> <!-- Upper plate region -->
+  <cell id="9"  material="4"  surfaces="6 -7 32 -38" /> <!-- Downcomer -->
+  <cell id="10" material="5"  surfaces="7 -8 31 -39" /> <!-- RPV -->
+  <cell id="11" material="6"  surfaces="5 -6 32 -34" /> <!-- Bottom of radial reflector -->
+  <cell id="12" material="7"  surfaces="5 -6 36 -38" /> <!-- Top of radial reflector -->
+
+  <!-- Fuel pin, cladding, cold water -->
+  <cell id="21" universe="1" material="1" surfaces="-1"   />
+  <cell id="22" universe="1" material="2" surfaces="1 -2" />
+  <cell id="23" universe="1" material="3" surfaces="2"    />
+
+  <!-- Instrumentation guide tube -->
+  <cell id="24" universe="2" material="3" surfaces="-3"   />
+  <cell id="25" universe="2" material="2" surfaces="3 -4" />
+  <cell id="26" universe="2" material="3" surfaces="4"    />
+
+  <!-- Fuel pin, cladding, hot water -->
+  <cell id="27" universe="3" material="1" surfaces="-1"   />
+  <cell id="28" universe="3" material="2" surfaces="1 -2" />
+  <cell id="29" universe="3" material="4" surfaces="2"    />
+
+  <!-- Instrumentation guide tube -->
+  <cell id="30" universe="4" material="4" surfaces="-3"   />
+  <cell id="31" universe="4" material="2" surfaces="3 -4" />
+  <cell id="32" universe="4" material="4" surfaces="4"    />
+
+  <!-- cell for water assembly (cold) -->
+  <cell id="50" universe="5" material="4" surfaces="34 -35" />
+
+  <!-- containing cell for fuel assembly -->
+  <cell id="60" universe="6" fill="100"   surfaces="34 -35" />
+
+  <!-- cell for water assembly (hot) -->
+  <cell id="70" universe="7" material="3" surfaces="35 -36" />
+
+  <!-- containing cell for fuel assembly -->
+  <cell id="80" universe="8" fill="101"   surfaces="35 -36" />
+
+  <!-- Fuel Assembly (Lower Half) -->
+  <lattice id="100">
+    <type>rectangular</type>
+    <dimension>17 17</dimension>
+    <lower_left>-10.71 -10.71</lower_left>
+    <width>1.26 1.26</width>
+    <universes>
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1
+      1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
+      1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+    </universes>
+  </lattice>
+
+  <!-- Fuel Assembly (Upper Half) -->
+  <lattice id="101">
+    <type>rectangular</type>
+    <dimension>17 17</dimension>
+    <lower_left>-10.71 -10.71</lower_left>
+    <width>1.26 1.26</width>
+    <universes>
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 4 3 3 4 3 3 4 3 3 3 3 3
+      3 3 3 4 3 3 3 3 3 3 3 3 3 4 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 4 3 3 3 3 3 3 3 3 3 4 3 3 3
+      3 3 3 3 3 4 3 3 4 3 3 4 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+    </universes>
+  </lattice>
+
+  <!-- Core Lattice (Lower Half) -->
+  <lattice id="200">
+    <type>rectangular</type>
+    <dimension>21 21</dimension>
+    <lower_left>-224.91 -224.91</lower_left>
+    <width>21.42 21.42</width>
+    <universes>
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 6 6 6 6 6 6 6 5 5 5 5 5 5 5
+      5 5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5
+      5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5
+      5 5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5
+      5 5 5 5 5 5 5 6 6 6 6 6 6 6 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+    </universes>
+  </lattice>
+
+  <!-- Core Lattice (Upper Half) -->
+  <lattice id="201">
+    <type>rectangular</type>
+    <dimension>21 21</dimension>
+    <lower_left>-224.91 -224.91</lower_left>
+    <width>21.42 21.42</width>
+    <universes>
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 8 8 8 8 8 8 8 7 7 7 7 7 7 7
+      7 7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7 7
+      7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7
+      7 7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7 7
+      7 7 7 7 7 7 7 8 8 8 8 8 8 8 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+    </universes>
+  </lattice>
+
+</geometry>
\ No newline at end of file
diff --git a/tests/test_score_nuscatter_yn/materials.xml b/tests/test_score_nuscatter_yn/materials.xml
new file mode 100644
index 0000000000..9c0b74f3f1
--- /dev/null
+++ b/tests/test_score_nuscatter_yn/materials.xml
@@ -0,0 +1,272 @@
+<?xml version="1.0"?>
+<materials>
+
+  <default_xs>71c</default_xs>
+
+  <!-- Fuel composition -->
+  <material id="1">
+    <density value="10.062" units="g/cm3" />
+    <nuclide name="U-234"   ao="4.9476e-6" />
+    <nuclide name="U-235"   ao="4.8218e-4" />
+    <nuclide name="U-236"   ao="9.0402e-5" />
+    <nuclide name="U-238"   ao="2.1504e-2" />
+    <nuclide name="Np-237"  ao="7.3733e-6" />
+    <nuclide name="Pu-238"  ao="1.5148e-6" />
+    <nuclide name="Pu-239"  ao="1.3955e-4" />
+    <nuclide name="Pu-240"  ao="3.4405e-5" />
+    <nuclide name="Pu-241"  ao="2.1439e-5" />
+    <nuclide name="Pu-242"  ao="3.7422e-6" />
+    <nuclide name="Am-241"  ao="4.5041e-7" />
+    <nuclide name="Am-242m" ao="9.2301e-9" />
+    <nuclide name="Am-243"  ao="4.7878e-7" />
+    <nuclide name="Cm-242"  ao="1.0485e-7" />
+    <nuclide name="Cm-243"  ao="1.4268e-9" />
+    <nuclide name="Cm-244"  ao="8.8756e-8" />
+    <nuclide name="Cm-245"  ao="3.5285e-9" />
+    <nuclide name="Mo-95"   ao="2.6497e-5" />
+    <nuclide name="Tc-99"   ao="3.2772e-5" />
+    <nuclide name="Ru-101"  ao="3.0742e-5" />
+    <nuclide name="Ru-103"  ao="2.3505e-6" />
+    <nuclide name="Ag-109"  ao="2.0009e-6" />
+    <nuclide name="Xe-135"  ao="1.0801e-8" />
+    <nuclide name="Cs-133"  ao="3.4612e-5" />
+    <nuclide name="Nd-143"  ao="2.6078e-5" />
+    <nuclide name="Nd-145"  ao="1.9898e-5" />
+    <nuclide name="Sm-147"  ao="1.6128e-6" />
+    <nuclide name="Sm-149"  ao="1.1627e-7" />
+    <nuclide name="Sm-150"  ao="7.1727e-6" />
+    <nuclide name="Sm-151"  ao="5.4947e-7" />
+    <nuclide name="Sm-152"  ao="3.0221e-6" />
+    <nuclide name="Eu-153"  ao="2.6209e-6" />
+    <nuclide name="Gd-155"  ao="1.5369e-9" />
+    <nuclide name="O-16"    ao="4.5737e-2" />
+  </material>
+
+  <!-- Cladding composition -->
+  <material id="2">
+    <density value="5.77" units="g/cm3" />
+    <nuclide name="Zr-90" ao="0.5145" />
+    <nuclide name="Zr-91" ao="0.1122" />
+    <nuclide name="Zr-92" ao="0.1715" />
+    <nuclide name="Zr-94" ao="0.1738" />
+    <nuclide name="Zr-96" ao="0.0280" />
+  </material>
+
+  <!-- Cold borated water -->
+  <material id="3">
+    <density value="0.07416" units="atom/b-cm" />
+    <nuclide name="H-1" ao="2.0" />
+    <nuclide name="O-16" ao="1.0" />
+    <nuclide name="B-10" ao="6.490e-4" />
+    <nuclide name="B-11" ao="2.689e-3" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Hot borated water -->
+  <material id="4">
+    <density value="0.06614" units="atom/b-cm" />
+    <nuclide name="H-1"  ao="2.0" />
+    <nuclide name="O-16" ao="1.0" />
+    <nuclide name="B-10" ao="6.490e-4" />
+    <nuclide name="B-11" ao="2.689e-3" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- RPV Composition -->
+  <material id="5">
+    <density value="7.9" units="g/cm3" />
+    <nuclide name="Fe-54"  wo="0.05437098" />
+    <nuclide name="Fe-56"  wo="0.88500663" />
+    <nuclide name="Fe-57"  wo="0.0208008" />
+    <nuclide name="Fe-58"  wo="0.00282159" />
+    <nuclide name="Ni-58"  wo="0.0067198" />
+    <nuclide name="Ni-60"  wo="0.0026776" />
+    <nuclide name="Ni-61"  wo="0.0001183" />
+    <nuclide name="Ni-62"  wo="0.0003835" />
+    <nuclide name="Ni-64"  wo="0.0001008" />
+    <nuclide name="Mn-55"  wo="0.01" />
+    <nuclide name="Mo-92"  wo="0.000849" />
+    <nuclide name="Mo-94"  wo="0.0005418" />
+    <nuclide name="Mo-95"  wo="0.0009438" />
+    <nuclide name="Mo-96"  wo="0.0010002" />
+    <nuclide name="Mo-97"  wo="0.0005796" />
+    <nuclide name="Mo-98"  wo="0.0014814" />
+    <nuclide name="Mo-100" wo="0.0006042" />
+    <nuclide name="Si-28"  wo="0.00367464" />
+    <nuclide name="Si-29"  wo="0.00019336" />
+    <nuclide name="Si-30"  wo="0.000132" />
+    <nuclide name="Cr-50"  wo="0.00010435" />
+    <nuclide name="Cr-52"  wo="0.002092475" />
+    <nuclide name="Cr-53"  wo="0.00024185" />
+    <nuclide name="Cr-54"  wo="6.1325e-05" />
+    <nuclide name="C-Nat"  wo="0.0025" />
+    <nuclide name="Cu-63"  wo="0.0013696" />
+    <nuclide name="Cu-65"  wo="0.0006304" />
+  </material>
+
+  <!-- Lower radial reflector -->
+  <material id="6">
+    <density value="4.32" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0095661" />
+    <nuclide name="O-16"  wo="0.0759107" />
+    <nuclide name="B-10"  wo="3.08409e-5" />
+    <nuclide name="B-11"  wo="1.40499e-4" />
+    <nuclide name="Fe-54" wo="0.035620772088" />
+    <nuclide name="Fe-56" wo="0.579805982228" />
+    <nuclide name="Fe-57" wo="0.01362750048" />
+    <nuclide name="Fe-58" wo="0.001848545204" />
+    <nuclide name="Ni-58" wo="0.055298376566" />
+    <nuclide name="Ni-60" wo="0.022034425592" />
+    <nuclide name="Ni-61" wo="0.000973510811" />
+    <nuclide name="Ni-62" wo="0.003155886695" />
+    <nuclide name="Ni-64" wo="0.000829500336" />
+    <nuclide name="Mn-55" wo="0.0182870" />
+    <nuclide name="Si-28" wo="0.00839976771" />
+    <nuclide name="Si-29" wo="0.00044199679" />
+    <nuclide name="Si-30" wo="0.0003017355" />
+    <nuclide name="Cr-50" wo="0.007251360806" />
+    <nuclide name="Cr-52" wo="0.145407678031" />
+    <nuclide name="Cr-53" wo="0.016806340306" />
+    <nuclide name="Cr-54" wo="0.004261520857" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Upper radial reflector / Top plate region -->
+  <material id="7">
+    <density value="4.28" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0086117" />
+    <nuclide name="O-16"  wo="0.0683369" />
+    <nuclide name="B-10"  wo="2.77638e-5" />
+    <nuclide name="B-11"  wo="1.26481e-4" />
+    <nuclide name="Fe-54" wo="0.035953677186" />
+    <nuclide name="Fe-56" wo="0.585224740891" />
+    <nuclide name="Fe-57" wo="0.01375486056" />
+    <nuclide name="Fe-58" wo="0.001865821363" />
+    <nuclide name="Ni-58" wo="0.055815129186" />
+    <nuclide name="Ni-60" wo="0.022240333032" />
+    <nuclide name="Ni-61" wo="0.000982608081" />
+    <nuclide name="Ni-62" wo="0.003185377845" />
+    <nuclide name="Ni-64" wo="0.000837251856" />
+    <nuclide name="Mn-55" wo="0.0184579" />
+    <nuclide name="Si-28" wo="0.00847831314" />
+    <nuclide name="Si-29" wo="0.00044612986" />
+    <nuclide name="Si-30" wo="0.000304557" />
+    <nuclide name="Cr-50" wo="0.00731912987" />
+    <nuclide name="Cr-52" wo="0.146766614995" />
+    <nuclide name="Cr-53" wo="0.01696340737" />
+    <nuclide name="Cr-54" wo="0.004301347765" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom plate region -->
+  <material id="8">
+    <density value="7.184" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0011505" />
+    <nuclide name="O-16"  wo="0.0091296" />
+    <nuclide name="B-10"  wo="3.70915e-6" />
+    <nuclide name="B-11"  wo="1.68974e-5" />
+    <nuclide name="Fe-54" wo="0.03855611055" />
+    <nuclide name="Fe-56" wo="0.627585036425" />
+    <nuclide name="Fe-57" wo="0.014750478" />
+    <nuclide name="Fe-58" wo="0.002000875025" />
+    <nuclide name="Ni-58" wo="0.059855207342" />
+    <nuclide name="Ni-60" wo="0.023850159704" />
+    <nuclide name="Ni-61" wo="0.001053732407" />
+    <nuclide name="Ni-62" wo="0.003415945715" />
+    <nuclide name="Ni-64" wo="0.000897854832" />
+    <nuclide name="Mn-55" wo="0.0197940" />
+    <nuclide name="Si-28" wo="0.00909197802" />
+    <nuclide name="Si-29" wo="0.00047842098" />
+    <nuclide name="Si-30" wo="0.000326601" />
+    <nuclide name="Cr-50" wo="0.007848910646" />
+    <nuclide name="Cr-52" wo="0.157390026871" />
+    <nuclide name="Cr-53" wo="0.018191270146" />
+    <nuclide name="Cr-54" wo="0.004612692337" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom nozzle region -->
+  <material id="9">
+    <density value="2.53" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0245014" />
+    <nuclide name="O-16"  wo="0.1944274" />
+    <nuclide name="B-10"  wo="7.89917e-5" />
+    <nuclide name="B-11"  wo="3.59854e-4" />
+    <nuclide name="Fe-54" wo="0.030411411144" />
+    <nuclide name="Fe-56" wo="0.495012237964" />
+    <nuclide name="Fe-57" wo="0.01163454624" />
+    <nuclide name="Fe-58" wo="0.001578204652" />
+    <nuclide name="Ni-58" wo="0.047211231662" />
+    <nuclide name="Ni-60" wo="0.018811987544" />
+    <nuclide name="Ni-61" wo="0.000831139127" />
+    <nuclide name="Ni-62" wo="0.002694352115" />
+    <nuclide name="Ni-64" wo="0.000708189552" />
+    <nuclide name="Mn-55" wo="0.0156126" />
+    <nuclide name="Si-28" wo="0.007171335558" />
+    <nuclide name="Si-29" wo="0.000377356542" />
+    <nuclide name="Si-30" wo="0.0002576079" />
+    <nuclide name="Cr-50" wo="0.006190885148" />
+    <nuclide name="Cr-52" wo="0.124142524198" />
+    <nuclide name="Cr-53" wo="0.014348496148" />
+    <nuclide name="Cr-54" wo="0.003638294506" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Top nozzle region -->
+  <material id="10">
+    <density value="1.746" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0358870" />
+    <nuclide name="O-16"  wo="0.2847761" />
+    <nuclide name="B-10"  wo="1.15699e-4" />
+    <nuclide name="B-11"  wo="5.27075e-4" />
+    <nuclide name="Fe-54" wo="0.02644016154" />
+    <nuclide name="Fe-56" wo="0.43037146399" />
+    <nuclide name="Fe-57" wo="0.0101152584" />
+    <nuclide name="Fe-58" wo="0.00137211607" />
+    <nuclide name="Ni-58" wo="0.04104621835" />
+    <nuclide name="Ni-60" wo="0.0163554502" />
+    <nuclide name="Ni-61" wo="0.000722605975" />
+    <nuclide name="Ni-62" wo="0.002342513875" />
+    <nuclide name="Ni-64" wo="0.0006157116" />
+    <nuclide name="Mn-55" wo="0.0135739" />
+    <nuclide name="Si-28" wo="0.006234853554" />
+    <nuclide name="Si-29" wo="0.000328078746" />
+    <nuclide name="Si-30" wo="0.0002239677" />
+    <nuclide name="Cr-50" wo="0.005382452306" />
+    <nuclide name="Cr-52" wo="0.107931450781" />
+    <nuclide name="Cr-53" wo="0.012474806806" />
+    <nuclide name="Cr-54" wo="0.003163190107" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Top of Fuel Assemblies -->
+  <material id="11">
+    <density value="3.044" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0162913" />
+    <nuclide name="O-16"  wo="0.1292776" />
+    <nuclide name="B-10"  wo="5.25228e-5" />
+    <nuclide name="B-11"  wo="2.39272e-4" />
+    <nuclide name="Zr-90" wo="0.43313403903" />
+    <nuclide name="Zr-91" wo="0.09549277374" />
+    <nuclide name="Zr-92" wo="0.14759527104" />
+    <nuclide name="Zr-94" wo="0.15280552077" />
+    <nuclide name="Zr-96" wo="0.02511169542" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom of Fuel Assemblies -->
+  <material id="12">
+    <density value="1.762" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0292856" />
+    <nuclide name="O-16"  wo="0.2323919" />
+    <nuclide name="B-10"  wo="9.44159e-5" />
+    <nuclide name="B-11"  wo="4.30120e-4" />
+    <nuclide name="Zr-90" wo="0.3741373658" />
+    <nuclide name="Zr-91" wo="0.0824858164" />
+    <nuclide name="Zr-92" wo="0.1274914944" />
+    <nuclide name="Zr-94" wo="0.1319920622" />
+    <nuclide name="Zr-96" wo="0.0216912612" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+</materials>
diff --git a/tests/test_score_nuscatter_yn/results.py b/tests/test_score_nuscatter_yn/results.py
new file mode 100644
index 0000000000..f4c32a893f
--- /dev/null
+++ b/tests/test_score_nuscatter_yn/results.py
@@ -0,0 +1,44 @@
+#!/usr/bin/env python
+
+import sys
+import numpy as np
+
+# import statepoint
+sys.path.append('../../src/utils')
+import statepoint
+
+# read in statepoint file
+if len(sys.argv) > 1:
+    sp = statepoint.StatePoint(sys.argv[1])
+else:
+    sp = statepoint.StatePoint('statepoint.10.binary')
+sp.read_results()
+
+# extract tally results and convert to vector
+results1 = sp.tallies[0].results
+shape1 = results1.shape
+size1 = (np.product(shape1))
+results1 = np.reshape(results1, size1)
+results2 = sp.tallies[1].results
+shape2 = results2.shape
+size2 = (np.product(shape2))
+results2 = np.reshape(results2, size2)
+
+# set up output string
+outstr = ''
+ 
+# write out k-combined
+outstr += 'k-combined:\n'
+outstr += "{0:12.6E} {1:12.6E}\n".format(sp.k_combined[0], sp.k_combined[1])
+
+# write out tally results
+outstr += 'tally 1:\n'
+for item in results1:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'tally 2:\n'
+for item in results2:
+  outstr += "{0:12.6E}\n".format(item)
+
+# write results to file
+with open('results_test.dat','w') as fh:
+    fh.write(outstr)
diff --git a/tests/test_score_nuscatter_yn/settings.xml b/tests/test_score_nuscatter_yn/settings.xml
new file mode 100644
index 0000000000..517637a59f
--- /dev/null
+++ b/tests/test_score_nuscatter_yn/settings.xml
@@ -0,0 +1,19 @@
+<?xml version="1.0"?>
+<settings>
+
+  <eigenvalue>
+    <batches>10</batches>
+    <inactive>5</inactive>
+    <particles>100</particles>
+  </eigenvalue>
+
+  <source>
+    <space type="box">
+      <parameters>
+        -160 -160 -183
+         160  160  183
+      </parameters>
+    </space>
+  </source>
+
+</settings>
diff --git a/tests/test_score_nuscatter_yn/tallies.xml b/tests/test_score_nuscatter_yn/tallies.xml
new file mode 100644
index 0000000000..b9f977f091
--- /dev/null
+++ b/tests/test_score_nuscatter_yn/tallies.xml
@@ -0,0 +1,14 @@
+<?xml version="1.0"?>
+<tallies>
+
+  <tally id="1">
+    <filter type="cell" bins="21" />
+    <scores>nuscatter-0</scores>
+  </tally>
+
+  <tally id="2">
+    <filter type="cell" bins="21" />
+    <scores>nuscatter-y4</scores>
+  </tally>
+
+</tallies>
diff --git a/tests/test_score_nuscatter_yn/test_score_nuscatter_yn.py b/tests/test_score_nuscatter_yn/test_score_nuscatter_yn.py
new file mode 100644
index 0000000000..6e54389fc1
--- /dev/null
+++ b/tests/test_score_nuscatter_yn/test_score_nuscatter_yn.py
@@ -0,0 +1,62 @@
+#!/usr/bin/env python
+
+import os
+from subprocess import Popen, STDOUT, PIPE, call
+import filecmp
+import glob
+from optparse import OptionParser
+
+parser = OptionParser()
+parser.add_option('--mpi_exec', dest='mpi_exec', default='')
+parser.add_option('--mpi_np', dest='mpi_np', default='3')
+parser.add_option('--exe', dest='exe')
+(opts, args) = parser.parse_args()
+cwd = os.getcwd()
+
+def test_run():
+    if opts.mpi_exec != '':
+        proc = Popen([opts.mpi_exec, '-np', opts.mpi_np, opts.exe, cwd],
+               stderr=STDOUT, stdout=PIPE)
+    else:
+        proc = Popen([opts.exe, cwd], stderr=STDOUT, stdout=PIPE)
+    print(proc.communicate()[0])
+    returncode = proc.returncode
+    assert returncode == 0, 'OpenMC did not exit successfully.'
+
+def test_created_statepoint():
+    statepoint = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    assert len(statepoint) == 1, 'Either multiple or no statepoint files exist.'
+    assert statepoint[0].endswith('binary') or statepoint[0].endswith('h5'),\
+        'Statepoint file is not a binary or hdf5 file.'
+
+def test_output_exists():
+    assert os.path.exists(os.path.join(cwd, 'tallies.out')), 'Tally output file does not exist.'
+
+def test_results():
+    statepoint = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    call(['python', 'results.py', statepoint[0]])
+    compare = filecmp.cmp('results_test.dat', 'results_true.dat')
+    if not compare:
+      os.rename('results_test.dat', 'results_error.dat')
+    assert compare, 'Results do not agree.'
+
+def teardown():
+    output = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    output.append(os.path.join(cwd, 'tallies.out'))
+    output.append(os.path.join(cwd, 'results_test.dat'))
+    for f in output:
+        if os.path.exists(f):
+            os.remove(f)
+
+if __name__ == '__main__':
+
+    # test for openmc executable
+    if opts.exe is None:
+        raise Exception('Must specify OpenMC executable from command line with --exe.')
+
+    # run tests
+    test_run()
+    test_created_statepoint()
+    test_output_exists()
+    test_results()
+    teardown()
diff --git a/tests/test_score_scatter_yn/geometry.xml b/tests/test_score_scatter_yn/geometry.xml
new file mode 100644
index 0000000000..2407bf74b3
--- /dev/null
+++ b/tests/test_score_scatter_yn/geometry.xml
@@ -0,0 +1,185 @@
+<?xml version="1.0"?>
+<geometry>
+
+  <surface id="1" type="z-cylinder" coeffs="0. 0. 0.41" />
+  <surface id="2" type="z-cylinder" coeffs="0. 0. 0.475" />
+  <surface id="3" type="z-cylinder" coeffs="0. 0. 0.56" />
+  <surface id="4" type="z-cylinder" coeffs="0. 0. 0.62" />
+  <surface id="5" type="z-cylinder" coeffs="0. 0. 187.6" />
+  <surface id="6" type="z-cylinder" coeffs="0. 0. 209.0" />
+  <surface id="7" type="z-cylinder" coeffs="0. 0. 229.0" />
+  <surface id="8" type="z-cylinder" coeffs="0. 0. 249.0" boundary="vacuum" />
+
+  <surface id="31" type="z-plane" coeffs="-229.0" boundary="vacuum" />
+  <surface id="32" type="z-plane" coeffs="-199.0" />
+  <surface id="33" type="z-plane" coeffs="-193.0" />
+  <surface id="34" type="z-plane" coeffs="-183.0" />
+  <surface id="35" type="z-plane" coeffs="0.0" />
+  <surface id="36" type="z-plane" coeffs="183.0" />
+  <surface id="37" type="z-plane" coeffs="203.0" />
+  <surface id="38" type="z-plane" coeffs="215.0" />
+  <surface id="39" type="z-plane" coeffs="223.0" boundary="vacuum" />
+
+  <!-- All geometry on base universe -->
+  <cell id="1"  fill="200"    surfaces="  -6 34 -35" /> <!-- Lower core -->
+  <cell id="2"  fill="201"    surfaces="  -6 35 -36" /> <!-- Upper core -->
+  <cell id="3"  material="8"  surfaces="  -7 31 -32" /> <!-- Lower core plate region -->
+  <cell id="4"  material="9"  surfaces="  -5 32 -33" /> <!-- Bottom nozzle region -->
+  <cell id="5"  material="12" surfaces="  -5 33 -34" /> <!-- Bottom FA region -->
+  <cell id="6"  material="11" surfaces="  -5 36 -37" /> <!-- Top FA region -->
+  <cell id="7"  material="10" surfaces="  -5 37 -38" /> <!-- Top nozzle region -->
+  <cell id="8"  material="7"  surfaces="  -7 38 -39" /> <!-- Upper plate region -->
+  <cell id="9"  material="4"  surfaces="6 -7 32 -38" /> <!-- Downcomer -->
+  <cell id="10" material="5"  surfaces="7 -8 31 -39" /> <!-- RPV -->
+  <cell id="11" material="6"  surfaces="5 -6 32 -34" /> <!-- Bottom of radial reflector -->
+  <cell id="12" material="7"  surfaces="5 -6 36 -38" /> <!-- Top of radial reflector -->
+
+  <!-- Fuel pin, cladding, cold water -->
+  <cell id="21" universe="1" material="1" surfaces="-1"   />
+  <cell id="22" universe="1" material="2" surfaces="1 -2" />
+  <cell id="23" universe="1" material="3" surfaces="2"    />
+
+  <!-- Instrumentation guide tube -->
+  <cell id="24" universe="2" material="3" surfaces="-3"   />
+  <cell id="25" universe="2" material="2" surfaces="3 -4" />
+  <cell id="26" universe="2" material="3" surfaces="4"    />
+
+  <!-- Fuel pin, cladding, hot water -->
+  <cell id="27" universe="3" material="1" surfaces="-1"   />
+  <cell id="28" universe="3" material="2" surfaces="1 -2" />
+  <cell id="29" universe="3" material="4" surfaces="2"    />
+
+  <!-- Instrumentation guide tube -->
+  <cell id="30" universe="4" material="4" surfaces="-3"   />
+  <cell id="31" universe="4" material="2" surfaces="3 -4" />
+  <cell id="32" universe="4" material="4" surfaces="4"    />
+
+  <!-- cell for water assembly (cold) -->
+  <cell id="50" universe="5" material="4" surfaces="34 -35" />
+
+  <!-- containing cell for fuel assembly -->
+  <cell id="60" universe="6" fill="100"   surfaces="34 -35" />
+
+  <!-- cell for water assembly (hot) -->
+  <cell id="70" universe="7" material="3" surfaces="35 -36" />
+
+  <!-- containing cell for fuel assembly -->
+  <cell id="80" universe="8" fill="101"   surfaces="35 -36" />
+
+  <!-- Fuel Assembly (Lower Half) -->
+  <lattice id="100">
+    <type>rectangular</type>
+    <dimension>17 17</dimension>
+    <lower_left>-10.71 -10.71</lower_left>
+    <width>1.26 1.26</width>
+    <universes>
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1
+      1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
+      1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+    </universes>
+  </lattice>
+
+  <!-- Fuel Assembly (Upper Half) -->
+  <lattice id="101">
+    <type>rectangular</type>
+    <dimension>17 17</dimension>
+    <lower_left>-10.71 -10.71</lower_left>
+    <width>1.26 1.26</width>
+    <universes>
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 4 3 3 4 3 3 4 3 3 3 3 3
+      3 3 3 4 3 3 3 3 3 3 3 3 3 4 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 4 3 3 3 3 3 3 3 3 3 4 3 3 3
+      3 3 3 3 3 4 3 3 4 3 3 4 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+    </universes>
+  </lattice>
+
+  <!-- Core Lattice (Lower Half) -->
+  <lattice id="200">
+    <type>rectangular</type>
+    <dimension>21 21</dimension>
+    <lower_left>-224.91 -224.91</lower_left>
+    <width>21.42 21.42</width>
+    <universes>
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 6 6 6 6 6 6 6 5 5 5 5 5 5 5
+      5 5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5
+      5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5
+      5 5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5
+      5 5 5 5 5 5 5 6 6 6 6 6 6 6 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+    </universes>
+  </lattice>
+
+  <!-- Core Lattice (Upper Half) -->
+  <lattice id="201">
+    <type>rectangular</type>
+    <dimension>21 21</dimension>
+    <lower_left>-224.91 -224.91</lower_left>
+    <width>21.42 21.42</width>
+    <universes>
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 8 8 8 8 8 8 8 7 7 7 7 7 7 7
+      7 7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7 7
+      7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7
+      7 7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7 7
+      7 7 7 7 7 7 7 8 8 8 8 8 8 8 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+    </universes>
+  </lattice>
+
+</geometry>
\ No newline at end of file
diff --git a/tests/test_score_scatter_yn/materials.xml b/tests/test_score_scatter_yn/materials.xml
new file mode 100644
index 0000000000..9c0b74f3f1
--- /dev/null
+++ b/tests/test_score_scatter_yn/materials.xml
@@ -0,0 +1,272 @@
+<?xml version="1.0"?>
+<materials>
+
+  <default_xs>71c</default_xs>
+
+  <!-- Fuel composition -->
+  <material id="1">
+    <density value="10.062" units="g/cm3" />
+    <nuclide name="U-234"   ao="4.9476e-6" />
+    <nuclide name="U-235"   ao="4.8218e-4" />
+    <nuclide name="U-236"   ao="9.0402e-5" />
+    <nuclide name="U-238"   ao="2.1504e-2" />
+    <nuclide name="Np-237"  ao="7.3733e-6" />
+    <nuclide name="Pu-238"  ao="1.5148e-6" />
+    <nuclide name="Pu-239"  ao="1.3955e-4" />
+    <nuclide name="Pu-240"  ao="3.4405e-5" />
+    <nuclide name="Pu-241"  ao="2.1439e-5" />
+    <nuclide name="Pu-242"  ao="3.7422e-6" />
+    <nuclide name="Am-241"  ao="4.5041e-7" />
+    <nuclide name="Am-242m" ao="9.2301e-9" />
+    <nuclide name="Am-243"  ao="4.7878e-7" />
+    <nuclide name="Cm-242"  ao="1.0485e-7" />
+    <nuclide name="Cm-243"  ao="1.4268e-9" />
+    <nuclide name="Cm-244"  ao="8.8756e-8" />
+    <nuclide name="Cm-245"  ao="3.5285e-9" />
+    <nuclide name="Mo-95"   ao="2.6497e-5" />
+    <nuclide name="Tc-99"   ao="3.2772e-5" />
+    <nuclide name="Ru-101"  ao="3.0742e-5" />
+    <nuclide name="Ru-103"  ao="2.3505e-6" />
+    <nuclide name="Ag-109"  ao="2.0009e-6" />
+    <nuclide name="Xe-135"  ao="1.0801e-8" />
+    <nuclide name="Cs-133"  ao="3.4612e-5" />
+    <nuclide name="Nd-143"  ao="2.6078e-5" />
+    <nuclide name="Nd-145"  ao="1.9898e-5" />
+    <nuclide name="Sm-147"  ao="1.6128e-6" />
+    <nuclide name="Sm-149"  ao="1.1627e-7" />
+    <nuclide name="Sm-150"  ao="7.1727e-6" />
+    <nuclide name="Sm-151"  ao="5.4947e-7" />
+    <nuclide name="Sm-152"  ao="3.0221e-6" />
+    <nuclide name="Eu-153"  ao="2.6209e-6" />
+    <nuclide name="Gd-155"  ao="1.5369e-9" />
+    <nuclide name="O-16"    ao="4.5737e-2" />
+  </material>
+
+  <!-- Cladding composition -->
+  <material id="2">
+    <density value="5.77" units="g/cm3" />
+    <nuclide name="Zr-90" ao="0.5145" />
+    <nuclide name="Zr-91" ao="0.1122" />
+    <nuclide name="Zr-92" ao="0.1715" />
+    <nuclide name="Zr-94" ao="0.1738" />
+    <nuclide name="Zr-96" ao="0.0280" />
+  </material>
+
+  <!-- Cold borated water -->
+  <material id="3">
+    <density value="0.07416" units="atom/b-cm" />
+    <nuclide name="H-1" ao="2.0" />
+    <nuclide name="O-16" ao="1.0" />
+    <nuclide name="B-10" ao="6.490e-4" />
+    <nuclide name="B-11" ao="2.689e-3" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Hot borated water -->
+  <material id="4">
+    <density value="0.06614" units="atom/b-cm" />
+    <nuclide name="H-1"  ao="2.0" />
+    <nuclide name="O-16" ao="1.0" />
+    <nuclide name="B-10" ao="6.490e-4" />
+    <nuclide name="B-11" ao="2.689e-3" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- RPV Composition -->
+  <material id="5">
+    <density value="7.9" units="g/cm3" />
+    <nuclide name="Fe-54"  wo="0.05437098" />
+    <nuclide name="Fe-56"  wo="0.88500663" />
+    <nuclide name="Fe-57"  wo="0.0208008" />
+    <nuclide name="Fe-58"  wo="0.00282159" />
+    <nuclide name="Ni-58"  wo="0.0067198" />
+    <nuclide name="Ni-60"  wo="0.0026776" />
+    <nuclide name="Ni-61"  wo="0.0001183" />
+    <nuclide name="Ni-62"  wo="0.0003835" />
+    <nuclide name="Ni-64"  wo="0.0001008" />
+    <nuclide name="Mn-55"  wo="0.01" />
+    <nuclide name="Mo-92"  wo="0.000849" />
+    <nuclide name="Mo-94"  wo="0.0005418" />
+    <nuclide name="Mo-95"  wo="0.0009438" />
+    <nuclide name="Mo-96"  wo="0.0010002" />
+    <nuclide name="Mo-97"  wo="0.0005796" />
+    <nuclide name="Mo-98"  wo="0.0014814" />
+    <nuclide name="Mo-100" wo="0.0006042" />
+    <nuclide name="Si-28"  wo="0.00367464" />
+    <nuclide name="Si-29"  wo="0.00019336" />
+    <nuclide name="Si-30"  wo="0.000132" />
+    <nuclide name="Cr-50"  wo="0.00010435" />
+    <nuclide name="Cr-52"  wo="0.002092475" />
+    <nuclide name="Cr-53"  wo="0.00024185" />
+    <nuclide name="Cr-54"  wo="6.1325e-05" />
+    <nuclide name="C-Nat"  wo="0.0025" />
+    <nuclide name="Cu-63"  wo="0.0013696" />
+    <nuclide name="Cu-65"  wo="0.0006304" />
+  </material>
+
+  <!-- Lower radial reflector -->
+  <material id="6">
+    <density value="4.32" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0095661" />
+    <nuclide name="O-16"  wo="0.0759107" />
+    <nuclide name="B-10"  wo="3.08409e-5" />
+    <nuclide name="B-11"  wo="1.40499e-4" />
+    <nuclide name="Fe-54" wo="0.035620772088" />
+    <nuclide name="Fe-56" wo="0.579805982228" />
+    <nuclide name="Fe-57" wo="0.01362750048" />
+    <nuclide name="Fe-58" wo="0.001848545204" />
+    <nuclide name="Ni-58" wo="0.055298376566" />
+    <nuclide name="Ni-60" wo="0.022034425592" />
+    <nuclide name="Ni-61" wo="0.000973510811" />
+    <nuclide name="Ni-62" wo="0.003155886695" />
+    <nuclide name="Ni-64" wo="0.000829500336" />
+    <nuclide name="Mn-55" wo="0.0182870" />
+    <nuclide name="Si-28" wo="0.00839976771" />
+    <nuclide name="Si-29" wo="0.00044199679" />
+    <nuclide name="Si-30" wo="0.0003017355" />
+    <nuclide name="Cr-50" wo="0.007251360806" />
+    <nuclide name="Cr-52" wo="0.145407678031" />
+    <nuclide name="Cr-53" wo="0.016806340306" />
+    <nuclide name="Cr-54" wo="0.004261520857" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Upper radial reflector / Top plate region -->
+  <material id="7">
+    <density value="4.28" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0086117" />
+    <nuclide name="O-16"  wo="0.0683369" />
+    <nuclide name="B-10"  wo="2.77638e-5" />
+    <nuclide name="B-11"  wo="1.26481e-4" />
+    <nuclide name="Fe-54" wo="0.035953677186" />
+    <nuclide name="Fe-56" wo="0.585224740891" />
+    <nuclide name="Fe-57" wo="0.01375486056" />
+    <nuclide name="Fe-58" wo="0.001865821363" />
+    <nuclide name="Ni-58" wo="0.055815129186" />
+    <nuclide name="Ni-60" wo="0.022240333032" />
+    <nuclide name="Ni-61" wo="0.000982608081" />
+    <nuclide name="Ni-62" wo="0.003185377845" />
+    <nuclide name="Ni-64" wo="0.000837251856" />
+    <nuclide name="Mn-55" wo="0.0184579" />
+    <nuclide name="Si-28" wo="0.00847831314" />
+    <nuclide name="Si-29" wo="0.00044612986" />
+    <nuclide name="Si-30" wo="0.000304557" />
+    <nuclide name="Cr-50" wo="0.00731912987" />
+    <nuclide name="Cr-52" wo="0.146766614995" />
+    <nuclide name="Cr-53" wo="0.01696340737" />
+    <nuclide name="Cr-54" wo="0.004301347765" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom plate region -->
+  <material id="8">
+    <density value="7.184" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0011505" />
+    <nuclide name="O-16"  wo="0.0091296" />
+    <nuclide name="B-10"  wo="3.70915e-6" />
+    <nuclide name="B-11"  wo="1.68974e-5" />
+    <nuclide name="Fe-54" wo="0.03855611055" />
+    <nuclide name="Fe-56" wo="0.627585036425" />
+    <nuclide name="Fe-57" wo="0.014750478" />
+    <nuclide name="Fe-58" wo="0.002000875025" />
+    <nuclide name="Ni-58" wo="0.059855207342" />
+    <nuclide name="Ni-60" wo="0.023850159704" />
+    <nuclide name="Ni-61" wo="0.001053732407" />
+    <nuclide name="Ni-62" wo="0.003415945715" />
+    <nuclide name="Ni-64" wo="0.000897854832" />
+    <nuclide name="Mn-55" wo="0.0197940" />
+    <nuclide name="Si-28" wo="0.00909197802" />
+    <nuclide name="Si-29" wo="0.00047842098" />
+    <nuclide name="Si-30" wo="0.000326601" />
+    <nuclide name="Cr-50" wo="0.007848910646" />
+    <nuclide name="Cr-52" wo="0.157390026871" />
+    <nuclide name="Cr-53" wo="0.018191270146" />
+    <nuclide name="Cr-54" wo="0.004612692337" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom nozzle region -->
+  <material id="9">
+    <density value="2.53" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0245014" />
+    <nuclide name="O-16"  wo="0.1944274" />
+    <nuclide name="B-10"  wo="7.89917e-5" />
+    <nuclide name="B-11"  wo="3.59854e-4" />
+    <nuclide name="Fe-54" wo="0.030411411144" />
+    <nuclide name="Fe-56" wo="0.495012237964" />
+    <nuclide name="Fe-57" wo="0.01163454624" />
+    <nuclide name="Fe-58" wo="0.001578204652" />
+    <nuclide name="Ni-58" wo="0.047211231662" />
+    <nuclide name="Ni-60" wo="0.018811987544" />
+    <nuclide name="Ni-61" wo="0.000831139127" />
+    <nuclide name="Ni-62" wo="0.002694352115" />
+    <nuclide name="Ni-64" wo="0.000708189552" />
+    <nuclide name="Mn-55" wo="0.0156126" />
+    <nuclide name="Si-28" wo="0.007171335558" />
+    <nuclide name="Si-29" wo="0.000377356542" />
+    <nuclide name="Si-30" wo="0.0002576079" />
+    <nuclide name="Cr-50" wo="0.006190885148" />
+    <nuclide name="Cr-52" wo="0.124142524198" />
+    <nuclide name="Cr-53" wo="0.014348496148" />
+    <nuclide name="Cr-54" wo="0.003638294506" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Top nozzle region -->
+  <material id="10">
+    <density value="1.746" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0358870" />
+    <nuclide name="O-16"  wo="0.2847761" />
+    <nuclide name="B-10"  wo="1.15699e-4" />
+    <nuclide name="B-11"  wo="5.27075e-4" />
+    <nuclide name="Fe-54" wo="0.02644016154" />
+    <nuclide name="Fe-56" wo="0.43037146399" />
+    <nuclide name="Fe-57" wo="0.0101152584" />
+    <nuclide name="Fe-58" wo="0.00137211607" />
+    <nuclide name="Ni-58" wo="0.04104621835" />
+    <nuclide name="Ni-60" wo="0.0163554502" />
+    <nuclide name="Ni-61" wo="0.000722605975" />
+    <nuclide name="Ni-62" wo="0.002342513875" />
+    <nuclide name="Ni-64" wo="0.0006157116" />
+    <nuclide name="Mn-55" wo="0.0135739" />
+    <nuclide name="Si-28" wo="0.006234853554" />
+    <nuclide name="Si-29" wo="0.000328078746" />
+    <nuclide name="Si-30" wo="0.0002239677" />
+    <nuclide name="Cr-50" wo="0.005382452306" />
+    <nuclide name="Cr-52" wo="0.107931450781" />
+    <nuclide name="Cr-53" wo="0.012474806806" />
+    <nuclide name="Cr-54" wo="0.003163190107" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Top of Fuel Assemblies -->
+  <material id="11">
+    <density value="3.044" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0162913" />
+    <nuclide name="O-16"  wo="0.1292776" />
+    <nuclide name="B-10"  wo="5.25228e-5" />
+    <nuclide name="B-11"  wo="2.39272e-4" />
+    <nuclide name="Zr-90" wo="0.43313403903" />
+    <nuclide name="Zr-91" wo="0.09549277374" />
+    <nuclide name="Zr-92" wo="0.14759527104" />
+    <nuclide name="Zr-94" wo="0.15280552077" />
+    <nuclide name="Zr-96" wo="0.02511169542" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom of Fuel Assemblies -->
+  <material id="12">
+    <density value="1.762" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0292856" />
+    <nuclide name="O-16"  wo="0.2323919" />
+    <nuclide name="B-10"  wo="9.44159e-5" />
+    <nuclide name="B-11"  wo="4.30120e-4" />
+    <nuclide name="Zr-90" wo="0.3741373658" />
+    <nuclide name="Zr-91" wo="0.0824858164" />
+    <nuclide name="Zr-92" wo="0.1274914944" />
+    <nuclide name="Zr-94" wo="0.1319920622" />
+    <nuclide name="Zr-96" wo="0.0216912612" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+</materials>
diff --git a/tests/test_score_scatter_yn/results.py b/tests/test_score_scatter_yn/results.py
new file mode 100644
index 0000000000..f4c32a893f
--- /dev/null
+++ b/tests/test_score_scatter_yn/results.py
@@ -0,0 +1,44 @@
+#!/usr/bin/env python
+
+import sys
+import numpy as np
+
+# import statepoint
+sys.path.append('../../src/utils')
+import statepoint
+
+# read in statepoint file
+if len(sys.argv) > 1:
+    sp = statepoint.StatePoint(sys.argv[1])
+else:
+    sp = statepoint.StatePoint('statepoint.10.binary')
+sp.read_results()
+
+# extract tally results and convert to vector
+results1 = sp.tallies[0].results
+shape1 = results1.shape
+size1 = (np.product(shape1))
+results1 = np.reshape(results1, size1)
+results2 = sp.tallies[1].results
+shape2 = results2.shape
+size2 = (np.product(shape2))
+results2 = np.reshape(results2, size2)
+
+# set up output string
+outstr = ''
+ 
+# write out k-combined
+outstr += 'k-combined:\n'
+outstr += "{0:12.6E} {1:12.6E}\n".format(sp.k_combined[0], sp.k_combined[1])
+
+# write out tally results
+outstr += 'tally 1:\n'
+for item in results1:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'tally 2:\n'
+for item in results2:
+  outstr += "{0:12.6E}\n".format(item)
+
+# write results to file
+with open('results_test.dat','w') as fh:
+    fh.write(outstr)
diff --git a/tests/test_score_scatter_yn/settings.xml b/tests/test_score_scatter_yn/settings.xml
new file mode 100644
index 0000000000..517637a59f
--- /dev/null
+++ b/tests/test_score_scatter_yn/settings.xml
@@ -0,0 +1,19 @@
+<?xml version="1.0"?>
+<settings>
+
+  <eigenvalue>
+    <batches>10</batches>
+    <inactive>5</inactive>
+    <particles>100</particles>
+  </eigenvalue>
+
+  <source>
+    <space type="box">
+      <parameters>
+        -160 -160 -183
+         160  160  183
+      </parameters>
+    </space>
+  </source>
+
+</settings>
diff --git a/tests/test_score_scatter_yn/tallies.xml b/tests/test_score_scatter_yn/tallies.xml
new file mode 100644
index 0000000000..41d61e9503
--- /dev/null
+++ b/tests/test_score_scatter_yn/tallies.xml
@@ -0,0 +1,14 @@
+<?xml version="1.0"?>
+<tallies>
+
+  <tally id="1">
+    <filter type="cell" bins="21" />
+    <scores>scatter-0</scores>
+  </tally>
+
+  <tally id="2">
+    <filter type="cell" bins="21" />
+    <scores>scatter-y4</scores>
+  </tally>
+
+</tallies>
diff --git a/tests/test_score_scatter_yn/test_score_scatter_yn.py b/tests/test_score_scatter_yn/test_score_scatter_yn.py
new file mode 100644
index 0000000000..6e54389fc1
--- /dev/null
+++ b/tests/test_score_scatter_yn/test_score_scatter_yn.py
@@ -0,0 +1,62 @@
+#!/usr/bin/env python
+
+import os
+from subprocess import Popen, STDOUT, PIPE, call
+import filecmp
+import glob
+from optparse import OptionParser
+
+parser = OptionParser()
+parser.add_option('--mpi_exec', dest='mpi_exec', default='')
+parser.add_option('--mpi_np', dest='mpi_np', default='3')
+parser.add_option('--exe', dest='exe')
+(opts, args) = parser.parse_args()
+cwd = os.getcwd()
+
+def test_run():
+    if opts.mpi_exec != '':
+        proc = Popen([opts.mpi_exec, '-np', opts.mpi_np, opts.exe, cwd],
+               stderr=STDOUT, stdout=PIPE)
+    else:
+        proc = Popen([opts.exe, cwd], stderr=STDOUT, stdout=PIPE)
+    print(proc.communicate()[0])
+    returncode = proc.returncode
+    assert returncode == 0, 'OpenMC did not exit successfully.'
+
+def test_created_statepoint():
+    statepoint = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    assert len(statepoint) == 1, 'Either multiple or no statepoint files exist.'
+    assert statepoint[0].endswith('binary') or statepoint[0].endswith('h5'),\
+        'Statepoint file is not a binary or hdf5 file.'
+
+def test_output_exists():
+    assert os.path.exists(os.path.join(cwd, 'tallies.out')), 'Tally output file does not exist.'
+
+def test_results():
+    statepoint = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    call(['python', 'results.py', statepoint[0]])
+    compare = filecmp.cmp('results_test.dat', 'results_true.dat')
+    if not compare:
+      os.rename('results_test.dat', 'results_error.dat')
+    assert compare, 'Results do not agree.'
+
+def teardown():
+    output = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    output.append(os.path.join(cwd, 'tallies.out'))
+    output.append(os.path.join(cwd, 'results_test.dat'))
+    for f in output:
+        if os.path.exists(f):
+            os.remove(f)
+
+if __name__ == '__main__':
+
+    # test for openmc executable
+    if opts.exe is None:
+        raise Exception('Must specify OpenMC executable from command line with --exe.')
+
+    # run tests
+    test_run()
+    test_created_statepoint()
+    test_output_exists()
+    test_results()
+    teardown()
diff --git a/tests/test_score_total_yn/geometry.xml b/tests/test_score_total_yn/geometry.xml
new file mode 100644
index 0000000000..2407bf74b3
--- /dev/null
+++ b/tests/test_score_total_yn/geometry.xml
@@ -0,0 +1,185 @@
+<?xml version="1.0"?>
+<geometry>
+
+  <surface id="1" type="z-cylinder" coeffs="0. 0. 0.41" />
+  <surface id="2" type="z-cylinder" coeffs="0. 0. 0.475" />
+  <surface id="3" type="z-cylinder" coeffs="0. 0. 0.56" />
+  <surface id="4" type="z-cylinder" coeffs="0. 0. 0.62" />
+  <surface id="5" type="z-cylinder" coeffs="0. 0. 187.6" />
+  <surface id="6" type="z-cylinder" coeffs="0. 0. 209.0" />
+  <surface id="7" type="z-cylinder" coeffs="0. 0. 229.0" />
+  <surface id="8" type="z-cylinder" coeffs="0. 0. 249.0" boundary="vacuum" />
+
+  <surface id="31" type="z-plane" coeffs="-229.0" boundary="vacuum" />
+  <surface id="32" type="z-plane" coeffs="-199.0" />
+  <surface id="33" type="z-plane" coeffs="-193.0" />
+  <surface id="34" type="z-plane" coeffs="-183.0" />
+  <surface id="35" type="z-plane" coeffs="0.0" />
+  <surface id="36" type="z-plane" coeffs="183.0" />
+  <surface id="37" type="z-plane" coeffs="203.0" />
+  <surface id="38" type="z-plane" coeffs="215.0" />
+  <surface id="39" type="z-plane" coeffs="223.0" boundary="vacuum" />
+
+  <!-- All geometry on base universe -->
+  <cell id="1"  fill="200"    surfaces="  -6 34 -35" /> <!-- Lower core -->
+  <cell id="2"  fill="201"    surfaces="  -6 35 -36" /> <!-- Upper core -->
+  <cell id="3"  material="8"  surfaces="  -7 31 -32" /> <!-- Lower core plate region -->
+  <cell id="4"  material="9"  surfaces="  -5 32 -33" /> <!-- Bottom nozzle region -->
+  <cell id="5"  material="12" surfaces="  -5 33 -34" /> <!-- Bottom FA region -->
+  <cell id="6"  material="11" surfaces="  -5 36 -37" /> <!-- Top FA region -->
+  <cell id="7"  material="10" surfaces="  -5 37 -38" /> <!-- Top nozzle region -->
+  <cell id="8"  material="7"  surfaces="  -7 38 -39" /> <!-- Upper plate region -->
+  <cell id="9"  material="4"  surfaces="6 -7 32 -38" /> <!-- Downcomer -->
+  <cell id="10" material="5"  surfaces="7 -8 31 -39" /> <!-- RPV -->
+  <cell id="11" material="6"  surfaces="5 -6 32 -34" /> <!-- Bottom of radial reflector -->
+  <cell id="12" material="7"  surfaces="5 -6 36 -38" /> <!-- Top of radial reflector -->
+
+  <!-- Fuel pin, cladding, cold water -->
+  <cell id="21" universe="1" material="1" surfaces="-1"   />
+  <cell id="22" universe="1" material="2" surfaces="1 -2" />
+  <cell id="23" universe="1" material="3" surfaces="2"    />
+
+  <!-- Instrumentation guide tube -->
+  <cell id="24" universe="2" material="3" surfaces="-3"   />
+  <cell id="25" universe="2" material="2" surfaces="3 -4" />
+  <cell id="26" universe="2" material="3" surfaces="4"    />
+
+  <!-- Fuel pin, cladding, hot water -->
+  <cell id="27" universe="3" material="1" surfaces="-1"   />
+  <cell id="28" universe="3" material="2" surfaces="1 -2" />
+  <cell id="29" universe="3" material="4" surfaces="2"    />
+
+  <!-- Instrumentation guide tube -->
+  <cell id="30" universe="4" material="4" surfaces="-3"   />
+  <cell id="31" universe="4" material="2" surfaces="3 -4" />
+  <cell id="32" universe="4" material="4" surfaces="4"    />
+
+  <!-- cell for water assembly (cold) -->
+  <cell id="50" universe="5" material="4" surfaces="34 -35" />
+
+  <!-- containing cell for fuel assembly -->
+  <cell id="60" universe="6" fill="100"   surfaces="34 -35" />
+
+  <!-- cell for water assembly (hot) -->
+  <cell id="70" universe="7" material="3" surfaces="35 -36" />
+
+  <!-- containing cell for fuel assembly -->
+  <cell id="80" universe="8" fill="101"   surfaces="35 -36" />
+
+  <!-- Fuel Assembly (Lower Half) -->
+  <lattice id="100">
+    <type>rectangular</type>
+    <dimension>17 17</dimension>
+    <lower_left>-10.71 -10.71</lower_left>
+    <width>1.26 1.26</width>
+    <universes>
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1
+      1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
+      1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+      1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+    </universes>
+  </lattice>
+
+  <!-- Fuel Assembly (Upper Half) -->
+  <lattice id="101">
+    <type>rectangular</type>
+    <dimension>17 17</dimension>
+    <lower_left>-10.71 -10.71</lower_left>
+    <width>1.26 1.26</width>
+    <universes>
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 4 3 3 4 3 3 4 3 3 3 3 3
+      3 3 3 4 3 3 3 3 3 3 3 3 3 4 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 4 3 3 3 3 3 3 3 3 3 4 3 3 3
+      3 3 3 3 3 4 3 3 4 3 3 4 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+      3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
+    </universes>
+  </lattice>
+
+  <!-- Core Lattice (Lower Half) -->
+  <lattice id="200">
+    <type>rectangular</type>
+    <dimension>21 21</dimension>
+    <lower_left>-224.91 -224.91</lower_left>
+    <width>21.42 21.42</width>
+    <universes>
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 6 6 6 6 6 6 6 5 5 5 5 5 5 5
+      5 5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5
+      5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5
+      5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5
+      5 5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5
+      5 5 5 5 5 5 5 6 6 6 6 6 6 6 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+      5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+    </universes>
+  </lattice>
+
+  <!-- Core Lattice (Upper Half) -->
+  <lattice id="201">
+    <type>rectangular</type>
+    <dimension>21 21</dimension>
+    <lower_left>-224.91 -224.91</lower_left>
+    <width>21.42 21.42</width>
+    <universes>
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 8 8 8 8 8 8 8 7 7 7 7 7 7 7
+      7 7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7 7
+      7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7
+      7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7
+      7 7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7 7
+      7 7 7 7 7 7 7 8 8 8 8 8 8 8 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+      7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7
+    </universes>
+  </lattice>
+
+</geometry>
\ No newline at end of file
diff --git a/tests/test_score_total_yn/materials.xml b/tests/test_score_total_yn/materials.xml
new file mode 100644
index 0000000000..9c0b74f3f1
--- /dev/null
+++ b/tests/test_score_total_yn/materials.xml
@@ -0,0 +1,272 @@
+<?xml version="1.0"?>
+<materials>
+
+  <default_xs>71c</default_xs>
+
+  <!-- Fuel composition -->
+  <material id="1">
+    <density value="10.062" units="g/cm3" />
+    <nuclide name="U-234"   ao="4.9476e-6" />
+    <nuclide name="U-235"   ao="4.8218e-4" />
+    <nuclide name="U-236"   ao="9.0402e-5" />
+    <nuclide name="U-238"   ao="2.1504e-2" />
+    <nuclide name="Np-237"  ao="7.3733e-6" />
+    <nuclide name="Pu-238"  ao="1.5148e-6" />
+    <nuclide name="Pu-239"  ao="1.3955e-4" />
+    <nuclide name="Pu-240"  ao="3.4405e-5" />
+    <nuclide name="Pu-241"  ao="2.1439e-5" />
+    <nuclide name="Pu-242"  ao="3.7422e-6" />
+    <nuclide name="Am-241"  ao="4.5041e-7" />
+    <nuclide name="Am-242m" ao="9.2301e-9" />
+    <nuclide name="Am-243"  ao="4.7878e-7" />
+    <nuclide name="Cm-242"  ao="1.0485e-7" />
+    <nuclide name="Cm-243"  ao="1.4268e-9" />
+    <nuclide name="Cm-244"  ao="8.8756e-8" />
+    <nuclide name="Cm-245"  ao="3.5285e-9" />
+    <nuclide name="Mo-95"   ao="2.6497e-5" />
+    <nuclide name="Tc-99"   ao="3.2772e-5" />
+    <nuclide name="Ru-101"  ao="3.0742e-5" />
+    <nuclide name="Ru-103"  ao="2.3505e-6" />
+    <nuclide name="Ag-109"  ao="2.0009e-6" />
+    <nuclide name="Xe-135"  ao="1.0801e-8" />
+    <nuclide name="Cs-133"  ao="3.4612e-5" />
+    <nuclide name="Nd-143"  ao="2.6078e-5" />
+    <nuclide name="Nd-145"  ao="1.9898e-5" />
+    <nuclide name="Sm-147"  ao="1.6128e-6" />
+    <nuclide name="Sm-149"  ao="1.1627e-7" />
+    <nuclide name="Sm-150"  ao="7.1727e-6" />
+    <nuclide name="Sm-151"  ao="5.4947e-7" />
+    <nuclide name="Sm-152"  ao="3.0221e-6" />
+    <nuclide name="Eu-153"  ao="2.6209e-6" />
+    <nuclide name="Gd-155"  ao="1.5369e-9" />
+    <nuclide name="O-16"    ao="4.5737e-2" />
+  </material>
+
+  <!-- Cladding composition -->
+  <material id="2">
+    <density value="5.77" units="g/cm3" />
+    <nuclide name="Zr-90" ao="0.5145" />
+    <nuclide name="Zr-91" ao="0.1122" />
+    <nuclide name="Zr-92" ao="0.1715" />
+    <nuclide name="Zr-94" ao="0.1738" />
+    <nuclide name="Zr-96" ao="0.0280" />
+  </material>
+
+  <!-- Cold borated water -->
+  <material id="3">
+    <density value="0.07416" units="atom/b-cm" />
+    <nuclide name="H-1" ao="2.0" />
+    <nuclide name="O-16" ao="1.0" />
+    <nuclide name="B-10" ao="6.490e-4" />
+    <nuclide name="B-11" ao="2.689e-3" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Hot borated water -->
+  <material id="4">
+    <density value="0.06614" units="atom/b-cm" />
+    <nuclide name="H-1"  ao="2.0" />
+    <nuclide name="O-16" ao="1.0" />
+    <nuclide name="B-10" ao="6.490e-4" />
+    <nuclide name="B-11" ao="2.689e-3" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- RPV Composition -->
+  <material id="5">
+    <density value="7.9" units="g/cm3" />
+    <nuclide name="Fe-54"  wo="0.05437098" />
+    <nuclide name="Fe-56"  wo="0.88500663" />
+    <nuclide name="Fe-57"  wo="0.0208008" />
+    <nuclide name="Fe-58"  wo="0.00282159" />
+    <nuclide name="Ni-58"  wo="0.0067198" />
+    <nuclide name="Ni-60"  wo="0.0026776" />
+    <nuclide name="Ni-61"  wo="0.0001183" />
+    <nuclide name="Ni-62"  wo="0.0003835" />
+    <nuclide name="Ni-64"  wo="0.0001008" />
+    <nuclide name="Mn-55"  wo="0.01" />
+    <nuclide name="Mo-92"  wo="0.000849" />
+    <nuclide name="Mo-94"  wo="0.0005418" />
+    <nuclide name="Mo-95"  wo="0.0009438" />
+    <nuclide name="Mo-96"  wo="0.0010002" />
+    <nuclide name="Mo-97"  wo="0.0005796" />
+    <nuclide name="Mo-98"  wo="0.0014814" />
+    <nuclide name="Mo-100" wo="0.0006042" />
+    <nuclide name="Si-28"  wo="0.00367464" />
+    <nuclide name="Si-29"  wo="0.00019336" />
+    <nuclide name="Si-30"  wo="0.000132" />
+    <nuclide name="Cr-50"  wo="0.00010435" />
+    <nuclide name="Cr-52"  wo="0.002092475" />
+    <nuclide name="Cr-53"  wo="0.00024185" />
+    <nuclide name="Cr-54"  wo="6.1325e-05" />
+    <nuclide name="C-Nat"  wo="0.0025" />
+    <nuclide name="Cu-63"  wo="0.0013696" />
+    <nuclide name="Cu-65"  wo="0.0006304" />
+  </material>
+
+  <!-- Lower radial reflector -->
+  <material id="6">
+    <density value="4.32" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0095661" />
+    <nuclide name="O-16"  wo="0.0759107" />
+    <nuclide name="B-10"  wo="3.08409e-5" />
+    <nuclide name="B-11"  wo="1.40499e-4" />
+    <nuclide name="Fe-54" wo="0.035620772088" />
+    <nuclide name="Fe-56" wo="0.579805982228" />
+    <nuclide name="Fe-57" wo="0.01362750048" />
+    <nuclide name="Fe-58" wo="0.001848545204" />
+    <nuclide name="Ni-58" wo="0.055298376566" />
+    <nuclide name="Ni-60" wo="0.022034425592" />
+    <nuclide name="Ni-61" wo="0.000973510811" />
+    <nuclide name="Ni-62" wo="0.003155886695" />
+    <nuclide name="Ni-64" wo="0.000829500336" />
+    <nuclide name="Mn-55" wo="0.0182870" />
+    <nuclide name="Si-28" wo="0.00839976771" />
+    <nuclide name="Si-29" wo="0.00044199679" />
+    <nuclide name="Si-30" wo="0.0003017355" />
+    <nuclide name="Cr-50" wo="0.007251360806" />
+    <nuclide name="Cr-52" wo="0.145407678031" />
+    <nuclide name="Cr-53" wo="0.016806340306" />
+    <nuclide name="Cr-54" wo="0.004261520857" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Upper radial reflector / Top plate region -->
+  <material id="7">
+    <density value="4.28" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0086117" />
+    <nuclide name="O-16"  wo="0.0683369" />
+    <nuclide name="B-10"  wo="2.77638e-5" />
+    <nuclide name="B-11"  wo="1.26481e-4" />
+    <nuclide name="Fe-54" wo="0.035953677186" />
+    <nuclide name="Fe-56" wo="0.585224740891" />
+    <nuclide name="Fe-57" wo="0.01375486056" />
+    <nuclide name="Fe-58" wo="0.001865821363" />
+    <nuclide name="Ni-58" wo="0.055815129186" />
+    <nuclide name="Ni-60" wo="0.022240333032" />
+    <nuclide name="Ni-61" wo="0.000982608081" />
+    <nuclide name="Ni-62" wo="0.003185377845" />
+    <nuclide name="Ni-64" wo="0.000837251856" />
+    <nuclide name="Mn-55" wo="0.0184579" />
+    <nuclide name="Si-28" wo="0.00847831314" />
+    <nuclide name="Si-29" wo="0.00044612986" />
+    <nuclide name="Si-30" wo="0.000304557" />
+    <nuclide name="Cr-50" wo="0.00731912987" />
+    <nuclide name="Cr-52" wo="0.146766614995" />
+    <nuclide name="Cr-53" wo="0.01696340737" />
+    <nuclide name="Cr-54" wo="0.004301347765" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom plate region -->
+  <material id="8">
+    <density value="7.184" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0011505" />
+    <nuclide name="O-16"  wo="0.0091296" />
+    <nuclide name="B-10"  wo="3.70915e-6" />
+    <nuclide name="B-11"  wo="1.68974e-5" />
+    <nuclide name="Fe-54" wo="0.03855611055" />
+    <nuclide name="Fe-56" wo="0.627585036425" />
+    <nuclide name="Fe-57" wo="0.014750478" />
+    <nuclide name="Fe-58" wo="0.002000875025" />
+    <nuclide name="Ni-58" wo="0.059855207342" />
+    <nuclide name="Ni-60" wo="0.023850159704" />
+    <nuclide name="Ni-61" wo="0.001053732407" />
+    <nuclide name="Ni-62" wo="0.003415945715" />
+    <nuclide name="Ni-64" wo="0.000897854832" />
+    <nuclide name="Mn-55" wo="0.0197940" />
+    <nuclide name="Si-28" wo="0.00909197802" />
+    <nuclide name="Si-29" wo="0.00047842098" />
+    <nuclide name="Si-30" wo="0.000326601" />
+    <nuclide name="Cr-50" wo="0.007848910646" />
+    <nuclide name="Cr-52" wo="0.157390026871" />
+    <nuclide name="Cr-53" wo="0.018191270146" />
+    <nuclide name="Cr-54" wo="0.004612692337" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom nozzle region -->
+  <material id="9">
+    <density value="2.53" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0245014" />
+    <nuclide name="O-16"  wo="0.1944274" />
+    <nuclide name="B-10"  wo="7.89917e-5" />
+    <nuclide name="B-11"  wo="3.59854e-4" />
+    <nuclide name="Fe-54" wo="0.030411411144" />
+    <nuclide name="Fe-56" wo="0.495012237964" />
+    <nuclide name="Fe-57" wo="0.01163454624" />
+    <nuclide name="Fe-58" wo="0.001578204652" />
+    <nuclide name="Ni-58" wo="0.047211231662" />
+    <nuclide name="Ni-60" wo="0.018811987544" />
+    <nuclide name="Ni-61" wo="0.000831139127" />
+    <nuclide name="Ni-62" wo="0.002694352115" />
+    <nuclide name="Ni-64" wo="0.000708189552" />
+    <nuclide name="Mn-55" wo="0.0156126" />
+    <nuclide name="Si-28" wo="0.007171335558" />
+    <nuclide name="Si-29" wo="0.000377356542" />
+    <nuclide name="Si-30" wo="0.0002576079" />
+    <nuclide name="Cr-50" wo="0.006190885148" />
+    <nuclide name="Cr-52" wo="0.124142524198" />
+    <nuclide name="Cr-53" wo="0.014348496148" />
+    <nuclide name="Cr-54" wo="0.003638294506" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Top nozzle region -->
+  <material id="10">
+    <density value="1.746" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0358870" />
+    <nuclide name="O-16"  wo="0.2847761" />
+    <nuclide name="B-10"  wo="1.15699e-4" />
+    <nuclide name="B-11"  wo="5.27075e-4" />
+    <nuclide name="Fe-54" wo="0.02644016154" />
+    <nuclide name="Fe-56" wo="0.43037146399" />
+    <nuclide name="Fe-57" wo="0.0101152584" />
+    <nuclide name="Fe-58" wo="0.00137211607" />
+    <nuclide name="Ni-58" wo="0.04104621835" />
+    <nuclide name="Ni-60" wo="0.0163554502" />
+    <nuclide name="Ni-61" wo="0.000722605975" />
+    <nuclide name="Ni-62" wo="0.002342513875" />
+    <nuclide name="Ni-64" wo="0.0006157116" />
+    <nuclide name="Mn-55" wo="0.0135739" />
+    <nuclide name="Si-28" wo="0.006234853554" />
+    <nuclide name="Si-29" wo="0.000328078746" />
+    <nuclide name="Si-30" wo="0.0002239677" />
+    <nuclide name="Cr-50" wo="0.005382452306" />
+    <nuclide name="Cr-52" wo="0.107931450781" />
+    <nuclide name="Cr-53" wo="0.012474806806" />
+    <nuclide name="Cr-54" wo="0.003163190107" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Top of Fuel Assemblies -->
+  <material id="11">
+    <density value="3.044" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0162913" />
+    <nuclide name="O-16"  wo="0.1292776" />
+    <nuclide name="B-10"  wo="5.25228e-5" />
+    <nuclide name="B-11"  wo="2.39272e-4" />
+    <nuclide name="Zr-90" wo="0.43313403903" />
+    <nuclide name="Zr-91" wo="0.09549277374" />
+    <nuclide name="Zr-92" wo="0.14759527104" />
+    <nuclide name="Zr-94" wo="0.15280552077" />
+    <nuclide name="Zr-96" wo="0.02511169542" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+  <!-- Bottom of Fuel Assemblies -->
+  <material id="12">
+    <density value="1.762" units="g/cm3" />
+    <nuclide name="H-1"   wo="0.0292856" />
+    <nuclide name="O-16"  wo="0.2323919" />
+    <nuclide name="B-10"  wo="9.44159e-5" />
+    <nuclide name="B-11"  wo="4.30120e-4" />
+    <nuclide name="Zr-90" wo="0.3741373658" />
+    <nuclide name="Zr-91" wo="0.0824858164" />
+    <nuclide name="Zr-92" wo="0.1274914944" />
+    <nuclide name="Zr-94" wo="0.1319920622" />
+    <nuclide name="Zr-96" wo="0.0216912612" />
+    <sab name="HH2O" xs="71t" />
+  </material>
+
+</materials>
diff --git a/tests/test_score_total_yn/results.py b/tests/test_score_total_yn/results.py
new file mode 100644
index 0000000000..f4c32a893f
--- /dev/null
+++ b/tests/test_score_total_yn/results.py
@@ -0,0 +1,44 @@
+#!/usr/bin/env python
+
+import sys
+import numpy as np
+
+# import statepoint
+sys.path.append('../../src/utils')
+import statepoint
+
+# read in statepoint file
+if len(sys.argv) > 1:
+    sp = statepoint.StatePoint(sys.argv[1])
+else:
+    sp = statepoint.StatePoint('statepoint.10.binary')
+sp.read_results()
+
+# extract tally results and convert to vector
+results1 = sp.tallies[0].results
+shape1 = results1.shape
+size1 = (np.product(shape1))
+results1 = np.reshape(results1, size1)
+results2 = sp.tallies[1].results
+shape2 = results2.shape
+size2 = (np.product(shape2))
+results2 = np.reshape(results2, size2)
+
+# set up output string
+outstr = ''
+ 
+# write out k-combined
+outstr += 'k-combined:\n'
+outstr += "{0:12.6E} {1:12.6E}\n".format(sp.k_combined[0], sp.k_combined[1])
+
+# write out tally results
+outstr += 'tally 1:\n'
+for item in results1:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'tally 2:\n'
+for item in results2:
+  outstr += "{0:12.6E}\n".format(item)
+
+# write results to file
+with open('results_test.dat','w') as fh:
+    fh.write(outstr)
diff --git a/tests/test_score_total_yn/settings.xml b/tests/test_score_total_yn/settings.xml
new file mode 100644
index 0000000000..517637a59f
--- /dev/null
+++ b/tests/test_score_total_yn/settings.xml
@@ -0,0 +1,19 @@
+<?xml version="1.0"?>
+<settings>
+
+  <eigenvalue>
+    <batches>10</batches>
+    <inactive>5</inactive>
+    <particles>100</particles>
+  </eigenvalue>
+
+  <source>
+    <space type="box">
+      <parameters>
+        -160 -160 -183
+         160  160  183
+      </parameters>
+    </space>
+  </source>
+
+</settings>
diff --git a/tests/test_score_total_yn/tallies.xml b/tests/test_score_total_yn/tallies.xml
new file mode 100644
index 0000000000..df6de1c134
--- /dev/null
+++ b/tests/test_score_total_yn/tallies.xml
@@ -0,0 +1,14 @@
+<?xml version="1.0"?>
+<tallies>
+
+  <tally id="1">
+    <filter type="cell" bins="10 21 22 23" />
+    <scores>total</scores>
+  </tally>
+  
+  <tally id="2">
+    <filter type="cell" bins="10 21 22 23" />
+    <scores>total-y4</scores>
+  </tally>
+
+</tallies>
diff --git a/tests/test_score_total_yn/test_score_total_yn.py b/tests/test_score_total_yn/test_score_total_yn.py
new file mode 100644
index 0000000000..6e54389fc1
--- /dev/null
+++ b/tests/test_score_total_yn/test_score_total_yn.py
@@ -0,0 +1,62 @@
+#!/usr/bin/env python
+
+import os
+from subprocess import Popen, STDOUT, PIPE, call
+import filecmp
+import glob
+from optparse import OptionParser
+
+parser = OptionParser()
+parser.add_option('--mpi_exec', dest='mpi_exec', default='')
+parser.add_option('--mpi_np', dest='mpi_np', default='3')
+parser.add_option('--exe', dest='exe')
+(opts, args) = parser.parse_args()
+cwd = os.getcwd()
+
+def test_run():
+    if opts.mpi_exec != '':
+        proc = Popen([opts.mpi_exec, '-np', opts.mpi_np, opts.exe, cwd],
+               stderr=STDOUT, stdout=PIPE)
+    else:
+        proc = Popen([opts.exe, cwd], stderr=STDOUT, stdout=PIPE)
+    print(proc.communicate()[0])
+    returncode = proc.returncode
+    assert returncode == 0, 'OpenMC did not exit successfully.'
+
+def test_created_statepoint():
+    statepoint = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    assert len(statepoint) == 1, 'Either multiple or no statepoint files exist.'
+    assert statepoint[0].endswith('binary') or statepoint[0].endswith('h5'),\
+        'Statepoint file is not a binary or hdf5 file.'
+
+def test_output_exists():
+    assert os.path.exists(os.path.join(cwd, 'tallies.out')), 'Tally output file does not exist.'
+
+def test_results():
+    statepoint = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    call(['python', 'results.py', statepoint[0]])
+    compare = filecmp.cmp('results_test.dat', 'results_true.dat')
+    if not compare:
+      os.rename('results_test.dat', 'results_error.dat')
+    assert compare, 'Results do not agree.'
+
+def teardown():
+    output = glob.glob(os.path.join(cwd, 'statepoint.10.*'))
+    output.append(os.path.join(cwd, 'tallies.out'))
+    output.append(os.path.join(cwd, 'results_test.dat'))
+    for f in output:
+        if os.path.exists(f):
+            os.remove(f)
+
+if __name__ == '__main__':
+
+    # test for openmc executable
+    if opts.exe is None:
+        raise Exception('Must specify OpenMC executable from command line with --exe.')
+
+    # run tests
+    test_run()
+    test_created_statepoint()
+    test_output_exists()
+    test_results()
+    teardown()