ERA

programme :ref:`readERA.m`:

file='CAT.r.2005.2007.6TU.nc';
z_deb=[ -30 -30 datenum(2005,1,1) 350];
z_fin=[ 40 40 datenum(2007,12,1) 350];

zone définie par les paramètres début et fin, pour le niveau, il fait la moyenne entre les deux niveaux définis

chemin='/homedata/eymard/OPERA_Trop/';
nom_var='r';
[mat_3D,lat ,lon,date,niveau]=readERA(file,z_deb, z_fin, nom_var,chemin);
[EOF2D,PC]=affiche_EOF_PC_3D(mat_3D,date,lon,lat);

options:choix des eof affichés: par défaut 1,2,3 sinon, ajouter [2 3 4] pour les trois suivants

pour changer d'autres paramètres (titre, par ex, regarder :ref:`affiche_eofpc_3Drev.m` pour identifier les variables à ajouter (le titre est le dixième ou à peu près)

Pour sauvegarder la figure en format png:

print -dpng figures/nomfig.png

AMSU

Avec :ref:`readAMSU.m`, on lit le fichier 2000-2007 par défaut, sinon déclarer le chemin.

on peut aussi définir z_deb et z_fin pour le choix de la zone et la période (en julien comme pour :ref:`readERA.m`):

z_deb=[ -30 -30 datenum(2005,1,1)];
z_fin=[ 40 40 datenum(2008,1,1)];
nom_var='b3';
file='b3.2000_2007.nc';
[mat_3D,jours,latitude ,longitude]=readAMSU(file,z_deb,z_fin);

puis lancer:

[EOF2D,PC]=affiche_EOF_PC_3D(mat_3D,jours,longitude,latitude);

Concernant b5, le problême vient du fait que dans ce canal, il y a des NaN persistants après bouchage des trous car le nombre de points manquants est important dans la région du golfe persique.

Il faut donc refaire une interpolation pour boucher les trous restants :

mat_3D_b=boucheNaN_3D(mat_3D);

puis refaire le calcul d'eof avec la nouvelle matrice:

[EOF2D,PC]=affiche_EOF_PC_3D(mat_3D_b,jours,longitude,latitude);

pour calculer des eofs sur une partie de l'année, chaque année, utiliser la fonction :ref:`choix_mois.m`:

[mat_3D_s, jours_s]=choix_mois(mat_3D,jours,mois_deb,mois_fin );

puis refaire tourner eof2d avec la bonne entrée.

Pour un fichier mat3D(lat, lon, jours), on sélectionne tous les mois compris entre mois_deb et mois_fin inclus.

ex: pour prendre d'avril à août, on choisit mois_deb=4,mois_fin=8 , de décembre a février: mois_deb=12,mois_fin=2, toute l'année (par défaut): mois_deb=1,mois_fin=12

Calcul des eof sans le cycle annuel : fonction:

[EOF2D,PC]=affiche_EOF_PC_3Dsscycle(mat3D_s,jours,longitude,latitude,titre_fig);