open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello Chi-Tuan,

I allowed myself to continue this post as another issue occurred. First I ran the T3D+THERMIC to calibrate the related coefficients and the modeled temperature matched the field data, but after the EUTRO module was included the temperature became much higher, from the figure attached we can see the values were somehow amplified...Do you have any idea why this happened? All the boundary and meteorological conditions were the same.

Best regards,
Yunhao

Hello Yunhao,

Which release of TELEMAC do you use? Can you upload your different steering files (T3D and WAQTEL steering files) for the 2 computations + Fortran files?
Maybe you can try your steering files for a simpler geometry file (e.g. a basin at rest) and see if the issue still occurs?

Running the heat_exchange example and adding the coupling with EUTRO (WATER QUALITY PROCESS = 55), the temperature of the test case seems similar with the original example. I had to decrease the time step + relax some solver parameters for the diffusion of tracers, and it seemed to work. What about the listing, are the convergence of solvers OK?

Chi-Tuan

Hello Chi-Tuan,

Thanks for your time, after a few tests I found that the problem was caused by setting the WATER QUALITY PROCESS = 330, with 55 the problem was gone. I use the latest v8p2.

There's another issue I encountered regarding those keywords with two parameters assigned to in the WAQTEL steering file, for example:

COEFFICIENTS OF ALGAL MORTALITY AT 20C = 0.08; 0.003

the latter one (0.003) was not actually read into the corresponding subroutine algae_death. If we print out these two parameters by simply adding the two lines below in algae_death.f,

WRITE(LU,*) 'M1= ',CMOR(1)
WRITE(LU,*) 'M2= ',CMOR(2)

we can find in the listing below the CMOR(2) was assigned with the wrong parameter.
M1= 8.0E-002
M2= 8.0E-002

This looks like a bug needs to be solved.

Best regards,
Yunhao Song

Hello Yunhao,

Good news that you have solved your issue.
The EUTRO module is a combo of BIOMASS and O2 modules, so no need to couple them with EUTRO.

Thanks for raising this bug, it is a problem when reading the keyword in lecdon_waqtel, there was the same problem for the keyword COEFFICIENTS OF AERATION FORMULA fixed in the trunk.
Anyway, checking other keywords, there is also a bug when reading the 2nd coefficient of the keyword ALGAL TOXICITY COEFFICIENTS that you may also encounter.

Please find enclosed a fix for this 3 bugs, it will be available in next release.

Hope this helps and thanks for feedback,

Chi-Tuan