Added related publications.
Updated notation for Darcy flow (piezometric height denoted by "p"). Improved proof of error estimate for flow.
Added some numerical results (w.i.p.).