Changeset 37


Ignore:
Timestamp:
01/27/09 17:14:49 (15 years ago)
Author:
pinsard
Message:

parametrisation of carte_eof234_jjas_....m

Location:
trunk/src/mode_sahelien
Files:
4 edited
1 copied

Legend:

Unmodified
Added
Removed

  • trunk/src/mode_sahelien/carte_eof234_jjas_1025.m

    r26 r37  
    11% 
     2 
     3%+ 
     4% 
     5% EVOLUTIONS 
     6% =========== 
     7% 
     8% $Id$ 
     9%  
     10% - fplod 2009-01-27T14:47:28Z aedon.locean-ipsl.upmc.fr (Darwin) 
     11% 
     12%   * replaced by carte_eof234_jjas_df.m 
     13% 
     14%- 
    215 
    316clear; 

  • trunk/src/mode_sahelien/carte_eof234_jjas_1030.m

    r30 r37  
    11% 
     2 
     3%+ 
     4% 
     5% EVOLUTIONS 
     6% =========== 
     7% 
     8% $Id$ 
     9% 
     10% - fplod 2009-01-27T14:47:28Z aedon.locean-ipsl.upmc.fr (Darwin) 
     11% 
     12%   * replaced by carte_eof234_jjas_df.m 
     13% 
     14%- 
     15 
    216clear; 
    317status=load('eof_1030_jjas.mat','C','L'); 

  • trunk/src/mode_sahelien/carte_eof234_jjas_df.m

    r26 r37  
    1 % 
    2  
    3 clear; 
    4 status=load('eof_1025_jjas.mat','C','L'); 
     1function carte_eof234_jjas_df(ndmin, ndmax) 
     2% carte_eof234_jjas_df ++explication between ndmin days and ndmax days 
     3 
     4%+ 
     5% module 
     6% ====== 
     7% 
     8% ++ 
     9% 
     10% DESCRIPTION 
     11% =========== 
     12% 
     13% plot ++ from file eof_10-30d_jjas.mat 
     14% 
     15% EXAMPLES 
     16% ======== 
     17% 
     18% Following line read eof_10-30d_jjas.mat, plot ++and produce a Postcript file 
     19% resultats_eof_jjas_2006_10-30d.eps  
     20% produce also carte_eofx_mode_sahel_10-30d.dat and  
     21% carte_eofx_mode_sahel_10-30d.ctl (GRADS files) and vecteurs_eof_10-30d.txt. 
     22% :: 
     23% 
     24% >> carte_eof234_jjas_df(10,30); 
     25% 
     26% SEE ALSO 
     27% ======== 
     28% 
     29% olr_eof_1025.m_ 
     30% 
     31% .. _olr_eof_1025.m : olr_eof_1025.m.html 
     32% 
     33% reconstitution_1025.m_ 
     34% 
     35% .. _reconstitution_1025.m : reconstitution_1025.m.html 
     36% 
     37% time_serie_jjas_10_25.m_ 
     38% 
     39% .. _time_serie_jjas_10_25.m : time_serie_jjas_10_25.m.html 
     40% 
     41% sauvegrads.m_ 
     42% 
     43% .. _sauvegrads.m : sauvegrads.m.html 
     44% 
     45% TODO 
     46% ==== 
     47% 
     48% improve description 
     49% 
     50% improve file pb 
     51% 
     52% use return function 
     53% 
     54% change filename of Postscript file because same as in carte_eof_ER_jjas_df.m 
     55% 
     56% directory for sauvegrads files 
     57% 
     58% control of sauvegrads call 
     59% 
     60% why ``-carte`` in contourf and sauvegrads instead of ``carte`` ? 
     61% 
     62% EVOLUTIONS 
     63% ========== 
     64% 
     65% $Id$ 
     66% 
     67% - fplod 2009-01-27T14:49:51Z aedon.locean-ipsl.upmc.fr (Darwin) 
     68% 
     69%   * created from carte_eof234_jjas_1025.m and carte_eof234_jjas_1030.m 
     70%     to replace both. 
     71% 
     72%     Days will be now written ``ndmin``\-``ndmaxd`` (ex: 10-30d) in filenames. 
     73% 
     74%   *  add sauvegrads on first EOF for homogenity 
     75% 
     76%- 
     77 
     78 
     79global IRCAAM_ID; 
     80global IRCAAM_OD; 
     81 
     82if nargin==2 
     83 % ndmin must be lower than ndmax 
     84 if (ndmin > ndmax) 
     85  disp(['ndmin = ', int2str(ndmin)]); 
     86  disp(['ndmax = ', int2str(ndmax)]); 
     87  error('eee : ndmin must be lower than ndmax'); 
     88 end 
     89else 
     90  usage='usage : carte_eof234_jjas_df(ndmin, ndmax)'; 
     91  disp(usage); 
     92  error('eee : wrong arguments numbers') 
     93end 
     94 
     95% ndmin and ndmax will be used in xx-yyd form in filenames 
     96df = [int2str(ndmin),'-',int2str(ndmax),'d']; 
     97clear ndmin; 
     98clear ndmax; 
     99 
     100 
     101fullfilename=[IRCAAM_ID,'eof_',df,'_jjas.mat']; 
     102status=load(fullfilename,'C','L'); 
     103clear fullfilename; 
     104 
    5105C=status.C; 
    6106L=status.L; 
    … …  
    34134 
    35135contourf(lon,lat,-carte,[-0.2:0.02:0.2],'LineStyle','none'); 
     136gradsname=['carte_eof1_mode_sahel_',df]; 
     137sauvegrads(gradsname,-carte,[-30 2.5 -10 2.5]); 
     138clear gradsname; 
     139 
    36140colorbar; 
    37141title('OLR EOF1','Fontsize',8); 
    … …  
    65169carte(1,2)=0.2001; 
    66170 
    67 contourf(lon,lat,carte,[-0.2:0.02:0.2],'LineStyle','none'); 
     171contourf(lon,lat,-carte,[-0.2:0.02:0.2],'LineStyle','none'); 
     172gradsname=['carte_eof2_mode_sahel_',df]; 
     173sauvegrads(gradsname,-carte,[-30 2.5 -10 2.5]); 
     174clear gradsname; 
    68175colorbar; 
    69176 
    … …  
    100207 
    101208contourf(lon,lat,-carte,[-0.2:0.02:0.2],'LineStyle','none'); 
     209gradsname=['carte_eof3_mode_sahel_',df]; 
     210sauvegrads(gradsname,-carte,[-30 2.5 -10 2.5]); 
     211clear gradsname; 
    102212colorbar; 
    103213title('OLR EOF3','Fontsize',8); 
    … …  
    133243 
    134244contourf(lon,lat,-carte,[-0.2:0.02:0.2],'LineStyle','none'); 
     245gradsname=['carte_eof4_mode_sahel_',df]; 
     246sauvegrads(gradsname,-carte,[-30 2.5 -10 2.5]); 
     247clear gradsname; 
    135248colorbar; 
    136249title('OLR EOF4','Fontsize',8); 
    … …  
    154267 
    155268bar(x,vecteurs,'LineStyle','none') 
     269fullfilename=[IRCAAM_OD,'vecteurs_eof_',df,'.txt']; 
     270save(fullfilename,'vecteurs','-ASCII'); 
     271clear fullfilename; 
     272 
    156273axis([0.5 20.5 0 25]); 
    157274set (gca,'YTick',[0:5:25],'YTickLabel',[0:5:25],'fontname','Arial','fontsize',6); 
    … …  
    159276set (gca,'tickDir','out'); 
    160277 
    161 print('-depsc2','resultats_eof_jjas_2006_1025.eps'); 
     278fullfilename=[IRCAAM_OD,'resultats_eof_jjas_2006_',df,'.eps']; 
     279print('-depsc2',fullfilename); 
     280clear fullfilename; 

  • trunk/src/mode_sahelien/carte_eof_jja_2006_1025.m

    r26 r37  
    11% 
     2 
     3%+ 
     4% 
     5% EVOLUTIONS 
     6% =========== 
     7% 
     8% $Id$ 
     9%  
     10% - fplod 2009-01-27T16:10:14Z aedon.locean-ipsl.upmc.fr (Darwin) 
     11% 
     12%   * replaced by carte_eof_ER_jjas_df.m 
     13% 
     14%- 
     15 
    216 
    317clear; 

  • trunk/src/mode_sahelien/carte_eof_jja_2006_2590.m

    r26 r37  
    11% 
     2 
     3%+ 
     4% 
     5% EVOLUTIONS 
     6% =========== 
     7% 
     8% $Id$ 
     9%  
     10% - fplod 2009-01-27T16:10:14Z aedon.locean-ipsl.upmc.fr (Darwin) 
     11% 
     12%   * replaced by carte_eof_ER_jjas_df.m 
     13% 
     14%- 
     15 
    216clear; 
    317close all; 
Note: See TracChangeset for help on using the changeset viewer.