I have defined a mesh in the same domain (in FreeFem) and interpolated the data from one mesh to the other. I found that the points near the surface of the square cylinder are not properly calculated. I have read in the documentation (page 196) that to speed-up the process the interior of the objects are also triangulated. That can be the reason of this variation?