Add subdomain id from libMesh to the vtk files

print out subdomain ids in vtk format.
the patch resembles #145