diff --git a/docs/julia/heater_3d.jl b/docs/julia/heater_3d.jl
index 5c331e42..90da2a00 100644
--- a/docs/julia/heater_3d.jl
+++ b/docs/julia/heater_3d.jl
@@ -29,7 +29,7 @@ using Femwell.Maxwell.Electrostatic
 using Femwell.Thermal
 
 dir = @__DIR__
-run(`python $dir/heater_3d_mesh.py > /dev/null`)
+run(`python $dir/heater_3d_mesh.py`)
 
 
 model = GmshDiscreteModel("mesh.msh")
@@ -136,6 +136,8 @@ GridapPETSc.with(args = split(options)) do
         solver = PETScLinearSolver(),
     )
 
+    stophere
+
     #createpvd("poisson_transient_solution") do pvd
     #    for (uₕ, t) in uₕₜ
     #        pvd[t] = createvtk(