Get almost all features working with MPI parallelization

Still to do:

+ Add an extra zone of padding to HDF5 files to account for shells if
  domain is non-periodic

+ Check in testing code from new m2tools Python module