open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi
I have some fortran code that forces the elevation on the open boundary using harmonics in formatted data file 2 with formatted data file 1 being a file that maps the global node number and the boundary number. The code works in scalar but falls over in parallel (get segmentation fault). I think the boundary numbers are not communicated properly between Sl and Bord in parallel.

I've tried following the advise in www.opentelemac.org/index.php/kunena/16-...ditions?limitstart=0 but to no avail.

Any comments would be greatly appreciated.

David

Hello,

In parallel you do not have all the points in a sub-domain, and their number is changed, with the property:

global number in original mesh = MESH%KNOLG%I(local number in subdomain)

Then in a subdomain you have :

local number in subdomain = MESH%NBOR%I(boundary number in sub-domain)

You also have :

original boundary number in sub-domain = BOUNDARY_COLOUR%I(local boundary number in sub-domain)


It is difficult to help you more since we do not see explanations on your algorithm, but as you never use any of these 3 arrays it is a hint that it cannot work in parallel.

With best regards,

Jean-Michel Hervouet