Including plotting of general scalar and vector fields on mesh

It would be nice to include a method which allows to plot some generic scalar or vector field on the mesh. Assume, somebody has set some anistropy and wants to look at it.

my_anis = ... # some vector field defined on mesh nodes
sample.plot(my_anis, *args)
sample.show()