Tighten tolerance for better fit with sim_simple.cpp.
Use harmonic averaging of mobility for pressure eqns.