Split off the plot functionality for the interface surface

[mattngc]

Currently there is one plot() utility which generates multiple utilities. It would be better to split this off into equlibrium and heatflux.