Use numdiff when available.
Merge branch 'master' of github.com:geodynamics/aspect into freesurface
Check that we are not using nonsense pressure normalization with free surface