open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

What is the best way to visualize the evolution of flow in a siphon ?

Regards

Hello,

Hum, I would take subroutine siphon.f in my Fortran file. You can see there the line that prints the discharge in the listing, it is the sequennce embedded into:

IF(ENTET) THEN
...
ENDIF

You can just after print the discharge on a file. Suppose that you want to fill a file with lines, every line with the time followed by the discharge.

In my example your file will be declared in the parameter file as :

FORMATTED RESULTS FILE (FICHIER DE RESULTATS FORMATE)

for example :

FORMATTED RESULTS FILE : siphon_discharge

Then in subroutine siphon.f :

Add after the line with USE_BIEF the line :

USE DECLARATIONS_TELEMAC2D, ONLY : T2D_FILES,AT

Just after the block IF(ENTET)... add the line :

WRITE(T2D_FILES(T2DRFO)%LU,*) AT,DSCE(I2)


T2D_FILES(T2DRFO)%LU is 29, except in case of coupling with Sisyphe or Tomawac, where it may change. You need not open or close the file, this is done elsewhere for all files.

Then you can use a simple program like maybe Excel to draw the curve after what is in your file..

I have not tested but it should work about this way,

With best regards,

Jean-Michel Hervouet

Thank you...

I will try to test this method

Regards

Thank you
The only problem with this solution is if you have more than one siphon, it write :
AT DSCE(spihon1)
AT DSCE(spihon2)
AT DSCE(spihon3)
AT DSCE(spihon4)
AT+dt DSCE(spihon1)
AT+dt DSCE(spihon2)
AT+dt DSCE(spihon3)
AT+dt DSCE(spihon4)
...

How can we have in each formatted results file only one flow of siphon ?

Thank for your help

Regards

Marine

Hello,

In this case you will have to deal yourself with your files, open and close them, keeping in mind that you must name them ../file1, ../file2, or possibly ..\file1 on Windows, because telemac-2D is executed in a temporary folder which is destroyed after job completion. Then you will have the problem of opening the files only once at the beginning, and you can maybe use the command INQUIRE to know if the file already exists.

Hope this helps,

Jean-Michel Hervouet

Hey,

Thank for your help !
The only problem is that this method works but just in scalar mode and I am in parallel mode

Regards

Marine

Then you should manage the processor number in order to avoid the attemps to use the same file from 2 different processors.
This is a little more complex but possible.

If it's just for some analysis, you could use the listing file in a text editor and select all lines that contain "DISCHARGE CLUVERT X" and copy to excel.

Not realy beautiful but easy to do idf you don't know how to program the write in parallel.

Hope this helps

:) I am sure with a basic programmation I can find a way to extract in different files the flow of each siphon.
My only problem with this solution is that exel is limited in number of line so I need to find a way to do this in another software like notepad

Thank for your help

Regards

Marine

Yes
With notepad it's quite easy
You select the string to search (in the mark tab) then you mark all the lines
You just have to select in the serch menu / bookmark copy all the mark lines


regards

Hello everyone,

I managed to generate the output for my culvert by implementing your suggestions into buse.f! I´m just wondering if the maximum output of parameters is limited to 6?
I would like to get the following parameters (attached my modified buse line):
AT,V2DPAR%R(I1),QMAX1,Q,DBUS(N),S1,S2,FTYP

The results file only prints until S1 and misses S2 and FTYP.
Is the file somehow limited in columns or any other reason?

Further the "VD2PAR%R" is confusing me. In scalar mode of my 1x1m regular grid it shows different areas, depending on the direction of the edges (see pictures).How is the area calculated exactly? In parallel mode it even drops to 6,7E-270. Any ideas, where this comes from?



Thanks
Best regards
Chris

ps: I´m running on V7p3r1, parallel