Loading xdmf data back in

Hi there,

Having successfully imported time averaged functions for the velocity and pressure fields, I was wondering if there was a way to interpolate these to a mixed element function space. The purpose of doing so is to extract the Jacobian matrix via the process outlined in:

Any help would again be much appreciated,
Thanks,
Ben