The source term so far did not take into account contributions from the
pore pressure gradient. This is fixed in this commit and makes the
discrete solution agree better with the exact solution, which is why the
reference solution is modified here.