Pressure is now nondimensionalized using the thermal diffusion time
instead of the rotation time.  This changes P/Ek in the equation set to
P.  and dPdr/Ek to dPdr.  This appears to stablize the band solves for
finite-differences.
Merge branch 'master' of https://github.com/geodynamics/rayleigh