If I have some code which produces MeshFunctions whilst running in parallel, can I save that all to one file?
MF = MeshFunction("size_t", Mesh, 1)
XFile = XDMFFile("test.xdmf")
XFile << MF
HFile =HDF5File("h5test.h5","w")
HFile.write(MF, "meshfunction/example")