open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

in a Telemac-2D model I have around 2000 boundaries.
In scalar mode the model works but in parallel mode the maximum number of boundaries is limited to 300 in the partel.f routine.

It is enough to change MAXFRO in the partel.f and to compile only the parallel directory?
Can I compile the parallel directory with Gfortran, replace those files at the cluster since I don't want to stress the cluster people and I don't have an Intel Fortran compiler.

For any hints I would be very glad!
Clemens

Hi Clemens,

You will need to change that value in declarations_telemac2d.f as well.

You need to use the same compiler and compile on the same machine as the rest of the code.

Hope it helps,
Yoann

Hi,
I have the same problem with a case that has about 600 solid boundaries. I changed everywhere in the sources files the MAXFRO variable to be: MAXFRO = 700 in :

partel.f
partel_para.f
partel_prelim.f
declarations_telemac2d.f

Then I recompiled the scalar and parallel mode using gfortran.

The case still cannot start if I run it in parallel mode. If I run it in scalar mode it runs fine. What could be the problem?

Thanks

Hello,

Changing partel.f and declarations_telemac2d.f should be enough, but if you have 600 solid boundaries, you perhaps have more than 700 boundaries as a whole ?

JMH

Thanks JMH for the reply. I have changed all the MAXFRO entries in the source code to MAXFRO = 3000. In my case I have 675 solid boundaries and 4 liquid boundaries. Even with the increase in MAXFRO I still have this error:

If I run the case in scalar mode TELEMAC2D does not complain at all. I checked the geo file and cli file to see if there was any incompatibility but so far everything looks good. Could this be a different problem?

Thanks,

Christian

Hello,

This is a problem in reading the geometry file. It could be due, not to the number of boundaries, but to another problem in the mesh, that is not seen in scalar mode. We would need to have the case to try it.

With best regards,

Jean-Michel Hervouet

Thanks Jean-Michel:

I tried to attach the case files but is exceeding the maximum size allowe in the forum is there any other way to send it to you?. It is a flooding case in which we have one inlet, one outlet and we want the water to flow free at the domain perimeter.

Christian

Hello,

for the print out statements in the if-conditions (v6p2r1) in partel.f line 2619 and in partel_para.f line 2010 MAXFRO is actually not used but an integer value:
----
IF(NILE.GT.300.OR.NFRSOL.GT.300.OR.NFRLIQ.GT.300) THEN
...
---
I hope this helps and maybe this could be changed for the next release?

Best regards,
Clemens

Hi Clemens,

You are right I didn't notice the line:
IF(NILE.GT.300.OR.NFRSOL.GT.300.OR.NFRLIQ.GT.300) THEN

I changed it to
IF(NILE.GT.3000.OR.NFRSOL.GT.3000.OR.NFRLIQ.GT.3000) THEN

in both partel.f line 2619 and in partel_para.f line 2010 and the case is running fine in parallel.

Thanks,

Christian

Well done !

We'll make sure that this is corrected in the next version.

JMH