InterpPointDataToFld failure with repeated points
When there are points with the same coordinates in the .pts file, which is common on the element interfaces if the points in .pts are generated from slices of DG solutions, the interpolation gives wrong results.
Add a checking subroutine to delete/average the points with repeated coordinates will solve this problem.