open TELEMAC-MASCARET The mathematically superior suite of solvers

hello everyone

I’m modelling heat exchange with atmosphere in a stratified lake with closed boundary and constant level planes (defined in “CONDIM” and “DOUBLE PRECISION FUNCTION TRA_PROF_Z”). I’m using Tel3Dv7p2r0 coupled with “WAQTEL” like in the “heat exchange” example.

For the BORD3D subroutine in my fortran file I used the one of the example “heat exchange”, but when starting the computation I get “GRACJG (BIEF) : EXCEEDING MAXIMUM ITERATIONS“ errors for every time step even with very low time steps (0.1 seconds).
Without the part “HEAT EXCHANGE WITH ATMOSPHERE” which calls “CALCS3D_THERMICS” the calculation runs (with a time step of 1 second) without error, and the influence of the weather seems to be taken into account.

So I have following questions:
1) What calculation is made if telemac3d is coupled with WAQTEL and “CALCS3D_THERMICS” is not called in the fortran file?
2) Why is the number of iteration exceeded when calling it? Is something wrong in my “calcs3d_thermic.f file? (I’ve never changed anything in it) Or is there an incompatibility in my fortran file?

I’ve another question concerning the literature:
In the Telemac manual part concerning the thermic module the following publication K. El-Kadi-Abderrezzak, 2012)„K. El-Kadi-Abderrezzak, L. M. (2012). Module de qualité d'eau TRACER - note de principe. EDF R&D“ is quiet often quoted but I can’t find it anywhere. I would be interested in reading it (even if it is in French), do someone know where I can find it?

Please find attached my fortran file, the CALCS3D_THERMICS.f file and my two .cas files.

Best regards

Peter

to rectify:

the guilty part is not the part calling “CALCS3D_THERMICS”
but this two lines:
ATABOS%ADR(IND_T)%P%TYPR='Q'
BTABOS%ADR(IND_T)%P%TYPR='Q

with following comments
! IMPORTANT:
! STATES THAT ATABOS AND BTABOS ARE NOT ZERO (SEE LIMI3D AND DIFF3D)
! OTHERWISE THEY WILL NOT BE CONSIDERED

I looked in limi3d.f and diff3d.f and couldn't find a description of "ATABOS" and "BTABOS".

Thank you in advance for any clarification.
Peter

To complement my previous posts:

The calculation without the two lines lead to ATABOS and BTABOS = 0 and the exchange with the atmosphere seems the not to be taken in account.

I finally succeeded to obtain a calculation running with this two lines activated, but only if there is no dry node. This seems restrictive and a bit strange to me, any advices? (I already tryed varying the turbulence model and time step as well as solver and advection scheme.)

I would be thankful for any help

Best regards

Peter

Did you solve your problems?
I have some similar situations