open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

I got the following error message when I have
- OPTION FOR LIQUID BOUNDARIES = 2;1
- PARALLEL PROCESSORS = 8

The model works fine when I use the default value instead of Thompson condition
- OPTION FOR LIQUID BOUNDARIES = 1;1
- PARALLEL PROCESSORS = 8
or when I use 6 processors.
- OPTION FOR LIQUID BOUNDARIES = 2;1
- PARALLEL PROCESSORS = 6




Anyone encountered similar issue? My guess is that this has to do with sub-domain without boundary node. Is there an easy fix to this problem?

Regards,

Edwin

Hello,

I think I see the point. It happens when a subdomain has no boundary points and it is due to the first 3 lines of thomps.f, where we use IT1(1:I) and IT1(I+1:2*I). The dimension of IT1 is declared as 2*NPTFR and it causes the problem, and in fact the real dimension of IT1 is NPOIN. So a solution is to take thomps.f in your Fortran file and give the real dimention of IT1 in the declaration : IT1(NPOIN) instead of IT1(2*NPTFR). We have already corrected a number of such cases and it is strange that this one did not pop up before.

We'll do this also in the next release.

With best regards,

Jean-Michel Hervouet

Thanks. That modification fixed the problem.

Cheers,

Edwin