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