Update node coordinates in paraview

You should use XDMFFile for this:

from dolfin import *

mesh = UnitSquareMesh(10,10)
V = FunctionSpace(mesh, "CG", 2)
Q = FunctionSpace(mesh, "CG", 1)
u = Function(V)
p = Function(Q)
outfile = XDMFFile("output.xdmf")
t = 0
outfile.write(u, t)
outfile.write(p, t)

and then follow the instructions about opening the file and extracting blocks in my previous post.

1 Like