gridder.interp returns raveled arrays instead of matrices

gridder.interp used to return 2D arrays with the interpolated data and x, y coordinates. This is inconsistent with the rest of fatiando, which uses 1D arrays. So I changed this and relevant code that uses the functionality. For most of it, the user is not affected (plotting functions work the same).