Note

Scripts disponibles sous /home/flamingo/vgarnier/TOOLS/VACUMM/PERSO/TTB

Validation de la marée : démarche générale

Composantes harmoniques

Pour lister les composantes harmoniques du potentiel générateur (par amplitude décroissante ou période décroissante):

ttb_showarg.py --amplitude --show

Sur une période de run limitée, quelles composantes harmoniques puis-je extraire?

ttb_get_wave_to_be_analyzed.py --run_duration 1704. --waves M2

>>>
Sorted by decreasing amplitudes, the list of waves is :  M2 K1 S2 O1 P1 N2 K2 Q1 M1 Nu2 J1 L2 T2 2N2 Mu2 OO1 Ro1 Pi1 2Q1 Phi1 E2 La2 KJ2 Sig1 Ki1 Psi1 Tta1 R2 MP1 MSqm KQ1 Mqm 3OK1 MS1 S1 SO1 2MN2S2 2NS2 3M2S2 ST1 OQ2 MNS2 MNuS2 ST2 ST3 2MK2 SNK2 NKM2 MSN2 2SM2 SKM2 2SN2 2SMu2 MQ3 2MK3 MO3 M3 SO3 MS3 MK3 SP3 S3 SK3 K3 2MNS4 N4 3MS4 MN4 MNu4 M4 SN4 ML4 NK4 MS4 MK4 2MSN4 S4 SK4 3MNK6 3MNS6 4MK6 3MNL6 4MS6 2MN6 2MNu6 3MSK6 M6 3MKS6 MSN6 2ML6 MSNu6 MNK6 2MS6 3MLN6 2MK6 MSL6 3MSN6 3MKN6 2SM6 MSK6 2M2N8 3MN8 M8 2MSN8 3ML8 3MS8 3MK8 2M2S8 2MSK8 2M2K8

Extraction de l atlas (1ère étape)

Dans un premier temps, extraire toutes les ondes proposées en sortie de ttb_get_list_wave.py:

comodo-detidor -l list_file_IROI.txt --variable-mask --only-atlases -a -p ./ALL -v XE -d M2 K1 S2 O1 P1 N2 K2 Q1 M1 Nu2 J1 L2 T2 2N2 Mu2 OO1 Ro1 Pi1 2Q1 Phi1 E2 La2 KJ2 Sig1 Ki1 Psi1 Tta1 R2 MP1 MSqm KQ1 Mqm 3OK1 MS1 S1 SO1 2MN2S2 2NS2 3M2S2 ST1 OQ2 MNS2 MNuS2 ST2 ST3 2MK2 SNK2 NKM2 MSN2 2SM2 SKM2 2SN2 2SMu2 MQ3 2MK3 MO3 M3 SO3 MS3 MK3 SP3 S3 SK3 K3 2MNS4 N4 3MS4 MN4 MNu4 M4 SN4 ML4 NK4 MS4 MK4 2MSN4 S4 SK4 3MNK6 3MNS6 4MK6 3MNL6 4MS6 2MN6 2MNu6 3MSK6 M6 3MKS6 MSN6 2ML6 MSNu6 MNK6 2MS6 3MLN6 2MK6 MSL6 3MSN6 3MKN6 2SM6 MSK6 2M2N8 3MN8 M8 2MSN8 3ML8 3MS8 3MK8 2M2S8 2MSK8 2M2K8

Analyse intermédiaire

Visualisation de l’amplitude de chaque composante à proximité des maragraphes disponibles dans la zone, et comparaison avec XXX:

XXXXX

Voir avec Florent le defaut de --wwave dans tidegauge_harmcomp.py

Extraction de l atlas (étape finale)

Création de l’atlas en éliminant les composantes d’amplitude trop faible, et en ordonnant les ondes par ordre d’amplitude décroissant. Finalement, on analyse les ondes astronomiques forcées aux conditions aux limites et les ondes issues d’intéractions non linéaires, plus quelques ondes radiationnelles.

comodo-detidor -l list_file_IROI.txt --variable-mask --only-atlases -a -p ./IROISE -v XE -d M2 K1 S2 O1 P1 N2 K2 Q1 M1 Nu2 J1 L2 T2 2N2 Mu2 OO1 Ro1 Pi1 2Q1 Phi1 E2 La2 KJ2 Sig1 Ki1 Psi1 Tta1 R2 MP1 MSqm KQ1 Mqm 3OK1 MS1 S1 SO1 2MN2S2 2NS2 3M2S2 ST1 OQ2 MNS2 MNuS2 ST2 ST3 2MK2 SNK2 NKM2 MSN2 2SM2 SKM2 2SN2 2SMu2 MQ3 2MK3 MO3 M3 SO3 MS3 MK3 SP3 S3 SK3 K3 2MNS4 M4 M6 N4 S4 M8

Visualisation des différentes composantes harmoniques

tidegauge_harmcomp.py --ogcmatlas /home1/caparmor/vgarnier/TOOLS/TTB/VALID_IROI/GCC/ALL/WAVE-XE-atlas.nc --tidegauge 'BREST' --print --fes2004 --fes2012 --cstfrance --waves 'M2,K1,S2,O1,P1,N2,K2,Q1,M1,Nu2,J1,L2,T2,2N2,Mu2,OO1,Ro1,Pi1,2Q1,Phi1,E2,La2,KJ2,Sig1,Ki1,Psi1,Tta1,R2,MP1,MSqm,KQ1,Mqm,3OK1,MS1,S1,SO1,2MN2S2,2NS2,3M2S2,ST1,OQ2,MNS2,MNuS2,ST2,ST3,2MK2,SNK2,NKM2,MSN2,2SM2,SKM2,2SN2,2SMu2,MQ3,2MK3,MO3,M3,SO3,MS3,MK3,SP3,S3,SK3,K3,2MNS4,M4,M6,N4,S4,M8'

Au vu de la figure, on ne garde que les amplitudes significatives (forcées aux frontières) et pour lesquelles on a une référence dans les atlas, soit:

M2 S2 N2 K2 M6 M4 N4 S4 M8 K1 O1 P1 Q1 M1 Nu2 J1 L2 T2 2N2 Mu2 OO1 2Q1 E2 La2 R2 MP1 MSqm KQ1 3OK1 2MN2S2 2NS2 2SM2 2SN2 2SMu2 MQ3 2MK3 M3 SO3 SP3 K3 2MNS4

et on les met dans un ordre d’amplitude décroissant

tidegauge_harmcomp.py --ogcmatlas /home1/caparmor/vgarnier/TOOLS/TTB/VALID_IROI/GCC/IROI/WAVE-XE-atlas.nc --tidegauge 'BREST' --print --fes2004 --fes2012 --cstfrance --waves M2,S2,N2,K2,M6,M4,N4,S4,M8,K1,O1,P1,Q1,M1,Nu2,J1,L2,T2,2N2,Mu2,OO1,2Q1,E2,La2,R2,MP1,MSqm,KQ1,3OK1,2MN2S2,2NS2,2SM2,2SN2,2SMu2,MQ3,2MK3,M3,SO3,SP3,K3,2MNS4

ou

tidegauge_harmcomp.py --ogcmatlas /home1/caparmor/vgarnier/TOOLS/TTB/VALID_IROI/GCC/IROI/WAVE-XE-atlas.nc --tidegauge 'BREST' --print --fes2004 --fes2012 --cstfrance --waves M2,S2,N2,K2,M6,M4,N4,S4,M8,K1,O1,P1,Q1,M1,Nu2,J1,L2,T2,2N2,Mu2,OO1

Pas de reference pour Ro1 Pi1 Phi1

1. spatiale
2. marégraphe

Séries temporelles, champ résiduel