[SD@TUB]

Hello FOAMers,

as long as sample utility and ParaFoam won't provide circular planes (cylindrical region) for post-processing issues, i want to put collected cells by routine

forAll(mesh.cells(),cellI)
{
// Loop over all the cells in the mesh
++numCells;

// Check containment of each cell center
const [*definition of sub domain*]

if [*cells within subdomain*]
{
++numInside;
totalVolume += mesh.V()[cellI]; //e.g. volume of collected cells
}
}

in a subset to grab them afterwards and reading data of interest. What would be a kind way of performing such iteration. I noticed that cellSet works via dictionary. Is it possible to adopt some functionality?


/Stefan

[SD@TUB]

One reply to myself:
The utility cellSet provide curved region via cellSetDict, e.g. cylindrical volumes which could be used to post-process selected cells afterwards.

Nevertheless, if anybody has a suggestion to perform such process within code, this would be very helpful! I tried adopting the appropriate code from cellSet.C without success.

~stefan