Changeset 42


Ignore:
Timestamp:
01/28/09 13:00:57 (15 years ago)
Author:
pinsard
Message:

parametrisation of composite_olr_phase....m; some unparametrized files in obsolete directory

Location:
trunk/src/mode_sahelien
Files:
11 edited
1 copied
8 moved

Legend:

Unmodified
Added
Removed

  • trunk/src/mode_sahelien/composite_olr_phase.m

    r24 r42  
     1function composite_olr_phase() 
     2% composite_olr_phase ++ 
     3%+ 
     4% module 
     5% ====== 
    16% 
     7% ++ 
     8% 
     9% DESCRIPTION 
     10% =========== 
     11% 
     12% From ``olr.day.mean.nc``, composite_olr_phase compute ++,  
     13% roduce ++ and plot ++ 
     14% 2 Postscript files time_serie_type_tps_phase14.eps and  
     15% time_serie_type_tps_phase58.eps are written. 
     16% 
     17% CAUTIONS 
     18% ======== 
     19% 
     20% By heritage of composite_olr_phase_num.m_ : 
     21% 
     22% Need Statistics_Toolbox because of nanmean 
     23% 
     24% can't work with octave because of ncload and nanmean 
     25% 
     26% 
     27% EXAMPLES 
     28% ======== 
     29% 
     30% :: 
     31% 
     32% >> composite_olr_phase(); 
     33% 
     34% SEE ALSO 
     35% ======== 
     36% 
     37% lagged_correlations_eof2_eof3_1030.m_ 
     38% 
     39% .. _lagged_correlations_eof2_eof3_1030.m : lagged_correlations_eof2_eof3_1030.m.html 
     40% 
     41% composite_olr_phase_num.m_ 
     42% 
     43% .. _composite_olr_phase_num.m : composite_olr_phase_num.m.html 
     44% 
     45% time_serie_histog_phase.m_ 
     46% 
     47% .. _time_serie_histog_phase.m : time_serie_histog_phase.m.html 
     48% 
     49% time_serie_histog_phase_num.m_ 
     50% 
     51% .. _time_serie_histog_phase_num.m : time_serie_histog_phase_num.m.html 
     52% 
     53% TODO 
     54% ==== 
     55% 
     56% improve description 
     57% 
     58% improve file pb 
     59% 
     60% use return function 
     61% 
     62% By heritage of composite_olr_phase_num.m_ : 
     63% 
     64% test at the beginning if Statistics_Toolbox  available 
     65% 
     66% test at the beginning if octave is running 
     67% 
     68% EVOLUTIONS 
     69% ========== 
     70% 
     71% $Id$ 
     72% 
     73% - fplod 2009-01-28T09:35:17Z aedon.locean-ipsl.upmc.fr (Darwin) 
     74%   
     75%   * replace calls to composite_olr_phase[1-8] by calls to 
     76%     composite_olr_phase_num 
     77%   * add comments in ReST 
     78% 
     79%- 
     80close all; 
    281 
    3 close all; 
     82global IRCAAM_OD; 
     83 
    484figure(1); 
    585orient('landscape'); 
    6 subplot(4,1,1); 
    7 composite_olr_phase1; 
    8 text(-240,10,'Phase1'); 
    9 subplot(4,1,2); 
    10 composite_olr_phase2; 
    11 text(-240,10,'Phase2'); 
    12 subplot(4,1,3); 
    13 composite_olr_phase3; 
    14 text(-240,10,'Phase3'); 
    15 subplot(4,1,4); 
    16 composite_olr_phase4; 
    17 text(-240,10,'Phase4'); 
     86phase_num=1; 
     87subplot(4,1,mod(phase_num-1,4)+1); 
     88composite_olr_phase_num(phase_num); 
     89title=['Phase', int2str(phase_num)]; 
     90text(-240,10,title); 
     91phase_num=phase_num+1; 
     92subplot(4,1,mod(phase_num-1,4)+1); 
     93composite_olr_phase_num(phase_num); 
     94title=['Phase', int2str(phase_num)]; 
     95text(-240,10,title); 
     96phase_num=phase_num+1; 
     97subplot(4,1,mod(phase_num-1,4)+1); 
     98composite_olr_phase_num(phase_num); 
     99title=['Phase', int2str(phase_num)]; 
     100text(-240,10,title); 
     101phase_num=phase_num+1; 
     102subplot(4,1,mod(phase_num-1,4)+1); 
     103composite_olr_phase_num(phase_num); 
     104title=['Phase', int2str(phase_num)]; 
     105text(-240,10,title); 
     106 
    18107 
    19108figure(2); 
    20109orient('landscape'); 
    21 subplot(4,1,1); 
    22 composite_olr_phase5; 
    23 text(-240,10,'Phase5'); 
    24 subplot(4,1,2); 
    25 composite_olr_phase6; 
    26 text(-240,10,'Phase6'); 
    27 subplot(4,1,3); 
    28 composite_olr_phase7; 
    29 text(-240,10,'Phase7'); 
    30 subplot(4,1,4); 
    31 composite_olr_phase8; 
    32 text(-240,10,'Phase8'); 
     110phase_num=phase_num+1; 
     111subplot(4,1,mod(phase_num-1,4)+1); 
     112composite_olr_phase_num(phase_num); 
     113title=['Phase', int2str(phase_num)]; 
     114text(-240,10,title); 
     115phase_num=phase_num+1; 
     116subplot(4,1,mod(phase_num-1,4)+1); 
     117composite_olr_phase_num(phase_num); 
     118title=['Phase', int2str(phase_num)]; 
     119text(-240,10,title); 
     120phase_num=phase_num+1; 
     121subplot(4,1,mod(phase_num-1,4)+1); 
     122composite_olr_phase_num(phase_num); 
     123title=['Phase', int2str(phase_num)]; 
     124text(-240,10,title); 
     125phase_num=phase_num+1; 
     126subplot(4,1,mod(phase_num-1,4)+1); 
     127composite_olr_phase_num(phase_num); 
     128title=['Phase', int2str(phase_num)]; 
     129text(-240,10,title); 
     130 
     131clear phase_num; 
    33132 
    34133figure(1); 
    35 print('-depsc2','composite_olr_phase_14.eps'); 
     134fullfilename=[IRCAAM_OD,'composite_olr_phase_14.eps']; 
     135print('-depsc2',fullfilename); 
     136clear fullfilename; 
     137 
    36138figure(2); 
    37 print('-depsc2','composite_olr_phase_58.eps'); 
     139fullfilename=[IRCAAM_OD,'composite_olr_phase_58.eps']; 
     140print('-depsc2',fullfilename); 
     141clear fullfilename; 

  • trunk/src/mode_sahelien/composite_olr_phase1.m

    r26 r42  
    11% 
     2%+ 
     3% 
     4% EVOLUTIONS 
     5% =========== 
     6% 
     7% $Id$ 
     8%  
     9% - fplod 2009-01-28T10:14:38Z aedon.locean-ipsl.upmc.fr (Darwin) 
     10% 
     11%   * replaced by composite_olr_phase_num.m 
     12% 
     13%- 
    214 
    315% initialisation 

  • trunk/src/mode_sahelien/composite_olr_phase2.m

    r26 r42  
    11% 
     2%+ 
     3% 
     4% EVOLUTIONS 
     5% =========== 
     6% 
     7% $Id$ 
     8% 
     9% - fplod 2009-01-28T10:14:38Z aedon.locean-ipsl.upmc.fr (Darwin) 
     10% 
     11%   * replaced by composite_olr_phase_num.m 
     12% 
     13%- 
    214 
    315% initialisation 

  • trunk/src/mode_sahelien/composite_olr_phase3.m

    r26 r42  
    11% 
     2%+ 
     3% 
     4% EVOLUTIONS 
     5% =========== 
     6% 
     7% $Id$ 
     8% 
     9% - fplod 2009-01-28T10:14:38Z aedon.locean-ipsl.upmc.fr (Darwin) 
     10% 
     11%   * replaced by composite_olr_phase_num.m 
     12% 
     13%- 
    214 
    315% initialisation 

  • trunk/src/mode_sahelien/composite_olr_phase4.m

    r26 r42  
    11% 
     2%+ 
     3% 
     4% EVOLUTIONS 
     5% =========== 
     6% 
     7% $Id$ 
     8% 
     9% - fplod 2009-01-28T10:14:38Z aedon.locean-ipsl.upmc.fr (Darwin) 
     10% 
     11%   * replaced by composite_olr_phase_num.m 
     12% 
     13%- 
    214 
    315% initialisation 

  • trunk/src/mode_sahelien/composite_olr_phase5.m

    r26 r42  
    11% 
     2%+ 
     3% 
     4% EVOLUTIONS 
     5% =========== 
     6% 
     7% $Id$ 
     8% 
     9% - fplod 2009-01-28T10:14:38Z aedon.locean-ipsl.upmc.fr (Darwin) 
     10% 
     11%   * replaced by composite_olr_phase_num.m 
     12% 
     13%- 
    214 
    315% initialisation 

  • trunk/src/mode_sahelien/composite_olr_phase6.m

    r26 r42  
    11% 
     2%+ 
     3% 
     4% EVOLUTIONS 
     5% =========== 
     6% 
     7% $Id$ 
     8% 
     9% - fplod 2009-01-28T10:14:38Z aedon.locean-ipsl.upmc.fr (Darwin) 
     10% 
     11%   * replaced by composite_olr_phase_num.m 
     12% 
     13%- 
    214 
    315% initialisation 

  • trunk/src/mode_sahelien/composite_olr_phase7.m

    r26 r42  
    11% 
     2%+ 
     3% 
     4% EVOLUTIONS 
     5% =========== 
     6% 
     7% $Id$ 
     8% 
     9% - fplod 2009-01-28T10:14:38Z aedon.locean-ipsl.upmc.fr (Darwin) 
     10% 
     11%   * replaced by composite_olr_phase_num.m 
     12% 
     13%- 
    214 
    315% initialisation 

  • trunk/src/mode_sahelien/composite_olr_phase8.m

    r26 r42  
    11% 
     2%+ 
     3% 
     4% EVOLUTIONS 
     5% =========== 
     6% 
     7% $Id$ 
     8% 
     9% - fplod 2009-01-28T10:14:38Z aedon.locean-ipsl.upmc.fr (Darwin) 
     10% 
     11%   * replaced by composite_olr_phase_num.m 
     12% 
     13%- 
    214 
    315% initialisation 

  • trunk/src/mode_sahelien/composite_olr_phase_num.m

    r38 r42  
    1 % 
    2  
    3 % initialisation 
    4 clear; 
    5  
    6 ncload('olr.day.mean.nc'); 
     1function composite_olr_phase_num(phase_num) 
     2% composite_olr_phase_num ++explication between 1 and 8 
     3 
     4%+ 
     5% module 
     6% ====== 
     7% 
     8% ++ 
     9% 
     10% DESCRIPTION 
     11% =========== 
     12% 
     13% From ``olr.day.mean.nc`` and phase\ *phase_num*\ .txt, composite_olr_phase_num  
     14% compute ++, plot ++ 
     15% 
     16% CAUTIONS 
     17% ======== 
     18% 
     19% Need Statistics_Toolbox because of nanmean 
     20% 
     21% can't work with octave because of ncload and nanmean 
     22% 
     23% EXAMPLES 
     24% ======== 
     25% 
     26% Following line read olr.day.mean.nc and phase\ *phase_num*\ .txt, compute ++, 
     27% plot ++  
     28% :: 
     29%  
     30% >> composite_olr_phase_num(1); 
     31% 
     32% SEE ALSO 
     33% ======== 
     34% 
     35% 
     36% lagged_correlations_eof2_eof3_1030.m_ 
     37% 
     38% .. _lagged_correlations_eof2_eof3_1030.m : lagged_correlations_eof2_eof3_1030.m.html 
     39% 
     40% composite_olr_phase.m_ 
     41% 
     42% .. _composite_olr_phase.m : composite_olr_phase.m.html 
     43% 
     44% time_serie_histog_phase.m_ 
     45% 
     46% .. _time_serie_histog_phase.m : time_serie_histog_phase.m.html 
     47% 
     48% time_serie_histog_phase_num.m_ 
     49% 
     50% .. _time_serie_histog_phase_num.m : time_serie_histog_phase_num.m.html 
     51% 
     52% TODO 
     53% ==== 
     54% 
     55% improve description 
     56% 
     57% improve file pb 
     58% 
     59% use return function 
     60% 
     61% externalize time selection on olr.day.mean.nc 
     62% 
     63% test at the beginning if Statistics_Toolbox  available 
     64% 
     65% test at the beginning if octave is running 
     66% 
     67% EVOLUTIONS 
     68% ========== 
     69% 
     70% $Id$ 
     71% 
     72% - fplod 2009-01-28T10:38:44Z aedon.locean-ipsl.upmc.fr (Darwin) 
     73% 
     74%   * created from composite_olr_phase8.m to replace it as well as 
     75%     composite_olr_phase[1-7].m 
     76% 
     77%- 
     78 
     79global IRCAAM_ID; 
     80global IRCAAM_OD; 
     81 
     82phase_num_min=1; 
     83phase_num_max=8; 
     84 
     85if nargin==1 
     86 % phase_num must be greater or equal to 1 and lower or equal to 8 
     87 if ((phase_num < phase_num_min) || (phase_num > phase_num_max)) 
     88  disp(['phase_num_min = ', int2str(phase_num_min)]); 
     89  disp(['phase_num_max = ', int2str(phase_num_max)]); 
     90  error('eee : phase_num must be between phase_num_min and phase_num_max'); 
     91 end 
     92else 
     93  usage='usage : composite_olr_phase_num(phase_num)'; 
     94  disp(usage); 
     95  error('eee : wrong arguments numbers') 
     96end 
     97 
     98fullfilename=[IRCAAM_ID,'olr.day.mean.nc']; 
     99ncload(fullfilename); 
     100clear fullfilename; 
    7101 
    8102% taille 12327x73x144 1er juin 1974-31 decembre2007 
    … …  
    151245clear olr2; 
    152246 
    153 status=load('phase8.txt'); 
    154 phase8=status; 
     247fullfilename=[IRCAAM_OD,'phase',int2str(phase_num),'.txt']; 
     248status=load(fullfilename); 
     249clear fullfilename; 
     250 
     251phase=status; 
    155252clear status; 
    156 phase8=reshape(phase8,122,28,3); 
     253phase=reshape(phase,122,28,3); 
    157254olr_filtre_compomax=zeros(90,21,144); 
    158255compteurmax=0; 
    … …  
    164261      % 1ere condition (laisser une marge de 10) 
    165262      if (b>11 && b<112); 
    166          if phase8(b,a,1)==1; 
     263         if phase(b,a,1)==1; 
    167264            compteurmax=compteurmax+1; 
    168265            olr_filtre_compomax(compteurmax,:,:)=olr1(b,a,:,:); 

  • trunk/src/mode_sahelien/time_serie_histog_phase.m

    r41 r42  
    1010% =========== 
    1111% 
    12 % From occu.txt and phase[1-8].txt, time_serie_histog_phase plots 
     12% From ``occu.txt`` and ``phase[1-8].txt``, time_serie_histog_phase plots 
    1313% for each WT and each phase distribution ++. 
    1414% 2 Postscript files time_serie_type_tps_phase14.eps and  
    … …  
    2525% ======== 
    2626% 
     27% lagged_correlations_eof2_eof3_1030.m_ 
    2728% 
    28 % composite_olr_phase1.m_ 
     29% .. _lagged_correlations_eof2_eof3_1030.m : lagged_correlations_eof2_eof3_1030.m.html 
    2930% 
    30 % .. _composite_olr_phase1.m : composite_olr_phase1.m.html 
     31% composite_olr_phase.m_ 
     32% 
     33% .. _composite_olr_phase.m : composite_olr_phase.m.html 
     34% 
     35% composite_olr_phase_num.m_ 
     36% 
     37% .. _composite_olr_phase_num.m : composite_olr_phase_num.m.html 
    3138% 
    3239% time_serie_histog_phase_num.m_ 
    3340% 
    3441% .. _time_serie_histog_phase_num.m : time_serie_histog_phase_num.m.html 
    35 % 
    36 % TODO 
    37 % ==== 
    3842% 
    3943% TODO 

  • trunk/src/mode_sahelien/time_serie_histog_phase_num.m

    r41 r42  
    1111% =========== 
    1212% 
    13 % From occu.txt and phase``phase_num``.txt, time_serie_histog_phase_num plots 
     13% From ``occu.txt`` and phase\ *phase_num*\ .txt, time_serie_histog_phase_num plots 
    1414% for each WT ++. 
    1515% 
    … …  
    2525% ======== 
    2626% 
    27 % 
    28 % composite_olr_phase1.m_ 
    29 % 
    30 % .. _composite_olr_phase1.m : composite_olr_phase1.m.html 
     27% lagged_correlations_eof2_eof3_1030.m_ 
     28% 
     29% .. _lagged_correlations_eof2_eof3_1030.m : lagged_correlations_eof2_eof3_1030.m.html 
     30% 
     31% composite_olr_phase_num.m_ 
     32% 
     33% .. _composite_olr_phase_num.m : composite_olr_phase_num.m.html 
     34% 
     35% composite_olr_phase.m_ 
     36% 
     37% .. _composite_olr_phase.m : composite_olr_phase.m.html 
    3138% 
    3239% time_serie_histog_phase.m_ 
Note: See TracChangeset for help on using the changeset viewer.