diff --git a/meshmode/discretization/visualization.py b/meshmode/discretization/visualization.py
index 03c17d93b94ab48f9f5038afcb48f8840ba5ec12..359acab6a73bc80e64c643134079d0430681ae59 100644
--- a/meshmode/discretization/visualization.py
+++ b/meshmode/discretization/visualization.py
@@ -138,6 +138,7 @@ class Visualizer(object):
                 nodes.append(0*nodes[0])
 
             args = tuple(nodes) + (vis_connectivity.reshape(-1, 3),)
+            kwargs["scalars"] = field
 
             mlab.triangular_mesh(*args, **kwargs)