open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi

We use a CORFON subroutine to affect the elevation of each points of our mesh, by reading a text file and with a loop like this :

Do I=1, NPOIN
read(11,*) a,b
ZF%R(a)=b
ENDDO

It works in no-parrallel mode.

But in parrallel mode, we use this loop :

DO I=1, 700657
read(11,*) a,b
ZF%R(MESH%KNOGL%I(a)),b
ENDDO

where we had to replace NPOIN by the value of the global number of nodes of our mesh, because NPOIN seems to bo the number of node of a subdomain in parrallel mode.

What's the name of the variable which give the global nimber of nodes of the mesh in parrallel mode ?

Thanks

Regards

Hi
Even you will have an access to the global number of nodes, your program couldn't work in parallel.
In each subdomain i, you had NPOIN(i).
So when you make a loop on your mesh in parallel, you work on your subdomain with a maximum of NPOIN(i) nodes.
In your program, your file had been created for the entire mesh. It's not comaptible.

for my opinion, it's easiest to manage the bathymetry before running your simulation (with Fudaa for example). Then you don't have to use corfon during the run.
Eventually, you make a first run of only one tiem step in scalar mode and you save only the bathymetry. After that you run your simulation in parallel with the previous result as geometry file.

regards