Moved all plotting functions into ionvis.py.
Initial attempt at using JSON files.

The simulation runs, but it appears that something goes wrong immediately.
This could be due to poor choices of parameters, since the current
default.json is using new stuff for Ca+ ions.