diff --git a/meshmode/mesh/visualization.py b/meshmode/mesh/visualization.py
index 601086c1570e1b2a87f57b68361010774a78aeb7..52711a0a8b3654a0ea89291c47a7c14ab19d8ed0 100644
--- a/meshmode/mesh/visualization.py
+++ b/meshmode/mesh/visualization.py
@@ -41,6 +41,11 @@ def draw_2d_mesh(mesh, draw_vertex_numbers=True, draw_element_numbers=True,
         for iel, el in enumerate(grp.vertex_indices):
             elverts = mesh.vertices[:, el]
 
+            from meshmode.mesh import TensorProductElementGroup
+            if isinstance(grp, TensorProductElementGroup) and grp.dim == 2:
+                elverts = elverts[:,
+                        np.array([0, 1, 3, 2])]
+
             pathdata = [
                 (Path.MOVETO, (elverts[0, 0], elverts[1, 0])),
                 ]