[sunagra27]

Dear all,


I want to extract faceLabels from a patch, so I wrote the code as shown below.

But, I get the patch with incomplete faces when visualized, although it shows the exact no of faces present in the patch.


1. I saved the faceLabels in a ".txt" file.
2. I used this facelabel in topoSet as "LabelToFace" and converted to VTK using "foamToVTK" for visualization.

3. What is wrong in the code, so that it shows incomplete patch faces although it shows correct amount of faces and also I believe that the faceLabel alignment is not proper.

4. Code

code
       #{

    
        const fvPatch& mesh = this->patch();
        const fvBoundaryMesh& boundaryMesh = mesh.boundaryMesh();
        
    
        for (int patch = 0; patch < 1; patch++)
        {
        
            const word& patchName = mesh.boundaryMesh()[patch].name();
            Info << patchName << endl;
            const fvPatch& cPatch = mesh.boundaryMesh()[patch];
            std::ofstream file;
            file.open("faceID.txt");

            forAll (cPatch, faceI)
            {
                label faceId = cPatch.start() + faceI;
                //Info << faceId << endl;
                file<< faceId << endl;
                file<< nl << endl;
            }
            file.close();
        }


        #};

5. Images

a. Patch I needed: Below is the patch that I needed. This image is taken after snappyHexMesh.

https://drive.google.com/file/d/1KUq...ew?usp=sharing

b. Patch I am obtaining: Below is the patch based on code. We can see it has incomplete faces.


https://drive.google.com/file/d/1Rdi...ew?usp=sharing

Any leads will be appreciated.


Thanks and Regards,
Sunag R A.