open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi again,

I have just transformed my geometry file and binary boundary file to spherical coordinates but I keep getting exactly the same error at line 286 of parres.f.

After checking the allocation of variables I think the error might be the following one:

In line 229 of parres.f IKLES is allocated as
ALLOCATE (IKLES(NELEM_GEO*NDP_GEO),STAT=IERR)

NDP_GEO=3 for my particular application (Triangles in the horizontal)

Then, in line 295 IKLE3D is allocated as

ALLOCATE (IKLE3D(NELEM_GEO*(NPLAN_INP-1)*NDP_INP),STAT=IERR)

with NDP_INP=4 when reading the binary boundary file (this is a 2D vertical file, the elements are kind of rectangles).

When IKE3D is filled, it reads from IKLES and it gets an error due to the different dimensions.

I wonder if the solution would just be replacing in line 229 of parres.f this
ALLOCATE (IKLES(NELEM_GEO*NDP_GEO),STAT=IERR)
by this
ALLOCATE (IKLES(NELEM_GEO*NDP_INP),STAT=IERR)

I am assuming that NDP_GEO==NDP_INP when reading the geometry file, but I could not check it.

Any ideas as welcomed.

Thanks a lot,
Luz

Hi Luz,

Ok thank you for the test in spherical coordinates. Now we suspect that T3DBND should not go through PARRES but through PARTEL itself. Maybe it's a mistake in the python scripts. Could you try to run partel.py on the file T3DBND in your temporary folder run?

Thanks,
Agnès

Hi Agnès,

Thank you for your reply. I have tried to run partel.py on T3DBND, and although I am not sure about the correct syntax, I have tried several possibilities and I always get the same error.

This is an example:

/gpfs/software/opentelemac/v7p2r2/scripts/python27/partel.py -c intpar --ncsize 24 --file T3DBND --bnd-file T3DCLI --input-file PARTEL.PAR

This is the error I get:

File "/gpfs/software/opentelemac/v7p2r2/scripts/python27/partel.py", line 154, in <module>
geoFile,geoFileFmt,partitionningMethod)
TypeError: runPARTEL() takes exactly 12 arguments (11 given)


I have been checking and, indeed, in partel.py function runPARTEL is called with 11 arguments, whereas in runcode.py, where runPARTEL is defined, it is indicated to use 12 arguments.

Could you please give me a hint on how to correct this and which would be the proper syntax for running partel.py??

Thanks a lot,
Luz

Hello,

For my part, I directly ran the partel executable in the temporary directory, after removing the existing T3DGEO* files. If I well understood, homere_partel first checks the existence of these partitioned T3DGEO* files (creating at the first call of partel during the partitioning of T3DGEO file), if they already exist, parres is called, if no it's partel that is called.
And this test still leads to the following error message with a segmentation fault.


+

---

+
PARTEL/PARRES: TELEMAC METISOLOGIC PARTITIONER

REBEKKA KOPMANN & JACEK A. JANKOWSKI (BAW)
JEAN-MICHEL HERVOUET (LNHE)
CHRISTOPHE DENIS (SINETICS)
YOANN AUDOUIN (LNHE)
PARTEL (C) COPYRIGHT 2000-2002
BUNDESANSTALT FUER WASSERBAU, KARLSRUHE

METIS 5.0.2 (C) COPYRIGHT 2012
REGENTS OF THE UNIVERSITY OF MINNESOTA

BIEF V7P2R0 (C) COPYRIGHT 2012 EDF
+

---

+


MAXIMUM NUMBER OF PARTITIONS: 100000

+

---

+

--INPUT FILE NAME <INPUT_NAME>:
T3DBND
INPUT: T3DBND
--INPUT FILE FORMAT <INPFORMAT> [MED,SERAFIN,SERAFIND]:
SERAFIN
INPUT: SERAFIN
--BOUNDARY CONDITIONS FILE NAME:
T3DCLI
INPUT: T3DCLI
--NUMBER OF PARTITIONS <NPARTS> [2 -100000]:
2
INPUT: 2
PARTITIONING METHOD <PMETHOD> [1 (METIS) OR 2 (SCOTCH)]:
1
--INPUT: 1
--CONTROL SECTIONS FILE NAME (OR RETURN) :

NO SECTIONS
--CONTROL ZONES FILE NAME (OR RETURN) :

NO ZONES
--WEIR FILE NAME (OR RETURN) :

NO WEIRS
--GEOMETRY FILE NAME <INPUT_NAME>:
T3DGEO
INPUT: T3DGEO
--GEOMETRY FILE FORMAT <GEOFORMAT> [MED,SERAFIN,SERAFIND]:
SERAFIN
INPUT: SERAFIN
+---- PARTEL: BEGINNING

---

+
NO BND ELEMENTS IN A 3D MESH


READ_MESH_INFO: TITLE=
NUMBER OF ELEMENTS: 549
NUMBER OF POINTS: 630

FORMAT NOT INDICATED IN TITLE


3D MESH DETECTED
NDP NODES PER ELEMENT: 4
NPLAN NUMBER OF MESH LEVELS: 10
NPOIN2 NUMBER OF 2D MESH NODES: 63
NPOIN NUMBER OF 3D MESH NODES: 630
NELEM2 NUMBER OF 2D MESH ELEMENTS: 61
NELEM NUMBER OF 3D MESH ELEMENTS: 549

THE INPUT FILE ASSUMED TO BE 3D
THERE ARE 2 TIME-DEPENDENT RECORDINGS
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
partel 00000000004E9869 Unknown Unknown Unknown
partel 00000000004E81E0 Unknown Unknown Unknown
partel 00000000004A4792 Unknown Unknown Unknown
partel 0000000000461263 Unknown Unknown Unknown
partel 0000000000465B1B Unknown Unknown Unknown
libpthread.so.0 0000003D57E0F7E0 Unknown Unknown Unknown
partel 00000000004452A2 Unknown Unknown Unknown
partel 0000000000409EA9 Unknown Unknown Unknown
partel 00000000004417C4 Unknown Unknown Unknown
partel 00000000004041F6 Unknown Unknown Unknown
libc.so.6 0000003D5761ED1D Unknown Unknown Unknown
partel 00000000004040E9 Unknown Unknown Unknown

Laurent

Hello,

Could you send me all the files of your case please? It will be easier to debug this.

Cheers,
Agnès

Hello Agnès,

All the files have already been posted by Luz in the first post of this topic, and a subset of the binary boundary file in the third post.
Let me know if I can do something to help.

Regards,
Laurent

Hello,

Yes, as Laurent says all the files are posted at the beginning of this forum thread.

Thank you for having a deeper look into this, it would be great to have this working soon!!

Cheers,
Luz

Hi,

I'm missing the tide files though, could you send them to me?

Cheers,
Agnès

Hi again,

These files are quite big to be posted here (47M and 27M, respectively). I could send them to you via Dropbox or something similar, but you can download them from

volkov.oce.orst.edu/tides/ES.html

If you click on the left (Downloads) in European Shelf 1/30, these file will be downloaded.

Cheers,
Luz

Or turn off tidal effect in the steering file, as it doesn't seem to be the problem !
Laurent