Add __name__ == '__main__' magic to test_upwind_flux
Add DG advection test (parallel too)
Run 2D and 3D steady advection tests in parallel

Now that we can do interior facet integrals in parallel, test steady
advection in parallel as well as serial.