open TELEMAC-MASCARET The mathematically superior suite of solvers

I am trying to output data through the routine 'UTIMP_TELEMAC2D'. The routine I have setup is very simple, but I cannot figure out why it does not work. Does anyone know what could be amiss? Custom routine is attached. I am running v6p3r2, serial mode, number of nodes is 2026. Much appreciated!

Hello,

Unit 1 is the geometry file, which is always open, so you have little chance. Try unit 99 instead, it will probably work better.

Regards,

Jean-Michel Hervouet

That's obvious... thanks!

Attached is fixed custom output routine for others.

Hello,

Hum... you have left the CLOSE(1) instead of CLOSE(99). I was just thinking that the geometry file could be closed after use, whereas we have a subroutine at the end which closes files. The lack of CLOSE(99) is apparently not a problem.

Regards,

JM

good catch!

I have another related question. I am trying to write tracer values to a text file. The problem is that tracer values appear to be wrong. The model is initialized with tracers equal to zero, I am monitoring tracers at the downstream boundary and tracers are being injected at the upstream boundary. The max tracer value possible is 1 because that is the highest concentration injected. The tracer outputs are always more than 1 (about 1.5), and at t=0 the tracer value is 8.0! Visualizing results in Rhino confirms these custom outputs are wrong. Any thoughts? Customization and case file are attached.

Much appreciated

Hello,

T1 is not the first tracer but a work array that may contain anything. Tracers are in block T. What you must print is :

T%ADR(1)%P%R(2306)

For tracer 2 it would be ADR(2) instead of ADR(1).

Regards,

JMH