diff --git a/src/utils/statepoint.py b/src/utils/statepoint.py
index 5698506c92..96b98ab671 100644
--- a/src/utils/statepoint.py
+++ b/src/utils/statepoint.py
@@ -199,8 +199,10 @@ class StatePoint(object):
                 self.cmfd_indices = self._get_int(4, path='cmfd/indicies')
                 self.k_cmfd = self._get_double(self.current_batch,
                               path='cmfd/k_cmfd')
-                self.cmfd_src = self._get_double(np.product(self.cmfd_indices),
+                self.cmfd_src = self._get_double_array(np.product(self.cmfd_indices),
                                 path='cmfd/cmfd_src')
+                self.cmfd_src = np.reshape(self.cmfd_src,
+                                tuple(self.cmfd_indices), order='F')
                 self.cmfd_entropy = self._get_double(self.current_batch,
                                     path='cmfd/cmfd_entropy')
                 self.cmfd_balance = self._get_double(self.current_batch,
@@ -598,7 +600,7 @@ class StatePoint(object):
         else:
             return [float(v) for v in self._get_data(n, 'd', 8)]
 
-    def _get_double(self, n=1, path=None):
+    def _get_double_array(self, n=1, path=None):
         if self._hdf5:
             return self._f[path].value
         else:
diff --git a/tests/test_cmfd_feed/results.py b/tests/test_cmfd_feed/results.py
index 533a624558..67bbe6ca4f 100644
--- a/tests/test_cmfd_feed/results.py
+++ b/tests/test_cmfd_feed/results.py
@@ -53,6 +53,30 @@ outstr += 'tally 4:\n'
 for item in results4:
   outstr += "{0:12.6E}\n".format(item)
 
+# write out cmfd answers
+outstr += 'cmfd indices\n'
+for item in sp.cmfd_indices:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'k cmfd\n'
+for item in sp.k_cmfd:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'cmfd entropy\n'
+for item in sp.cmfd_entropy:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'cmfd balance\n'
+for item in sp.cmfd_balance:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'cmfd dominance ratio\n'
+for item in sp.cmfd_dominance:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'cmfd openmc source comparison\n'
+for item in sp.cmfd_srccmp:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'cmfd source\n'
+cmfdsrc = np.reshape(sp.cmfd_src, np.product(sp.cmfd_indices), order='F')
+for item in cmfdsrc:
+  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_cmfd_feed/results_true.dat b/tests/test_cmfd_feed/results_true.dat
index ad4c85c761..20ed778cd2 100644
--- a/tests/test_cmfd_feed/results_true.dat
+++ b/tests/test_cmfd_feed/results_true.dat
@@ -1,14 +1,14 @@
 k-combined:
-1.167124E+00 1.217343E-02
+1.167124E+00 1.217344E-02
 tally 1:
 1.126891E+01
-1.275755E+01
-2.086556E+01
-4.369501E+01
+1.275756E+01
+2.086557E+01
+4.369504E+01
 2.853629E+01
-8.169559E+01
-3.404506E+01
-1.165019E+02
+8.169563E+01
+3.404507E+01
+1.165020E+02
 3.723908E+01
 1.389881E+02
 3.760048E+01
@@ -16,86 +16,86 @@ tally 1:
 3.455469E+01
 1.197186E+02
 2.837503E+01
-8.086291E+01
-2.151789E+01
-4.648315E+01
+8.086287E+01
+2.151788E+01
+4.648312E+01
 1.194541E+01
-1.432331E+01
+1.432330E+01
 tally 2:
-2.292173E+01
-2.650792E+01
-1.598183E+01
-1.288692E+01
-2.194143E+00
-2.510557E-01
-4.217816E+01
-8.938583E+01
+2.292174E+01
+2.650794E+01
+1.598184E+01
+1.288693E+01
+2.194144E+00
+2.510560E-01
+4.217817E+01
+8.938588E+01
 2.974466E+01
-4.447038E+01
-3.939269E+00
-7.882413E-01
-5.742289E+01
+4.447040E+01
+3.939271E+00
+7.882419E-01
+5.742290E+01
 1.654599E+02
-4.067581E+01
-8.308454E+01
-5.566955E+00
-1.570072E+00
-6.810586E+01
+4.067582E+01
+8.308457E+01
+5.566956E+00
+1.570073E+00
+6.810587E+01
 2.331019E+02
-4.831543E+01
+4.831544E+01
 1.174673E+02
-6.244999E+00
-1.967938E+00
+6.245000E+00
+1.967939E+00
 7.258131E+01
 2.646756E+02
 5.185146E+01
-1.351425E+02
-6.622488E+00
+1.351426E+02
+6.622489E+00
 2.231695E+00
 7.246115E+01
 2.641295E+02
-5.159542E+01
+5.159541E+01
 1.339441E+02
 6.702239E+00
 2.277825E+00
-6.716397E+01
+6.716396E+01
 2.266815E+02
-4.765063E+01
-1.141627E+02
-6.384098E+00
+4.765062E+01
+1.141626E+02
+6.384097E+00
 2.063760E+00
-5.549542E+01
-1.545361E+02
+5.549541E+01
+1.545360E+02
 3.940819E+01
-7.794880E+01
-5.183375E+00
-1.369483E+00
-4.155685E+01
-8.679439E+01
-2.926744E+01
-4.307145E+01
-3.938852E+00
-7.877949E-01
-2.335180E+01
-2.758241E+01
-1.622683E+01
-1.331996E+01
-2.248921E+00
-2.594934E-01
+7.794877E+01
+5.183373E+00
+1.369482E+00
+4.155684E+01
+8.679434E+01
+2.926743E+01
+4.307143E+01
+3.938851E+00
+7.877944E-01
+2.335179E+01
+2.758239E+01
+1.622682E+01
+1.331995E+01
+2.248920E+00
+2.594933E-01
 tally 3:
 1.537807E+01
 1.193746E+01
 1.044189E+00
 5.627797E-02
 2.861614E+01
-4.118322E+01
+4.118324E+01
 1.846208E+00
 1.729007E-01
-3.921462E+01
-7.723576E+01
+3.921463E+01
+7.723579E+01
 2.427723E+00
 2.983153E-01
-4.653207E+01
+4.653208E+01
 1.089818E+02
 3.128394E+00
 4.945409E-01
@@ -111,15 +111,15 @@ tally 3:
 1.061019E+02
 2.936424E+00
 4.356307E-01
-3.783939E+01
-7.188489E+01
+3.783938E+01
+7.188486E+01
 2.485680E+00
 3.117618E-01
 2.819358E+01
-3.997588E+01
+3.997585E+01
 1.875908E+00
 1.772000E-01
-1.562421E+01
+1.562420E+01
 1.236020E+01
 1.013140E+00
 5.260420E-02
@@ -160,8 +160,8 @@ tally 4:
 0.000000E+00
 0.000000E+00
 0.000000E+00
-3.126365E+00
-4.930682E-01
+3.126366E+00
+4.930686E-01
 0.000000E+00
 0.000000E+00
 0.000000E+00
@@ -208,10 +208,10 @@ tally 4:
 0.000000E+00
 0.000000E+00
 0.000000E+00
-5.456410E+00
-1.499416E+00
-2.727755E+00
-3.788867E-01
+5.456412E+00
+1.499417E+00
+2.727756E+00
+3.788869E-01
 0.000000E+00
 0.000000E+00
 0.000000E+00
@@ -256,10 +256,10 @@ tally 4:
 0.000000E+00
 0.000000E+00
 0.000000E+00
-7.453480E+00
-2.790159E+00
-5.240103E+00
-1.386895E+00
+7.453481E+00
+2.790160E+00
+5.240104E+00
+1.386896E+00
 0.000000E+00
 0.000000E+00
 0.000000E+00
@@ -304,10 +304,10 @@ tally 4:
 0.000000E+00
 0.000000E+00
 0.000000E+00
-8.741563E+00
-3.830225E+00
-7.137968E+00
-2.555042E+00
+8.741565E+00
+3.830226E+00
+7.137969E+00
+2.555043E+00
 0.000000E+00
 0.000000E+00
 0.000000E+00
@@ -354,8 +354,8 @@ tally 4:
 0.000000E+00
 9.221731E+00
 4.266573E+00
-8.474116E+00
-3.604647E+00
+8.474118E+00
+3.604648E+00
 0.000000E+00
 0.000000E+00
 0.000000E+00
@@ -401,9 +401,9 @@ tally 4:
 0.000000E+00
 0.000000E+00
 9.054937E+00
-4.124537E+00
-9.066119E+00
-4.133833E+00
+4.124536E+00
+9.066120E+00
+4.133834E+00
 0.000000E+00
 0.000000E+00
 0.000000E+00
@@ -448,8 +448,8 @@ tally 4:
 0.000000E+00
 0.000000E+00
 0.000000E+00
-8.362253E+00
-3.517245E+00
+8.362251E+00
+3.517244E+00
 9.086548E+00
 4.156118E+00
 0.000000E+00
@@ -496,10 +496,10 @@ tally 4:
 0.000000E+00
 0.000000E+00
 0.000000E+00
-6.871573E+00
-2.369190E+00
-8.399584E+00
-3.549784E+00
+6.871571E+00
+2.369189E+00
+8.399583E+00
+3.549783E+00
 0.000000E+00
 0.000000E+00
 0.000000E+00
@@ -544,10 +544,10 @@ tally 4:
 0.000000E+00
 0.000000E+00
 0.000000E+00
-4.998103E+00
-1.258194E+00
-7.145296E+00
-2.568115E+00
+4.998101E+00
+1.258193E+00
+7.145294E+00
+2.568114E+00
 0.000000E+00
 0.000000E+00
 0.000000E+00
@@ -592,10 +592,10 @@ tally 4:
 0.000000E+00
 0.000000E+00
 0.000000E+00
-2.716076E+00
-3.753568E-01
-5.427593E+00
-1.484154E+00
+2.716075E+00
+3.753566E-01
+5.427591E+00
+1.484153E+00
 0.000000E+00
 0.000000E+00
 0.000000E+00
@@ -642,8 +642,8 @@ tally 4:
 0.000000E+00
 0.000000E+00
 0.000000E+00
-3.056645E+00
-4.702349E-01
+3.056644E+00
+4.702345E-01
 0.000000E+00
 0.000000E+00
 0.000000E+00
@@ -652,3 +652,124 @@ tally 4:
 0.000000E+00
 0.000000E+00
 0.000000E+00
+cmfd indices
+1.000000E+01
+1.000000E+00
+1.000000E+00
+1.000000E+00
+k cmfd
+0.000000E+00
+0.000000E+00
+0.000000E+00
+0.000000E+00
+1.158665E+00
+1.173348E+00
+1.180514E+00
+1.168646E+00
+1.152485E+00
+1.154989E+00
+1.152285E+00
+1.152275E+00
+1.149476E+00
+1.154912E+00
+1.164393E+00
+1.170145E+00
+1.169161E+00
+1.167966E+00
+1.170153E+00
+1.170989E+00
+cmfd entropy
+0.000000E+00
+0.000000E+00
+0.000000E+00
+0.000000E+00
+3.229801E+00
+3.226763E+00
+3.225011E+00
+3.230316E+00
+3.232043E+00
+3.228854E+00
+3.227439E+00
+3.228206E+00
+3.226749E+00
+3.222635E+00
+3.218882E+00
+3.219196E+00
+3.218057E+00
+3.220324E+00
+3.219099E+00
+3.219847E+00
+cmfd balance
+0.000000E+00
+0.000000E+00
+0.000000E+00
+0.000000E+00
+6.244292E-03
+5.210436E-03
+4.212224E-03
+3.321578E-03
+3.067254E-03
+3.076634E-03
+2.601637E-03
+2.666441E-03
+2.490817E-03
+2.571411E-03
+3.013819E-03
+2.802836E-03
+2.725356E-03
+2.475281E-03
+2.279123E-03
+2.110816E-03
+cmfd dominance ratio
+0.000000E+00
+0.000000E+00
+0.000000E+00
+0.000000E+00
+5.509393E-01
+5.472954E-01
+5.461377E-01
+5.518823E-01
+5.540420E-01
+5.513822E-01
+5.505071E-01
+5.498933E-01
+5.486126E-01
+5.450413E-01
+5.437663E-01
+5.421995E-01
+5.417191E-01
+5.437232E-01
+5.429525E-01
+5.428436E-01
+cmfd openmc source comparison
+0.000000E+00
+0.000000E+00
+0.000000E+00
+0.000000E+00
+1.136193E-02
+1.059214E-02
+9.675100E-03
+7.296572E-03
+5.812484E-03
+6.138879E-03
+6.070297E-03
+5.791570E-03
+6.020032E-03
+7.287183E-03
+7.068423E-03
+6.380521E-03
+5.828010E-03
+4.741191E-03
+4.848114E-03
+4.686772E-03
+cmfd source
+4.367444E-02
+7.866919E-02
+1.050799E-01
+1.366581E-01
+1.421824E-01
+1.393453E-01
+1.269612E-01
+1.060419E-01
+7.920472E-02
+4.218277E-02
diff --git a/tests/test_cmfd_nofeed/results.py b/tests/test_cmfd_nofeed/results.py
index 533a624558..67bbe6ca4f 100644
--- a/tests/test_cmfd_nofeed/results.py
+++ b/tests/test_cmfd_nofeed/results.py
@@ -53,6 +53,30 @@ outstr += 'tally 4:\n'
 for item in results4:
   outstr += "{0:12.6E}\n".format(item)
 
+# write out cmfd answers
+outstr += 'cmfd indices\n'
+for item in sp.cmfd_indices:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'k cmfd\n'
+for item in sp.k_cmfd:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'cmfd entropy\n'
+for item in sp.cmfd_entropy:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'cmfd balance\n'
+for item in sp.cmfd_balance:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'cmfd dominance ratio\n'
+for item in sp.cmfd_dominance:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'cmfd openmc source comparison\n'
+for item in sp.cmfd_srccmp:
+  outstr += "{0:12.6E}\n".format(item)
+outstr += 'cmfd source\n'
+cmfdsrc = np.reshape(sp.cmfd_src, np.product(sp.cmfd_indices), order='F')
+for item in cmfdsrc:
+  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_cmfd_nofeed/results_true.dat b/tests/test_cmfd_nofeed/results_true.dat
index 6f00290374..be594a3fba 100644
--- a/tests/test_cmfd_nofeed/results_true.dat
+++ b/tests/test_cmfd_nofeed/results_true.dat
@@ -652,3 +652,124 @@ tally 4:
 0.000000E+00
 0.000000E+00
 0.000000E+00
+cmfd indices
+1.000000E+01
+1.000000E+00
+1.000000E+00
+1.000000E+00
+k cmfd
+0.000000E+00
+0.000000E+00
+0.000000E+00
+0.000000E+00
+1.158665E+00
+1.177323E+00
+1.179341E+00
+1.169660E+00
+1.179827E+00
+1.180220E+00
+1.177699E+00
+1.179321E+00
+1.180728E+00
+1.180996E+00
+1.186680E+00
+1.181455E+00
+1.178393E+00
+1.176470E+00
+1.172593E+00
+1.167042E+00
+cmfd entropy
+0.000000E+00
+0.000000E+00
+0.000000E+00
+0.000000E+00
+3.229801E+00
+3.227279E+00
+3.224715E+00
+3.226705E+00
+3.226404E+00
+3.227329E+00
+3.227407E+00
+3.228601E+00
+3.230423E+00
+3.228093E+00
+3.226358E+00
+3.226226E+00
+3.227324E+00
+3.228594E+00
+3.229513E+00
+3.229820E+00
+cmfd balance
+0.000000E+00
+0.000000E+00
+0.000000E+00
+0.000000E+00
+6.244292E-03
+5.300885E-03
+5.198276E-03
+3.759109E-03
+3.994559E-03
+3.064086E-03
+2.749389E-03
+2.388483E-03
+2.304428E-03
+2.188627E-03
+2.269076E-03
+2.016776E-03
+1.622043E-03
+1.470984E-03
+1.367144E-03
+1.402525E-03
+cmfd dominance ratio
+0.000000E+00
+0.000000E+00
+0.000000E+00
+0.000000E+00
+5.509393E-01
+5.473805E-01
+5.434034E-01
+5.469278E-01
+5.460554E-01
+5.474063E-01
+5.486271E-01
+5.485095E-01
+5.490619E-01
+5.473524E-01
+5.466605E-01
+5.471154E-01
+5.480465E-01
+5.485550E-01
+5.500141E-01
+5.502154E-01
+cmfd openmc source comparison
+0.000000E+00
+0.000000E+00
+0.000000E+00
+0.000000E+00
+1.136193E-02
+1.007349E-02
+8.776995E-03
+6.729971E-03
+6.921055E-03
+5.867589E-03
+4.366276E-03
+4.051899E-03
+3.406846E-03
+4.043983E-03
+4.864903E-03
+3.711687E-03
+3.124413E-03
+2.675579E-03
+2.523549E-03
+2.724884E-03
+cmfd source
+4.534769E-02
+8.593793E-02
+1.059054E-01
+1.317554E-01
+1.358070E-01
+1.352224E-01
+1.290227E-01
+1.093443E-01
+7.943056E-02
+4.222668E-02