- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/SIMULS_IRCAAM/composite_olr_AS_eof1_30100.m
r18 r26 5 5 close all; 6 6 olr=zeros(122,29,64,128); 7 cd('INDE') 7 8 8 ncload('rlut_d.AsNQIVIV.197106-09.nc','rlut'); 9 9 olr(:,1,:,:)=rlut; clear rlut; … … 93 93 olr(:,29,:,:)=rlut; clear rlut; 94 94 95 cd('..');96 95 % selction du domaine geographique 97 96 olr1=zeros(122,29,64,128); … … 100 99 clear olr; 101 100 102 load eof1_AS_30100.txt 101 status=load('eof1_AS_30100.txt'); 102 eof1_AS_30100=status; 103 clear status; 103 104 104 105 % on met sous la forme annee mois … … 135 136 if (delta(b+1)<0 && delta(b)>0 && olr_eof(b,a)>(olrbis)); 136 137 compteurmax=compteurmax+1; 137 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:);138 olr_filtre_compomax(compteurmax,1,:,:)=olr1(b-20,a,:,:); 138 139 olr_filtre_compomax(compteurmax,2,:,:)=olr1(b-19,a,:,:); 139 140 olr_filtre_compomax(compteurmax,3,:,:)=olr1(b-18,a,:,:); … … 218 219 y=olr_filtre_compomax(1:nb_value,a,b,c); 219 220 student1(a,b,c)=ttest(x,y,.2); 220 end 221 end 222 end 221 end; 222 end; 223 end; 223 224 224 225 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences … … 259 260 end; 260 261 end; 261 end 262 end; 262 263 subplot(7,1,8-e); 263 264 palette … … 273 274 else 274 275 set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 275 end 276 end; 276 277 set (gca,'YTick',[-10:10:30],'YTickLabel',['10S';' 0';'10N';'20N';'30N'],'fontname','Arial','fontsize',6); 277 278 if e==4; … … 285 286 axis([-180 180 -20 40]); 286 287 if e==7; 287 title('OLR composite EOF1 30100j simulation AS','fontsize',8) ;288 end 288 title('OLR composite EOF1 30100j simulation AS','fontsize',8) ; 289 end; 289 290 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 290 291 pos_vert=e/10; … … 293 294 294 295 figure(1); 295 print -depsc2 composite_olr_EOF1_30100j_simulation_AS_part1.eps;296 print('-depsc2','composite_olr_EOF1_30100j_simulation_AS_part1.eps'); 296 297 297 298 figure(2); … … 321 322 end; 322 323 end; 323 end 324 end; 324 325 subplot(7,1,8-e); 325 326 palette … … 335 336 else 336 337 set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 337 end 338 end; 338 339 set (gca,'YTick',[-10:10:30],'YTickLabel',['10S';' 0';'10N';'20N';'30N'],'fontname','Arial','fontsize',6); 339 340 if e==4; … … 347 348 axis([-180 180 -20 40]); 348 349 if e==7; 349 350 end 350 title('OLR composite EOF1 30100j simulation AS','fontsize',8) ; 351 end; 351 352 text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 352 353 pos_vert=e/10; … … 355 356 356 357 figure(2); 357 print -depsc2 composite_olr_EOF1_30100j_simulation_AS_part2.eps;358 print('-depsc2','composite_olr_EOF1_30100j_simulation_AS_part2.eps');
Note: See TracChangeset
for help on using the changeset viewer.