WIP. Added Dirichlet BCs and new plotting. Also normalizing PSF to 1. Introducing boundary rate Gamma.