Poisson equation with two subdomains and Newton cooling type discontinuity at the interface

For anyone interested, see solution at Imposing a discontinuity at interface using DG method - #19 by smesc