Coefficients de marée.

oui, j'ai toujours entendu dire que c'est franco-français !
pas bête d'ailleurs le coup du coeff ...
par contre pourquoi de 0 à 120 ?...

oups ... :oups: oui, coquille !

coefficient-coefficient oui c'est seulement français, les "autres" utilisent les termes springs/neaps qui correspondent à Vives Eaux/Mortes Eaux, mais il manque la caractérisation numérique d'une marée en particulier. On dit "springs" quand en France on pourrait dire 95-90-100 etc
.
Pour faire un exemple, pour les interpolations pour les courants de marée, en France on a le graphique dans lequel on entre avec le coefficient et on interpole entre les deux valeurs pour 45 e 95, dans les autres pays on utilise l'amplitude du jour par rapport à l'amplitude standard de springs et neaps.
.
On remarque par contre parfois une utilisation pas correcte du coefficient: dans un certain nombre de ports on vous demande le tirant d'eau puis on vous dit "aujourd'hui on a 90 donc là bas vous n'allez pas toucher" ce qui n'est pas sûr à 100% car on peut avoir "coeff 90", disons 4m de marnage, soit avec PM 6.0m BM 2m, soit avec PM 5.5m et BM 1.5m, ce qui pour la semelle de la quille n'est pas la même chose..

Oui.
Si j'ai bien capté, le marnage étant la hauteur d'eau calculée en référence à la BM précédente ?

@ Yves (...) : Mais pas le même marnage

@ Yves (...) : le coefficient entre 2 PM ? Tu additionnes les 2 coeff et tu divises par 2
:-D

@Yves :

En fait pourquoi 20 à 120 plutôt que 0 à 100.

Si j'ai bien compris, cela date d'il y a longtemps. D'une époque où on ne savait calculer avec certitude ni la plus basse marée basse, ni la plus haute marée haute. Parce que les calculs était basé sur des mesures faites à Brest, si mes souvenirs sont bons.

Alors, par précaution, on a décalé l'échelle de 20 vers le haut. Comme ça, peut être un jour y aura t il un coefficient de 19 ou un coefficient de 121 ... Qui sait.

Effectivement, ici on n'utilise pas le coefficient de marée.

Les tables donnent la hauteur et l'heure pour les ports de références. Quant aux ports secondaires elles donnent la différence en hauteur et en temps par rapport au port de référence.

C'est là que ça se complique dans ma tête.

Intuitivement, je n'ai jamais été attiré par le Golfe du Mexique et c'est sans doute un début d'explication. :reflechi:

Dans mon bled, les coefficients ne me sont pas d une grande utilité. Pour y naviguer, ou plutôt pour pouvoir entre et sortir, ce qui m intéresse c est la hauteur d eau et rien que la hauteur d eau. Le coef, j ai beau chercher, pour mon "utilisation" je vois pas ce que je peux en faire sauf la lecture des cartes de courant simplifiée.

Galveston !
Un effroyable épisode de submersion pendant un ouragan en 1900, bien pire que la Nouvelle Orléans pendant Katrina
fr.wikipedia.org[...]lveston

7mètres??? Hubert, c est ridicule. Viens me voir a la maison. Tu boiras un bon coup et tu verras une VRAIE amplitude de maree ;-)

@ Hubert, de Cherbourg : mais voir surtout la hauteur de marnage

Même chose: valeurs pour les 'neaps' (mortes-eaux) et 'springs' (vive-eaux) dans les pays anglo-saxons. De toute façon, l'extrapolation est peu fiable, juste indicative...

Veronique: Je navigue en permanence avec trois cartes de courants différentes et, en cas de doute, je consulte une autre carte. C'est pourquoi, j'affirme qu'il existe des différences importantes. Consulte mon blog et tu comprendras: karibario.blogspot.com .

Cela dit, le meilleur indicateur reste les bouées de casiers que les pécheurs mettent sur notre chemin...

En croisière, je suis nettement plus cool et me contente de jeter un oeil sur les tendances du courant et les passages à niveau .

@ viking35 : imprécises ? peut être mais suffisamment précises pour nous indiquer où on va; indispensables pour faire une estime qui se révèle souvent pas trop mauvaise donc cartes de courant pas si imprécises que ça

45 = mortes eaux moyennes = neap; 95 = vives eaux moyennes = spring

Moi je m'en fous j'ai tout qui est indiqué directement sur mes cartes Navionics sous Android la la lère. Il n'y a qu'à lire et c'est pratique quand je dois prendre un mouillage j'ai sous les yeux la hauteur mini maxi et la hauteur actuellement.
Il faut vivre avec son temps. :langue2:

peut être l’inventeur des coefficients

fr.wikipedia.org[...]Laplace

www.meteolafleche.com[...]15.html

:litjournal:

Oui, le vent et la pression atmosphérique, sans oublier la configuration des lieux.

Tout à fait d'accord, mais le vent peut avoir deux conséquences :
- S'il souffle dans une direction constante pendant plusieurs jours, il pourra renforcer au affaiblir le courant.
- sur une durée éventuellement plus courte, le vent qui souffle du large vers la côte fera monter le niveau de la mer et l'inverse. Si on ajoute "algébriquement" l'effet de la pression atmosphérique, ce sont autant de corrections à apporter à la hauteur d'eau calculée. Chez nous, cela joue également pour les heures d'ouverture/fermeture de portes.

Je vais essayer d'être un peu didactique au sujet de cette histoire d'onde harmoniques.

Pour calculer la marée, on empile plusieurs courbes sinusoïdales. C'est un peu comme le timbre qui caractérise un son. La fréquence la plus basse du set de courbes sinusoïdales va définir la note, par exemple un do. Les autres ondes de fréquence plus élevées vont caractériser le timbre, c'es en gros ce qui permet de distinguer le son d'une flute de celui d'une trompette.

Pour la marée, on utilise un peu le même type d'analyse. On aura un set de fonctions sinusoïdales variant une fois par jour, ce seront les ondes diurnes et elles sont prépondérantes dans les ports à marée diurnes. D'autres ondes auront un variation semi diurnes. Ce sont les plus importantes sur le littoral Français métropolitain Manche/Atlantique. Elles déreminent des marées semi diurnes avec deux PM et deux BM par jour. Pour avoir unrésultat précis, on ajoute aussi des ondes quart diurnes, sixième diurnes, etc.

Pour mémoire je met en image le marégramme pour Lorient ce jour calculé avec Marées dans le Monde, ainsi que les coubes obtenues en ne tenant compte que des ondes diurnes, semi diurnes et quart diurnes.

Comme on le voit, les ondes semi diurnes sont de loin prépondérantes, alors que les ondes diurnes sont très faibles, c'est bien sur la caractéristique des ports bretons, mais là, au voisinage de l'équinoxe et de ses grandes marées, le phénomèe est amplifié. C'est pourquoi les grandes marées d'équinoxe sont grandes en Europe, mais dans les régions à marées diurnes, c'est plutôt au voisinage des solstices qu'on observe de grands marnages.

Puisque chez nous ce sont les ondes semi diurnes qui comptent le plus, ce sont elles qui décrivent le mieux l'importance de la marée, d'ou la notion de coefficient qui varie de 20 à 120 quand la marée va de son minimum à son maximum théorique.

Bien sur ce coef n'a un sens que si on est dans un port à marée semi diurne et le SHOM ne le calcule que pour les ports atlantiques de la France métropolitaine. C'est une approximation, mais elle est bien pratique. Mon chantier sait bien qu'il peut gruter mon bateau si le coef est supérieur à 70. Sinon, il n'y aura pas assez d'eau sur la cale.

Une remarque : le coef décrit l'amplitude de la PLEINE MER réduite aux ondes semi diurnes. Il ne s'applique pas aux basses mer.

Enfin, on sait que les vives eaux sont plus importantes quand la Lune et le Soleil génèrent des ondes de marée synchrone, en gros autour de la nouvelle et de la pleine Lune. Il ne faut pas confondre avec le phénomène des marées d'équinoxe. Ce qui fait l'importance des marées d'équinoxe, c'est qu'à ce moment le Soleil est sur le plan de l'équateur. Les astronomes disent que sa déclinaison est nulle. Un des aspects extraordinaires de la marée du siècle du 21 mars prochain, c'est que la Lune elle même va traverser le plan de l'équateur dans les jours précédant les grandes marées.

Et pour moi, ca flirte avec mes limites intellectuelles! (Facilement atteintes, je le concède....)

Ci dessous, une petite animation flash pour les gars "limités", comme Yves !
.
La plupart des phénomènes périodiques sont le fruit de la superposition d'ondes de fréquences multiples de la fondamentale, appelée harmonique de rang 1 (somme le son, décrit par gonio-coque)
La décomposition spectrale dite de Fourrier, permet de déterminer l'amplitude et le déphasage de chaque harmonique.
.
Sur cette animation simple (sans déphasage), amusez-vous à faire varier le curseur d'amplitude des harmoniques et regardez l'onde obtenue.
.
Pour la marée en Bretagne, baissez l’harmonique de rang 1 (marée diurne) à une valeur de l'ordre de 0,2 et mettez celle de rang 2 (semi diurne) à 1. Vous retrouverez l'allure de la courbe de marée de gonio-coque.
.
PS : le premier qui réussit à faire un carré et un triangle sans utiliser la solution décrite à droite, a gagné toute ma considération !
.
www.ostralo.net[...]ues.swf

PS 1 - On n'utilise pas les séries de Fourier pour calculer les coefficients des marées mais un raisonnement bien plus astronomique incluant les sommes et les différences des harmoniques des 'pulsations' lunaires et solaires (qui sont nombreuses)
En voici une toute petite collection alors que le SHOM en utilise maintenant 120 !
PS 2 Un véritable jeux avec les harmoniques devraient inclure les phases (qui sont toutes égales à zéro dans le lien ci-dessus)

Ch'u allé ici,
www.cpge.eu[...]ier.pdf

puis là !
f5zv.pagesperso-orange.fr[...]B04.htm

En avant première,

Je ne peux résister à l'envie de vous faire partager le fruit de mes élucubrations de l'année passée. J'ai passé en temps cumulé 13,5 jours à bord sur mon catway en 2014 à Locmiquélic en rade de Lorient. J'en ai profité pour monitorer la profondeur en branchant le sondeur sur l'ordi via un port COM.

J'ai extrait les constantes harmoniques en utilisant la méthode de calcul exploitée par Marées dans le Monde. Bien sur je néglige l'onde annuelle Sa, la seule à longue période prise en compte par le soft, voici ce que cela donne. Et on ne peut pas moyenner les effets météorologiques sur une aussi courte observation. La courbe en rouge est obenue à partir des données du SHOM, la bleue c'est la mienne.

Le décalage de 30 minutes, c'est ce que j'ai tout le temps avec cette méthode. pour mémoire, j'obtient une valeur de Zo de 304 cm, quand le SHOM la fixe à 309 cm. Ca fait une différence de 5 cm sur le zéro des cartes.

Je ne prends en compte qu'une vingtaine d'harmoniques. avec 120 ondes élémentaires, il faudrait au moins 3 mois pour sortir une solution.

L'idée, ce serait de diffuser ce module de calcul avec la prochaine versiion de Marées dans le Monde afin de constituer une base de donnée contributive de constantes harmoniques libre de droit, mais je ne sais pas trop si ce projet est vraiment réaliste.

Répondre à ta question n'est pas évident parce que j'utilise une formule de calcul simplifiée éditée par le SHOM dans les années 70 pour les calculatrices programmables.

En réalité, on ne fournit que 10 ondes élémentaires : Sa, Q1, O1, K1, N2, M2, S2, Mn4, M4 et Ms4. 11 autres sont déduites des premières. Je calcule une approximation pour 2N2, Nu2, Mu2, L2, T2, K2, P1. Je calcule aussi les corrections nodales pour M2, O2, K2 et K1.

Par exemple T2 est obtenu en prenant une amplitude de1/17 X celle de S2 et une situation à 273° de celle de S2. Donc si on bidouille S2, on agit aussi sur T2.

Mais bon, si tu veux faire joujou, tu n'as qu'a rentrer dans Marées dans le Monde les valeurs que je te mets en PJ, ce sont celles que j'ai calculées après deux semaines d'observation en rade de Lorient. Tu pouras alors les modifier pour voir ce que cela donne.

Ces données sont un peu différentes des vieilles données numérisées par Jean-Pierre LaPointe à partir de la 'Table des Marées des Grands Ports du Monde', Ouvrage No. 540 - Service Hydrographique et Océanographique de la Marine...

rel_pedigree: 1250000000 (Via LaPointe)

!units: meters

!longitude: -3.3500

!latitude: 47.7500

Lorient, France
01:00 :Europe/Paris
3.0000 meters
J1 0.0050 161.00
K1 0.0500 103.00
K2 0.1490 165.70
L2 0.0450 151.00
M1 0.0050 44.50
M2 1.5210 130.00
M4 0.0790 87.00
N2 0.3140 109.00
2N2 0.0420 88.00
O1 0.0650 346.00
OO1 0.0030 220.00
P1 0.0170 94.20
Q1 0.0220 298.00
2Q1 0.0020 -130.10
R2 0.0040 164.30
S2 0.5470 163.00
T2 0.0320 161.70
LDA2 0.0110 145.30
MU2 0.0370 97.00
NU2 0.0610 111.80
RHO1 0.0020 -64.20
MN4 0.0400 36.00
MS4 0.0220 183.00
SA 0.0460 209.00

Jamais de FFT (transformée de Fourier Rapide) sur les données des marées. Beaucoup de composantes ont des périodes tellement proches qu'il faudrait des enregistrement TTTRRRREEEESSS longs pour les séparer. La bonne vieille méthode des astronomes depuis Laplace est la meilleure, avec ajustement aux moindres carrées. Mais les calculs restent LOURDS !

Voici les résultats pours ondes suivantes.

On n'a rien avant 3,5 jours, car le soft ne fournit pas de solution tant que l'acquisition est insuffisante.

Intuitivement, il me semble que pour améliorer un peu la précision au delà de 15jours d'observation, il faudrait énormément prolonger la période de mesure, mais ça je le vérifierai en avril au moment de la remise à l'eau.

Des pollutions de fil comme ça, on en voudrait tous les jours !

Merci pour toutes ces explications détaillées et les liens que le profane que je suis va "essayer" de digérer.

Pour François,

Les données de Mérak semblent provenir d'un crack du logiciel "Marées". Elles sont assez ancienne quand même, environ 20 ans mais je n'ose pas les inclure dans Marées dans le Monde car j'ai un doute quant à la légalité de la source. C'est dommage car elles sont géniales pour la France Métro et l'Europe.

Je te mets en PJ les hamoniques que j'ai extraite à partir des hauteurs d'eau publiées par le SHOM (un an d'observation, pas de 1 heure. Comme tu vois, y a pas épais comme différence.

Pour Vincent,

Pour le décalage, je n'ai pas d'explication.
Ca ne semble pas tenir du lieu puisque les marées calculées par le SHOM dans le coin que ce soit à Hennebont au fond du Blavet, à Lorient, à Port Louis ou à Groix sont toutes synchrones. Dans ces conditions, je ne vois pas pourquoi le resultat serait différent à Locmiquélic.
Je me suis demandé si c'est du aux frictions sur le catway, mais là non plus ça ne tient pas, puisque je suis en avance et non pas en retard par rapport au SHOM.
Ce ne semble pas non plus être lié aux algorithmes de calcul parce que si j'observe les étales en suivant mon sondeur mes calculs semblent coller avec ce que je vois.

Il me semble que tu bosse avec les solutions FES si tu exploite 16 ondes. Marées dans le Monde est plus rusé. La formule de calcul est une cuisine utilisant des approximations qui au final avec seulement 10 ondes en approche 10 autres pour un résultat correct quel que soit le type de marée. C'est du au cerveau des gars du SOM et je trouve que c'est brillant.

La suite de mon post est un copié collé de la routine de calcul que j'utilise, essaye de voir si tu peux récupérer ça pour ton usage perso. Tu verra que c'est assez différent des formulations plus orthodoxes.

Sub CalculeHauteur()

' ***** L'ENTREE SE FAIT DANS LE FUSEAU DES TABLES

‘JourJulienTable désigne le jour julien compté à partir du premier janvier 1900.5 exprimé dans le fuseau des tables de marées
‘T désigne le temps en siècle juliens depuis le premier janvier 1900
‘TempsParHeure indique l’heure du calcul exprimé dans le fuseau des tables
‘AML = Anomalie moyenne de la Lune
‘LML = Longitude moyenne de la Lune
‘LMS = Longitude moyenne du Soleil
‘AMS = Anomalie moyenne du Soleil
‘LNL = Longitude du nœud ascendant de la Lune
‘LPS = Longitude du périgée du Soleil
‘LPL = Longitude du périgée de la Lune
‘HauteurdEau = Hauteur d’eau exprimée en mm
‘Zo = Zo en cm
‘Asa, Ak1, Ao1, Aq1, Am2, As2, Aq2, Amn4 etc désignent l’amplitude en mm des ondes Sa, K1, O1, Q1, S2, MN4, etc.
‘De même Gsa, Gk1 Go1, etc désignent la phase des ondes Sa, K1, O1, etc exprimées en degrés – une heure de temps équivaut à 15°

Dim n(1 To 6) As Integer

'ANOMALIES, ELONGATIONS, NOEUDS
'Longitude moyenne de la Lune
LML = 270.434164 + 481267.8831 * T - 0.001133 * T * T + 0.0000019 * T * T * T

'Anomalie Moyenne du Soleil
AML = 296.104608 + 477198.8491 * T + 0.009192 * T * T + 0.0000144 * T * T * T

'Longitude du noeud ascendant de la Lune
LNL = 259.183275 - 1934.142 * T + 0.002078 * T * T + 0.0000022 * T * T * T

'Longitude moyenne du Soleil
LMS = 279.69668 + 36000.76892 * T + 0.0003025 * T * T

'Anomalie moyenne du Soleil
AMS = 358.47583 + 35999.04975 * T - 0.00015 * T * T - 0.0000033 * T * T * T

'Corrections de ces valeurs, termes additifs
LML = LML + 0.000233 * Sin((51.2 + 20.2 * T) * Pi / 180)
AMS = AMS - 0.001778 * Sin((51.2 + 20.2 * T) * Pi / 180)
AML = AML + 0.000817 * Sin((51.2 + 20.2 * T) * Pi / 180)
Double1 = 0.003964 * Sin((346.56 + 132.87 * T - 0.0091731 * T * T) * Pi / 180)
LML = LML + Double1
AML = AML + Double1
LML = LML + 0.001964 * Sin(LNL * Pi / 180)
AML = AML + 0.002541 * Sin(LNL * Pi / 180)

'Calcul de la longitude des périgées
LPL = LML - AML
LPS = LMS - AMS

'Ramener les résultats entre 0 et 360°
LML = LML - 360 * (LML \ 360)
LNL = -LNL + 360 * (LNL \ 360)
LMS = LMS - 360 * (LMS \ 360)
LPL = LPL - 360 * (LPL \ 360)
LPS = LPS - 360 * (LPS \ 360)

'transformation en radiants
LML = LML * Pi / 180
LMS = LMS * Pi / 180
LNL = LNL * Pi / 180
LPL = LPL * Pi / 180
LPS = LPS * Pi / 180

' ***** Calcul de la hauteur d'eau *****

'**niveau de référence**

HauteurdEau = Zo * 10

'**Sa**

j = 0
n(1) = 0
n(2) = 1
n(3) = 0
n(4) = 0
n(5) = 0
n(6) = 0

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + Asa * Cos(v - Gsa * Pi / 180)

' ***** K1 *****

j = 1
n(1) = 0
n(2) = 1
n(3) = 0
n(4) = 0
n(5) = 0
n(6) = 1

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + Ak1 * Cos(v - Gk1 * Pi / 180)

' ***** O1 *****

j = 1
n(1) = -2
n(2) = 1
n(3) = 0
n(4) = 0
n(5) = 0
n(6) = -1

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + Ao1 * Cos(v - Go1 * Pi / 180)

' ***** Q1 *****

j = 1
n(1) = -3
n(2) = 1
n(3) = 1
n(4) = 0
n(5) = 0
n(6) = -1

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + Aq1 * Cos(v - Gq1 * Pi / 180)

' ***** P1 *****

j = 1
n(1) = 0
n(2) = -1
n(3) = 0
n(4) = 0
n(5) = 0
n(6) = 1

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau - 1 / 3 * Ak1 * Cos(v - Gk1 * Pi / 180)

' ***** o1 (correction nodale de O1) *****

j = 1
n(1) = -2
n(2) = 1
n(3) = 0
n(4) = -1
n(5) = 0
n(6) = -1

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + 1 / 5.3 * Ao1 * Cos(v - Go1 * Pi / 180)

' ***** k1 (correction nodale de K1) *****

j = 1
n(1) = 0
n(2) = 1
n(3) = 0
n(4) = 1
n(5) = 0
n(6) = 1

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + 1 / 7.4 * Ak1 * Cos(v - Gk1 * Pi / 180)

' ***** M2 *****

j = 2
n(1) = -2
n(2) = 2
n(3) = 0
n(4) = 0
n(5) = 0
n(6) = 0

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + Am2 * Cos(v - Gm2 * Pi / 180)

' ***** N2 *****

j = 2
n(1) = -3
n(2) = 2
n(3) = 1
n(4) = 0
n(5) = 0
n(6) = 0

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + An2 * Cos(v - Gn2 * Pi / 180)

' ***** S2 *****

j = 2
n(1) = 0
n(2) = 0
n(3) = 0
n(4) = 0
n(5) = 0
n(6) = 0

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + As2 * Cos(v - Gs2 * Pi / 180)

' ***** 2 N2 *****

j = 2
n(1) = -4
n(2) = 2
n(3) = 2
n(4) = 0
n(5) = 0
n(6) = 0

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + 1 / 7.6 * An2 * Cos(v - Gn2 * Pi / 180)

' ***** mu2 *****

j = 2
n(1) = -4
n(2) = 4
n(3) = 0
n(4) = 0
n(5) = 0
n(6) = 0

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + 1 / 6.3 * An2 * Cos(v - Gn2 * Pi / 180)

' ***** v2 *****

j = 2
n(1) = -3
n(2) = 4
n(3) = -1
n(4) = 0
n(5) = 0
n(6) = 0

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + 1 / 5.3 * An2 * Cos(v - Gn2 * Pi / 180)

' ***** L2 *****

j = 2
n(1) = -1
n(2) = 2
n(3) = -1
n(4) = 0
n(5) = 0
n(6) = 0

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau - 1 / 35 * Am2 * Cos(v - Gm2 * Pi / 180)

' ***** K2 *****

j = 2
n(1) = 0
n(2) = 2
n(3) = 0
n(4) = 0
n(5) = 0
n(6) = 0

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + 1 / 3.7 * As2 * Cos(v - Gs2 * Pi / 180)

' ***** T2 *****

j = 2
n(1) = 0
n(2) = -1
n(3) = 0
n(4) = 0
n(5) = 0
n(6) = 0

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + 1 / 17 * As2 * Cos(v - Gs2 * Pi / 180 + 283 * Pi / 180)

' ***** m2 (correction nodale de M2) *****

j = 2
n(1) = -2
n(2) = 2
n(3) = 0
n(4) = -1
n(5) = 0
n(6) = 0

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau - 1 / 27 * Am2 * Cos(v - Gm2 * Pi / 180)

' ***** k2 (Correction nodale de K2) *****

j = 2
n(1) = 0
n(2) = 2
n(3) = 0
n(4) = 1
n(5) = 0
n(6) = 0

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + 1 / 12 * As2 * Cos(v - Gs2 * Pi / 180)

' ***** MN4 *****

j = 4
n(1) = -5
n(2) = 4
n(3) = 1
n(4) = 0
n(5) = 0
n(6) = 0

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + Amn4 * Cos(v - Gmn4 * Pi / 180)

' ***** M4 *****

j = 4
n(1) = -4
n(2) = 4
n(3) = 0
n(4) = 0
n(5) = 0
n(6) = 0

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + Am4 * Cos(v - Gm4 * Pi / 180)

' ***** MS4 *****

j = 4
n(1) = -2
n(2) = 2
n(3) = 0
n(4) = 0
n(5) = 0
n(6) = 0

v = 15 * j * TempsParHeure * Pi / 180
v = v + n(1) * LML + n(2) * LMS + n(3) * LPL
v = v + n(4) * LNL + n(5) * LPS + n(6) * Pi / 2

HauteurdEau = HauteurdEau + Ams4 * Cos(v - Gms4 * Pi / 180)

End Sub

Pour ce qui est des harmoniques du plug in de OpenCpn, je me tâte vraiment pour les utiliser.. Après tout OpenCpn est en open source, donc pourquoi ne pas les récupérer? j'hésite encore.

Je ne crois pas avoir commis d'erreur dans l'écriture des algorithme d'extraction. Pour preuve, voici ce que j'obtiens pour la courbe de ce jour en exploitant les harmoniques calculées à partir des hauteur d'eau fournies par le SHOM (à une correction minime de Zo près). On est très près de ce que fournissent OpenCpn et le SHOM sur son site.