incorrect diffusion flux discretization
The current implementation of the diffusion fluxes is based on slope reconstructions. the resulting scheme depends a lot on the limiter that is used in the reconstruction which is likely a numerical artifact. We should change to a standard discretization of the diffusion terms.