1 | % |
---|
2 | |
---|
3 | clear; |
---|
4 | close all; |
---|
5 | |
---|
6 | status=load('eof_olr_annuel.mat','C','olr'); |
---|
7 | C=status.C; |
---|
8 | olr=status.olr; |
---|
9 | clear status; |
---|
10 | |
---|
11 | olr=olr'; |
---|
12 | PC1=olr*C; |
---|
13 | |
---|
14 | %%%%%%% cronique de la VP1 |
---|
15 | time_serie_olr1=PC1(:,425); |
---|
16 | time_serie_olr2=PC1(:,424); |
---|
17 | time_serie_olr3=PC1(:,423); |
---|
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; |
---|
33 | end; |
---|
34 | |
---|
35 | status=save('eof1_noaa_annuel.txt','time_serie_olr1','-ASCII'); |
---|
36 | status=save('eof2_noaa_annuel.txt','time_serie_olr2','-ASCII'); |
---|
37 | status=save('eof3_noaa_annuel.txt','time_serie_olr3','-ASCII'); |
---|
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'); |
---|
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'); |
---|
60 | print('-depsc2','time_serie_eof_noaa_annuel.eps'); |
---|