Changeset 15 for trunk/src/mode_sahelien/composite_olr_EOF12_ER1030.m
- Timestamp:
- 12/23/08 10:36:27 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/mode_sahelien/composite_olr_EOF12_ER1030.m
r10 r15 153 153 k=122*28; 154 154 olrbis=reshape(olr_eof,k,1); 155 olr_std=std(olrbis) 156 olr_moy=mean(olrbis) 155 olr_std=std(olrbis); 156 olr_moy=mean(olrbis); 157 157 for a=1:3416 158 158 olrbis(a)=(olrbis(a)-olr_moy)/olr_std; … … 176 176 end; 177 177 for b=1:122; % boucle sur les jours 178 if (b>20 & b<112); % 1ere condition (laisser une marge de 10)178 if (b>20 && b<112); % 1ere condition (laisser une marge de 10) 179 179 % cas 1 : max deltas apres negetifs et deltas avant postifs 180 if (delta(b+1)<0 & delta(b)>0& olr_eof(b,a)>1);181 % if (delta(b+1)<0 & delta(b)>0& olr_eof(b,a)>olr_std);180 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>1); 181 % if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>olr_std); 182 182 compteurmax=compteurmax+1; 183 183 format short g; … … 217 217 olr_filtre_compomax(compteurmax,31,:,:)=olr1(b+10,a,:,:); 218 218 % cas 2 : min deltas apres positifs et deltas avant negatifs 219 elseif (delta(b+1)>0 & delta(b)<0& olr_eof(b,a)<-1);220 % elseif (delta(b+1)>0 & delta(b)<0& olr_eof(b,a)<(-olr_std));219 elseif (delta(b+1)>0 && delta(b)<0 && olr_eof(b,a)<-1); 220 % elseif (delta(b+1)>0 && delta(b)<0 && olr_eof(b,a)<(-olr_std)); 221 221 compteurmin=compteurmin+1; 222 222 format short g … … 259 259 end; 260 260 end; 261 compteurmin_eof12_er1030=compteurmin 262 compteurmax_eof12_er1030=compteurmax 261 compteurmin_eof12_er1030=compteurmin; 262 compteurmax_eof12_er1030=compteurmax; 263 263 264 264 % on fait le test de Student … … 297 297 if student1(f,b,c)==1; 298 298 delta_olr_filtre_min(b,c)=delta_olr_filtre_min(b,c); 299 else ;299 else 300 300 % delta_olr_filtre_min(b,c)=NaN; 301 301 end; … … 328 328 if e==4; 329 329 co=colorbar; 330 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6);330 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 331 331 end; 332 332 hold on; cartemonde1; … … 339 339 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 340 340 pos_vert=e/10; 341 set (gca,'position',[0.2 pos_vert 0.5 0.1])341 set (gca,'position',[0.2 pos_vert 0.5 0.1]) 342 342 end; 343 343 … … 357 357 if student1(f,b,c)==1; 358 358 delta_olr_filtre_min(b,c)=delta_olr_filtre_min(b,c); 359 else ;359 else 360 360 % delta_olr_filtre_min(b,c)=NaN; 361 361 end; … … 388 388 if e==4; 389 389 co=colorbar; 390 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6);390 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 391 391 end; 392 392 hold on; cartemonde1; … … 399 399 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 400 400 pos_vert=e/10; 401 set (gca,'position',[0.2 pos_vert 0.5 0.1])401 set (gca,'position',[0.2 pos_vert 0.5 0.1]) 402 402 end; 403 403
Note: See TracChangeset
for help on using the changeset viewer.