Thank you very much for sharing this piece of code, this is very helpful. I'm a bit confused with the way to set the source: - how should I add another dirac point source somewhere on the grid ? - how should I make the point source time dependent (e.g. a sinusoidal function) ?