|
Welcome,
Guest
|
|
TOPIC: error in the simulation telemac2D
error in the simulation telemac2D 9 years 3 months ago #17796
|
Hi
I run a simulation using new formula of sediment transport that I added in subroutine qsforme , I added this new subroutine to my fortran file ( with other subroutine that I used ) But the simulation indicates the following error: C:\Users\Desktop\Baptist2>telemac2D.py cas.txt
Loading Options and Configurations
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
... parsing configuration file: C:\opentelemac-mascaret\v6p2\config\systel.cfg
Running your CAS file for:
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+> configuration: win7gfor
+> root: C:\opentelemac-mascaret\v6p2
+> version v6p2
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Running cas.txt with telemac2d under C:\Users\Desktop\Baptist2
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
... reading module dictionary
... simulation en Francais avec cas.txt
... simulation en Francais avec transport.txt
copying: conlim C:\Users\morri\Desktop\Baptist2\cas.txt_2015-07-24-14h37min
13s\T2DCLI
copying: flume2014.f C:\Users\morri\Desktop\Baptist2\cas.txt_2015-07-24-14h
37min13s\t2dfort.f
re-writing: C:\Users\morri\Desktop\Baptist2\cas.txt_2015-07-24-14h37min13s\T2D
CAS
copying: geo C:\Users\morri\Desktop\Baptist2\cas.txt_2015-07-24-14h37min13s
\T2DGEO
copying: telemac2dv6p2.dico C:\Users\morri\Desktop\Baptist2\cas.txt_2015-07
-24-14h37min13s\T2DDICO
re-writing: C:\Users\morri\Desktop\Baptist2\cas.txt_2015-07-24-14h37min13s\SIS
CAS
copying: geo C:\Users\morri\Desktop\Baptist2\cas.txt_2015-07-24-14h37min13s
\SISGEO
copying: sisyphev6p2.dico C:\Users\morri\Desktop\Baptist2\cas.txt_2015-07-2
4-14h37min13s\SISDICO
copying: conlim C:\Users\morri\Desktop\Baptist2\cas.txt_2015-07-24-14h37min
13s\SISCLI
Running your simulation :
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C:\Users\\Desktop\Baptist2\cas.txt_2015-07-24-14h37min13s\out_flume2014.exe
LISTING DE TELEMAC-2D ---------------------------------------------------------
---------------------
TTTTT EEEEE L EEEEE M M AAAAA CCCCC
T E L E MM MM A A C
T EEE L EEE M M M AAAAA C
T E L E M M A A C
T EEEEE LLLLL EEEEE M M A A CCCCC
2D VERSION 6.2 FORTRAN 90
WITH SEVERAL TRACERS
COUPLED WITH SISYPHE AND TOMAWAC
********************************************
* LECDON: *
* APRES APPEL DE DAMOCLES *
* VERIFICATION DES DONNEES LUES *
* SUR LE FICHIER DES PARAMETRES *
********************************************
SORTIE DE LECDON. TITRE DE L'ETUDE :
Couplage avec sisyphe
OUVERTURE DES FICHIERS POUR TELEMAC2D
*****************************
* ALLOCATION DE LA MEMOIRE *
*****************************
READGEO1 : TITRE= C:\Users\KOUKA\Desktop\Baptist2\/canal.MAT
NOMBRE D'ELEMENTS: 3236
NOMBRE REEL DE POINTS: 1786
MXPTEL (BIEF) : NOMBRE MAXIMUM D'ELEMENTS VOISINS D'
UN POINT : 8
NOMBRE MAXIMUM DE POINTS VOISINS D'UN
POINT : 8
CORRXY (BIEF) : PAS DE MODIFICATION DES COORDONNEES
MAILLAGE : MESH ALLOUE
****************************************
* FIN DE L'ALLOCATION DE LA MEMOIRE : *
****************************************
LISTING DE SISYPHE AVEC COUPLAGE-----------------------------------------------
-------------------------------
SSSS I SSSS Y Y PPPP H H EEEEE
S I S Y Y P P H H E
SSS I SSS Y PPPP HHHHH EEEE
S I S Y P H H E
SSSS I SSSS Y P H H EEEEE
VERSION 6.2
COUPLED WITH TELEMAC-2D INTERNALLY
ATTENTION ! A LA LIGNE 2554 DU DICTIONNAIRE :
LE NOMBRE DE VALEURS PAR DEFAUT 10 EST DIFFERENT DE LA TAILLE ANNONC
EE 2
ATTENTION ! A LA LIGNE 2637 DU DICTIONNAIRE :
LE NOMBRE DE VALEURS PAR DEFAUT 10 EST DIFFERENT DE LA TAILLE ANNONC
EE 2
ATTENTION ! A LA LIGNE 2675 DU DICTIONNAIRE :
LE NOMBRE DE VALEURS PAR DEFAUT 10 EST DIFFERENT DE LA TAILLE ANNONC
EE 2
ATTENTION ! A LA LIGNE 3357 DU DICTIONNAIRE :
LE NOMBRE DE VALEURS PAR DEFAUT 13 EST DIFFERENT DE LA TAILLE ANNONC
EE 1
********************************************
* LECDON: *
* APRES APPEL DE DAMOCLES *
* VERIFICATION DES DONNEES LUES *
* SUR LE FICHIER DES PARAMETRES *
********************************************
OUVERTURE DES FICHIERS POUR SISYPHE
*******************************
* ALLOCATION DE LA MEMOIRE *
*******************************
READGEO1 : TITRE= C:\Users\KOUKA\Desktop\Baptist2\/canal.MAT
NOMBRE D'ELEMENTS: 3236
NOMBRE REEL DE POINTS: 1786
MXPTEL (BIEF) : NOMBRE MAXIMUM D'ELEMENTS VOISINS D'
UN POINT : 8
NOMBRE MAXIMUM DE POINTS VOISINS D'UN
POINT : 8
CORRXY (BIEF) : PAS DE MODIFICATION DES COORDONNEES
MAILLAGE : MESH_S ALLOUE
****************************************
* FIN DE L'ALLOCATION DE LA MEMOIRE : *
****************************************
INBIEF (BIEF) : MACHINE NON VECTORIELLE (SELON VOS DONNEES)
STRCHE (BIEF) : PAS DE MODIFICATION DU FROTTEMENT
IL Y A 2 FRONTIERE(S) LIQUIDE(S) :
FRONTIERE 1 :
DEBUT AU POINT DE BORD 160 , DE NUMERO GLOBAL 11
ET DE COORDONNEES : 15.85000 0.000000
FIN AU POINT DE BORD 168 , DE NUMERO GLOBAL 13
ET DE COORDONNEES : 15.85000 0.8000000
FRONTIERE 2 :
DEBUT AU POINT DE BORD 327 , DE NUMERO GLOBAL 15
ET DE COORDONNEES : 0.000000 0.8000000
FIN AU POINT DE BORD 1 , DE NUMERO GLOBAL 2
ET DE COORDONNEES : 0.000000 0.000000
IL Y A 2 FRONTIERE(S) SOLIDE(S) :
FRONTIERE 1 :
DEBUT AU POINT DE BORD 1 , DE NUMERO GLOBAL 2
ET DE COORDONNEES : 0.000000 0.000000
FIN AU POINT DE BORD 160 , DE NUMERO GLOBAL 11
ET DE COORDONNEES : 15.85000 0.000000
FRONTIERE 2 :
DEBUT AU POINT DE BORD 168 , DE NUMERO GLOBAL 13
ET DE COORDONNEES : 15.85000 0.8000000
FIN AU POINT DE BORD 327 , DE NUMERO GLOBAL 15
ET DE COORDONNEES : 0.000000 0.8000000
CORFON (TELEMAC2D) : PAS DE MODIFICATION DU FOND
================================================================================
ITERATION 0 TEMPS : 0.0000 S
TELEMAC2D COUPLE AVEC : INTER-SISYPHE
TELEMAC-2D : COUPLAGE INTERNE AVEC SISYPHE
PART 0 : INITIALISING SISYPHE
INBIEF (BIEF) : MACHINE NON VECTORIELLE (SELON VOS DONNEES)
IL Y A 2 FRONTIERE(S) LIQUIDE(S) :
FRONTIERE 1 :
DEBUT AU POINT DE BORD 160 , DE NUMERO GLOBAL 11
ET DE COORDONNEES : 15.85000 0.000000
FIN AU POINT DE BORD 168 , DE NUMERO GLOBAL 13
ET DE COORDONNEES : 15.85000 0.8000000
FRONTIERE 2 :
DEBUT AU POINT DE BORD 327 , DE NUMERO GLOBAL 15
ET DE COORDONNEES : 0.000000 0.8000000
FIN AU POINT DE BORD 1 , DE NUMERO GLOBAL 2
ET DE COORDONNEES : 0.000000 0.000000
IL Y A 2 FRONTIERE(S) SOLIDE(S) :
FRONTIERE 1 :
DEBUT AU POINT DE BORD 1 , DE NUMERO GLOBAL 2
ET DE COORDONNEES : 0.000000 0.000000
FIN AU POINT DE BORD 160 , DE NUMERO GLOBAL 11
ET DE COORDONNEES : 15.85000 0.000000
FRONTIERE 2 :
DEBUT AU POINT DE BORD 168 , DE NUMERO GLOBAL 13
ET DE COORDONNEES : 15.85000 0.8000000
FIN AU POINT DE BORD 327 , DE NUMERO GLOBAL 15
ET DE COORDONNEES : 0.000000 0.8000000
STRCHE (BIEF) : PAS DE MODIFICATION DU FROTTEMENT
INITIALISATION EN CAS DE COUPLAGE : PART= 0
STOP :
LE TAUX DE TRANSPORT DOIT ETRE CALCULE DANS QSFORM
PLANTE : ARRET DU PROGRAMME APRES ERREUR
_____________
runcode::main:
C:\Users\Desktop\Baptist2:
|runCAS: fail to run
| C:\Users\Desktop\Baptist2\cas.txt_2015-07-24-14h37min13s\out_flu
me2014.exeI didn’t understand this error to correct Thank you in advance |
|
The administrator has disabled public write access.
|
error in the simulation telemac2D 9 years 3 months ago #17797
|
Hi
the message seems quite clear! You probably made a mistake in your adapted version of qsform. Sisyphe indicate that the transport rate ("Taux de transport") was not computed in the subroutine. Maybe you should join the subroutine to allow us to well understand what you've done regards |
|
Christophe
The administrator has disabled public write access.
|
error in the simulation telemac2D 9 years 3 months ago #17798
|
Hello
thank you for your reply and this the subroutine that I added I changed some parameter thank you in advance
Attachments:
|
|
The administrator has disabled public write access.
|
error in the simulation telemac2D 9 years 3 months ago #17799
|
Hello,
You can see in the original subroutine that you can cancel the last part that stops you once you have implemented your formula. However the example shows that QSC and QSS are computed, whereas you have only done QSC (but I am not sure about this). Other remarks: Rather write DM**3 than DM**3.0D0 Rather write 0.1D0 than (1.D0/10.D0) With best regards, Jean-Michel Hervouet |
|
The administrator has disabled public write access.
|
error in the simulation telemac2D 9 years 3 months ago #17801
|
Hello
thank you JMH for your reply, that's help me!! just there is another error that still and that's what there isn't à calculation !!!! (ATTENTION ! A LA LIGNE 2554 DU DICTIONNAIRE :
LE NOMBRE DE VALEURS PAR DEFAUT 10 EST DIFFERENT DE LA TAILLE ANNONC
EE 2)
STING DE TELEMAC-2D ---------------------------------------------------------
---------------------
TTTTT EEEEE L EEEEE M M AAAAA CCCCC
T E L E MM MM A A C
T EEE L EEE M M M AAAAA C
T E L E M M A A C
T EEEEE LLLLL EEEEE M M A A CCCCC
2D VERSION 6.2 FORTRAN 90
WITH SEVERAL TRACERS
COUPLED WITH SISYPHE AND TOMAWAC
********************************************
* LECDON: *
* APRES APPEL DE DAMOCLES *
* VERIFICATION DES DONNEES LUES *
* SUR LE FICHIER DES PARAMETRES *
********************************************
SORTIE DE LECDON. TITRE DE L'ETUDE :
Couplage avec sisyphe
OUVERTURE DES FICHIERS POUR TELEMAC2D
*****************************
* ALLOCATION DE LA MEMOIRE *
*****************************
READGEO1 : TITRE= C:\Users\KOUKA\Desktop\Baptist2\/canal.MAT
NOMBRE D'ELEMENTS: 3236
NOMBRE REEL DE POINTS: 1786
MXPTEL (BIEF) : NOMBRE MAXIMUM D'ELEMENTS VOISINS D'
UN POINT : 8
NOMBRE MAXIMUM DE POINTS VOISINS D'UN
POINT : 8
CORRXY (BIEF) : PAS DE MODIFICATION DES COORDONNEES
MAILLAGE : MESH ALLOUE
****************************************
* FIN DE L'ALLOCATION DE LA MEMOIRE : *
****************************************
LISTING DE SISYPHE AVEC COUPLAGE-----------------------------------------------
-------------------------------
SSSS I SSSS Y Y PPPP H H EEEEE
S I S Y Y P P H H E
SSS I SSS Y PPPP HHHHH EEEE
S I S Y P H H E
SSSS I SSSS Y P H H EEEEE
VERSION 6.2
COUPLED WITH TELEMAC-2D INTERNALLY
ATTENTION ! A LA LIGNE 2554 DU DICTIONNAIRE :
LE NOMBRE DE VALEURS PAR DEFAUT 10 EST DIFFERENT DE LA TAILLE ANNONC
EE 2
ATTENTION ! A LA LIGNE 2637 DU DICTIONNAIRE :
LE NOMBRE DE VALEURS PAR DEFAUT 10 EST DIFFERENT DE LA TAILLE ANNONC
EE 2
ATTENTION ! A LA LIGNE 2675 DU DICTIONNAIRE :
LE NOMBRE DE VALEURS PAR DEFAUT 10 EST DIFFERENT DE LA TAILLE ANNONC
EE 2
ATTENTION ! A LA LIGNE 3357 DU DICTIONNAIRE :
LE NOMBRE DE VALEURS PAR DEFAUT 13 EST DIFFERENT DE LA TAILLE ANNONC
EE 1
********************************************
* LECDON: *
* APRES APPEL DE DAMOCLES *
* VERIFICATION DES DONNEES LUES *
* SUR LE FICHIER DES PARAMETRES *
********************************************
OUVERTURE DES FICHIERS POUR SISYPHE
*******************************
* ALLOCATION DE LA MEMOIRE *
*******************************
READGEO1 : TITRE= C:\Users\KOUKA\Desktop\Baptist2\/canal.MAT
NOMBRE D'ELEMENTS: 3236
NOMBRE REEL DE POINTS: 1786
MXPTEL (BIEF) : NOMBRE MAXIMUM D'ELEMENTS VOISINS D'
UN POINT : 8
NOMBRE MAXIMUM DE POINTS VOISINS D'UN
POINT : 8
CORRXY (BIEF) : PAS DE MODIFICATION DES COORDONNEES
MAILLAGE : MESH_S ALLOUE
****************************************
* FIN DE L'ALLOCATION DE LA MEMOIRE : *
****************************************
INBIEF (BIEF) : MACHINE NON VECTORIELLE (SELON VOS DONNEES)
STRCHE (BIEF) : PAS DE MODIFICATION DU FROTTEMENT
IL Y A 2 FRONTIERE(S) LIQUIDE(S) :
FRONTIERE 1 :
DEBUT AU POINT DE BORD 160 , DE NUMERO GLOBAL 11
ET DE COORDONNEES : 15.85000 0.000000
FIN AU POINT DE BORD 168 , DE NUMERO GLOBAL 13
ET DE COORDONNEES : 15.85000 0.8000000
FRONTIERE 2 :
DEBUT AU POINT DE BORD 327 , DE NUMERO GLOBAL 15
ET DE COORDONNEES : 0.000000 0.8000000
FIN AU POINT DE BORD 1 , DE NUMERO GLOBAL 2
ET DE COORDONNEES : 0.000000 0.000000
IL Y A 2 FRONTIERE(S) SOLIDE(S) :
FRONTIERE 1 :
DEBUT AU POINT DE BORD 1 , DE NUMERO GLOBAL 2
ET DE COORDONNEES : 0.000000 0.000000
FIN AU POINT DE BORD 160 , DE NUMERO GLOBAL 11
ET DE COORDONNEES : 15.85000 0.000000
FRONTIERE 2 :
DEBUT AU POINT DE BORD 168 , DE NUMERO GLOBAL 13
ET DE COORDONNEES : 15.85000 0.8000000
FIN AU POINT DE BORD 327 , DE NUMERO GLOBAL 15
ET DE COORDONNEES : 0.000000 0.8000000
CORFON (TELEMAC2D) : PAS DE MODIFICATION DU FOND
================================================================================
ITERATION 0 TEMPS : 0.0000 S
TELEMAC2D COUPLE AVEC : INTER-SISYPHE
TELEMAC-2D : COUPLAGE INTERNE AVEC SISYPHE
PART 0 : INITIALISING SISYPHE
INBIEF (BIEF) : MACHINE NON VECTORIELLE (SELON VOS DONNEES)
IL Y A 2 FRONTIERE(S) LIQUIDE(S) :
FRONTIERE 1 :
DEBUT AU POINT DE BORD 160 , DE NUMERO GLOBAL 11
ET DE COORDONNEES : 15.85000 0.000000
FIN AU POINT DE BORD 168 , DE NUMERO GLOBAL 13
ET DE COORDONNEES : 15.85000 0.8000000
FRONTIERE 2 :
DEBUT AU POINT DE BORD 327 , DE NUMERO GLOBAL 15
ET DE COORDONNEES : 0.000000 0.8000000
FIN AU POINT DE BORD 1 , DE NUMERO GLOBAL 2
ET DE COORDONNEES : 0.000000 0.000000
IL Y A 2 FRONTIERE(S) SOLIDE(S) :
FRONTIERE 1 :
DEBUT AU POINT DE BORD 1 , DE NUMERO GLOBAL 2
ET DE COORDONNEES : 0.000000 0.000000
FIN AU POINT DE BORD 160 , DE NUMERO GLOBAL 11
ET DE COORDONNEES : 15.85000 0.000000
FRONTIERE 2 :
DEBUT AU POINT DE BORD 168 , DE NUMERO GLOBAL 13
ET DE COORDONNEES : 15.85000 0.8000000
FIN AU POINT DE BORD 327 , DE NUMERO GLOBAL 15
ET DE COORDONNEES : 0.000000 0.8000000
STRCHE (BIEF) : PAS DE MODIFICATION DU FROTTEMENT
INITIALISATION EN CAS DE COUPLAGE : PART= 0
192.46068004249150 6.5424903256652688E-007 6.9481909461161256E-002
13.372535545545920 4.5329686090869722E-002 0.15626573596224702 1
.5328125000000001
QSCL REINITIALISE DANS INIT_TRANSPORT
POUR LA CLASSE 1
QS REINITIALISE DANS INIT_TRANSPORT
================================================================================
ITERATION 0 TEMPS : 0.0000 S
SISYPHE COUPLE AVEC : TELEMAC2D
================================================================================
ITERATION 1 TEMPS : 0.0500 S
USING STREAMLINE VERSION 6.2 FOR CHARACTERISTICS
--------------------------------------------------------------------------------
ETAPE DE CONVECTION
--------------------------------------------------------------------------------
ETAPE DE DIFFUSION-PROPAGATION
GRACJG (BIEF) : 7 ITERATIONS, PRECISION ABSOLUE : 0.6287347E-06
--------------------------------------------------------------------------------
BILAN DE VOLUME D'EAU
VOLUME DANS LE DOMAINE : 3.894912 M3
FLUX FRONTIERE 1 : 0.4245799E-01 M3/S ( >0 : ENTRANT <0 : SORTANT )
FLUX FRONTIERE 2 : -0.000000 M3/S ( >0 : ENTRANT <0 : SORTANT )
ERREUR RELATIVE EN VOLUME A T = 0.5000E-01 S : -0.7512451E-05
================================================================================
ITERATION 1 TEMPS : 0.0500 S
192.46068004249150 6.5424903256652688E-007 6.9481892528737932E-002
13.372532286720212 4.5329675044232694E-002 0.15626565577189824 1
.5328125000000001
EQUATION DE CHARRIAGE RESOLUE EN 1 ITERATIONS
BILAN DE MASSE (EN VOLUME, VIDES INCLUS) :
SOMME DES EVOLUTIONS : NaN M3
FRONTIERE 1 FLUX EN CHARRIAGE = NaN M3/S >0=ENTRANT, AVEC VIDES
FRONTIERE 2 FLUX EN CHARRIAGE = NaN M3/S >0=ENTRANT, AVEC VIDES
TOTAL FLUX EN CHARRIAGE = NaN M3/S >0=ENTRANT, AVEC VIDES
VOLUME PERDU : NaN M3 ( SI <0 VOLUME SORTI )
EVOLUTION MAXIMUM : Infinity NOEUD : 344
EVOLUTION MINIMUM : -100.0000 NOEUD : 25
EVOLUTION MAXIMUM TOTALE : Infinity NOEUD : 344
EVOLUTION MINIMUM TOTALE : -100.0000 NOEUD : 25
================================================================================
ITERATION 2 TEMPS : 0.1000 S
--------------------------------------------------------------------------------
ETAPE DE CONVECTION
--------------------------------------------------------------------------------
ETAPE DE DIFFUSION-PROPAGATION
GRACJG (BIEF) : MAX D'ITERATIONS ATTEINT: 150 PRECISION RELATIVE:
NaN
--------------------------------------------------------------------------------
BILAN DE VOLUME D'EAU
VOLUME DANS LE DOMAINE : NaN M3
FLUX FRONTIERE 1 : 0.9286401E-01 M3/S ( >0 : ENTRANT <0 : SORTANT )
FLUX FRONTIERE 2 : -0.5994257E-03 M3/S ( >0 : ENTRANT <0 : SORTANT )
ERREUR RELATIVE EN VOLUME A T = 0.1000 S : NaN
================================================================================
ITERATION 2 TEMPS : 0.1000 S
192.46068004249150 6.5424903256652688E-007 NaN
NaN NaN NaN 1
.5328125000000001
EQUATION DE CHARRIAGE RESOLUE EN 1 ITERATIONS
BILAN DE MASSE (EN VOLUME, VIDES INCLUS) :
SOMME DES EVOLUTIONS : NaN M3
FRONTIERE 1 FLUX EN CHARRIAGE = NaN M3/S >0=ENTRANT, AVEC VIDES
FRONTIERE 2 FLUX EN CHARRIAGE = NaN M3/S >0=ENTRANT, AVEC VIDES
TOTAL FLUX EN CHARRIAGE = NaN M3/S >0=ENTRANT, AVEC VIDES
VOLUME PERDU : NaN M3 ( SI <0 VOLUME SORTI )
EVOLUTION MAXIMUM : 0.000000 NOEUD : 1
EVOLUTION MINIMUM : 0.000000 NOEUD : 1
EVOLUTION MAXIMUM TOTALE : 0.000000 NOEUD : 1
EVOLUTION MINIMUM TOTALE : -100.0000 NOEUD : 25
================================================================================
ITERATION 3 TEMPS : 0.1500 S
--------------------------------------------------------------------------------
ETAPE DE CONVECTION
INTERP: POINT 1 NOT IN ELEMENT 2003
SHP(1,IP)= NaN
SHP(2,IP)= NaN
SHP(3,IP)= NaN
EPSILO= 9.9999999999999995E-007 IPID= 0
PLANTE : ARRET DU PROGRAMME APRES ERREUR
_____________
runcode::main:
C:\Users\Desktop\Baptist2:
|runCAS: fail to run
| C:\Users\Desktop\Baptist2\cas.txt_2015-07-28-07h51min22s\out_qsf
orm.exethank you in advance |
|
The administrator has disabled public write access.
|
error in the simulation telemac2D 9 years 3 months ago #17803
|
Hi
You probably have to check the computed value of your qsform formulae as there was NaN and infinity at the end of the bedload equation resolution. hope this helps |
|
Christophe
The administrator has disabled public write access.
|
error in the simulation telemac2D 9 years 3 months ago #17807
|
Hello,
Yes, try to add QSS=0.D0 just as you did QSC=your formula, because QSS is also done in the example, but an insight of sediment modelers would be useful here... Regards, JMH |
|
The administrator has disabled public write access.
|
error in the simulation telemac2D 9 years 3 months ago #17816
|
HELLO JMH and C.Coulet
thank you , for your useful sugestions!! I verified my formula and I made some changes ( here in my case I added a new varaibles sisyphe, I tested the calculation of these two varaibles, and I see the new varaibles calculated in resulats, so I used these varaibles in the formula ( I used prive to define the varaible in the new formula ), I added also as proposed JMH ( Qss= 0.D0 ) but when I run the model I find these error : INITIALISATION EN CAS DE COUPLAGE : PART= 0 maybe because I used the new variable in the formula I find this type of problem ?, or these varaible as added but I if i'd like to use them in the calculation shall I added them in interface sisyphe subroutine !!!! any suggestion please to correct this error ??? thank you in advance win7gfor\telemac2dv6p2.lib C:\opentelemac-mascaret\v6p2\telemac2d\tel2d_v6p2\win
7gfor\sisyphev6p2.lib C:\opentelemac-mascaret\v6p2\telemac2d\tel2d_v6p2\win7gfor
\tomawacv6p2.lib C:\opentelemac-mascaret\v6p2\telemac2d\tel2d_v6p2\win7gfor\bief
v6p2.lib C:\opentelemac-mascaret\v6p2\telemac2d\tel2d_v6p2\win7gfor\damoclesv6p2
.lib C:\opentelemac-mascaret\v6p2\telemac2d\tel2d_v6p2\win7gfor\mumpsvoidv6p2.li
b C:\opentelemac-mascaret\v6p2\telemac2d\tel2d_v6p2\win7gfor\paravoidv6p2.lib C:
\opentelemac-mascaret\v6p2\telemac2d\tel2d_v6p2\win7gfor\specialv6p2.lib -l gfor
tran
Built by Equation Solution <http://www.Equation.com>.
Using built-in specs.
COLLECT_GCC=gfortran
COLLECT_LTO_WRAPPER=c:/opentelemac-mascaret/gcc/bin/../libexec/gcc/x86_64-w64-mi
ngw32/4.7.1/lto-wrapper.exe
Target: x86_64-w64-mingw32
Configured with: ../gcc-4.7.1-mingw/configure --host=x86_64-w64-mingw32 --build=
x86_64-unknown-linux-gnu --target=x86_64-w64-mingw32 --prefix=/home/gfortran/gcc
-home/binary/mingw32/native/x86_64/gcc/4.7.1 --with-sysroot=/home/gfortran/gcc-h
ome/binary/mingw32/cross/x86_64/gcc/4.7.1 --with-gcc --with-gnu-ld --with-gnu-as
--with-gmp=/home/gfortran/gcc-home/binary/mingw32/native/x86_64/gmp --with-mpfr
=/home/gfortran/gcc-home/binary/mingw32/native/x86_64/mpfr --with-mpc=/home/gfor
tran/gcc-home/binary/mingw32/native/x86_64/mpc --with-ppl=/home/gfortran/gcc-hom
e/binary/mingw32/native/x86_64/ppl --with-cloog=/home/gfortran/gcc-home/binary/m
ingw32/native/x86_64/cloog --with-host-libstdcxx='-lstdc++ -lsupc++ -lm' --enabl
e-targets=i686-w64-mingw32,x86_64-w64-mingw32 --enable-cloog-backend=ppl --enabl
e-lto --enable-languages=c,c++,fortran --enable-libgomp --enable-threads=win32 -
-enable-static --enable-shared=lto-plugin --enable-plugins --enable-ld=yes --ena
ble-libquadmath --enable-libquadmath-support --disable-nls --disable-tls --disab
le-win32-registry
Thread model: win32
gcc version 4.7.1 (GCC)
Reading specs from c:/opentelemac-mascaret/gcc/bin/../lib/gcc/x86_64-w64-mingw32
/4.7.1/../../../../lib/libgfortran.spec
rename spec lib to liborig
COLLECT_GCC_OPTIONS='-fconvert=big-endian' '-frecord-marker=4' '-v' '-o' 'flume2
015.exe' '-mtune=generic' '-march=x86-64'
COMPILER_PATH=c:/opentelemac-mascaret/gcc/bin/../libexec/gcc/x86_64-w64-mingw32/
4.7.1/;c:/opentelemac-mascaret/gcc/bin/../libexec/gcc/;c:/opentelemac-mascaret/g
cc/bin/../lib/gcc/x86_64-w64-mingw32/4.7.1/../../../../x86_64-w64-mingw32/bin/
EQ_LIBRARY_PATH=c:/opentelemac-mascaret/gcc/bin/../lib/gcc/x86_64-w64-mingw32/4.
7.1/;c:/opentelemac-mascaret/gcc/bin/../lib/gcc/;c:/opentelemac-mascaret/gcc/bin
/../lib/gcc/x86_64-w64-mingw32/4.7.1/../../../../x86_64-w64-mingw32/lib/../lib/;
c:/opentelemac-mascaret/gcc/bin/../lib/gcc/x86_64-w64-mingw32/4.7.1/../../../../
lib/;c:/opentelemac-mascaret/gcc/bin/../lib/gcc/x86_64-w64-mingw32/4.7.1/../../.
./../x86_64-w64-mingw32/lib/;c:/opentelemac-mascaret/gcc/bin/../lib/gcc/x86_64-w
64-mingw32/4.7.1/../../../
COLLECT_GCC_OPTIONS='-fconvert=big-endian' '-frecord-marker=4' '-v' '-o' 'flume2
015.exe' '-mtune=generic' '-march=x86-64'
c:/opentelemac-mascaret/gcc/bin/../libexec/gcc/x86_64-w64-mingw32/4.7.1/collect
2.exe --sysroot=/home/gfortran/gcc-home/binary/mingw32/cross/x86_64/gcc/4.7.1 -m
i386pep -Bdynamic -o flume2015.exe c:/opentelemac-mascaret/gcc/bin/../lib/gcc/x
86_64-w64-mingw32/4.7.1/../../../../x86_64-w64-mingw32/lib/../lib/crt2.o c:/open
telemac-mascaret/gcc/bin/../lib/gcc/x86_64-w64-mingw32/4.7.1/../../../../x86_64-
w64-mingw32/lib/../lib/crtbegin.o -Lc:/opentelemac-mascaret/gcc/bin/../lib/gcc/x
86_64-w64-mingw32/4.7.1 -Lc:/opentelemac-mascaret/gcc/bin/../lib/gcc -Lc:/opente
lemac-mascaret/gcc/bin/../lib/gcc/x86_64-w64-mingw32/4.7.1/../../../../x86_64-w6
4-mingw32/lib/../lib -Lc:/opentelemac-mascaret/gcc/bin/../lib/gcc/x86_64-w64-min
gw32/4.7.1/../../../../lib -Lc:/opentelemac-mascaret/gcc/bin/../lib/gcc/x86_64-w
64-mingw32/4.7.1/../../../../x86_64-w64-mingw32/lib -Lc:/opentelemac-mascaret/gc
c/bin/../lib/gcc/x86_64-w64-mingw32/4.7.1/../../.. t2dfort.o --start-group C:\op
entelemac-mascaret\v6p2\telemac2d\tel2d_v6p2\win7gfor\telemac2dv6p2.lib C:\opent
elemac-mascaret\v6p2\telemac2d\tel2d_v6p2\win7gfor\sisyphev6p2.lib C:\opentelema
c-mascaret\v6p2\telemac2d\tel2d_v6p2\win7gfor\tomawacv6p2.lib C:\opentelemac-mas
caret\v6p2\telemac2d\tel2d_v6p2\win7gfor\biefv6p2.lib C:\opentelemac-mascaret\v6
p2\telemac2d\tel2d_v6p2\win7gfor\damoclesv6p2.lib C:\opentelemac-mascaret\v6p2\t
elemac2d\tel2d_v6p2\win7gfor\mumpsvoidv6p2.lib C:\opentelemac-mascaret\v6p2\tele
mac2d\tel2d_v6p2\win7gfor\paravoidv6p2.lib C:\opentelemac-mascaret\v6p2\telemac2
d\tel2d_v6p2\win7gfor\specialv6p2.lib -lgfortran -lmingw32 -lgcc -lmoldname -lmi
ngwex -lmsvcrt -lquadmath -lm -lmingw32 -lgcc -lmoldname -lmingwex -lmsvcrt -lad
vapi32 -lshell32 -luser32 -lkernel32 -lmingw32 -lgcc -lmoldname -lmingwex -lmsvc
rt c:/opentelemac-mascaret/gcc/bin/../lib/gcc/x86_64-w64-mingw32/4.7.1/../../../
../x86_64-w64-mingw32/lib/../lib/crtend.o
Running your simulation :
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C:\Users\Desktop\Baptist2\cas.txt_2015-07-30-06h46min27s\out_flume2015.exe
LISTING DE TELEMAC-2D ---------------------------------------------------------
---------------------
TTTTT EEEEE L EEEEE M M AAAAA CCCCC
T E L E MM MM A A C
T EEE L EEE M M M AAAAA C
T E L E M M A A C
T EEEEE LLLLL EEEEE M M A A CCCCC
2D VERSION 6.2 FORTRAN 90
WITH SEVERAL TRACERS
COUPLED WITH SISYPHE AND TOMAWAC
********************************************
* LECDON: *
* APRES APPEL DE DAMOCLES *
* VERIFICATION DES DONNEES LUES *
* SUR LE FICHIER DES PARAMETRES *
********************************************
SORTIE DE LECDON. TITRE DE L'ETUDE :
Couplage avec sisyphe
OUVERTURE DES FICHIERS POUR TELEMAC2D
*****************************
* ALLOCATION DE LA MEMOIRE *
*****************************
READGEO1 : TITRE= C:\Users\KOUKA\Desktop\Baptist2\/canal.MAT
NOMBRE D'ELEMENTS: 3236
NOMBRE REEL DE POINTS: 1786
MXPTEL (BIEF) : NOMBRE MAXIMUM D'ELEMENTS VOISINS D'
UN POINT : 8
NOMBRE MAXIMUM DE POINTS VOISINS D'UN
POINT : 8
CORRXY (BIEF) : PAS DE MODIFICATION DES COORDONNEES
MAILLAGE : MESH ALLOUE
****************************************
* FIN DE L'ALLOCATION DE LA MEMOIRE : *
****************************************
LISTING DE SISYPHE AVEC COUPLAGE-----------------------------------------------
-------------------------------
SSSS I SSSS Y Y PPPP H H EEEEE
S I S Y Y P P H H E
SSS I SSS Y PPPP HHHHH EEEE
S I S Y P H H E
SSSS I SSSS Y P H H EEEEE
VERSION 6.2
COUPLED WITH TELEMAC-2D INTERNALLY
ATTENTION ! A LA LIGNE 2554 DU DICTIONNAIRE :
LE NOMBRE DE VALEURS PAR DEFAUT 10 EST DIFFERENT DE LA TAILLE ANNONC
EE 2
ATTENTION ! A LA LIGNE 2637 DU DICTIONNAIRE :
LE NOMBRE DE VALEURS PAR DEFAUT 10 EST DIFFERENT DE LA TAILLE ANNONC
EE 2
ATTENTION ! A LA LIGNE 2675 DU DICTIONNAIRE :
LE NOMBRE DE VALEURS PAR DEFAUT 10 EST DIFFERENT DE LA TAILLE ANNONC
EE 2
ATTENTION ! A LA LIGNE 3357 DU DICTIONNAIRE :
LE NOMBRE DE VALEURS PAR DEFAUT 13 EST DIFFERENT DE LA TAILLE ANNONC
EE 1
********************************************
* LECDON: *
* APRES APPEL DE DAMOCLES *
* VERIFICATION DES DONNEES LUES *
* SUR LE FICHIER DES PARAMETRES *
********************************************
OUVERTURE DES FICHIERS POUR SISYPHE
*******************************
* ALLOCATION DE LA MEMOIRE *
*******************************
READGEO1 : TITRE= C:\Users\KOUKA\Desktop\Baptist2\/canal.MAT
NOMBRE D'ELEMENTS: 3236
NOMBRE REEL DE POINTS: 1786
MXPTEL (BIEF) : NOMBRE MAXIMUM D'ELEMENTS VOISINS D'
UN POINT : 8
NOMBRE MAXIMUM DE POINTS VOISINS D'UN
POINT : 8
CORRXY (BIEF) : PAS DE MODIFICATION DES COORDONNEES
MAILLAGE : MESH_S ALLOUE
****************************************
* FIN DE L'ALLOCATION DE LA MEMOIRE : *
****************************************
CALLING LECLIM
BACK FROM LECLIM
CALLING INBIEF
INBIEF (BIEF) : MACHINE NON VECTORIELLE (SELON VOS DONNEES)
BACK FROM INBIEF
STRCHE (BIEF) : PAS DE MODIFICATION DU FROTTEMENT
CALLING FRICTION_CHOICE
BACK FROM FRICTION_CHOICE
CALLING FRONT2
IL Y A 2 FRONTIERE(S) LIQUIDE(S) :
FRONTIERE 1 :
DEBUT AU POINT DE BORD 160 , DE NUMERO GLOBAL 11
ET DE COORDONNEES : 15.85000 0.000000
FIN AU POINT DE BORD 168 , DE NUMERO GLOBAL 13
ET DE COORDONNEES : 15.85000 0.8000000
FRONTIERE 2 :
DEBUT AU POINT DE BORD 327 , DE NUMERO GLOBAL 15
ET DE COORDONNEES : 0.000000 0.8000000
FIN AU POINT DE BORD 1 , DE NUMERO GLOBAL 2
ET DE COORDONNEES : 0.000000 0.000000
IL Y A 2 FRONTIERE(S) SOLIDE(S) :
FRONTIERE 1 :
DEBUT AU POINT DE BORD 1 , DE NUMERO GLOBAL 2
ET DE COORDONNEES : 0.000000 0.000000
FIN AU POINT DE BORD 160 , DE NUMERO GLOBAL 11
ET DE COORDONNEES : 15.85000 0.000000
FRONTIERE 2 :
DEBUT AU POINT DE BORD 168 , DE NUMERO GLOBAL 13
ET DE COORDONNEES : 15.85000 0.8000000
FIN AU POINT DE BORD 327 , DE NUMERO GLOBAL 15
ET DE COORDONNEES : 0.000000 0.8000000
BACK FROM FRONT2
CALLING CORNOR
BACK FROM CORNOR
CALLING MASBAS2D
BACK FROM MASBAS2D
CALLING CORFON
CORFON (TELEMAC2D) : PAS DE MODIFICATION DU FOND
BACK FROM CORFON
CALLING CONDIN
BACK FROM CONDIN
CALLING PROPIN
BACK FROM PROPIN
CALLING FRICTION_CHOICE
BACK FROM FRICTION_CHOICE
================================================================================
ITERATION 0 TEMPS : 0.0000 S
CALLING PRERES_TELEMAC2D
BACK FROM PRERES_TELEMAC2D
CALLING DESIMP
BACK FROM DESIMP
CALLING HPROPA
BACK FROM HPROPA
CALLING DE CHPCON
BACK FROM CHPCON
CALLING BILAN
BACK FROM BILAN
TELEMAC2D COUPLE AVEC : INTER-SISYPHE
TELEMAC-2D : COUPLAGE INTERNE AVEC SISYPHE
PREMIER APPEL DE SISYPHE
PART 0 : INITIALISING SISYPHE
INBIEF (BIEF) : MACHINE NON VECTORIELLE (SELON VOS DONNEES)
IL Y A 2 FRONTIERE(S) LIQUIDE(S) :
FRONTIERE 1 :
DEBUT AU POINT DE BORD 160 , DE NUMERO GLOBAL 11
ET DE COORDONNEES : 15.85000 0.000000
FIN AU POINT DE BORD 168 , DE NUMERO GLOBAL 13
ET DE COORDONNEES : 15.85000 0.8000000
FRONTIERE 2 :
DEBUT AU POINT DE BORD 327 , DE NUMERO GLOBAL 15
ET DE COORDONNEES : 0.000000 0.8000000
FIN AU POINT DE BORD 1 , DE NUMERO GLOBAL 2
ET DE COORDONNEES : 0.000000 0.000000
IL Y A 2 FRONTIERE(S) SOLIDE(S) :
FRONTIERE 1 :
DEBUT AU POINT DE BORD 1 , DE NUMERO GLOBAL 2
ET DE COORDONNEES : 0.000000 0.000000
FIN AU POINT DE BORD 160 , DE NUMERO GLOBAL 11
ET DE COORDONNEES : 15.85000 0.000000
FRONTIERE 2 :
DEBUT AU POINT DE BORD 168 , DE NUMERO GLOBAL 13
ET DE COORDONNEES : 15.85000 0.8000000
FIN AU POINT DE BORD 327 , DE NUMERO GLOBAL 15
ET DE COORDONNEES : 0.000000 0.8000000
STRCHE (BIEF) : PAS DE MODIFICATION DU FROTTEMENT
INITIALISATION EN CAS DE COUPLAGE : PART= 0
_____________
runcode::main:
C:\Users\morri\Desktop\Baptist2:
|runCAS: fail to run
| C:\Users\Desktop\Baptist2\cas.txt_2015-07-30-06h46min27s\out_flu
me2015.exe |
|
The administrator has disabled public write access.
|
error in the simulation telemac2D 9 years 3 months ago #17817
|
Hello,
This last line in the listing is not an error, just a message, the problem is that probably the program gets stuck somewhere after this, so please post your modified subroutine so that we can see if their is a mistake. Regards, JMH |
|
The administrator has disabled public write access.
|
|
Moderators: pham
