Changeset 16 for trunk/src/SIMULS_IRCAAM/composite_olr_AFR_eof1_30100.m
- Timestamp:
- 01/06/09 10:53:18 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/SIMULS_IRCAAM/composite_olr_AFR_eof1_30100.m
r15 r16 1 % 2 1 3 % initialisation 2 clear; close all; 4 clear; 5 close all; 3 6 olr=zeros(122,29,64,128); 4 7 … … 95 98 olr1(:,:,:,65:128)=olr(:,:,:,1:64); 96 99 clear olr; 97 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%98 100 99 101 load eof1_AFR_30100.txt 100 102 101 103 % on met sous la forme annee mois 102 olr_eof=reshape(eof1_AFR_30100,122,29); clear eof1_AFR_30100; 104 olr_eof=reshape(eof1_AFR_30100,122,29); 105 clear eof1_AFR_30100; 103 106 104 107 % on calcule l ecart type sur la serie temporelle de l OLR 105 k=122*29; 108 k=122*29; 106 109 olrbis=reshape(olr_eof,k,1); 107 olrbis=std(olrbis); 110 olrbis=std(olrbis); 108 111 109 112 % intialisation des tableaux 110 olr_filtre_compomax=NaN*ones(50,31,64,128); 113 olr_filtre_compomax=NaN*ones(50,31,64,128); 111 114 olr_filtre_compomin=NaN*ones(50,31,64,128); 112 115 113 compteurmax=0; 116 compteurmax=0; 114 117 compteurmin=0; 115 118 116 for a=1:29; % boucle sur les annees 117 clear delta; 119 % boucle sur les annees 120 for a=1:29; 121 clear delta; 118 122 delta=zeros(122,1); 119 for b=1:122; % boucle sur les jours : on prepare avant 123 % boucle sur les jours : on prepare avant 124 for b=1:122; 120 125 if b>1; 121 126 delta(b)=olr_eof(b,a)-olr_eof(b-1,a); 122 127 end; 123 128 end; 124 for b=1:122; % boucle sur les jours 125 if (b>20 && b<112); % 1ere condition (laisser une marge de 10) 126 % cas 1 : max deltas apres negetifs et deltas avant postifs 127 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 129 % boucle sur les jours 130 for b=1:122; 131 % 1ere condition (laisser une marge de 10) 132 if (b>20 && b<112); 133 % cas 1 : max deltas apres negetifs et deltas avant postifs 134 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 128 135 compteurmax=compteurmax+1; 129 136 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); … … 160 167 olr_filtre_compomax(compteurmax,30,:,:)=olr1(b+9,a,:,:); 161 168 olr_filtre_compomax(compteurmax,31,:,:)=olr1(b+10,a,:,:); 162 % cas 2 : min deltas apres positifs et deltas avant negatifs 169 % cas 2 : min deltas apres positifs et deltas avant negatifs 163 170 elseif (delta(b+1)>0 && delta(b)<0 && olr_eof(b,a)<(-olrbis)); 164 171 compteurmin=compteurmin+1; … … 199 206 end; 200 207 end; 201 end; 208 end; 202 209 end; 203 210 % on fait le test de Student … … 217 224 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 218 225 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 219 olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); clear olr_filtre_compomin; 220 olr_filtre_compomax1(:,:,:)=nanmean(olr_filtre_compomax); clear olr_filtre_compomax; 221 222 x=[-10:1:10]'; %' 223 224 figure(1); orient('landscape'); 225 226 lat=lat'; 226 olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); 227 clear olr_filtre_compomin; 228 olr_filtre_compomax1(:,:,:)=nanmean(olr_filtre_compomax); 229 clear olr_filtre_compomax; 230 231 x=[-10:1:10]'; 232 233 figure(1); 234 orient('landscape'); 235 236 lat=lat'; 227 237 lon=lon-180; 228 238 for e=1:7; … … 255 265 delta_olr_filtre_min(1,2)=50.0001; 256 266 fin=-30+24*2.5; 257 267 258 268 clear fin 259 269 fin=-10+16*2.5; … … 263 273 else 264 274 set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 265 end 275 end 266 276 set (gca,'YTick',[-10:10:30],'YTickLabel',['10S';' 0';'10N';'20N';'30N'],'fontname','Arial','fontsize',6); 267 277 if e==4; 268 co=colorbar; 278 co=colorbar; 269 279 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 270 280 end; 271 hold on; cartemonde1; 281 hold on; 282 cartemonde1; 272 283 set (gca,'tickDir','out'); 273 284 axis equal 274 285 axis([-180 180 -20 40]); 275 if e==7; 286 if e==7; 276 287 title('OLR composite EOF1 30100j simulation AFR ','fontsize',8) ; 277 288 end 278 289 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 279 290 pos_vert=e/10; 280 set (gca,'position',[0.2 pos_vert 0.5 0.1]) 291 set (gca,'position',[0.2 pos_vert 0.5 0.1]) 281 292 end; 282 293 283 figure(1); print -depsc2 composite_olr_EOF1_30100j_simulation_AFR_part1.eps; 284 285 figure(2); orient('landscape'); 294 figure(1); 295 print -depsc2 composite_olr_EOF1_30100j_simulation_AFR_part1.eps; 296 297 figure(2); 298 orient('landscape'); 286 299 287 300 for e=1:7; … … 314 327 delta_olr_filtre_min(1,2)=50.0001; 315 328 fin=-30+24*2.5; 316 329 317 330 clear fin 318 331 fin=-10+16*2.5; … … 322 335 else 323 336 set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 324 end 337 end 325 338 set (gca,'YTick',[-10:10:30],'YTickLabel',['10S';' 0';'10N';'20N';'30N'],'fontname','Arial','fontsize',6); 326 339 if e==4; 327 co=colorbar; 340 co=colorbar; 328 341 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 329 342 end; 330 hold on; cartemonde1; 343 hold on; 344 cartemonde1; 331 345 set (gca,'tickDir','out'); 332 346 axis equal … … 337 351 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 338 352 pos_vert=e/10; 339 set (gca,'position',[0.2 pos_vert 0.5 0.1]) 353 set (gca,'position',[0.2 pos_vert 0.5 0.1]) 340 354 end; 341 355 342 figure(2); print -depsc2 composite_olr_EOF1_30100j_simulation_AFR_part2.eps; 343 344 345 356 figure(2); 357 print -depsc2 composite_olr_EOF1_30100j_simulation_AFR_part2.eps;
Note: See TracChangeset
for help on using the changeset viewer.