Improve output of results and add checkpointing
The current output scheme doesn't stitch the MPI subdomains together during visualization (at least with VisIt). Ideally, we want something that can output 3D fields to a common format like VTK that works well in parallel. The file format should be directly openable in VisIt and/or ParaView.
Related or not, we should have the ability to save/load checkpoints.
ADIOS is one option here, but there might be others.