open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello everyone.
I'm trying to use Telemac3D in combination with AED2. The cas file I use is similar to the waq3d_aed2_flume example. With 02 open boundaries including 01 Q boundary and 01 SL boundary. However, despite trying many ways, the model still stops at this time step. I don't know the cause or how to fix it. I look forward to receiving help on solutions from everyone.
Thank you very much.
Jason Ha.

================================================================ ============================================
ITERATION 4366 TIME 0 D 6 H 3 MN 50.0000 S ( 21830.0000 S)
================================================================ ============================================
ADVERTISING STEP
DIFFUSION OF WN
GRACJG (BIEF): 22 ITERATIONS, ABSOLUTE PRECISION: 0.7015249E-08
PROPAGATION AND DIFFUSION WITH WAVE EQUATION
GMRES (BIEF): 16 ITERATIONS, ABSOLUTE PRECISION: 0.5252798E-08
DYNAMIC PRESSURE STAGE
GMRES (BIEF): 41 ITERATIONS, ABSOLUTE PRECISION: 0.9532586E-08
VELOCITY PROJECTION STEP
NAN AT I = 1 LEV = 2
H_AED2(LEV) = 2.9419250979767497E-008 FLUX(I,LEV) = NaN
TOP OF COLUMN @ 1 BOTTOM OF COLUMN @ 10
NAN VALUE
Note: The following floating-point exceptions are signaling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_OVERFLOW_FLAG IEEE_UNDERFLOW_FLAG IEEE_DENORMAL

================================================================ ==========================================

Hello Jason,

As told a few times in other posts, if a segmentation fault or a NaN (Not a Number) occurs, you should use a debug configuration with debug options to investigate.
See e.g. the S11.gfortran.debug configuration in the $HOMETEL/configs/systel.edf.cfg configuration file, in particular the flag fflags_debug_gfo for gfortran compiler.
fflags_debug_gfo: -g -Wall -fcheck=all -fbacktrace -fbounds-check -finit-integer=-1 -finit-real=nan -ffpe-trap=invalid,zero,overflow

It will show in which subroutine the issue/nan/segfault occurs, the first suspicious line and may help you to change something in your computation.

Anyway, can you also provide your steering files + potention Fortran files?

Do you have dry areas in your domain? That may explain this issue. I am not sure the coupling TELEMAC-3D with AED2 works correctly with dry areas.

What about the results of your TELEMAC-3D model without coupling with AED2, does it work?
If yes, does it work coupling with WAQTEL/THERMIC (WATER QUALITY PROCESS = 11)? Does it work with ATMOSPHERE-WATER EXCHANGE MODEL = 2 and/or 1? You may get errors with option 2 as it does work well with dry areas (if you have some).

Hope this helps,

Chi-Tuan

Hello, thank you very much for your answer.
Telemac3D works very well.
Telemac3D coupling with WAQTEL / THERMIC and ATMOSPHERE-WATER EXCHANGE MODEL = 1 works well when the parameter EVAPORATION = NO.
Telemac3D coupling with AED2 still does not work, even though the EVAPORATION = NO parameter. The AED2 cas file I use is similar to the cas file of the example \examples\waqtel\waq3d_aed2_flume.
Jason Ha.

Hello Jason,

You have not answered every question:
- can you give the listing of the execution with debug options, hard to help without seeing the first error,
- can you upload the TWO steering files (TELEMAC-3D + WAQTEL) + potential Fortran files, even though you think the WAQTEL steering file is similar to an existing example,
- "Do you have dry areas in your domain?"
- Does it work with ATMOSPHERE-WATER EXCHANGE MODEL = 2?

Chi-Tuan

Hello
I still have a problem related to Postel 3D temple. I followed the instructions in the example but did not cut the result successfully. I have attached point and cross section photos.
Looking forward to your guidance and answers.
Thank you.
Jason Ha.

Hello,

Can you upload the POSTEL-3D steering file + one sample of result (you can decrease the size to the last writing e.g.: use run_telfile.py alter -f -1 r3d.slf r3d_last.slf)?

In addition, can you upload the listing of the POSTEL-3D run please? Hard to help without knowing what you consider as an unsuccessful cut.

Chi-Tuan