Changeset 31 for trunk


Ignore:
Timestamp:
01/27/09 12:44:36 (15 years ago)
Author:
pinsard
Message:

first parametrisation with days filter

Location:
trunk/src
Files:
3 edited
1 copied

Legend:

Unmodified
Added
Removed

  • trunk/src/ircaam_startup.m

    r24 r31  
    1313% Define : 
    1414%  - path for matlab/octave commands used in IRCAAM project 
    15 %  - IO directories 
     15%  - IO directories in global variables 
    1616% 
    1717% This file can be either run classicaly when currect directory 
    … …  
    4040% $Id$ 
    4141% 
     42% - fplod 2009-01-27T11:25:50Z aedon.locean-ipsl.upmc.fr (Darwin) 
     43% 
     44%   * define global variables 
     45% 
    4246% - fplod 2009-01-07T11:12:31Z aedon.locean-ipsl.upmc.fr (Darwin) 
    4347% 
    … …  
    6468% 
    6569% test if IRCAAM environment set 
    66 ircaam=getenv('IRCAAM'); 
    67 if ( isempty(ircaam) ) 
     70global IRCAAM; 
     71IRCAAM=getenv('IRCAAM'); 
     72if ( isempty(IRCAAM) ) 
    6873 error('eee : IRCAAM environment not set'); 
    6974else 
    70  disp(['iii : IRCAAM =' ircaam]) 
     75 disp(['iii : IRCAAM =' IRCAAM]) 
    7176end 
    7277% 
    7378% add path 
    74 addpath(ircaam); 
    75 addpath([ircaam '/mode_sahelien/']); 
    76 addpath([ircaam '/SIMULS_IRCAAM/']); 
    77 addpath([ircaam '/indice_flore/']); 
     79addpath(IRCAAM); 
     80addpath([IRCAAM '/mode_sahelien/']); 
     81addpath([IRCAAM '/SIMULS_IRCAAM/']); 
     82addpath([IRCAAM '/indice_flore/']); 
    7883% 
    7984% IO directories 
    8085% test if IRCAAM_ID environment set 
    81 ircaam_id=getenv('IRCAAM_ID'); 
    82 if ( isempty(ircaam_id) ) 
     86 
     87global IRCAAM_ID; 
     88IRCAAM_ID=getenv('IRCAAM_ID'); 
     89if ( isempty(IRCAAM_ID) ) 
    8390 error('eee : IRCAAM_ID environment not set'); 
    8491end 
    8592% 
    8693% test if IRCAAM_OD environment set 
    87 ircaam_od=getenv('IRCAAM_OD'); 
    88 if ( isempty(ircaam_od) ) 
     94global IRCAAM_OD; 
     95IRCAAM_OD=getenv('IRCAAM_OD'); 
     96if ( isempty(IRCAAM_OD) ) 
    8997 error('eee : IRCAAM_OD environment not set'); 
    9098end 

  • trunk/src/mode_sahelien/time_serie_ER_rossby_1030.m

    r30 r31  
    11% 
     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%- 
    215 
    316clear; 

  • trunk/src/mode_sahelien/time_serie_ER_rossby_30100.m

    r30 r31  
    11% 
     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 
    216 
    317clear; 

  • trunk/src/mode_sahelien/time_serie_ER_rossby_df.m

    r30 r31  
     1function 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% ====== 
    17% 
     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%- 
    257 
    3 clear; 
     58global IRCAAM_ID; 
    459 
    5 status=load('eof_rossby_jjas_2006_1030.mat','C','olr'); 
     60if 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 
     67else 
     68  usage='usage : time_serie_ER_rossby_df(ndmin, ndmax)'; 
     69  disp(usage); 
     70  error('eee : wrong arguments numbers') 
     71end 
     72 
     73% ndmin and ndmax will be used in xx-yyd form in filenames 
     74df = [int2str(ndmin),'-',int2str(ndmax),'d']; 
     75clear ndmin; 
     76clear ndmax; 
     77 
     78fullfilename=[IRCAAM_ID,'eof_rossby_jjas_2006_',df,'.mat']; 
     79status=load(fullfilename,'C','olr'); 
     80clear fullfilename; 
    681C=status.C; 
    782olr=status.olr; 
    … …  
    29104day=reshape(day,3416,1); 
    30105 
    31 table_ER_1030=zeros(3416,4); 
     106table_ER=zeros(3416,4); 
    32107for a=1:3416; 
    33108     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); 
    38113     end; 
    39114end; 
    40115 
    41 save('table_ER_1030.txt','table_ER_1030','-ASCII'); 
     116fullfilename=[IRCAAM_ID,'table_ER_',df,'.txt']; 
     117save(fullfilename,'table_ER','-ASCII'); 
     118clear fullfilename; 
    42119 
    43120a=1; 
    … …  
    52129bb=plot(x,-koro2,'k'); 
    53130cc=legend('EOF1','EOF2'); 
    54 set (cc,'fontsize',6) 
     131set (cc,'fontsize',6); 
    55132legend boxoff 
    56133y=zeros(122,1); 
Note: See TracChangeset for help on using the changeset viewer.