open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

I'm calibrating a large estuary model. Is it possible to save the results of certain grid points (calibration points) every timestep and save the results of all other points for example every hour?

regards,

Sven

Hi Sven

This is possible and very easy to do.
You just have to use the KEYWORD "LIST OF POINTS" which is associate with "NAME OF POINTS"
With this keywords, you will have the result written at each time step in your listing file.
If you prefer to send this result in a specific file, I think you will be able to modify the preres_telemac2d.f subroutine in you personal Fortran file to change the output.

Regards
Christophe

Hi Christophe,

I used the keywords "LIST OF POINTS" and "NAME OF POINTS" but I didn't get results in my listing file. I get my mass balance results right in the listing but nothing about my points. Is there another keyword required?

Cheers,
Simon

Hi Simon
Those 2 keywords are sufficient but you could only print the value of MAXZ, TMXZ, MAXV, TMXV so if you dont ask for these variables in your steering file, you have no print in your listing.
If you want to print other value, you should modify the subroutine preres_telemac2d.f

Hope this helps

Hello,

There is a user subroutine which is called at every time step, where you can write any output that you fancy. The name of the subroutine is "utimp" or "utimp_telemac2d" if I remember well. In this subroutine all variables are seen thanks to a line "USE DECLARATIONS_TELEMAC2D", you only need to be cautious with these declarations, most variables like the depth H being declared as BIEF_OBF structure, so that their double precision data are in component H%R.

With best regards,

Jean-Michel Hervouet

Hello,
I would like to write results in defined output points too.
As Jean-Michel said, I use the utimp subroutine.
The problem is that I run my simulation in parallel mode and I think that I can only write the PID 0 results..
Is there a solution to this problem ?
Thanks

Youen

Hello,

A tentative code sample below (not tested), suppose that you want to write the velocity U of point 123 (original mesh):

DOUBLE PRECISION P_DMAX,P_DMIN
EXTERNAL P_DMAX,P_DMIN
DOUBLE PRECISION U123

U123=0.D0
DO I=1,NPOIN
IF(MESH%KNOLG%I(I).EQ.123) U123=U%R(I)
ENDDO
U123=P_DMAX(U123)+P_DMIN(U123)
WRITE(LU,*) 'U OF POINT 123=',U123


Note that: PRINT*,'U OF POINT 123=',U123

would cause all the processors to write on the screen. With the WRITE, LU of processors other than 0 is redirected to log files.

IF(IPID.EQ.0) PRINT*,'U OF POINT 123=',U123

also works.

With best regards,

Jean-Michel Hervouet