open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello

I want to run a model in parallel mode, using the MASKOB routine :

MASKEL(I)=0.D0, where I the the index of the element taht I want to mask.


How can I transform the global index of my elements to the local index in each sub-mesh ?

For the point, I konw the MESH%KNOLG%I(I) instruction, but what is the good one for element ?

Regards

Hello,

You can use the following integer arrays:
MESH%KNOLG%I : gives the original global number of nodes
MESH%KNOGL%I : gives the local number of node as a function of their original number
This mens that MESH%KNOGL%I(MESH%KNOLG%I(J))=J

You can see the description of these arrays in the commented lines of routine paragl.f (bief)

With my best regards.

Riadh

Thanks....

MESH%KNOGL%I gives the local number of node as a function of their original number

But I need a function that gives thne local number of ELEMENT as a function of their original number, since i want to use the MASKEL .....

Hello,

There is no such information for elements, I'm afraid you will have to do a loop on all elements, to find the one that has the correct 3 global numbers of points, and there you will have to use KNOLG. Depending on where you do this, you can use the connectivity IKLE as IKLE(NELEM,3), or if not available as such the 3 global numbers of points of the element will be MESH%IKLE%I(IELEM),MESH%IKLE%I(IELEM+NELEM),MESH%IKLE%I(IELEM+2*NELEM), this for the local numbering, and then apply KNOLG to find the original numbering.

I hope this helps, it is rather cumbersome but actually the KNOLG for elements is never used except for your case.

R>Egards,

Jean-Michel Hervouet