Use interpolation rather than projection for reconstruction.
Try to make the reconstruction function compatible with MPI : use of interpolate_nonmatching_mesh_any from fenicstools