diff --git a/applications/stokes_dg/stokes_dg.ufl b/applications/stokes_dg/stokes_dg.ufl
index 1347c0d05037206f03f5ac89dfc1803a72d1033b..8c2ac036b6bb8961a312037b88866e5a1aff4ccb 100644
--- a/applications/stokes_dg/stokes_dg.ufl
+++ b/applications/stokes_dg/stokes_dg.ufl
@@ -16,20 +16,21 @@ ds = ds(subdomain_id=1, subdomain_data=bctype)
 n = FacetNormal(cell)('+')
 eps = -1.0
 sigma = 1.0
+h_e = Min(CellVolume(cell)('+'), CellVolume(cell)('-')) / FacetArea(cell)
 
 r = inner(grad(u), grad(v))*dx \
   - p*div(v)*dx \
   - q*div(u)*dx \
   + inner(avg(grad(u))*n, jump(v))*dS \
-  + sigma * inner(jump(u), jump(v))*dS \
+  + sigma / h_e * inner(jump(u), jump(v))*dS \
   - eps * inner(avg(grad(v))*n, jump(u))*dS \
   - avg(p)*inner(jump(v), n)*dS \
   - avg(q)*inner(jump(u), n)*dS \
   - inner(grad(u)*n, v)*ds \
-  + p*inner(v, n)*ds \
-  + q*inner(u, n)*ds \
+  + sigma / h_e * inner(u-g_v, v)*ds \
   + eps * inner(grad(v)*n, u-g_v)*ds \
-  + sigma * inner(u-g_v, v)*ds \
-  - q*inner(g_v, n)*ds
+  + p*inner(v, n)*ds \
+  + q*inner(u-g_v, n)*ds
 
 forms = [r]
+exact_solution = g_v, 8*(1.-x[0])
\ No newline at end of file