open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

I am trying to run a 3d case. The 2d case worked pretty well.(both with v6p0)

After scanning the boundaries, I have the error below and I don't know where it comes from, do you have any idea ?


CPSTMT : FORBIDDEN CASE FOR X AND Y:
X=M
ELEMENTS 41 AND 41
Y=DIFF
ELEMENTS 41 AND 41
Y IS SMALLER THAN X
X HAS A DIAGONAL OF TYPE Q
WITH A SIZE OF 804030
AND OFF-DIAGONAL TERMS OF TYPE S
WITH A SIZE OF ******
Y HAS A DIAGONAL OF TYPE ?
AND A MAXIMUM SIZE OF 1
AND OFF-DIAGONAL TERMS OF TYPE ?
AND A MAXIMUM SIZE OF 1

PLANTE: PROGRAM STOPPED AFTER AN ERROR
RETURNING EXIT CODE:

Thanks

Juliette

Hello Juliette,

It looks very much like an old bug in Telemac-3D that has now been corrected in version 6.2, so the best would be to upgrade ?

Regards,

Jean-Michel

Jean-Michel,

I have done what you said but I still have the same error with 3D v6p2. The only difference is the line "WITH A SIZE OF ******" became "WITH A SIZE OF 18411660".

Best Regards,

Juliette

Hello Juliette,

Can you attach your steering file, and Fortran file if you have changed something please?

Best regards,

Chi-Tuan

Hello,

Yes, we would need even the whole case to see what option leads to this.

Regards,

Jean-Michel

Jean-Michel, Chi-Tuan,

Here are my steering file and my fortran file. If you need other files (.geo, .cli, .lqd, ini), let me know and I will sen you that personnaly.

Thank,
Regards

Juliette

Hi all,

Just to say that with the option NON-HYDROSTATIC VERSION = YES, the error disappears. I will investigate to know why this error raises when using the hydrostatic hypothesis.

Anyway, the non-hydrostatic version is strongly recommended to all TELEMAC-3D users (please see some test cases like lock exchange).

Regards,

Chi-Tuan

Chi-Tuan,

I have the same error of Juliette. As you say, i put the option NON-HYDROSTATIC VERSION = YES. I want to know if you have a solution for the hydrostatic hypothesis?

Thanks
Florian

Hello,

Sorry I had found a solution that seemed to work for Juliette but have forgotten to post the solution.

In the enclosed archive, you have two POINT_TELEMAC3D subroutines: the original one and one corrected to enable the hydrostatic pressure hypothesis with SCHEME FOR DIFFUSION OF VELOCITIES = 0 called point_telemac3d_corr_bug_hydrostat_et_diff.f file.





The lines create the error are:
TYPDIA = '0'
IF (DIF(1) .OR. NONHYD) TYPDIA = 'Q'
TYPEXT = '0'
IF (DIF(1) .OR. NONHYD) TYPEXT = 'S'
CALL BIEF_ALLMAT (MDIFF, 'DIFF ',
* IELM3, IELM3, CFG, TYPDIA, TYPEXT,MESH3D)
and should be changed by:
TYPDIA = 'Q'
TYPEXT = 'S'
CALL BIEF_ALLMAT (MDIFF, 'DIFF ',
* IELM3, IELM3, CFG, TYPDIA, TYPEXT,MESH3D)
like in point_telemac3d_corr_bug_hydrostat_et_diff.f

Nevertheless, the non-hydrostatic version is strongly advised.

Hope this helps,

Chi-Tuan

Chi-Tuan,


Thanks for your post.
I work with a non-hydrostatic version but i wanted to understand why the simulation fails.

Thanks a lot
Regards
Florian