Add DG local operator to the SimulationTransport class
Add support for discontinuous Galerkin methods in the transport model.
The main idea is to adapt the current Transport model to accept two local operators in a similar way as it is done now for the Richards model. The first local operator is naturally the finite volume local operator which is currently working, and the second operator would be a newly created local operator for the diffusion transport equations with for the discontinuous Galerkin method.
How to test the implementation?
Test for the finite volume is extensible for the DG cases, and they should work out of the box with higher precision.
See #72 (closed)