Output VTK file from domain in mesh example
Add VTK output of the mesh
Write one VTK file for requested subdomains
Give vtk output binary and compression options
Add subdomain id and mpi rank to vtk output