Changeset 16 for trunk/src/SIMULS_IRCAAM/composite_olr_CTL_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_CTL_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('CTL') … … 90 93 ncload('rlut_d.CtIV.200006-09.nc','rlut','lat','lon'); 91 94 olr(:,29,:,:)=rlut; clear rlut; 92 93 95 94 96 cd('..'); … … 98 100 olr1(:,:,:,65:128)=olr(:,:,:,1:64); 99 101 clear olr; 100 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%101 102 102 103 load eof2_CTL_30100.txt 103 104 104 105 % on met sous la forme annee mois 105 olr_eof=reshape(eof2_CTL_30100,122,29); clear eof2_CTL_30100; 106 olr_eof=reshape(eof2_CTL_30100,122,29); 107 clear eof2_CTL_30100; 106 108 107 109 % on calcule l ecart type sur la serie temporelle de l OLR 108 k=122*29; 110 k=122*29; 109 111 olrbis=reshape(olr_eof,k,1); 110 olrbis=std(olrbis); 112 olrbis=std(olrbis); 111 113 112 114 % intialisation des tableaux 113 olr_filtre_compomax=NaN*ones(50,31,64,128); 115 olr_filtre_compomax=NaN*ones(50,31,64,128); 114 116 olr_filtre_compomin=NaN*ones(50,31,64,128); 115 117 116 compteurmax=0; 118 compteurmax=0; 117 119 compteurmin=0; 118 120 119 for a=1:29; % boucle sur les annees 120 clear delta; 121 % boucle sur les annees 122 for a=1:29; 123 clear delta; 121 124 delta=zeros(122,1); 122 for b=1:122; % boucle sur les jours : on prepare avant 125 % boucle sur les jours : on prepare avant 126 for b=1:122; 123 127 if b>1; 124 128 delta(b)=olr_eof(b,a)-olr_eof(b-1,a); 125 129 end; 126 130 end; 127 for b=1:122; % boucle sur les jours 128 if (b>20 && b<112); % 1ere condition (laisser une marge de 10) 129 130 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 131 % boucle sur les jours 132 for b=1:122; 133 % 1ere condition (laisser une marge de 10) 134 if (b>20 && b<112); 135 136 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 131 137 compteurmax=compteurmax+1; 132 138 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); … … 163 169 olr_filtre_compomax(compteurmax,30,:,:)=olr1(b+9,a,:,:); 164 170 olr_filtre_compomax(compteurmax,31,:,:)=olr1(b+10,a,:,:); 165 171 166 172 elseif (delta(b+1)>0 && delta(b)<0 && olr_eof(b,a)<(-olrbis)); 167 173 compteurmin=compteurmin+1; … … 202 208 end; 203 209 end; 204 end; 210 end; 205 211 end; 206 212 % on fait le test de Student … … 220 226 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 221 227 % on obtent un tableau de la forme (jours-avant-apres,lon,lat) 222 olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); clear olr_filtre_compomin; 223 olr_filtre_compomax1(:,:,:)=nanmean(olr_filtre_compomax); clear olr_filtre_compomax; 224 225 x=[-10:1:10]'; %' 226 227 figure(1); orient('landscape'); 228 229 lat=lat'; 228 olr_filtre_compomin1(:,:,:)=nanmean(olr_filtre_compomin); 229 clear olr_filtre_compomin; 230 olr_filtre_compomax1(:,:,:)=nanmean(olr_filtre_compomax); 231 clear olr_filtre_compomax; 232 233 x=[-10:1:10]'; 234 235 figure(1); 236 orient('landscape'); 237 238 lat=lat'; 230 239 lon=lon-180; 231 240 for e=1:7; … … 258 267 delta_olr_filtre_min(1,2)=50.0001; 259 268 fin=-30+24*2.5; 260 269 261 270 clear fin 262 271 fin=-10+16*2.5; … … 266 275 else 267 276 set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 268 end 277 end 269 278 set (gca,'YTick',[-10:10:30],'YTickLabel',['10S';' 0';'10N';'20N';'30N'],'fontname','Arial','fontsize',6); 270 279 if e==4; 271 co=colorbar; 280 co=colorbar; 272 281 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 273 282 end; 274 hold on; cartemonde1; 283 hold on; 284 cartemonde1; 275 285 set (gca,'tickDir','out'); 276 286 axis equal 277 287 axis([-180 180 -20 40]); 278 if e==7; 279 title('OLR composite EOF2 30100j simulation CTL ','fontsize',8) ;288 if e==7; 289 title('OLR composite EOF2 30100j simulation CTL ','fontsize',8) ; 280 290 end 281 291 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 282 292 pos_vert=e/10; 283 set (gca,'position',[0.2 pos_vert 0.5 0.1]) 293 set (gca,'position',[0.2 pos_vert 0.5 0.1]) 284 294 end; 285 295 286 figure(1); print -depsc2 composite_olr_EOF2_30100j_simulation_CTL_part1.eps; 287 288 figure(2); orient('landscape'); 296 figure(1); 297 print -depsc2 composite_olr_EOF2_30100j_simulation_CTL_part1.eps; 298 299 figure(2); 300 orient('landscape'); 289 301 290 302 for e=1:7; … … 317 329 delta_olr_filtre_min(1,2)=50.0001; 318 330 fin=-30+24*2.5; 319 331 320 332 clear fin 321 333 fin=-10+16*2.5; … … 325 337 else 326 338 set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 327 end 339 end 328 340 set (gca,'YTick',[-10:10:30],'YTickLabel',['10S';' 0';'10N';'20N';'30N'],'fontname','Arial','fontsize',6); 329 341 if e==4; 330 co=colorbar; 342 co=colorbar; 331 343 set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 332 344 end; 333 hold on; cartemonde1; 345 hold on; 346 cartemonde1; 334 347 set (gca,'tickDir','out'); 335 348 axis equal … … 340 353 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 341 354 pos_vert=e/10; 342 set (gca,'position',[0.2 pos_vert 0.5 0.1]) 355 set (gca,'position',[0.2 pos_vert 0.5 0.1]) 343 356 end; 344 357 345 figure(2); print -depsc2 composite_olr_EOF2_30100j_simulation_CTL_part2.eps; 346 347 348 358 figure(2); 359 print -depsc2 composite_olr_EOF2_30100j_simulation_CTL_part2.eps;
Note: See TracChangeset
for help on using the changeset viewer.