[jhint9]

Is there any way to extract the coordinates of surface cells? As in the 8 corners of each individual cell on the surface of an object once it is meshed. It's not imperative, but I have been tasked with trying. Thanks for any info!

[T.D.]

Hi,

Try this:

forAll(mesh.C(),i)
{
    const unallocLabelList& ne = mesh.cellPoints()[i];  
    forAll(ne, neI)
    {
        const unallocLabelList& ne_adj = mesh.pointCells()[ne[neI]]; 
        forAll(ne_adj, ne_adjI)
        {
            
        vector Coord = mesh.C()[ne_adj[ne_adjI]];

        // then Do what you want with the vector Coord
        }
    }
}

Reagrds,
T.D.

Originally Posted by jhint9

Is there any way to extract the coordinates of surface cells? As in the 8 corners of each individual cell on the surface of an object once it is meshed. It's not imperative, but I have been tasked with trying. Thanks for any info!

[jhint9]

Thanks for the info!
I saved the code in its own file name "coor"
I tried it and I'm getting a syntax error that reads:

./coor: line 1: syntax error near unexpected token `mesh.C'
./coor: line 1: `forAll(mesh.C(),i)'

I tried removing some parenthesis but then it didn't like the "i" argument.
Any ideas? Thanks again!

- Josh

[J-P]

Originally Posted by jhint9

Thanks for the info!
I saved the code in its own file name "coor"
I tried it and I'm getting a syntax error that reads:

./coor: line 1: syntax error near unexpected token `mesh.C'
./coor: line 1: `forAll(mesh.C(),i)'

I tried removing some parenthesis but then it didn't like the "i" argument.
Any ideas? Thanks again!

- Josh

To me it seems that you tried executing the code as a shell script, which is incorrect. What you should do is to include the given piece of code to your solver or utility, compile it, and then run it with the case or mesh you are working with. The OF user's guide will tell you how to get started with making your own version of a solver or an utility.

J-P

[jhint9]

That probably explains it. I'll spend some time looking into that. Thanks!

- Josh

[Hxchange]

Originally Posted by J-P

To me it seems that you tried executing the code as a shell script, which is incorrect. What you should do is to include the given piece of code to your solver or utility, compile it, and then run it with the case or mesh you are working with. The OF user's guide will tell you how to get started with making your own version of a solver or an utility.

J-P

Hey, I have the same question ,and I want to know which file should I add your code into? Sorry I dont konw much about OpenFOAM.
I used snappyMesh and chtMultiRegionSimpleFoam ,where should I put your code into?
I want to get the cooridates of all the cells ,8 points for each cells, instead of only the cell centre points.
Thank you!