diff --git a/openpnm/network/CubicDual.py b/openpnm/network/CubicDual.py
index <HASH>..<HASH> 100644
--- a/openpnm/network/CubicDual.py
+++ b/openpnm/network/CubicDual.py
@@ -89,7 +89,7 @@ class CubicDual(GenericNetwork):
         single_dim = shape == 1
         shape[single_dim] = 2
         dual = Cubic(shape=shape-1, spacing=1)
-        faces = [['front', 'back'], ['left', 'right'], ['top', 'bottom']]
+        faces = [['left', 'right'], ['front', 'back'], ['top', 'bottom']]
         faces = [faces[i] for i in np.where(~single_dim)[0]]
         faces = np.array(faces).flatten().tolist()
         dual.add_boundary_pores(faces)
changed the order of the boundary faces in the faces array to the new state