Unfortunately I can’t run your example on my current machine. But it looks like in your formulation for a there’s a restriction missing. It might be the variable h here,
+ gamma/h * inner(jump(u_trial, n), jump(v, n)) * dS
You could try something like h("+").
See also 1d dG advection/diffusion probem