@mary_h I believe that when exporting data to vtk files, dolfin takes the nodal dofs of the output if the function belongs to a CG space.
To get a function for visualization with element-wise P3 polynomials you should try the write_checkpoint method of the XDMFFile class (see the docken answer in this post for instance).