Changeset 16 for trunk/src/mode_sahelien/composite_olr_ER1030.m
- Timestamp:
- 01/06/09 10:53:18 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/mode_sahelien/composite_olr_ER1030.m
r15 r16 1 % 2 1 3 % initialisation 2 clear; close all; 4 clear; 5 close all; 3 6 4 7 cd('..'); … … 25 28 olr1=zeros(siz,21,144); 26 29 olr1(:,:,1:72)=olr(:,:,73:144); 27 olr1(:,:,73:144)=olr(:,:,1:72); clear olr; 30 olr1(:,:,73:144)=olr(:,:,1:72); 31 clear olr; 28 32 29 33 % on permture le temps et l'espce … … 37 41 vecteur=zeros(siz,1); 38 42 for c=1:siz; 39 vecteur(c)=olr1(a,b,c); 43 vecteur(c)=olr1(a,b,c); 40 44 end 41 45 std_c=nanmean(vecteur); 42 46 for c=1:siz; 43 olr1(a,b,c)=olr1(a,b,c)-std_c; 44 end 47 olr1(a,b,c)=olr1(a,b,c)-std_c; 48 end 45 49 end 46 end; 47 48 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 50 end; 51 49 52 % on met sous la forme longitude latitude jour 50 51 53 52 54 % on selectionne les mois d ete (attention aux annees bixestiles!) … … 146 148 olr1=permute(olr,[3 4 1 2]); 147 149 clear olr; 148 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 150 149 151 load eof12_rossby_1030.txt 150 152 151 152 153 % on met sous la forme annee mois 153 olr_eof=reshape(eof12_rossby_1030,122,28); clear eof12_rossby_1030; 154 olr_eof=reshape(eof12_rossby_1030,122,28); 155 clear eof12_rossby_1030; 154 156 155 157 % on calcule l ecart type sur la serie temporelle de l OLR 156 k=122*28; 158 k=122*28; 157 159 olrbis=reshape(olr_eof,k,1); 158 olrbis=std(olrbis); 160 olrbis=std(olrbis); 159 161 160 162 % intialisation des tableaux 161 olr_filtre_compomax=NaN*ones(50,31,21,144); 163 olr_filtre_compomax=NaN*ones(50,31,21,144); 162 164 olr_filtre_compomin=NaN*ones(50,31,21,144); 163 165 164 compteurmax=0; 166 compteurmax=0; 165 167 compteurmin=0; 166 168 167 for a=1:28; % boucle sur les annees 168 clear delta; 169 % boucle sur les annees 170 for a=1:28; 171 clear delta; 169 172 delta=zeros(122,1); 170 for b=1:122; % boucle sur les jours : on prepare avant 173 % boucle sur les jours : on prepare avant 174 for b=1:122; 171 175 if b>1; 172 176 delta(b)=olr_eof(b,a)-olr_eof(b-1,a); 173 177 end; 174 178 end; 175 for b=1:122; % boucle sur les jours 176 if (b>20 && b<112); % 1ere condition (laisser une marge de 10) 177 % cas 1 : max deltas apres negetifs et deltas avant postifs 178 if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && olr_eof(b,a)>olrbis); 179 % boucle sur les jours 180 for b=1:122; 181 % 1ere condition (laisser une marge de 10) 182 if (b>20 && b<112); 183 % cas 1 : max deltas apres negetifs et deltas avant postifs 184 if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && olr_eof(b,a)>olrbis); 179 185 compteurmax=compteurmax+1; 180 186 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); … … 211 217 olr_filtre_compomax(compteurmax,30,:,:)=olr1(b+9,a,:,:); 212 218 olr_filtre_compomax(compteurmax,31,:,:)=olr1(b+10,a,:,:); 213 % cas 2 : min deltas apres positifs et deltas avant negatifs 219 % cas 2 : min deltas apres positifs et deltas avant negatifs 214 220 elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && olr_eof(b,a)<(-olrbis)); 215 221 compteurmin=compteurmin+1; … … 250 256 end; 251 257 end; 252 end; 258 end; 253 259 end; 254 260 % on fait le test de Student … … 268 274 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 269 275 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 270 olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); clear olr_filtre_compomin; 271 olr_filtre_compomax1(:,:,:)=nanmean(olr_filtre_compomax); clear olr_filtre_compomax; 272 273 x=[-10:1:10]'; %' 274 275 figure(1); orient('landscape'); 276 277 lat=lat(21:41)'; 276 olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); 277 clear olr_filtre_compomin; 278 olr_filtre_compomax1(:,:,:)=nanmean(olr_filtre_compomax); 279 clear olr_filtre_compomax; 280 281 x=[-10:1:10]'; 282 283 figure(1); 284 orient('landscape'); 285 286 lat=lat(21:41)'; 278 287 lon=lon-180; 279 288 for e=1:7; … … 306 315 delta_olr_filtre_min(1,2)=50.0001; 307 316 fin=-30+24*2.5; 308 317 309 318 clear fin 310 319 fin=-10+16*2.5 … … 331 340 else 332 341 set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 333 end 342 end 334 343 set (gca,'YTick',[-10:10:40],'YTickLabel',[' ';' 0';'10N';'20N';'30N';' '],'fontname','Arial','fontsize',6); 335 344 if e==4; 336 co=colorbar; 345 co=colorbar; 337 346 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 338 347 end; 339 hold on; cartemonde1; 348 hold on; 349 cartemonde1; 340 350 set (gca,'tickDir','out'); 341 351 axis equal … … 346 356 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 347 357 pos_vert=e/10; 348 set (gca,'position',[0.2 pos_vert 0.5 0.1]) 358 set (gca,'position',[0.2 pos_vert 0.5 0.1]) 349 359 end; 350 360 351 figure(1); print -depsc2 composite_olr_jjas_rossby_1030_eof2a.eps; 352 353 figure(2); orient('landscape'); 361 figure(1); 362 print -depsc2 composite_olr_jjas_rossby_1030_eof2a.eps; 363 364 figure(2); 365 orient('landscape'); 354 366 355 367 for e=1:7; … … 382 394 delta_olr_filtre_min(1,2)=50.0001; 383 395 fin=-30+24*2.5; 384 396 385 397 clear fin 386 398 fin=-10+16*2.5; … … 407 419 else 408 420 set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 409 end 421 end 410 422 set (gca,'YTick',[-10:10:40],'YTickLabel',[' ';' 0';'10N';'20N';'30N';' '],'fontname','Arial','fontsize',6); 411 423 if e==4; 412 co=colorbar; 424 co=colorbar; 413 425 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 414 426 end; 415 hold on; cartemonde1; 427 hold on; 428 cartemonde1; 416 429 set (gca,'tickDir','out'); 417 430 axis equal … … 422 435 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 423 436 pos_vert=e/10; 424 set (gca,'position',[0.2 pos_vert 0.5 0.1]) 437 set (gca,'position',[0.2 pos_vert 0.5 0.1]) 425 438 end; 426 439 427 figure(2); print -depsc2 composite_olr_jjas_rossby_1030_eof2b.eps; 428 429 440 figure(2); 441 print -depsc2 composite_olr_jjas_rossby_1030_eof2b.eps;
Note: See TracChangeset
for help on using the changeset viewer.