Change tolerance and plot errors.

* Use a tolerance of 1e-8 for the solves.
* Use BiCGStab for the velocity solve.
* Create a directory for each convergence test.
* Create nested directories if required.
* Generate a plot of the error at each grid point.