- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/mode_sahelien/composite_olr_sahel_1025.m
r18 r26 5 5 close all; 6 6 7 cd('..');8 7 ncload('olr.day.mean.nc'); 9 cd('mode_sahelien'); 8 10 9 % taille 12327x73x144 1er juin 1974-31 decembre2007 11 10 % on enleve les annees 1974 a 1978 nombre de jours a enlever au debut 214+365+366+365+365 … … 42 41 for c=1:siz; 43 42 vecteur(c)=olr1(a,b,c); 44 end 43 end; 45 44 std_c=nanmean(vecteur); 46 45 for c=1:siz; 47 46 olr1(a,b,c)=olr1(a,b,c)-std_c; 48 end 49 end 47 end; 48 end; 50 49 end; 51 50 … … 146 145 clear olr; 147 146 148 %load ind_olr_10_25_eof234.txt 149 load eof24_1025.txt 147 status=load('eof24_1025.txt'); 148 eof24_1025=status; 149 clear status; 150 150 % on met sous la forme annee mois 151 151 olr_eof=reshape(eof24_1025,122,28); … … 264 264 y=olr_filtre_compomax(1:nb_value,a,b,c); 265 265 student1(a,b,c)=ttest(x,y,.2); 266 end 267 end 268 end 266 end; 267 end; 268 end; 269 269 270 270 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences … … 305 305 end; 306 306 end; 307 end 307 end; 308 308 subplot(7,1,8-e); 309 309 palette; … … 319 319 else 320 320 set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 321 end 321 end; 322 322 set (gca,'YTick',[-10:10:40],'YTickLabel',[' ';' 0';'10N';'20N';'30N';' '],'fontname','Arial','fontsize',6); 323 323 if e==4; … … 332 332 if e==7; 333 333 title('composites OLR filtre 10-25 jours EOF234','fontsize',8) ; 334 end 334 end; 335 335 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 336 336 pos_vert=e/10; … … 339 339 340 340 figure(1); 341 print -depsc2 composite_olr_jjas_sahel_10_25_eof234_part1.eps;341 print('-depsc2','composite_olr_jjas_sahel_10_25_eof234_part1.eps'); 342 342 343 343 figure(2); … … 367 367 end; 368 368 end; 369 end 369 end; 370 370 subplot(7,1,8-e); 371 371 palette … … 381 381 else 382 382 set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 383 end 383 end; 384 384 set (gca,'YTick',[-10:10:40],'YTickLabel',[' ';' 0';'10N';'20N';'30N';' '],'fontname','Arial','fontsize',6); 385 385 if e==4; … … 394 394 if e==7; 395 395 title('composites OLR filtre 10-25 jours EOF234','fontsize',8) ; 396 end 396 end; 397 397 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 398 398 pos_vert=e/10; … … 401 401 402 402 figure(2); 403 print -depsc2 composite_olr_jjas_sahel_10_25_eof234_part2.eps;403 print('-depsc2','composite_olr_jjas_sahel_10_25_eof234_part2.eps');
Note: See TracChangeset
for help on using the changeset viewer.