TELEMAC-MASCARET Forum: Telemac 2D coupling with sisyphe (1/1)

Welcome, Guest

TOPIC: Telemac 2D coupling with sisyphe

Telemac 2D coupling with sisyphe 12 years 5 months ago #4795

nguyenthong_BK

Dear all,
I have a steering file (run.cas) for hydrodynamic model as following and program works very well in my case:

/---------------------------------------------------------------------
/ TELEMAC2D Version v6p0 Apr 9, 2012
/ TELEMAC 2D : CUA SONG SOAI RAP
/---------------------------------------------------------------------
/ ENTREES-SORTIES, FICHIERS
/---------------------------------------------------------------------
FICHIER DES RESULTATS              ='KetQua.res'
FICHIER DES CONDITIONS AUX LIMITES ='Conlim.cli'
FICHIER DES FRONTIERES LIQUIDES    ='QH.qsl'
FICHIER FORTRAN                    ='t2d_Fortran.f'
FICHIER DES FONDS                  ='Land_100_River.xyz'
FICHIER DE GEOMETRIE               ='Geo.slf'
FICHIER DES PARAMETRES             ='Run.cas'
FICHIER DE DONNEES FORMATE 2       ='Table_Trieu.txt'
TITRE ='TELEMAC 2D : CUA SONG SOAI RAP'
/---------------------------------------------------------------------
/ ENTREES-SORTIES, GRAPHIQUES ET LISTING
/---------------------------------------------------------------------
PERIODE DE SORTIE LISTING             =360
VARIABLES POUR LES SORTIES GRAPHIQUES ='H,U,V,S,B,W,Q,X,Y,P,M'
PERIODE POUR LES SORTIES GRAPHIQUES   =60
BILAN DE MASSE                        =OUI

SECTIONS DE CONTROLE           =35394;31602  ;  47271;37043
/---------------------------------------------------------------------
/ EQUATIONS
/---------------------------------------------------------------------
COEFFICIENT D'INFLUENCE DU VENT       =0.6D-03
MODELE DE TURBULENCE                  =1
COEFFICIENT DE FROTTEMENT             =35.D0
OPTION POUR LA DIFFUSION DES VITESSES =2
COEFFICIENT DE DIFFUSION DES VITESSES =1.D-04
LOI DE FROTTEMENT SUR LE FOND         =2
/---------------------------------------------------------------------
/ EQUATIONS, CONDITIONS INITIALES
/---------------------------------------------------------------------
CONDITIONS INITIALES ='COTE CONSTANTE'
COTE INITIALE        =20.
/---------------------------------------------------------------------
/ EQUATIONS, CONDITIONS LIMITES
/---------------------------------------------------------------------
DEBITS IMPOSES                      =0;0;0;0
PROFILS DE VITESSE                  =1;1;1;1
COTES IMPOSEES                      =0;0;0;0
OPTION POUR LES FRONTIERES LIQUIDES =2;1;1;2;2
/---------------------------------------------------------------------
/ PARAMETRES NUMERIQUES
/---------------------------------------------------------------------
PRECONDITIONNEMENT             =2
TRAITEMENT DU SYSTEME LINEAIRE =2
/-------------------
/ INITIALISATION 
PAS DE TEMPS                   =30.
PROCESSEURS PARALLELES         =0
NOMBRE DE PAS DE TEMPS         =23040
STOCKAGE DES MATRICES          =3
FORME DE LA CONVECTION         =1;5
/---------------------------------------------------------------------
/ PARAMETRES NUMERIQUES, SOLVEUR
/---------------------------------------------------------------------
OPTION DU SOLVEUR                    =3
SOLVEUR                              =1
PRECISION DU SOLVEUR                 =1.E-4
MAXIMUM D'ITERATIONS POUR LE SOLVEUR =200
/---------------------------------------------------------------------
/ PARAMETRES NUMERIQUES, VITESSE-CELERITE-HAUTEUR
/---------------------------------------------------------------------
IMPLICITATION POUR LA HAUTEUR =0.55
MASS-LUMPING SUR H            =1.
ORDRE DU TIR INITIAL POUR H   =1
IMPLICITATION POUR LA VITESSE =0.55

However, when I would like coupling with sisyphe by to add in the steering file as:

FICHIER DES PARAMETRES DE SISYPHE = 'sis_run.cas'
COUPLAGE AVEC = 'INTER-SISYPHE'
PERIODE DE COUPLAGE = 5

Le program has an error: ……file not found, unit 27 ,…… when it was executed. This is the file declared by : FICHIER DE DONNEES FORMATE 2 ='Table_Trieu.txt' in steering file (run.cas). Attention that this is the file where I give the boundary condition of hydrodynamic problem and the program normally read this file (I’m always in the same directory).
I would like to receive your suggestion.
Thanks you so much.
Thong NGUYEN
HCMUT, Viet Nam

The administrator has disabled public write access.

Re: Telemac 2D coupling with sisyphe 12 years 5 months ago #4796

jmhervouet

Hello Thong,

This is normal, you cannot use hardcoded unit numbers like 27 in coupling (there is a competition for unit numbers between the coupled programs). For overcoming this problem, the unit numbers are redistributed, so you must use the official parameterisation which is, instead of 27:

T2D_FILES(T2DFO1)%LU

T2DF01 is an integer equal to 27, T2D_FILES is an array of structures BIEF_FILE, each structure containing components which give information on a file, and the component LU is the logical unit. Without coupling you have T2D_FILES(T2DFO1)%LU=27, in coupling it will change.

If you do not have T2D_FILES and T2DF01 in your subroutine, add at the beginning:

USE DECLARATIONS_TELEMAC2D, ONLY : T2D_FILES,T2DFO1

I'm sure this will help...

Wth best regards,

Jean-Michel Hervouet

The administrator has disabled public write access.

Re: Telemac 2D coupling with sisyphe 12 years 5 months ago #4797

nguyenthong_BK	Dear Mr. Jean-Michel Hervouet, Thank you very much for your help. My problem was solved as your suggestion. Best regards Thong NGUYEN, HCMUT, Viet Nam
	The administrator has disabled public write access.