diff --git a/docs/julia/heater_3d.jl b/docs/julia/heater_3d.jl
index 4686296f..b9df4a20 100644
--- a/docs/julia/heater_3d.jl
+++ b/docs/julia/heater_3d.jl
@@ -125,6 +125,8 @@ GridapPETSc.with(args = split(options)) do
         Dict(get_tag_from_name(labels, u) => v for (u, v) in thermal_diffisitivities)
     ϵ_diffisitivities(tag) = thermal_diffisitivities[tag]
 
+    dontplottransient
+
     uₕₜ = calculate_temperature_transient(
         ϵ_conductivities ∘ τ,
         ϵ_diffisitivities ∘ τ,
@@ -147,7 +149,7 @@ GridapPETSc.with(args = split(options)) do
     #end
     sums = [(t, ∑(∫(u)dΩ_w) / ∑(∫(1)dΩ_w)) for (u, t) in uₕₜ]
 
-    dontplottransient
+    #
 
     figure = Figure()
     ax = Axis(figure[1, 1], ylabel = "Temperature / K", xlabel = "time / ms")
diff --git a/docs/julia/heater_3d_mesh.py b/docs/julia/heater_3d_mesh.py
index c5413fb8..95068f17 100644
--- a/docs/julia/heater_3d_mesh.py
+++ b/docs/julia/heater_3d_mesh.py
@@ -37,8 +37,8 @@
 
 resolutions = {
     "core": {"resolution": 0.3},
-    "via2": {"resolution": 0.2},
-    "via1": {"resolution": 0.2},
+    "via2": {"resolution": 0.3},
+    "via1": {"resolution": 0.3},
     "heater": {"resolution": 0.5},
 }
 geometry = get_mesh(