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