Simulations Numériques

1. Présentation du logiciel TELEMAC

2. Maillage utilisé

C-----------------------------------------------------------------------
C variable declarees localement par l'utilisateur
C lecture des donnees du SHOM
rewind(26) DO 111 i=1,114
DO 112 j=0,10
READ(26,150) ch(i,j,1),ch(i,j,2)
112 CONTINUE
111 CONTINUE
150 format(f8.1,2X,f6.2)
c------------------------------------------------------
c CONDITIONS LIMITES EN NIVEAU
c calcul simplifié pour M2,N2,S2,M4
c le temps T0 correspond à la PM pour une marée de vive eau exeptionnelle
t0=17100.
c------------------------------------------------------
c calcul pour tous les points frontiere de xe
c ( surface libre / au 0 cote marine )
c------------------------------------------------------
DO 113 i=1,114
c Niveau moyen pris constant sur la zone (NM=2.9m)
haut=2900.
c Onde N2: ch(i,5,1et2)
haut=haut+ch(i,5,1)*cosd(0.0078999*(TEMPS+t0)-ch(i,5,2))
c Onde M2: ch(i,6,1et2)
haut=haut+ch(i,6,1)*cosd(0.0080512*(TEMPS+t0)-ch(i,6,2))
c Onde S2: ch(i,7,1et2)
haut=haut+ch(i,7,1)*cosd(0.0083333*(TEMPS+t0)-ch(i,7,2))
c Onde M4: ch(i,9,1et2)
haut=haut+ch(i,9,1)*cosd(0.0161024*(TEMPS+t0)-ch(i,9,2))
xe(i)=haut/1000.
113 CONTINUE
C BOUCLE SUR TOUS LES POINTS FRONTIERE
DO 5 K=1,NPTFR
C Test concernant le type de frontiére liquide. On verifie que
C l'on a bien une condition de type Hauteur imposée
C ATTENTION C les points xe(1 à 35) correspondent à la frontiere Sud de L'Ouest vers l'Est
C soit aux points de num locale 2 à 36 C les points xe(36 à 113) correspondent à la frontiere Ouest du Nord au Sud
Csoit aux points de num locale 450 à 527
C le point xe(114) correspond au coin SO de num locale 1
C
IF(K.GE.2 .AND. K.LE.36) THEN
LIHBOR(K)=5
HBOR(K) = max(0.D0 , (xe(K-1)-ZF(NBOR(K))))
H(NBOR(K))=HBOR(K)
ENDIF
IF(K.GE.2 .AND. K.LE.36) THEN
LIHBOR(K)=5
HBOR(K) = max(0.D0 , (xe(K-1)-ZF(K)))
H(K)=HBOR(K)
ENDIF
IF(K.GE.450 .AND. K.LE.527) THEN
LIHBOR(K)=5
HBOR(K) = max(0.D0 , (xe(K-414)-ZF(K)))
H(K)=HBOR(K)
ENDIF
IF(K.EQ.1) THEN
LIHBOR(K)=5
HBOR(K) = max(0.D0 , (xe(114)-ZF(K)))
H(K)=HBOR(K)
ENDIF