Changeset 27
- Timestamp:
- 08/26/09 12:46:40 (15 years ago)
- Location:
- trunk/RESULTATS
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/RESULTATS/carteetstd_SST.m
- Property svn:keywords set to Id
r25 r27 1 1 %CARTEETSTD_SST permet de visualiser les cartes mensuelles moyennes 2 2 % de temperature de surface 3 % ainsi que les cartes annuelles moyennes et les ecart type3 % ainsi que les cartes annuelles moyennes et les écart type 4 4 % correspondants 5 5 … … 19 19 % ``carteetstd_SST`` permet de visualiser les cartes mensuelles moyennes 20 20 % de temperature de surface 21 % ainsi que les cartes annuelles moyennes et les ecart type21 % ainsi que les cartes annuelles moyennes et les écart type 22 22 % correspondants. 23 23 % … … 36 36 % (figure 3). 37 37 % 38 % ``carteetstd_SST`` sauve ces images dans ``./carteetstd_SST_[123]. eps``.38 % ``carteetstd_SST`` sauve ces images dans ``./carteetstd_SST_[123].ps``. 39 39 % 40 40 % EXAMPLES … … 70 70 % .. _contlab.m : contlab.m.html 71 71 % 72 % 73 % KNOWN PROBLEMS 74 % ============== 75 % 76 % pb avec les diacritiques à l'impression de .ps dans le titres des figures 77 % sous octave : é se transforme en o barré. 78 % 72 79 % TODO 73 80 % ==== … … 80 87 % $Id$ 81 88 % 89 % - fplod 2009-08-26T09:56:51Z aedon.locean-ipsl.upmc.fr (Darwin) 90 % 91 % * nom des mois à la place du numéro dans les titres 92 % * remplacement de ``eps`` par ``ps` pour pouvoir imprimer les figures 93 % produites par ``octave`` 94 % 95 % 96 % - jmignot 2009-08-25 97 % 98 % * suppression d'appel à contlab et colorbartype sous octave 99 % 82 100 % - fplod 2009-08-25T10:40:16Z aedon.locean-ipsl.upmc.fr (Darwin) 83 101 % 84 % 102 % * modif commentaires liée au dernier chgt de Juliette 85 103 % 86 104 % - jmignot 2009-08-25 … … 142 160 % MANIPULATIONS TEMPORELLES: MOYENNES 143 161 %--------------------------------------------------- 144 tab=reshape(tab,size(tab,1),size(tab,2),12,nannees); %organisation du tableau sous la forme [lat lon mois ann ees]162 tab=reshape(tab,size(tab,1),size(tab,2),12,nannees); %organisation du tableau sous la forme [lat lon mois années] 145 163 146 164 % Construction de moyennes annuelles … … 152 170 tab_mens=squeeze(mean(tab,4)); 153 171 154 %calcul des ecart-types172 %calcul des écart-types 155 173 %------------------------ 156 % ecart-type mensuel moyen:174 %écart-type mensuel moyen: 157 175 std_mens=squeeze(std(tab_mens,0,3)); 158 176 159 % ecart-type interannuel moyen:177 %écart-type interannuel moyen: 160 178 std_an=squeeze(std(tab_an,0,3)); 161 179 … … 212 230 213 231 %titre: nom de chaque mois 214 title(['mois ' num2str(iplot)]) 232 cmois = datestr(datenum([2000 iplot 1 00 00 00]), 'mmm'); 233 title(cmois); 234 clear cmois; 215 235 216 236 end; … … 227 247 228 248 % sauvegarde de la figure 229 printer=' eps';249 printer='ps'; 230 250 print_printer=['-d', printer]; 231 251 fullfilename=['./',mfilename,'_',num2str(ifigure),'.',printer]; … … 279 299 end; 280 300 281 %titre: num ero de l'annee301 %titre: numéro de l'année 282 302 title(['annee ' num2str(iplot+1948-1)]) 283 303 end; … … 294 314 295 315 % sauvegarde de la figure 296 printer=' eps';316 printer='ps'; 297 317 print_printer=['-d', printer]; 298 318 fullfilename=['./',mfilename,'_',num2str(ifigure),'.',printer]; … … 303 323 304 324 305 % ecart type325 % écart type 306 326 %----------- 307 327 … … 313 333 initfig; 314 334 315 % ecart-type mensuel moyen335 %écart-type mensuel moyen 316 336 317 337 %niveaux de couleurs … … 355 375 356 376 357 % ecart-type interannuel moyen377 %écart-type interannuel moyen 358 378 359 379 %niveaux de couleurs … … 379 399 ylabel('latitude') 380 400 381 title(' ecart-type interannuel')401 title('écart-type interannuel') 382 402 %barre de couleur 383 403 if (run_octave == 0) … … 393 413 394 414 % sauvegarde de la figure 395 printer=' eps';415 printer='ps'; 396 416 print_printer=['-d', printer]; 397 417 fullfilename=['./',mfilename,'_',num2str(ifigure),'.',printer]; -
trunk/RESULTATS/eof_NCEP_an.m
- Property svn:keywords set to Id
r25 r27 8 8 % 9 9 % le programme permet egalement de visualiser les resultats. Voir le doc 10 % .pdf joint pour l'interpretation de cette visualisation. 10 % .pdf joint pour l'interpretation de cette visualisation. 11 11 % 12 12 % Date: Avril-Mai 2008 … … 18 18 % ====== 19 19 % 20 % ``eof_NCEP_an`` 20 % ``eof_NCEP_an`` 21 21 % 22 22 % DESCRIPTION 23 23 % =========== 24 24 % 25 % ``eof_NCEP_an`` lit les fichiers 25 % ``eof_NCEP_an`` lit les fichiers 26 26 % ``../DONNEES/lsmask_tropatl_PacandMedblancs.mat`` 27 27 % et ``../DONNEES/skt.mon.tropatl.mat``. … … 34 34 % 35 35 % ``eof_NCEP_an`` affiche trois figures : 36 % 36 37 % - ++ 38 % 37 39 % - ++ 40 % 38 41 % - valeurs propres de la matrice de covariance 39 42 % 40 % ``eof_NCEP_an`` sauve ces images dans ``./eof_SSTan_tropAlt30N-20S ++_[123].eps``.43 % ``eof_NCEP_an`` sauve ces images dans ``./eof_SSTan_tropAlt30N-20S_[123].ps``. 41 44 % 42 45 % ++ … … 77 80 % tester la branche NetCDF 78 81 % 79 % sous octave, triangle bizarre dans la dernière barre du graphique valeurs 82 % sous octave, triangle bizarre dans la dernière barre du graphique valeurs 80 83 % propres. 81 84 % … … 84 87 % 85 88 % $Id$ 89 % 90 % - fplod 2009-08-26T10:26:35Z aedon.locean-ipsl.upmc.fr (Darwin) 91 % 92 % * remplacement de ``eps`` par ``ps` pour pouvoir imprimer les figures 93 % produites par ``octave`` 94 % 95 % - jmignot 2009-08-25 96 % 97 % * suppression d'appel à contlab et colorbartype sous octave 86 98 % 87 99 % - fplod 2009-08-25T10:40:16Z aedon.locean-ipsl.upmc.fr (Darwin) … … 98 110 % - fplod 2009-08-21T15:04:16Z aedon.locean-ipsl.upmc.fr (Darwin) 99 111 % 100 % * as octave save default format is ASCII and no default extension 101 % while matlab save default format is MAT v5 mat-file (little endian) and 112 % * as octave save default format is ASCII and no default extension 113 % while matlab save default format is MAT v5 mat-file (little endian) and 102 114 % default extension is ``.mat``, write more precise save instruction 103 115 % * add header … … 130 142 131 143 nannees=60; 132 144 133 145 repertoire=[ '..' filesep 'DONNEES' filesep ] % a adapter a chaque cas: chemin d acces aux donnees 134 146 … … 137 149 file='lsmask_tropatl_PacandMedblancs.mat' 138 150 load([repertoire file]); 139 ny=size(mask,1); nx=size(mask,2); % dimensions spatiales des tableaux 151 ny=size(mask,1); nx=size(mask,2); % dimensions spatiales des tableaux 140 152 141 153 % Data … … 146 158 load([repertoire2 file]); % chargement des donnees 147 159 148 160 149 161 %SI NETCDF 150 162 %repertoire2=[repertoire 'NETCDF' filesep]% a adapter a chaque cas: chemin d acces aux donnees … … 152 164 % f=netcdf([repertoire file],'nowrite'); 153 165 % tab=f{'skt'}(:,:,:); 154 166 155 167 156 168 tab=permute(tab,[2 3 1]);% permutation des dimensions pour organiser le tableau sous la forme [lat lon temps] … … 179 191 %---------------- 180 192 %Si on veut éventuellement attribuer un poids différent aux différents 181 %points spatiaux, on le fait ici. 193 %points spatiaux, on le fait ici. 182 194 183 195 % Préparation du zoom géographique: sur le masque, j'identifie les points … … 187 199 limite_latS=-20; %limite du domaine d'etuade en latitude 188 200 limite_latN=30; %limite du domaine d'etuade en latitude 189 lat2D=repmat(lat,[1 size(mask,2)]); 201 lat2D=repmat(lat,[1 size(mask,2)]); 190 202 mask(lat2D>limite_latN&mask==0)=4; 191 203 mask(lat2D<limite_latS&mask==0)=4; … … 219 231 % matrice de covariance: 220 232 C = Z * Z'; 221 % eval(['save ' repertoire 'MatCov_SST_tropAtl' num2str(limite_latN) 'N' num2str(limite_latS) 'S.asc C -ascii']) 222 223 % diagonalisation 224 [E,Lambda] = eig(C); 225 % Lambda est la matrice diagonale des valeurs propres et E la matrice dont les colonnes sont les vecteurs propres correspondants. 233 % eval(['save ' repertoire 'MatCov_SST_tropAtl' num2str(limite_latN) 'N' num2str(limite_latS) 'S.asc C -ascii']) 234 235 % diagonalisation 236 [E,Lambda] = eig(C); 237 % Lambda est la matrice diagonale des valeurs propres et E la matrice dont les colonnes sont les vecteurs propres correspondants. 226 238 % Soit Z*E = E*Lambda 227 % Les vecteurs propres sont appelés les EOF. Chaque EOF peut être considéré comme une carte de l'espace physique. 228 229 230 %% 239 % Les vecteurs propres sont appelés les EOF. Chaque EOF peut être considéré comme une carte de l'espace physique. 240 241 242 %% 231 243 %%%% calcul de la matrice donnant les pourcentages de variance 232 244 Lambda = diag(Lambda); 233 245 pctg_var = round(Lambda / sum(Lambda) * 1000) / 10; 234 246 235 247 236 248 %On definit les composantes principales de Q comme la decomposition du signal decrit par Q sur les EOF E. 237 249 238 250 P = Z' * E; 239 240 %% 251 252 %% 241 253 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 242 254 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% … … 288 300 %niveaux de couleurs 289 301 lev=[-2:.02:2]; 290 clip_lev=[-0.06 0.06] 302 clip_lev=[-0.06 0.06]; 291 303 292 304 figure(1); 293 305 initfig; 294 figure(2); 306 figure(2); 295 307 initfig; 296 308 297 309 for i=1:4 298 310 299 311 %subplot(4,2,2*i-1) 300 312 ifigure=1; … … 309 321 z=squeeze(TABEOF(:,:,i)); 310 322 contourf(x,y,z,lev); 311 caxis(clip_lev) 323 caxis(clip_lev); 312 324 hold on 313 325 whorldmap … … 329 341 if (run_octave == 0) 330 342 hpal=jet(100); %carte des couleurs: du bleu au rouge par defaut. 331 pos1=get(gca,'Position') 343 pos1=get(gca,'Position'); 332 344 pos=[pos1(1) pos1(2)-.075 pos1(3) .013]; 333 345 colorbartype(pos,lev,1,clip_lev,hpal,0); … … 336 348 colorbar('East') 337 349 end; 338 350 339 351 % sauvegarde de la figure 340 printer=' eps';352 printer='ps'; 341 353 print_printer=['-d', printer]; 342 354 fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer]; … … 344 356 clear fullfilename; 345 357 clear printer; 346 358 347 359 %subplot(4,2,2*i) 348 360 ifigure=2; … … 357 369 358 370 % sauvegarde de la figure 359 printer=' eps';371 printer='ps'; 360 372 print_printer=['-d', printer]; 361 373 fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer]; … … 369 381 % Graphe pourcentages de variance expliqué 370 382 ifigure=3; 371 figure(ifigure); 383 figure(ifigure); 372 384 bar(flipud(pctg_var)) 373 385 xlim([0 10]) … … 377 389 378 390 % sauvegarde de la figure 379 printer=' eps';391 printer='ps'; 380 392 print_printer=['-d', printer]; 381 393 fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer]; -
trunk/RESULTATS/regandsignif_seriestemporelles.m
- Property svn:keywords set to Id
r24 r27 28 28 % 29 29 % ``regandsignif_seriestemporelles`` sauve cette image dans 30 % ``./eof_SSTan_tropAtl30N-20S.mat_regandsignif_seriestemporelles_1. eps``.30 % ``./eof_SSTan_tropAtl30N-20S.mat_regandsignif_seriestemporelles_1.ps``. 31 31 % 32 32 % EXAMPLES … … 233 233 234 234 % sauvegarde de la figure 235 printer=' eps';235 printer='ps'; 236 236 print_printer=['-d', printer]; 237 237 fullfilename=['.' filesep file '_' mfilename '_' num2str(ifigure) '.' printer];
Note: See TracChangeset
for help on using the changeset viewer.