Mise en œuvre

Le programme permettant de filtrer les perturbations de la SST effectue donc les opérations suivantes :

• lecture des paramètres, avec entres autres l'année initiale, le nombre d'enregistrements temporels, la valeur de $\tau$, le nombre de variables à traiter et le nombre de fichiers annuels ;
• Lecture du mask et des facteurs d'échelle ;
• Pour chaque fichier, lecture des données ;
• Si demandé, test de l'adjoint du filtre récursif :
• Si demandé, comparaison du filtre récursif avec la valeur analytique de la fonction SOAR ;
• Calcul de la moyenne et de l'écart-type du champs initial pour chaque date ;
• Recentrage du champs initial pour avoir une moyenne nulle à chaque date ;
• Filtrage du champs initial ;
• Calcul de la moyenne et de l'écart-type du champs filtré pour chaque date ;
• Filtrage du rapport des variabilités ( $\sigma_{ini} / \sigma_{filtered}$) ;
• Multiplication du champs filtré par le rapport des variabilités filtrées ;
• Ecrétage du champs obtenu aux valeurs limites du champs initial ;
• Calcul de la moyenne du champs initial pour chaque date ;
• Recentrage du champs obtenu pour avoir une moyenne nulle à chaque date ;
• Calcul de la moyenne et de l'écart-type du champs final pour chaque date ;
• Ecriture du champs filtré dans des fichiers annuels.

Le champs obtenu a donc l'auto-corrélation temporelle voulue, ainsi que la même variabilité que le champs initial.