Changeset 31 for trunk/src/mode_sahelien
- Timestamp:
- 01/27/09 12:44:36 (15 years ago)
- Location:
- trunk/src/mode_sahelien
- Files:
-
- 2 edited
- 1 copied
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/mode_sahelien/time_serie_ER_rossby_1030.m
r30 r31 1 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.m 13 % 14 %- 2 15 3 16 clear; -
trunk/src/mode_sahelien/time_serie_ER_rossby_30100.m
r30 r31 1 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.m 13 % 14 %- 15 2 16 3 17 clear; -
trunk/src/mode_sahelien/time_serie_ER_rossby_df.m
r30 r31 1 function time_serie_ER_rossby_df(ndmin, ndmax) 2 % time_serie_ER_rossby_df ++explication between ndmin days and ndmax days 3 4 %+ 5 % module 6 % ====== 1 7 % 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 table_ER_10-30d.txt file 19 % :: 20 % 21 % >> time_serie_ER_rossby_df(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 % 34 % carte_eof_ER_jjas_1030.m_ 35 % 36 % .. _carte_eof_ER_jjas_1030.m : carte_eof_ER_jjas_1030.m.html 37 % 38 % TODO 39 % ==== 40 % 41 % improve description 42 % improve file pb 43 % use return function 44 % 45 % EVOLUTIONS 46 % ========== 47 % 48 % $Id$ 49 % 50 % - fplod 2009-01-27T10:07:29Z aedon.locean-ipsl.upmc.fr (Darwin) 51 % 52 % * created from time_serie_ER_rossby_1030.m to replace it as well as 53 % time_serie_ER_rossby_30100.m 54 % days will be now written ``ndmin``\-``ndmaxd`` (ex: 10-30d) in filenames 55 % 56 %- 2 57 3 clear;58 global IRCAAM_ID; 4 59 5 status=load('eof_rossby_jjas_2006_1030.mat','C','olr'); 60 if nargin==2 61 % ndmin must be lower than ndmax 62 if (ndmin > ndmax) 63 disp(['ndmin = ', int2str(ndmin)]); 64 disp(['ndmax = ', int2str(ndmax)]); 65 error('eee : ndmin must be lower than ndmax'); 66 end 67 else 68 usage='usage : time_serie_ER_rossby_df(ndmin, ndmax)'; 69 disp(usage); 70 error('eee : wrong arguments numbers') 71 end 72 73 % ndmin and ndmax will be used in xx-yyd form in filenames 74 df = [int2str(ndmin),'-',int2str(ndmax),'d']; 75 clear ndmin; 76 clear ndmax; 77 78 fullfilename=[IRCAAM_ID,'eof_rossby_jjas_2006_',df,'.mat']; 79 status=load(fullfilename,'C','olr'); 80 clear fullfilename; 6 81 C=status.C; 7 82 olr=status.olr; … … 29 104 day=reshape(day,3416,1); 30 105 31 table_ER _1030=zeros(3416,4);106 table_ER=zeros(3416,4); 32 107 for a=1:3416; 33 108 for b=1:4; 34 table_ER _1030(a,1)=year(a);35 table_ER _1030(a,2)=day(a);36 table_ER _1030(a,3)=time_serie_olr1(a);37 table_ER _1030(a,4)=time_serie_olr2(a);109 table_ER(a,1)=year(a); 110 table_ER(a,2)=day(a); 111 table_ER(a,3)=time_serie_olr1(a); 112 table_ER(a,4)=time_serie_olr2(a); 38 113 end; 39 114 end; 40 115 41 save('table_ER_1030.txt','table_ER_1030','-ASCII'); 116 fullfilename=[IRCAAM_ID,'table_ER_',df,'.txt']; 117 save(fullfilename,'table_ER','-ASCII'); 118 clear fullfilename; 42 119 43 120 a=1; … … 52 129 bb=plot(x,-koro2,'k'); 53 130 cc=legend('EOF1','EOF2'); 54 set (cc,'fontsize',6) 131 set (cc,'fontsize',6); 55 132 legend boxoff 56 133 y=zeros(122,1);
Note: See TracChangeset
for help on using the changeset viewer.