[pad]

Hi all,

I am experiencing a strange problem for which I would like your input if possible.

I have a volute geometry from which I would like to extract the min / max radius of the inflow patch (it's either cylindrical or conical in shape). For this I use the following approach written as a python script.

1) Identify startFace and nFaces for the inflow patch in the boundary file
2) Go through faces file starting from ID startFace and ending at startFace+nFaces to extract point IDs
3) Take all point IDs from 2) and find the corresponding point coordinates from points file.

Now judging from the pointsOK.png picture where the extrated points are displayed this approach is usable. However, when I increase the number of layers in the snappyHexMesh I suddenly get the points extraction shown in pointsNotOK.png.

How can this be? The inflow patch shows up OK in paraview in both cases.

Any comments / ideas are highly appreciated.

Update: problem still persists but I found a workaround using foamToVTK to extract points on patch. Still a strange problem though :-)

Kind regards,
Per