Dirac delta in boundary conditions for MethodOfLines.jl

Then maybe try the WENO scheme?