diff --git a/examples/serialized_custom_source/show_flux.py b/examples/serialized_custom_source/show_flux.py
new file mode 100644
index 0000000000..6f54943018
--- /dev/null
+++ b/examples/serialized_custom_source/show_flux.py
@@ -0,0 +1,14 @@
+import matplotlib.pyplot as plt
+import openmc
+
+# Get the flux from the statepoint
+with openmc.StatePoint('statepoint.10.h5') as sp:
+    flux = sp.tallies[1].mean
+    flux.shape = (50, 50)
+
+# Plot the flux
+fig, ax = plt.subplots()
+ax.imshow(flux, origin='lower', extent=(-5.0, 5.0, -5.0, 5.0))
+ax.set_xlabel('x [cm]')
+ax.set_ylabel('y [cm]')
+plt.show()