Generalise error estimates from velocity module
In !35 (merged) we introduced a method for an error estimate introduced in a paper by B. Hess. It would be nice if we could generalise this and move the algorithm from the transport module into utils.py.