[17] | 1 | % |
---|
[21] | 2 | |
---|
[17] | 3 | clear; |
---|
| 4 | close all; |
---|
[25] | 5 | |
---|
[24] | 6 | status=load('eof_olr_noaa_annuel_jas_1020.mat','C','olr'); |
---|
| 7 | C=status.C; |
---|
| 8 | olr=status.olr; |
---|
| 9 | clear status; |
---|
[17] | 10 | |
---|
| 11 | olr=olr'; |
---|
| 12 | PC1=olr*C; |
---|
| 13 | |
---|
| 14 | %%%%%%% cronique de la VP1 |
---|
| 15 | time_serie_olr1=PC1(:,325); |
---|
| 16 | time_serie_olr2=PC1(:,324); |
---|
| 17 | time_serie_olr3=PC1(:,323); |
---|
| 18 | |
---|
| 19 | dt_olr1=std(time_serie_olr1); |
---|
| 20 | dt_olr2=std(time_serie_olr2); |
---|
| 21 | dt_olr3=std(time_serie_olr3); |
---|
| 22 | |
---|
| 23 | moy_olr1=mean(time_serie_olr1); |
---|
| 24 | moy_olr2=mean(time_serie_olr2); |
---|
| 25 | moy_olr3=mean(time_serie_olr3); |
---|
| 26 | |
---|
| 27 | orient('landscape') |
---|
| 28 | subplot(2,2,1); |
---|
| 29 | for a=1:28; |
---|
| 30 | time_serie_olr1(a)=(time_serie_olr1(a)-moy_olr1)/dt_olr1; |
---|
| 31 | time_serie_olr2(a)=(time_serie_olr2(a)-moy_olr2)/dt_olr2; |
---|
| 32 | time_serie_olr3(a)=(time_serie_olr3(a)-moy_olr3)/dt_olr3; |
---|
[25] | 33 | end; |
---|
[17] | 34 | |
---|
[23] | 35 | status=save('eof1_noaa_annuel_20.txt','time_serie_olr1','-ASCII'); |
---|
| 36 | status=save('eof2_noaa_annuel_20.txt','time_serie_olr2','-ASCII'); |
---|
| 37 | status=save('eof3_noaa_annuel_20.txt','time_serie_olr3','-ASCII'); |
---|
[17] | 38 | |
---|
| 39 | x=[1:28]'; |
---|
| 40 | koro1=time_serie_olr1; |
---|
| 41 | koro2=time_serie_olr2; |
---|
| 42 | koro3=time_serie_olr3; |
---|
| 43 | |
---|
| 44 | clear x; |
---|
| 45 | x=[1979:2006]'; |
---|
| 46 | aa=plot(x,-koro1,'r'); |
---|
| 47 | hold on |
---|
| 48 | bb=plot(x,-koro2,'k'); |
---|
| 49 | cc=plot(x,-koro3,'g'); |
---|
| 50 | dd=legend('eof1','eof2','eof3'); |
---|
| 51 | legend boxoff |
---|
| 52 | set (dd,'Fontsize',8); |
---|
| 53 | y=zeros(28,1); |
---|
| 54 | plot(x,y,'k:'); |
---|
| 55 | title('serie temporelle eof NOAA annuel domaine 10S-20N'); |
---|
| 56 | axis([1979 2006 -3 3]); |
---|
| 57 | set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); |
---|
| 58 | set (gca,'XTick',[1980:5:2005],'XTickLabel',[1980:5:2005],'fontname','Arial','fontsize',6); |
---|
| 59 | set (gca,'tickDir','out'); |
---|
[24] | 60 | print('-depsc2','time_serie_eof_noaa_annuel_20.eps'); |
---|