Changeset 15 for trunk/src/mode_sahelien/composite_wind925_EOF1_30100.m
- Timestamp:
- 12/23/08 10:36:27 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/mode_sahelien/composite_wind925_EOF1_30100.m
r10 r15 74 74 end; 75 75 for b=1:122; % boucle sur les jours 76 if (b>20 & b<112); % 1ere condition (laisser une marge de 10)76 if (b>20 && b<112); % 1ere condition (laisser une marge de 10) 77 77 % cas 1 : max deltas apres negetifs et deltas avant postifs 78 if (delta(b+1)<0 & delta(b)<0 & delta(b-1)>0 & delta(b-2)>0& ind_olr(b,a)>ind_olr2);78 if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2); 79 79 compteurmax=compteurmax+1; 80 80 uwnd925_filtre_compomax(compteurmax,1,:,:)=uwnd(b-20,a,:,:); … … 112 112 uwnd925_filtre_compomax(compteurmax,31,:,:)=uwnd(b+10,a,:,:); 113 113 % cas 2 : min deltas apres positifs et deltas avant negatifs 114 elseif (delta(b+1)>0 & delta(b)>0 & delta(b-1)<0 & delta(b-2)<0& ind_olr(b,a)<(-ind_olr2));114 elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && ind_olr(b,a)<(-ind_olr2)); 115 115 compteurmin=compteurmin+1; 116 116 … … 231 231 end; 232 232 for b=1:122; % boucle sur les jours 233 if (b>20 & b<112); % 1ere condition (laisser une marge de 10)233 if (b>20 && b<112); % 1ere condition (laisser une marge de 10) 234 234 % cas 1 : max deltas apres negetifs et deltas avant postifs 235 if (delta(b+1)<0 & delta(b)<0 & delta(b-1)>0 & delta(b-2)>0& ind_olr(b,a)>ind_olr2);235 if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2); 236 236 compteurmax=compteurmax+1; 237 237 vwnd925_filtre_compomax(compteurmax,1,:,:)=vwnd(b-20,a,:,:); … … 269 269 vwnd925_filtre_compomax(compteurmax,31,:,:)=vwnd(b+10,a,:,:); 270 270 % cas 2 : min deltas apres positifs et deltas avant negatifs 271 elseif (delta(b+1)>0 & delta(b)>0 & delta(b-1)<0 & delta(b-2)<0& ind_olr(b,a)<(-ind_olr2));271 elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && ind_olr(b,a)<(-ind_olr2)); 272 272 compteurmin=compteurmin+1; 273 273 … … 352 352 if e==1; 353 353 set (gca,'XTick',[-60:20:120],'XTickLabel',[' 60W';' 40W';' 20W';' 0 ';' 20E';' 40E';' 60E'; ' 80E';'100E';'120E'],'fontsize',6); 354 else ;354 else 355 355 set (gca,'XTick',[-60:20:120],'XTickLabel',[' ';' ';' ';' ';' ';' ';' ';' ';' ';' '],'fontsize',6); 356 356 end … … 366 366 text(-150,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 367 367 pos_vert=e/6.1; 368 % set (gca,'position',[0.2 pos_vert 0.3 1/6.1])369 set (gca,'position',[0.2 pos_vert 0.4 1/6.1])368 % set (gca,'position',[0.2 pos_vert 0.3 1/6.1]) 369 set (gca,'position',[0.2 pos_vert 0.4 1/6.1]) 370 370 end; 371 371 … … 411 411 if e==1; 412 412 set (gca,'XTick',[-60:20:120],'XTickLabel',[' 60W';' 40W';' 20W';' 0 ';' 20E';' 40E';' 60E'; ' 80E';'100E';'120E'],'fontsize',6); 413 else ;413 else 414 414 set (gca,'XTick',[-60:20:120],'XTickLabel',[' ';' ';' ';' ';' ';' ';' ';' ';' ';' '],'fontsize',6); 415 415 end … … 425 425 text(-150,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 426 426 pos_vert=e/6.1; 427 set (gca,'position',[0.2 pos_vert 0.4 1/6.1])427 set (gca,'position',[0.2 pos_vert 0.4 1/6.1]) 428 428 end; 429 429 … … 465 465 cartev=interp2(lon-180,lat,delta_vwnd925_filtre,xi,yi); 466 466 aa=quiver(xi,yi,carteu,cartev,'k'); 467 467 if e==1; 468 468 set (gca,'XTick',[-60:20:120],'XTickLabel',[' 60W';' 40W';' 20W';' 0 ';' 20E';' 40E';' 60E'; ' 80E';'100E';'120E'],'fontsize',6); 469 else ;469 else 470 470 set (gca,'XTick',[-60:20:120],'XTickLabel',[' ';' ';' ';' ';' ';' ';' ';' ';' ';' '],'fontsize',6); 471 471 end … … 481 481 text(-150,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 482 482 pos_vert=e/6.1; 483 set (gca,'position',[0.2 pos_vert 0.4 1/6.1])483 set (gca,'position',[0.2 pos_vert 0.4 1/6.1]) 484 484 end; 485 485 … … 521 521 cartev=interp2(lon-180,lat,delta_vwnd925_filtre,xi,yi); 522 522 aa=quiver(xi,yi,carteu,cartev,'k'); 523 523 if e==1; 524 524 set (gca,'XTick',[-60:20:120],'XTickLabel',[' 60W';' 40W';' 20W';' 0 ';' 20E';' 40E';' 60E'; ' 80E';'100E';'120E'],'fontsize',6); 525 else ;525 else 526 526 set (gca,'XTick',[-60:20:120],'XTickLabel',[' ';' ';' ';' ';' ';' ';' ';' ';' ';' '],'fontsize',6); 527 527 end … … 537 537 text(-150,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 538 538 pos_vert=e/6.1; 539 set (gca,'position',[0.2 pos_vert 0.4 1/6.1])539 set (gca,'position',[0.2 pos_vert 0.4 1/6.1]) 540 540 end; 541 541
Note: See TracChangeset
for help on using the changeset viewer.