[phani45]

Hi foamers,

I tried to print out the face centroid values on a plane with fixed X value and lying on Y-axis.
I used the code snippet:

{


label patchID=mesh.boundaryMesh().findPatchID("left_wall ");

const labelList& own = mesh.faceOwner();


forAll(mesh.boundaryMesh()[patchID].faceCells(),cellf)
{

Info<<mesh.Cf()[own[cellf]]<<endl;
}

}

I'm used 10 cells in Y-direction . The length in Y-Direction is 1.0 starting from origin.

I'm getting following output:
(0.05 0.1 0.05)
(0.05 0.1 0.05)
(0.025 0.2 0.05)
(0.025 0.2 0.05)
(0.1 0.1 0.05)
(0.1 0.1 0.05)
(0.075 0.2 0.05)
(0.075 0.2 0.05)
(0.15 0.1 0.05)
(0.15 0.1 0.05)

I don't see the X value being same (which it is supposed to be 0.0) neither the Y -value varying from 0 to 1 as 0,0.1,0.2 etc....

Can anyone of you tell me where I went wrong ?

Regards,
Phanindra

[hjasak]

Think about it. Cf is a surfaceVector Field, meaning that it is defined on FACES. You are asking the patch for CELLS that are next to the patch faces:

const labelList& own = mesh.faceOwner();

and then addressing owner (which gives you a CELL index) into a list which is defined in FACES.

Enough said.

[phani45]

I was trying with many alternatives but couldn't find one. Can you tell me how to solve this problem please.

[theDoomGuy]

It should be like this for someone who has run into the same problem:

const surfaceVectorField Cfs = mesh.cf();

forAll(Cfs, faceI)
{
    vector Cf = Cfs[faceI];
    Info << "centroid of face " << faceI << " = " << Cf << nl << endl;
}

This works for me. I am also a beginner so any type of criticism (constructive) is welcome. The above code provides the centroid of internal faces but, I assume you could modify the code for extracting information on the baoundary cells.