For my project, I need to identify the element which are in a level set fonction and recreate the boundary with the exterior element edge.
I didn't fine such subroutine an LevelSet module.
So I make my own subroutine, but I have some difficulties this the EdgeIndexes.
Code: Select all
DO i=1,nbrEltSolid
EltSolidB = ElementSolid(i)
Element =>Solver % Mesh % Elements(EltSolidB)
EdgeIndexes => Element % EdgeIndexes
DO j=1,Element % Type % NumberOfEdges ! boule sur j, sur le nombre d'edge
write(*,*) "element n°", EltSolidB, "edge n°", j
k = Element % EdgeIndexes(j)
END DO ! fin boule sur j, sur le nombre d'edge
END DO ! fin boule sur i, sur les élements solides
Thanks,
Antoine