open TELEMAC-MASCARET The mathematically superior suite of solvers

Good afternoon,
I am trying to carry out a telemac 2d simulation with a lagrangian drift.
Without any lagrangian drit, I have no trouble to run telemac and the results are sensible.
However, adding a lagrangian drift in the cas file make the simulation fail at the time corresponding to the one of the start of lagnragian drift (DEBLAG in the fortran file).
As suggested in the telemac manual, I added in the cas file the variable A and G for print out and I integrated and modified the LAGRAN fortran script in mine. I did the tests in both scalar and paralel mode (telemac v6p2 and v6p3)

Does someone have any idea on the origin of the problem and a potential solution to it ?

Thanks in advance for any answer
Best regards
Renaud

Hello,

I checked and we have no test case for Lagrangian drifts, so probably something has been broken and we were not aware of it.

If you can provide your case so that I can work on it it would be nice. If not tell me and I'll set up another test case.

If your case can be given as an example to other users it would be even better.

With best regards,

Jean-Michel Hervouet

Hello,
Thanks for your quick reaction

What I found so far is that the lagragian drift works if
FINLAG= DEBLAG+1 in my fortan file
otherwise it fails.
please find in attachment the files I am working with (I run telemac 2d v6p3 on a windows 7 machine).

As you will see, the run fails with the current version of the fortran file. But if you uncomment line 366 and comment line 367, then you impose FINLAG= DEBLAG+1 and everything works.

Thanks for your help
best regards
Renaud Hostache

No files attached, maybe in a zip it will be OK ?

JMH

I do not now what went wrong

I got the files, thank you. When trying to run the case it immediately stops on file par.txt, which is a problem, because you open it on unit 1, which is the unit of the geometry file. Moreover this file starts with two integers and continues with reals whereas you try to read only two double precision numbers.

JMH

No, this was due to a backslash in the path to par.txt, while all other paths of your files have slashes. No I get probably the error that you have with Lagrangian drifts.

JMH

Thanks for your answers, I modified things as you suggested:
I read par in unit 10 now,
I replace the the backslash by a slash and
I converted the two first lines into doubles in the par.txt file.
But the issue with lagrangian drifts still remain.
Thanks again a lot for your help
Best regards
Renaud

Hello,

A first problem was that the drifts are written in the result file at iteration 0 before they are initialised, but there is something else more difficult to find after... patience...

Regards,

JMH

Hello,

Here is a correction (subroutines preres_telemac2d and derlag), I hope it works for you. However it will not work in parallel, that would be a rather long development that needs to be decided and scheduled.

With best regards,

Jean-Michel Hervouet