source: trunk/src/mode_sahelien/obsolete/time_serie_ER_rossby_30100_complet.m @ 42

Last change on this file since 42 was 35, checked in by pinsard, 15 years ago

unparametrized files in obsolete directory

File size: 1.6 KB
Line 
1%
2
3%+
4%
5% EVOLUTIONS
6% ===========
7%
8% $Id$
9%
10% - fplod 2009-01-27T11:17:40Z aedon.locean-ipsl.upmc.fr (Darwin)
11%
12%   * replaced by time_serie_ER_rossby_df_complet.m
13%
14%-
15
16clear;
17close all;
18
19figure(1);
20orient('landscape')
21figure(2);
22orient('landscape')
23figure(3);
24orient('landscape')
25
26status=load('eof_rossby_jjas_2006_30100.mat','C','olr');
27C=status.C;
28olr=status.olr;
29clear status;
30
31olr=olr';
32PC1=olr*C;
33
34time_serie_olr1=PC1(:,425);
35time_serie_olr2=PC1(:,424);
36
37for a=1:28;
38   if a<12;
39      figure(1);
40      subplot(4,3,a);
41   elseif a<23;
42      figure(2);
43      subplot(4,3,a-11);
44   else
45      figure(3);
46      subplot(4,3,a-22);
47   end;
48   clear koro1
49   debut=(a-1)*122+1;
50   fin=a*122;
51   koro1=time_serie_olr1(debut:fin);
52   koro2=time_serie_olr2(debut:fin);
53   x=[1:122]';
54   aa=plot(x,-koro1,'r');
55   hold on;
56   bb=plot(x,-koro2,'k');
57%   cc=legend('EOF1','EOF2')
58%   set (cc,'fontsize',6)
59%   legend boxoff;
60   y=zeros(122,1);
61   plot(x,y,'k:');
62   if a==2;
63       text(0,183,'time series EOF1 and EOF2 OLR ER 30-100d');
64   elseif a==13;
65       text(0,183,'time series EOF1 and EOF2 OLR ER 30-100d');
66   elseif a==24;
67       text(0,183,'time series EOF1 and EOF2 OLR ER 30-100d');
68   end;
69   title(a+1978,'FontSize',8);
70   set (gca,'tickDir','out');
71axis([1 122 -100 100]);
72
73set (gca,'YTick',[-100:20:100],'YTickLabel',[-100:20:100],'fontname','Arial','fontsize',6);
74set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6);
75end;
76
77figure(1);
78print('-depsc2','time_serie_eof_ER_30-100j_part1.eps');
79figure(2);
80print('-depsc2','time_serie_eof_ER_30-100j_part2.eps');
81figure(3);
82print('-depsc2','time_serie_eof_ER_30-100j_part3.eps');
Note: See TracBrowser for help on using the repository browser.