How to run:

mpiexec -np <number of processes> python TwoDAdvection_setup.py

With visualisation add -d and with HDF5 output add -f filename.hdf
