Programming custom sample utility
 

Hello all,

I was trying to create a postprocessing utility to calculate boundary layer parameters. What I want to do is:

- iterate over each faces of a patch,
- create a line normal to each face,
- sample velocities over the line,
- Do some calculations to compute the parameters (just some integrals along the line)
- Save the values on the boundary face as a new field.

My problem is that I can't find a valid function to do the sampling of the velocity over the line.
Could it be possible to use the sample utility of OpenFOAM from inside another utility?
Does anyone has an idea on how to do this? :confused:

Thank You!
Raul