diff --git a/test/test_grudge.py b/test/test_grudge.py
index 7d7ffbcc2e07e379542ed4db652e92a313ab4af7..fbc1f50f66eae8e03fbd74fde592e11137652b77 100644
--- a/test/test_grudge.py
+++ b/test/test_grudge.py
@@ -673,16 +673,13 @@ def test_surface_divergence_theorem(ctx_factory, mesh_name, visualize=False):
         eoc_local.add_data_point(h_max, err_local)
 
         if visualize:
-            r = cl.array.to_device(queue, op_local)
-            r = cl.clmath.log10(cl.clmath.fabs(r) + 1.0e-16)
+            from grudge.shortcuts import make_visualizer
+            vis = make_visualizer(discr, vis_order=builder.order)
 
-            from meshmode.discretization.visualization import make_visualizer
-            vis = make_visualizer(queue, discr, vis_order=builder.order)
-
-            filename = "test_surface_divergence_theorem_error_{:04d}".format(i)
+            filename = f"surface_divergence_theorem_{mesh_name}_{i:04d}.vtu"
             vis.write_vtk_file(filename, [
-                ("r", r)
-                ], overwrite=True, legend=False)
+                ("r", actx.np.log10(op_local))
+                ], overwrite=True)
 
     # }}}