Changeset 633 for trunk/docs/docs_dev/source/guides/read_me.rst
- Timestamp:
- 02/11/14 18:12:25 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/docs/docs_dev/source/guides/read_me.rst
r443 r633 38 38 programme :ref:`readERA.m`: 39 39 40 :: 40 .. code-block:: matlab 41 41 42 >>file='CAT.r.2005.2007.6TU.nc';43 >>z_deb=[ -30 -30 datenum(2005,1,1) 350];44 >>z_fin=[ 40 40 datenum(2007,12,1) 350];42 file='CAT.r.2005.2007.6TU.nc'; 43 z_deb=[ -30 -30 datenum(2005,1,1) 350]; 44 z_fin=[ 40 40 datenum(2007,12,1) 350]; 45 45 46 46 zone définie par les paramètres début et fin, pour le niveau, il fait la moyenne entre les deux niveaux définis 47 47 48 :: 48 .. code-block:: matlab 49 49 50 >>chemin='/homedata/eymard/OPERA_Trop/';51 >>nom_var='r';52 >>[mat_3D,lat ,lon,date,niveau]=readERA(file,z_deb, z_fin, nom_var,chemin);53 >>[EOF2D,PC]=affiche_EOF_PC_3D(mat_3D,date,lon,lat);50 chemin='/homedata/eymard/OPERA_Trop/'; 51 nom_var='r'; 52 [mat_3D,lat ,lon,date,niveau]=readERA(file,z_deb, z_fin, nom_var,chemin); 53 [EOF2D,PC]=affiche_EOF_PC_3D(mat_3D,date,lon,lat); 54 54 55 55 options:choix des eof affichés: par défaut 1,2,3 sinon, ajouter [2 3 4] pour les trois suivants … … 57 57 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) 58 58 59 Pour sauvegarder la figure en format png: :59 Pour sauvegarder la figure en format png: 60 60 61 >> print -dpng figures/nomfig.png 61 .. code-block:: matlab 62 63 print -dpng figures/nomfig.png 62 64 63 65 AMSU 64 66 ---- 65 67 66 Avec :ref:`readAMSU.m`, on lit le fichier 2000-2007 par défaut, sinon 68 Avec :ref:`readAMSU.m`, on lit le fichier 2000-2007 par défaut, sinon 67 69 déclarer le chemin. 68 70 69 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`): :71 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`): 70 72 71 >> z_deb=[ -30 -30 datenum(2005,1,1)]; 72 >> z_fin=[ 40 40 datenum(2008,1,1)]; 73 >> nom_var='b3'; 74 >> file='b3.2000_2007.nc'; 75 >> [mat_3D,jours,latitude ,longitude]=readAMSU(file,z_deb,z_fin); 73 .. code-block:: matlab 76 74 77 puis lancer:: 75 z_deb=[ -30 -30 datenum(2005,1,1)]; 76 z_fin=[ 40 40 datenum(2008,1,1)]; 77 nom_var='b3'; 78 file='b3.2000_2007.nc'; 79 [mat_3D,jours,latitude ,longitude]=readAMSU(file,z_deb,z_fin); 78 80 79 >> [EOF2D,PC]=affiche_EOF_PC_3D(mat_3D,jours,longitude,latitude); 81 puis lancer: 82 83 .. code-block:: matlab 84 85 [EOF2D,PC]=affiche_EOF_PC_3D(mat_3D,jours,longitude,latitude); 80 86 81 87 Concernant b5, le problême vient du fait que dans ce canal, il y a des NaN … … 83 89 important dans la région du golfe persique. 84 90 85 Il faut donc refaire une interpolation pour boucher les trous restants : :91 Il faut donc refaire une interpolation pour boucher les trous restants : 86 92 87 >> mat_3D_b=boucheNaN_3D(mat_3D); 93 .. code-block:: matlab 88 94 89 puis refaire le calcul d'eof avec la nouvelle matrice:: 95 mat_3D_b=boucheNaN_3D(mat_3D); 90 96 91 >> [EOF2D,PC]=affiche_EOF_PC_3D(mat_3D_b,jours,longitude,latitude); 97 puis refaire le calcul d'eof avec la nouvelle matrice: 92 98 93 pour calculer des eofs sur une partie de l'année, chaque année, utiliser la fonction :ref:`choix_mois.m`:: 99 .. code-block:: matlab 94 100 95 >> [mat_3D_s, jours_s]=choix_mois(mat_3D,jours,mois_deb,mois_fin ); 101 [EOF2D,PC]=affiche_EOF_PC_3D(mat_3D_b,jours,longitude,latitude); 102 103 pour calculer des eofs sur une partie de l'année, chaque année, utiliser la fonction :ref:`choix_mois.m`: 104 105 .. code-block:: matlab 106 107 [mat_3D_s, jours_s]=choix_mois(mat_3D,jours,mois_deb,mois_fin ); 96 108 97 109 puis refaire tourner eof2d avec la bonne entrée. … … 102 114 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 103 115 104 Calcul des eof sans le cycle annuel : fonction: :116 Calcul des eof sans le cycle annuel : fonction: 105 117 106 >> [EOF2D,PC]=affiche_EOF_PC_3Dsscycle(mat3D_s,jours,longitude,latitude,titre_fig); 118 .. code-block:: matlab 107 119 120 [EOF2D,PC]=affiche_EOF_PC_3Dsscycle(mat3D_s,jours,longitude,latitude,titre_fig); 121
Note: See TracChangeset
for help on using the changeset viewer.