[32] | 1 | function time_serie_ER_rossby_df_complet(ndmin, ndmax) |
---|
| 2 | % time_serie_ER_rossby_df_complet ++explication between ndmin days and ndmax days |
---|
| 3 | |
---|
| 4 | %+ |
---|
| 5 | % module |
---|
| 6 | % ====== |
---|
[17] | 7 | % |
---|
[32] | 8 | % ++ |
---|
| 9 | % |
---|
| 10 | % DESCRIPTION |
---|
| 11 | % =========== |
---|
| 12 | % |
---|
| 13 | % plot ++ from file eof_rossby_jjas_2006_10-30d.mat |
---|
| 14 | % |
---|
| 15 | % EXAMPLES |
---|
| 16 | % ======== |
---|
| 17 | % |
---|
| 18 | % Following line read eof_rossby_jjas_2006_10-30d.mat, plot ++and produce Postscript files time_serie_eof_ER_10-30d_part*.eps. |
---|
| 19 | % :: |
---|
| 20 | % |
---|
| 21 | % >> time_serie_ER_rossby_df_complet(10,30); |
---|
| 22 | % |
---|
| 23 | % SEE ALSO |
---|
| 24 | % ======== |
---|
| 25 | % |
---|
| 26 | % olr_eof_jjas2006_rossby_1030.m_ |
---|
| 27 | % |
---|
| 28 | % .. _olr_eof_jjas2006_rossby_1030.m : olr_eof_jjas2006_rossby_1030.m.html |
---|
| 29 | % |
---|
| 30 | % carte_eofER_1030_article.m_ |
---|
| 31 | % |
---|
| 32 | % .. _carte_eofER_1030_article.m : carte_eofER_1030_article.m.html |
---|
| 33 | % |
---|
[33] | 34 | % carte_eof_ER_jjas_df.m_ |
---|
[32] | 35 | % |
---|
[33] | 36 | % .. _carte_eof_ER_jjas_df.m : carte_eof_ER_jjas_df.m.html |
---|
[32] | 37 | % |
---|
| 38 | % time_serie_ER_rossby_df.m_ |
---|
| 39 | % |
---|
| 40 | % .. _time_serie_ER_rossby_df.m : time_serie_ER_rossby_df.m.html |
---|
| 41 | % |
---|
| 42 | % TODO |
---|
| 43 | % ==== |
---|
| 44 | % |
---|
| 45 | % improve description |
---|
[33] | 46 | % |
---|
[32] | 47 | % improve file pb |
---|
[33] | 48 | % |
---|
[32] | 49 | % use return function |
---|
| 50 | % |
---|
| 51 | % EVOLUTIONS |
---|
| 52 | % ========== |
---|
| 53 | % |
---|
| 54 | % - fplod 2009-01-27T11:51:31Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
| 55 | % |
---|
| 56 | % * created from time_serie_ER_rossby_1030_complet.m to replace it as well as |
---|
| 57 | % time_serie_ER_rossby_30100_complet.m |
---|
| 58 | % |
---|
[33] | 59 | % Days will be now written ``ndmin``\-``ndmaxd`` (ex: 10-30d) in filenames. |
---|
| 60 | % |
---|
[32] | 61 | %- |
---|
[17] | 62 | |
---|
[32] | 63 | global IRCAAM_ID; |
---|
| 64 | global IRCAAM_OD; |
---|
| 65 | |
---|
| 66 | if nargin==2 |
---|
| 67 | % ndmin must be lower than ndmax |
---|
| 68 | if (ndmin > ndmax) |
---|
| 69 | disp(['ndmin = ', int2str(ndmin)]); |
---|
| 70 | disp(['ndmax = ', int2str(ndmax)]); |
---|
| 71 | error('eee : ndmin must be lower than ndmax'); |
---|
| 72 | end |
---|
| 73 | else |
---|
| 74 | usage='usage : time_serie_ER_rossby_df(ndmin, ndmax)'; |
---|
| 75 | disp(usage); |
---|
| 76 | error('eee : wrong arguments numbers') |
---|
| 77 | end |
---|
| 78 | |
---|
| 79 | % ndmin and ndmax will be used in xx-yyd form in filenames |
---|
| 80 | df = [int2str(ndmin),'-',int2str(ndmax),'d']; |
---|
| 81 | clear ndmin; |
---|
| 82 | clear ndmax; |
---|
| 83 | |
---|
[17] | 84 | close all; |
---|
| 85 | |
---|
| 86 | figure(1); |
---|
| 87 | orient('landscape') |
---|
| 88 | figure(2); |
---|
| 89 | orient('landscape') |
---|
| 90 | figure(3); |
---|
| 91 | orient('landscape') |
---|
| 92 | |
---|
[32] | 93 | fullfilename=[IRCAAM_ID,'eof_rossby_jjas_2006_',df,'.mat']; |
---|
| 94 | status=load(fullfilename,'C','olr'); |
---|
| 95 | clear fullfilename; |
---|
[24] | 96 | C=status.C; |
---|
| 97 | olr=status.olr; |
---|
| 98 | clear status; |
---|
| 99 | |
---|
[17] | 100 | olr=olr'; |
---|
| 101 | PC1=olr*C; |
---|
| 102 | |
---|
| 103 | time_serie_olr1=PC1(:,425); |
---|
| 104 | time_serie_olr2=PC1(:,424); |
---|
| 105 | |
---|
| 106 | for a=1:28; |
---|
| 107 | if a<12; |
---|
| 108 | figure(1); |
---|
| 109 | subplot(4,3,a); |
---|
| 110 | elseif a<23; |
---|
| 111 | figure(2); |
---|
| 112 | subplot(4,3,a-11); |
---|
| 113 | else |
---|
| 114 | figure(3); |
---|
| 115 | subplot(4,3,a-22); |
---|
[21] | 116 | end; |
---|
[17] | 117 | clear koro1 |
---|
| 118 | debut=(a-1)*122+1; |
---|
| 119 | fin=a*122; |
---|
| 120 | koro1=time_serie_olr1(debut:fin); |
---|
| 121 | koro2=time_serie_olr2(debut:fin); |
---|
| 122 | x=[1:122]'; |
---|
| 123 | aa=plot(x,-koro1,'r'); |
---|
| 124 | hold on; |
---|
| 125 | bb=plot(x,-koro2,'k'); |
---|
| 126 | % cc=legend('EOF1','EOF2') |
---|
| 127 | % set (cc,'fontsize',6) |
---|
| 128 | % legend boxoff; |
---|
| 129 | y=zeros(122,1); |
---|
| 130 | plot(x,y,'k:'); |
---|
| 131 | if a==2; |
---|
[32] | 132 | text(0,183,['time series EOF1 and EOF2 OLR ER ' df]); |
---|
[25] | 133 | elseif a==13; |
---|
[32] | 134 | text(0,183,['time series EOF1 and EOF2 OLR ER ' df]); |
---|
[17] | 135 | elseif a==24; |
---|
[32] | 136 | text(0,183,['time series EOF1 and EOF2 OLR ER ' df]); |
---|
[21] | 137 | end; |
---|
[17] | 138 | title(a+1978,'FontSize',8); |
---|
| 139 | set (gca,'tickDir','out'); |
---|
| 140 | axis([1 122 -100 100]); |
---|
| 141 | |
---|
| 142 | set (gca,'YTick',[-100:20:100],'YTickLabel',[-100:20:100],'fontname','Arial','fontsize',6); |
---|
| 143 | set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); |
---|
[21] | 144 | end; |
---|
[17] | 145 | |
---|
| 146 | figure(1); |
---|
[32] | 147 | fullfilename=[IRCAAM_OD,'time_serie_eof_ER_',df,'_part1.eps']; |
---|
| 148 | print('-depsc2',fullfilename); |
---|
[17] | 149 | figure(2); |
---|
[32] | 150 | fullfilename=[IRCAAM_OD,'time_serie_eof_ER_',df,'_part2.eps']; |
---|
| 151 | print('-depsc2',fullfilename); |
---|
[17] | 152 | figure(3); |
---|
[32] | 153 | fullfilename=[IRCAAM_OD,'time_serie_eof_ER_',df,'_part3.eps']; |
---|
| 154 | print('-depsc2',fullfilename); |
---|