Color plot: better shading options

[wfranssen]

Make default 'shading=auto', to have a square at each point, and not between points as is now probably the case.

Perhaps also allow for Gouraud shading option.

See https://matplotlib.org/stable/gallery/images_contours_and_fields/pcolormesh_grids.html