Changeset 15 for trunk/src/mode_sahelien/composite_hgt500_eof23_30100.m
- Timestamp:
- 12/23/08 10:36:27 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/mode_sahelien/composite_hgt500_eof23_30100.m
r10 r15 106 106 end; 107 107 for b=1:122; % boucle sur les jours 108 if (b>20 & b<112); % 1ere condition (laisser une marge de 10)108 if (b>20 && b<112); % 1ere condition (laisser une marge de 10) 109 109 % cas 1 : max deltas apres negetifs et deltas avant postifs 110 if (delta(b+1)<0 & delta(b)<0 & delta(b-1)>0 & delta(b-2)>0& ind_olr(b,a)>ind_olr2);110 if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2); 111 111 compteurmax=compteurmax+1; 112 112 hgt500_filtre_compomax(compteurmax,1,:,:)=hgt(b-20,a,:,:); … … 144 144 hgt500_filtre_compomax(compteurmax,31,:,:)=hgt(b+10,a,:,:); 145 145 % cas 2 : min deltas apres positifs et deltas avant negatifs 146 elseif (delta(b+1)>0 & delta(b)>0 & delta(b-1)<0 & delta(b-2)<0& ind_olr(b,a)<(-ind_olr2));146 elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && ind_olr(b,a)<(-ind_olr2)); 147 147 compteurmin=compteurmin+1; 148 148 … … 184 184 end; 185 185 end; 186 compteurmin 187 compteurmax 186 disp(['iii : compteurmin =', int2str(compteurmin)]); 187 disp(['iii : compteurmax =', int2str(compteurmax)]); 188 188 % on fait le test de Student 189 189 nb_value=min(compteurmin,compteurmax); … … 223 223 if student1(f,b,c)==1; 224 224 delta_hgt500_filtre_min(b,c)=delta_hgt500_filtre_min(b,c); 225 else ;225 else 226 226 delta_hgt500_filtre_min(b,c)=NaN; 227 227 end; … … 246 246 if e==1; 247 247 set (gca,'XTick',[-60:30:60],'XTickLabel',['60W';'30W';' 0 ';'30W';'60W'],'fontsize',6); 248 else ;248 else 249 249 set (gca,'XTick',[-60:30:60],'XTickLabel',[' ';' ';' ';' ';' ']); 250 250 end; … … 252 252 if e==4; 253 253 co=colorbar; 254 set (co,'ytick',[-40:20:40],'yticklabel',[-40:20:40],'fontname','Arial','fontsize',6);254 set (co,'ytick',[-40:20:40],'yticklabel',[-40:20:40],'fontname','Arial','fontsize',6); 255 255 end; 256 256 hold on; cartemonde1; … … 263 263 text(-270,50,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 264 264 pos_vert=e/10; 265 set (gca,'position',[0.2 pos_vert 0.5 0.1]);265 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 266 266 end; 267 267 … … 281 281 if student1(f,b,c)==1; 282 282 delta_hgt500_filtre_min(b,c)=delta_hgt500_filtre_min(b,c); 283 else ;283 else 284 284 delta_hgt500_filtre_min(b,c)=NaN; 285 285 end; … … 305 305 if e==1; 306 306 set (gca,'XTick',[-60:30:60],'XTickLabel',['60W';'30W';' 0 ';'30W';'60W'],'fontsize',6); 307 else ;307 else 308 308 set (gca,'XTick',[-60:30:60],'XTickLabel',[' ';' ';' ';' ';' ']); 309 309 end; … … 311 311 if e==4; 312 312 co=colorbar; 313 set (co,'ytick',[-40:20:40],'yticklabel',[-40:20:40],'fontname','Arial','fontsize',6);313 set (co,'ytick',[-40:20:40],'yticklabel',[-40:20:40],'fontname','Arial','fontsize',6); 314 314 end; 315 315 hold on; cartemonde1; … … 322 322 text(-270,50,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 323 323 pos_vert=e/10; 324 set (gca,'position',[0.2 pos_vert 0.5 0.1]);324 set (gca,'position',[0.2 pos_vert 0.5 0.1]); 325 325 end; 326 326 print -depsc2 compo_hgt500_eof23_30100b.eps;
Note: See TracChangeset
for help on using the changeset viewer.