Changeset 15 for trunk/src/mode_sahelien/composite_uwnd925_EOF24_1030.m
- Timestamp:
- 12/23/08 10:36:27 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/mode_sahelien/composite_uwnd925_EOF24_1030.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 vwnd925_filtre_compomax(compteurmax,1,:,:)=vwnd(b-20,a,:,:); … … 112 112 vwnd925_filtre_compomax(compteurmax,31,:,:)=vwnd(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 … … 208 208 elseif e==6; 209 209 sauvegrads('composite_vwnd925_mode_sahel_1030_moins16',delta_vwnd925,[-180 2.5 -10 2.5]); 210 else ;210 else 211 211 sauvegrads('composite_vwnd925_mode_sahel_1030_moins18',delta_vwnd925,[-180 2.5 -10 2.5]); 212 212 end … … 215 215 if e==1; 216 216 set (gca,'XTick',[-180:60:180],'XTickLabel',[' 180';'120W';' 60W';' 0';' 60E';'120E';' 180'],'fontname','Arial','fontsize',6); 217 else ;217 else 218 218 set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 219 219 end; … … 221 221 if e==4; 222 222 co=colorbar; 223 set (co,'ytick',[-5:1:5],'yticklabel',[-5:1:5],'fontname','Arial','fontsize',6);223 set (co,'ytick',[-5:1:5],'yticklabel',[-5:1:5],'fontname','Arial','fontsize',6); 224 224 end; 225 225 hold on; cartemonde1; … … 232 232 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 233 233 pos_vert=e/10; 234 set (gca,'position',[0.2 pos_vert 0.5 0.1]);234 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 235 235 end; 236 236 … … 281 281 elseif e==6; 282 282 sauvegrads('composite_vwnd925_mode_sahel_1030_moins2',delta_vwnd925,[-180 2.5 -10 2.5]); 283 else ;283 else 284 284 sauvegrads('composite_vwnd925_mode_sahel_1030_moins4',delta_vwnd925,[-180 2.5 -10 2.5]); 285 285 end; … … 290 290 if e==1; 291 291 set (gca,'XTick',[-180:60:180],'XTickLabel',[' 180';'120W';' 60W';' 0';' 60E';'120E';' 180'],'fontname','Arial','fontsize',6); 292 else ;292 else 293 293 set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 294 294 end; … … 296 296 if e==4; 297 297 co=colorbar; 298 set (co,'ytick',[-5:1:5],'yticklabel',[-5:1:5],'fontname','Arial','fontsize',6);298 set (co,'ytick',[-5:1:5],'yticklabel',[-5:1:5],'fontname','Arial','fontsize',6); 299 299 end; 300 300 hold on; cartemonde1; … … 307 307 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 308 308 pos_vert=e/10; 309 set (gca,'position',[0.2 pos_vert 0.5 0.1]);309 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 310 310 end; 311 311
Note: See TracChangeset
for help on using the changeset viewer.