open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I apologise in advance in case this question was already addressed somewhere else in the forum (I couldn't find it).

Working with user fortran files for TELEMAC2D (but also TELEMAC3D) parallel calculations, I would need to know the global (2D) node number from the local node index and processor number. I know the BOUNDARY_COLOUR variable helps for boundary nodes, but what about interior nodes?

More concretely, say the global mesh contains 10000 nodes and TELEMAC is run on 12 procs, considering a local loop index I=1,NPOIN(local,<10000) running on the subdomain related to processor J, I would need to know the global node index K=f(I,J) described by index I in subdomain J, with K ranging from 1 to NPOIN(global,=10000) and f is the function I am looking for.

Is that possible?

Thank you in advance for your help.
JBS.

Hello,

Every telemac module has a variable named MESH (or something like that).

This variable is a user-defined type with other variables inside it. The ones you are looking for are KNOGL (Global to Local) and KNOLG (Local to Global).

A Global node number is the number of your entire mesh, as you created it. A Local node number is the number of the node in each parallel mesh.

So, in order to convert a global mesh node to a local one (and vice versa) you use these two variables.

I don't quite remember which variable converts global to local, but the usage is the same, so you only have to test which one is the one you need.

Best Regards,

Phelype

Thank you Phelype. It seems to work!

In the user subroutine (meteo.f here), I introduced the variable MESH (in 'USE DECLARATIONS_TELEMAC2D,ONLY : ...') and used :


The results look good so I guess it did the job!

Best regards.
JBS.