diff --git a/fenicsR13/fenicsR13.py b/fenicsR13/fenicsR13.py
index 3cb37b67f6c46ed2eb58be1dc0d948c00a1ca27c..c56ac8370282f75526f9ae9406f6399ac72437f9 100755
--- a/fenicsR13/fenicsR13.py
+++ b/fenicsR13/fenicsR13.py
@@ -132,7 +132,10 @@ def main():
                 if p == len(mesh_names) - 1:  # after last mesh
                     postp = Postprocessor(data, params["output_folder"])
                     postp.write_errors()
-                    postp.plot_errors(show_plot)
+                    if len(mesh_names) != 1:
+                        postp.plot_errors(show_plot)
+                    else:
+                        print("Skip convergence plot (only 1 mesh)")
 
             solver = None
             gc.collect()