open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello all,

I added in my .cas file the two following key words : "FICHIER DES SECTIONS DE CONTROLE" and "FICHIER DE SORTIE DES SECTIONS DE CONTROLE" and added in the same folder as my .cas file the following "sect_controle.txt" and "sorties_sect_controle.txt" files so that my .cas file looks like :

FICHIER DES SECTIONS DE CONTROLE : sect_controle.txt
FICHIER DE SORTIE DES SECTIONS DE CONTROLE : sorties_sect_controle.txt

But when running Telemac I've got the following error message :

"... reading module dictionary
... simulation en Francais avec TelemacBayonne.cas
re-writing: /home/globc/barthelemy/HYDRO/MASCARET_DA/MascaretDA/v6p2r1/telemac2d/Bayonne/TelemacBayonne.cas_2014-03-06-11h31min53s/T2DCAS
copying: TelemacBayonne_ConditionInitiale.slf /home/globc/barthelemy/HYDRO/MASCARET_DA/MascaretDA/v6p2r1/telemac2d/Bayonne/TelemacBayonne.cas_2014-03-06-11h31min53s/T2DPRE
copying: telemac2dv6p2.dico /home/globc/barthelemy/HYDRO/MASCARET_DA/MascaretDA/v6p2r1/telemac2d/Bayonne/TelemacBayonne.cas_2014-03-06-11h31min53s/T2DDICO
copying: TelemacBayonne_ConditionsLimites.liq /home/globc/barthelemy/HYDRO/MASCARET_DA/MascaretDA/v6p2r1/telemac2d/Bayonne/TelemacBayonne.cas_2014-03-06-11h31min53s/T2DIMP
copying: TelemacBayonne_TypeDeConditionsLimites.cli /home/globc/barthelemy/HYDRO/MASCARET_DA/MascaretDA/v6p2r1/telemac2d/Bayonne/TelemacBayonne.cas_2014-03-06-11h31min53s/T2DCLI
copying: TelemacBayonne_Geometrie_Ks.slf /home/globc/barthelemy/HYDRO/MASCARET_DA/MascaretDA/v6p2r1/telemac2d/Bayonne/TelemacBayonne.cas_2014-03-06-11h31min53s/T2DGEO
runCAS:
|processLIT: file does not exist: sect_controle.txt"

Do you know what is the problem ? The "sect_controle.txt" file should be in an other folder ?

Hello

This looks strange. At a first glance, we can not conclude about the problem.
Try first to put the names between quotes "sect_controle.txt" and verify that line length does not exceed 73 characters.
otherwise, send us your files and I will diagnose the problem.
With my best regards

Riadh ATA

The problem is that I made an error in the file name. But now when I run Telemac it begins normally and then it fails to run with the following error message :

*****************************
* ALLOCATION DE LA MEMOIRE *
*****************************

READGEO1 : TITRE= newSelafin

NOMBRE D'ELEMENTS: 8388
NOMBRE REEL DE POINTS: 4276
MXPTEL (BIEF) : NOMBRE MAXIMUM D'ELEMENTS VOISINS D' UN POINT : 9
NOMBRE MAXIMUM DE POINTS VOISINS D'UN POINT : 10
SEGBOR (BIEF) : NOMBRE DE SEGMENTS DE BORD = 162
EN COMPTANT CEUX DUS A LA DECOMPOSITION DE DOMAINE
CORRXY (BIEF) : PAS DE MODIFICATION DES COORDONNEES

MAILLAGE : MESH ALLOUE
POINT_TELEMAC2D: SECTIONS DEFINED IN THE SECTIONS INPUT FILE
rank 0 in job 1 guinea.cerfacs.fr_44272 caused collective abort of all ranks
exit status of rank 0: return code 2
_____________
runcode::main:
/home/globc/barthelemy/HYDRO/MASCARET_DA/MascaretDA/v6p2r1/telemac2d/Bayonne:
|runCAS: fail to run
| /usr/lib64/mpich2/bin/mpirun -wdir /home/globc/barthelemy/HYDRO/MASCARET_DA/MascaretDA/v6p2r1/telemac2d/Bayonne/TelemacBayonne.cas_2014-03-06-16h11min08s -n 8 /home/globc/barthelemy/HYDRO/MASCARET_DA/MascaretDA/v6p2r1/telemac2d/Bayonne/TelemacBayonne.cas_2014-03-06-16h11min08s/out_telemac2dv6p2

and my sect_controle.txt is as follows :

# Control sections definition
1 -1
essai_Pont_Blanc
31945 27801

I have a last question : is it possible to save in a txt file only some variables such a U, V, H ? If so how to do it ?

I mean some variables but at only one point.

Hello

yes sure you can, you have to write some Fortran lines in the your fortran file. However, the easiest way to have these data is to do it by post-precessing the result file.
Actually, with BlueKenue or Fudaa,(or Tecplot) you can export any variable you wish to txt file. This is easy and free, because when you ask Telemac to write these data directly, it will cost CPU and memory and may introduce errors.

With my best regards,
Riadh ATA

this option is also achieved by post-processing tools.

Thank you Riadh, that's what I finally did, I used BlueKenue to compute the velocity at a given node.

I have an other question I want to compute some quantities that Telemac2d doesn't compute by itself such as the scalar product of the velocity with a given vector what is the "best" subroutine fortran to do this ? I mean : where is it the most suitable ?

Hello

Look at routine p_dot.f (library bief, in folder Sources/utils).
It contains a function that computes dot product between two vectors X and Y and it takes into account also parallel aspects.

With my best regards

Riadh ATA