open TELEMAC-MASCARET The mathematically superior suite of solvers

hi,all
I want to set 59994 sources in my domain,and I tried to do it via the program <source_telemac2d.f>.Attached my fortran file below. But the program seemed wrong,because the value of tracer worked out is zero all the time.
Can somebody help me look at my program?

Hello,

Looking at your Fortran file, I see things that cannot work in the following:

OPEN(1,FILE='xy.txt') !read the CADXY file
OPEN(2,FILE='TSCE.txt') !read the tracer values file
DO I=1,NREJET
READ(1,*) XSCE(I), YSCE(I) !ABSCISSAE OF SOURCES,ORDINATES OF SOURCES
READ(2,*) TSCE(I,1) !tracer values
DSCE(I)=0.1 !water discharge
ENDDO
CLOSE(1)
CLOSE(2)

There are 2 problems: the executable is in a temporary folder, so seen from this folder your files xy.txt and TSCE.txt should be '../xy.txt' and '../TSCE.txt' on Linux and possibly \ instead of / on Windows. I do not understand why this did not crash in your run (you should have a message file xy.txt not found).

The other problem is that channel 1 is taken for the geometry file, so instead of channels 1 and 2 just choose 97 and 98, which are free.

Another thing that I do not understand is why you have so many sources. If it is a source per point in the mesh we could do it more simply, for example by using the rain which can be made variable in space in subroutine prosou.f.

With best regards,

Jean-Michel Hervouet

Thank you. Yeah, It is a source per point in the mesh.

I have check the prosou.f,this program is PREPARES THE SOURCE TERMS IN THE CONTINUITY EQUATION AND IN THE DYNAMIC EQUATIONS.

But I want to PREPARES THE SOURCE TERMS IN THE TRACER EQUATION

How can I do?

Hello,

I would try with the rain. If you condider your sources as rain variable in space all the complicated work of integrating the sources in the various continuiy equations will be already done. For the tracers there is a keyword "VALUES OF TRACERS IN THE RAIN" that might do it, if the values depend on the points we probably can program it easily by replacing the values of the keyword by a local value.

With best regards,

Jean-Michel Hervouet