diff --git a/examples/python/pincell_multigroup/build-xml.py b/examples/python/pincell_multigroup/build-xml.py
index ee74233c1..cc57b61e2 100644
--- a/examples/python/pincell_multigroup/build-xml.py
+++ b/examples/python/pincell_multigroup/build-xml.py
@@ -1,3 +1,5 @@
+import numpy as np
+
 import openmc
 import openmc.mgxs
 
@@ -26,7 +28,7 @@ uo2_xsdata.set_total(
      0.5644058])
 uo2_xsdata.set_absorption([8.0248E-03, 3.7174E-03, 2.6769E-02, 9.6236E-02,
                            3.0020E-02, 1.1126E-01, 2.8278E-01])
-uo2_xsdata.set_scatter_matrix(
+scatter_matrix = np.array(
     [[[0.1275370, 0.0423780, 0.0000094, 0.0000000, 0.0000000, 0.0000000, 0.0000000],
       [0.0000000, 0.3244560, 0.0016314, 0.0000000, 0.0000000, 0.0000000, 0.0000000],
       [0.0000000, 0.0000000, 0.4509400, 0.0026792, 0.0000000, 0.0000000, 0.0000000],
@@ -34,6 +36,8 @@ uo2_xsdata.set_scatter_matrix(
       [0.0000000, 0.0000000, 0.0000000, 0.0001253, 0.2714010, 0.0102550, 0.0000000],
       [0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0012968, 0.2658020, 0.0168090],
       [0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0085458, 0.2730800]]])
+scatter_matrix = np.swapaxes(np.swapaxes(scatter_matrix, 0, 1), 1, 2)
+uo2_xsdata.set_scatter_matrix(scatter_matrix)
 uo2_xsdata.set_fission([7.21206E-03, 8.19301E-04, 6.45320E-03,
                         1.85648E-02, 1.78084E-02, 8.30348E-02,
                         2.16004E-01])
@@ -50,7 +54,7 @@ h2o_xsdata.set_total([0.15920605, 0.412969593, 0.59030986, 0.58435,
 h2o_xsdata.set_absorption([6.0105E-04, 1.5793E-05, 3.3716E-04,
                            1.9406E-03, 5.7416E-03, 1.5001E-02,
                            3.7239E-02])
-h2o_xsdata.set_scatter_matrix(
+scatter_matrix = np.array(
     [[[0.0444777, 0.1134000, 0.0007235, 0.0000037, 0.0000001, 0.0000000, 0.0000000],
       [0.0000000, 0.2823340, 0.1299400, 0.0006234, 0.0000480, 0.0000074, 0.0000010],
       [0.0000000, 0.0000000, 0.3452560, 0.2245700, 0.0169990, 0.0026443, 0.0005034],
@@ -58,6 +62,8 @@ h2o_xsdata.set_scatter_matrix(
       [0.0000000, 0.0000000, 0.0000000, 0.0000714, 0.1391380, 0.5118200, 0.0612290],
       [0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0022157, 0.6999130, 0.5373200],
       [0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.1324400, 2.4807000]]])
+scatter_matrix = np.swapaxes(np.swapaxes(scatter_matrix, 0, 1), 1, 2)
+h2o_xsdata.set_scatter_matrix(scatter_matrix)
 
 mg_cross_sections_file = openmc.MGXSLibrary(groups)
 mg_cross_sections_file.add_xsdatas([uo2_xsdata, h2o_xsdata])
diff --git a/examples/xml/pincell_multigroup/mgxs.h5 b/examples/xml/pincell_multigroup/mgxs.h5
index 1d5561b00..ae3b17e90 100644
Binary files a/examples/xml/pincell_multigroup/mgxs.h5 and b/examples/xml/pincell_multigroup/mgxs.h5 differ