diff --git a/tests/test_tally_nuclides/results.py b/tests/test_tally_nuclides/results.py
index 93371960a..3aa0946a5 100644
--- a/tests/test_tally_nuclides/results.py
+++ b/tests/test_tally_nuclides/results.py
@@ -15,10 +15,9 @@ else:
 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)
+results0 = sp.tallies[0].results
+results1 = sp.tallies[1].results
+results = np.concatenate((results0.flatten(), results1.flatten()))
 
 # set up output string
 outstr = ''
diff --git a/tests/test_tally_nuclides/results_true.dat b/tests/test_tally_nuclides/results_true.dat
index ce7400500..e429381dd 100644
--- a/tests/test_tally_nuclides/results_true.dat
+++ b/tests/test_tally_nuclides/results_true.dat
@@ -17,3 +17,11 @@ tallies:
 5.127366E-01
 5.602079E+00
 6.310407E+00
+7.251620E+00
+1.056750E+01
+1.649542E+00
+5.462124E-01
+1.598276E+00
+5.127366E-01
+5.602079E+00
+6.310407E+00
diff --git a/tests/test_tally_nuclides/tallies.xml b/tests/test_tally_nuclides/tallies.xml
index d8c294de1..cf20668c8 100644
--- a/tests/test_tally_nuclides/tallies.xml
+++ b/tests/test_tally_nuclides/tallies.xml
@@ -6,4 +6,9 @@
     <scores>total absorption fission scatter</scores>
   </tally>
 
+  <tally id="2">
+    <nuclides>Pu-239</nuclides>
+    <scores>total absorption fission scatter</scores>
+  </tally>
+
 </tallies>