diff --git a/scripts/run_experiments.jl b/scripts/run_experiments.jl
index ce1cc2c32f665048f3085de0787c1b9331da851a..62ecfe455752ed4fa329443ac9c4ac81d89c3775 100644
--- a/scripts/run_experiments.jl
+++ b/scripts/run_experiments.jl
@@ -619,7 +619,8 @@ function estimate_res!(st::L1L2TVState)
     cellf(hcell; g, u, nablau, p1, p2, tdata) =
          st.alpha2 *
             st.T(adjoint, tdata, st.alpha2 * st.T(tdata, u) - g) +
-            st.T(adjoint, tdata, p1)
+            st.T(adjoint, tdata, p1) +
+            (isSnabla ? zero(p1) : st.beta * u)
     facetf(hfacet, n; g, u, nablau, p1, p2, tdata) =
         #display(n' * p2)
         # norm2 is calculated later after both facet contributions are