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