Current plotting

@cwg's current plotting code with some tests added and some cosmetic changes.

Still to do:

+ [ ] Change from Gaussian smearing to polynomial with finite support (e.g 1 - r**2 for r < 1 and 0 for r >= 1)

EDIT: due to numpy.where using a polynomial smearing is actually no more efficient