open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello all,

For the purpose of some development with telemac2D run in parallel I would like to access the NACHB array (described in the "Data structure specific to parallelism" section of the telemac guide for programming) in the time loop of the subroutine telemac2D.f. Can you tell me which module I should use so that I can access this array ? I hope I am clear enough.

Regards

Looks like NACHB is a sub_object of the BIEF_MESH type and its type is BIEF_OBJ, so I can access the values of NACHB in telemac2d.f subroutine like this : "MESH%NACHB%I(23)" for example.

What I don't get is :

1) why MESH%NACHB is a vector whereas in the telemac guide for programming it is said that NACHB is an array with 5 rows ;
2) why SIZE(MESH%NACHB%I(:)) can be greater than NPTIR on a given proc. For example in my case SIZE(MESH%NACHB%I(:))=230 and NPTIR = 23 ;
3) why some values of MESH%NACHB%I(:) can be greater than the number of proc I use to run TELEMAC in parallel. For example I run TELEMAC on 4 procs and MESH%NACHB%I(21) = 9.

Could you please help me ?

Sorry for pultiple posts, I must precise that what I'm looking for is a way to compute in telemac2d.f, for a given node at the boundary of a sub-domain, the number of sub-domains this node belongs to.

Regards

Hello

Yes you can use NACHB structure directly, however be careful in its use, because it is not obvious. You can look at subroutine voisin.f in bief library to have an idea how it is used. you can look at subroutine parini.f to see how it is built.
Depending on your algorithm, you can obtain this information without a use of nachb.

with my best regards

Riadh ATA

Thank you Riadh, yes the use of NACHB is far from being obvious ! However, if there exists another solution to know to how many sub-domain a node belongs to I am likely to learn it. Could you please telle me what is this solution ?

Regards

OK, I finally understood how the NACHB structure is composed, I wrote a little routine that compute the global number of node that belong to several sub-domain and it gives me the same results.

Regards