added test suite for surface tally, with partial currents both ways (fuel to water and water to fuel) and couple of filters relevant to crossing surfaces

src/tally.F90

@@ -3138,7 +3138,7 @@ contains
   end subroutine score_collision_tally
 
 !===============================================================================
-! score_surface_tally is called at every surface crossing and can be used to 
+! score_surface_tally is called at every surface crossing and can be used to
 ! tally partial currents between two cells
 !===============================================================================
 

tests/test_surface_tally/geometry.xml

@@ -0,0 +1,11 @@
+<?xml version='1.0' encoding='utf-8'?>
+<geometry>
+    <cell id="10000" material="10000" name="fuel" region="-10000" universe="10000" />
+    <cell id="10001" material="10001" name="water" region="10000 -10002 10003 -10004" universe="10000" />
+    <cell fill="10000" id="10002" name="root cell" region="10001 -10002 10003 -10004" universe="0" />
+    <surface coeffs="0 0 0.4" id="10000" name="Fuel OR" type="z-cylinder" />
+    <surface boundary="reflective" coeffs="-0.62992" id="10001" name="left" type="x-plane" />
+    <surface boundary="reflective" coeffs="0.62992" id="10002" name="right" type="x-plane" />
+    <surface boundary="reflective" coeffs="-0.62992" id="10003" name="bottom" type="y-plane" />
+    <surface boundary="reflective" coeffs="0.62992" id="10004" name="top" type="y-plane" />
+</geometry>

tests/test_surface_tally/materials.xml

@@ -0,0 +1,21 @@
+<?xml version='1.0' encoding='utf-8'?>
+<materials>
+    <material id="10000" name="UO2 fuel at 2.4% wt enrichment">
+        <density units="g/cm3" value="10.29769" />
+        <nuclide ao="0.00019522896365555328" name="U234" scattering="iso-in-lab" />
+        <nuclide ao="0.024299485497345803" name="U235" scattering="iso-in-lab" />
+        <nuclide ao="0.9755052855389987" name="U238" scattering="iso-in-lab" />
+        <nuclide ao="1.999242" name="O16" scattering="iso-in-lab" />
+        <nuclide ao="0.000758" name="O17" scattering="iso-in-lab" />
+    </material>
+    <material id="10001" name="Borated water">
+        <density units="g/cm3" value="0.740582" />
+        <nuclide ao="7.928e-06" name="B10" scattering="iso-in-lab" />
+        <nuclide ao="3.2072e-05" name="B11" scattering="iso-in-lab" />
+        <nuclide ao="0.04999221300000001" name="H1" scattering="iso-in-lab" />
+        <nuclide ao="7.787e-06" name="H2" scattering="iso-in-lab" />
+        <nuclide ao="0.023990904" name="O16" scattering="iso-in-lab" />
+        <nuclide ao="9.096e-06" name="O17" scattering="iso-in-lab" />
+        <sab name="c_H_in_H2O" />
+    </material>
+</materials>

tests/test_surface_tally/results_true.dat

@@ -0,0 +1,36 @@
+k-combined:
+1.177844E+00 1.713522E-02
+tally 1:
+9.400000E-01
+1.778300E-01
+2.561000E+00
+1.315759E+00
+8.618000E+00
+1.487929E+01
+2.559600E+01
+1.311144E+02
+4.700000E-02
+4.870000E-04
+1.010000E-01
+2.207000E-03
+3.780000E-01
+2.889000E-02
+1.089000E+00
+2.432570E-01
+tally 2:
+9.070000E-01
+1.648450E-01
+2.495000E+00
+1.247251E+00
+8.515000E+00
+1.452994E+01
+2.539500E+01
+1.290580E+02
+3.800000E-02
+3.220000E-04
+8.400000E-02
+1.480000E-03
+3.540000E-01
+2.528400E-02
+9.770000E-01
+1.957910E-01

tests/test_surface_tally/settings.xml

@@ -0,0 +1,12 @@
+<?xml version='1.0' encoding='utf-8'?>
+<settings>
+    <run_mode>eigenvalue</run_mode>
+    <particles>1000</particles>
+    <batches>10</batches>
+    <inactive>5</inactive>
+    <source strength="1.0">
+        <space type="fission">
+            <parameters>-0.62992 -0.62992 -1 0.62992 0.62992 1</parameters>
+        </space>
+    </source>
+</settings>

tests/test_surface_tally/tallies.xml

@@ -0,0 +1,34 @@
+<?xml version='1.0' encoding='utf-8'?>
+<tallies>
+    <filter id="10003" type="cellfrom">
+        <bins>10000</bins>
+    </filter>
+    <filter id="10004" type="cell">
+        <bins>10001</bins>
+    </filter>
+    <filter id="10000" type="energy">
+        <bins>0.0 4000000.0 20000000.0</bins>
+    </filter>
+    <filter id="10001" type="polar">
+        <bins>0.0 0.785398163397 3.14159265359</bins>
+    </filter>
+    <filter id="10002" type="azimuthal">
+        <bins>0.0 0.785398163397 3.14159265359</bins>
+    </filter>
+    <filter id="10005" type="cellfrom">
+        <bins>10001</bins>
+    </filter>
+    <filter id="10006" type="cell">
+        <bins>10000</bins>
+    </filter>
+    <tally id="0" name="10000-10001">
+        <filters>10003 10004 10000 10001 10002</filters>
+        <scores>partial_current</scores>
+        <estimator>analog</estimator>
+    </tally>
+    <tally id="1" name="10001-10000">
+        <filters>10005 10006 10000 10001 10002</filters>
+        <scores>partial_current</scores>
+        <estimator>analog</estimator>
+    </tally>
+</tallies>

tests/test_surface_tally/test_surface_tally.py

@@ -0,0 +1,11 @@
+#!/usr/bin/env python
+
+import os
+import sys
+sys.path.insert(0, os.pardir)
+from testing_harness import TestHarness
+
+
+if __name__ == '__main__':
+    harness = TestHarness('statepoint.10.h5')
+    harness.main()