IceGrid::kBelowHeight() should call MPI_Abort().

This should be ported to 0.6.
Better error messages. Port to 0.6.

Now PISM always checks if re-gridded ice thickness exceeds the height of
the computational domain and stops if it does.