Thanks for the hints on loading data from "csv" file with spline function.
I was able to generate the profile and test run using Navier Stokes equation, with turbulent boundary layer (Re=1E6). I have attached some images.


Going forward, I would want to check drag and lift coefficient with respect to available data from airfoil database. I used the tutorial on the flow around the cylinder, with the calculation of the lift coefficient.
Equation for the local drag force is set to : nx*p+miu*(-2*nx*ux-ny*(uy+vx))
Is there any reference documentation, apart from the tutorial, to explain the construction of this parametric equation from the integral equation ?
I set the equation for the local lift force to : ny*p+miu*(-2*ny*uy-nx*(ux+vy)) but want to check the correctness.
The combination with Matlab capacities makes CFD / FETools very powerfull.
Thanks for the support