diff --git a/examples/amoc.py b/examples/amoc.py
index cdba29b..0cb9e22 100644
--- a/examples/amoc.py
+++ b/examples/amoc.py
@@ -53,7 +53,6 @@ def main():
 
     nx = 60
     ny = 30
-    nz = 1
 
     # Define the problem
     parameters = {'Problem Type': 'AMOC',
@@ -66,7 +65,7 @@ def main():
                   # Give back extra output (this is also more expensive)
                   'Verbose': False}
 
-    interface = Interface(parameters, nx, ny, nz)
+    interface = Interface(parameters, nx, ny)
     continuation = Continuation(interface, newton_tolerance=1e-6)
 
     # First increase the temperature forcing to the desired value
diff --git a/examples/ldc2.py b/examples/ldc2.py
index 76a3151..2d63356 100644
--- a/examples/ldc2.py
+++ b/examples/ldc2.py
@@ -14,7 +14,6 @@ def main():
     ''' An example of performing a continuation for a 2D lid-driven cavity and computing eigenvalues along the way'''
     nx = 16
     ny = nx
-    nz = 1
 
     # Define the problem
     parameters = {'Problem Type': 'Lid-driven Cavity',
@@ -22,7 +21,7 @@ def main():
                   'Reynolds Number': 1,
                   'Lid Velocity': 0}
 
-    interface = Interface(parameters, nx, ny, nz)
+    interface = Interface(parameters, nx, ny)
     continuation = Continuation(interface)
 
     # Compute an initial guess
diff --git a/examples/ldc3.py b/examples/ldc3.py
index 8cc9347..fed22ac 100644
--- a/examples/ldc3.py
+++ b/examples/ldc3.py
@@ -24,7 +24,6 @@ def main():
     ''' An example of performing a "poor man's continuation" for a 2D lid-driven cavity using time integration'''
     nx = 16
     ny = nx
-    nz = 1
 
     # Define the problem
     parameters = {'Problem Type': 'Lid-driven Cavity',
@@ -37,12 +36,12 @@ def main():
                   'Verbose': False,
                   'Theta': 1}
 
-    interface = Interface(parameters, nx, ny, nz)
+    interface = Interface(parameters, nx, ny)
 
     # Store data for computing the bifurcation diagram using postprocessing
     data = Data()
 
-    n = interface.discretization.dof * nx * ny * nz
+    n = interface.discretization.dof * nx * ny
     x = numpy.random.random(n)
 
     mu_list = []
diff --git a/examples/qg.py b/examples/qg.py
index 3e21639..006a420 100644
--- a/examples/qg.py
+++ b/examples/qg.py
@@ -28,7 +28,6 @@ def main():
 
     nx = 32
     ny = nx
-    nz = 1
 
     # Define the problem
     parameters = {'Problem Type': 'Double Gyre',
@@ -39,7 +38,7 @@ def main():
                   # Give back extra output (this is also more expensive)
                   'Verbose': False}
 
-    interface = Interface(parameters, nx, ny, nz)
+    interface = Interface(parameters, nx, ny)
     continuation = Continuation(interface)
 
     # First activate the wind stress