Changeset 42 for trunk/src/mode_sahelien/composite_olr_phase_num.m

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

File:

• trunk/src/mode_sahelien/composite_olr_phase_num.m (copied) (copied from trunk/src/mode_sahelien/composite_olr_phase8.m) (3 diffs)

trunk/src/mode_sahelien/composite_olr_phase_num.m

@@ -1 @@
-%
-
-% initialisation
-clear;
-
-ncload('olr.day.mean.nc');
+function composite_olr_phase_num(phase_num)
+% composite_olr_phase_num ++explication between 1 and 8
+
+%+
+% module
+% ======
+%
+% ++
+%
+% DESCRIPTION
+% ===========
+%
+% From ``olr.day.mean.nc`` and phase\ *phase_num*\ .txt, composite_olr_phase_num 
+% compute ++, plot ++
+%
+% CAUTIONS
+% ========
+%
+% Need Statistics_Toolbox because of nanmean
+%
+% can't work with octave because of ncload and nanmean
+%
+% EXAMPLES
+% ========
+%
+% Following line read olr.day.mean.nc and phase\ *phase_num*\ .txt, compute ++,
+% plot ++ 
+% ::
+% 
+% >> composite_olr_phase_num(1);
+%
+% SEE ALSO
+% ========
+%
+%
+% lagged_correlations_eof2_eof3_1030.m_
+%
+% .. _lagged_correlations_eof2_eof3_1030.m : lagged_correlations_eof2_eof3_1030.m.html
+%
+% composite_olr_phase.m_
+%
+% .. _composite_olr_phase.m : composite_olr_phase.m.html
+%
+% time_serie_histog_phase.m_
+%
+% .. _time_serie_histog_phase.m : time_serie_histog_phase.m.html
+%
+% time_serie_histog_phase_num.m_
+%
+% .. _time_serie_histog_phase_num.m : time_serie_histog_phase_num.m.html
+%
+% TODO
+% ====
+%
+% improve description
+%
+% improve file pb
+%
+% use return function
+%
+% externalize time selection on olr.day.mean.nc
+%
+% test at the beginning if Statistics_Toolbox  available
+%
+% test at the beginning if octave is running
+%
+% EVOLUTIONS
+% ==========
+%
+% $Id$
+%
+% - fplod 2009-01-28T10:38:44Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * created from composite_olr_phase8.m to replace it as well as
+%     composite_olr_phase[1-7].m
+%
+%-
+
+global IRCAAM_ID;
+global IRCAAM_OD;
+
+phase_num_min=1;
+phase_num_max=8;
+
+if nargin==1
+ % phase_num must be greater or equal to 1 and lower or equal to 8
+ if ((phase_num < phase_num_min) || (phase_num > phase_num_max))
+  disp(['phase_num_min = ', int2str(phase_num_min)]);
+  disp(['phase_num_max = ', int2str(phase_num_max)]);
+  error('eee : phase_num must be between phase_num_min and phase_num_max');
+ end
+else
+  usage='usage : composite_olr_phase_num(phase_num)';
+  disp(usage);
+  error('eee : wrong arguments numbers')
+end
+
+fullfilename=[IRCAAM_ID,'olr.day.mean.nc'];
+ncload(fullfilename);
+clear fullfilename;
 
 % taille 12327x73x144 1er juin 1974-31 decembre2007
@@ -151 +245 @@
 clear olr2;
 
-status=load('phase8.txt');
-phase8=status;
+fullfilename=[IRCAAM_OD,'phase',int2str(phase_num),'.txt'];
+status=load(fullfilename);
+clear fullfilename;
+
+phase=status;
 clear status;
-phase8=reshape(phase8,122,28,3);
+phase=reshape(phase,122,28,3);
 olr_filtre_compomax=zeros(90,21,144);
 compteurmax=0;
@@ -164 +261 @@
       % 1ere condition (laisser une marge de 10)
       if (b>11 && b<112);
-         if phase8(b,a,1)==1;
+         if phase(b,a,1)==1;
             compteurmax=compteurmax+1;
             olr_filtre_compomax(compteurmax,:,:)=olr1(b,a,:,:);