Index: unk/PROGRAMMES/Reg_seriestemporelles.m
===================================================================
--- /trunk/PROGRAMMES/Reg_seriestemporelles.m	(revision 20)
+++ 	(revision )
@@ -1,211 +1,0 @@
-%REG_SERIESTEMPORELLES calcule la regression lineaire d'une serie temporelle issue du champ de donnees initial
-% (selectionne en un point d'espace particulier) et une composante principale choisie      
-%
-% Date: Avril-Mai 2008
-% Auteur: juliette.mignot@locean-ipsl.upmc.fr
-% Contexte Ecole d'ete ENVI-STAT 2008
-%%
-
-%+
-%
-% module
-% ======
-%
-% ``Reg_seriestemporelles``
-%
-% DESCRIPTION
-% ===========
-%
-% ``Reg_seriestemporelles`` lit les fichiers fichiers 
-% ``../DONNEES/lsmask_tropatl_PacandMedblancs.mat``,
-% ``../DONNEES/skt.mon.tropatl.mat``.
-%
-% ``Reg_seriestemporelles`` lit aussi le fichier 
-% ``../DONNEES//EOF_SSTmens_tropAtl30N-20S.mat``
-% produit par eof_NCEP_sais.m_.
-%
-% ``Reg_seriestemporelles`` calcule ++
-%
-% ``Reg_seriestemporelles`` affiche ++
-%
-% ``Reg_seriestemporelles`` sauve cette image dans 
-% ``./EOF_SSTmens_tropAtl30N-20S.mat_Reg_seriestemporelles_1.eps``.
-%
-% EXAMPLES
-% ========
-%
-% ::
-%
-% >> tpacpandreg_startup
-% >> Reg_seriestemporelles
-%
-%
-% SEE ALSO
-% ========
-%
-% tpacpandreg_startup.m_
-%
-% .. _tpacpandreg_startup.m : tpacpandreg_startup.m.html
-%
-% eof_NCEP_sais.m_
-%
-% .. _eof_NCEP_sais.m : eof_NCEP_sais.m.html
-%
-% TODO
-% ====
-%
-% improve description
-%
-% pourquoi ``return`` à la fin
-%
-% ajout xlabel et ylabel
-%
-% EVOLUTIONS
-% ==========
-%
-% $Id$
-%
-% - fplod 2009-08-24T08:46:31Z aedon.locean-ipsl.upmc.fr (Darwin)
-%
-%   * add header
-%   * répertoire de données relatifs au répertoire courant
-%   * suppression de ``clear all`` présent dans ``tpacpandreg_startup.m``
-%   * as octave save default format is ASCII and no default extension 
-%     while matlab save default format is MAT v5 mat-file (little endian) and 
-%     default extension is ``.mat``, write more precise save instruction
-%   * save the figure
-%
-%-
-
-%% PARAMETRES DU PROGRAMME
-% selection de la composante principale a considerer
-ipc=1; %rang de la composante principale que l'on va considerer ici
-
-% les coordonnees du point d'espace auquel on veut comparer cette
-% composante principale doit etre choisi plus bas, apres chargement des
-% longitudes et latitudes
-%%
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%%%%              CHARGEMENT ET PREPARATION DES DONNEES
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-nannees=60;
-
-repertoire=[ '..' filesep 'DONNEES' filesep ]; % a adapter a chaque cas: chemin d acces aux donnees
-
-% Mask
-%-----
-   file='lsmask_tropatl_PacandMedblancs'
-   load([repertoire file]);
-% Data
-%--------
-  repertoire2=[repertoire 'MAT' filesep] % a adapter a chaque cas: chemin d acces aux donnees
-  file=['skt.mon.tropatl.mat']
-  load([repertoire2 file]); % chargement des donnees
-  tab=permute(tab,[2 3 1]); % permutation des dimensions pour organiser le tableau sous la forme [lat lon temps]
-
-% Application du masque
-%----------------------
-mask3D=repmat(mask,[1 1 size(tab,3)]);
-tab(mask3D~=0)=NaN;
-
-
-%  MANIPULATIONS TEMPORELLES: MOYENNES 
-
-% Construction de moyennes mensuelles saisonnieres
-%---------------------------------------------------
-tab=reshape(tab,size(tab,1),size(tab,2),12,nannees);
-tab=squeeze(mean(tab,4));
-
-%Calcul des anomalies
-%tab_mean=squeeze(mean(tab,3));
-%tab_mean=repmat(tab_mean,[1 1 size(tab,3)]);
-%anom=tab-tab_mean;
-
-%%
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%%%% selection d'un point d'espace
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-% le point d'espace est selectionne en reperant son indice i,j dans le
-% tableau de longitude et de latitude
-
-%i=find(lon>-96 & lon<-94);
-%j=find(lat>19 & lat<21);
-i=find(lon>-96 & lon<-94);
-j=find(lat>19 & lat<21);
-%i=find(lon>-21 & lon<-19);
-%j=find(lat>19 & lat<21);
-%i=find(lon>4 & lon<6)
-%j=find(lat>-11 & lat<-9)
-
-SST=squeeze(tab(j,i,:));
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%%%%% chargement des composantes principales
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-file=['EOF_SSTmens_tropAtl30N-20S.mat']; % nom du fichier sous lequel on avait enregistre le resultat de l'ACP en moyenne annuelle
-titledataeof='SST' % je definis ce nom pour pouvoir l'inclure dans le graphique plus bas. 
-if (run_octave == 0)
- load([repertoire file]);
-else
- load('-v7',[repertoire file]);
-end
-PC=PC(:,ipc);
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%%%%               NUAGE DE POINTS
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-ifigure=1;
-figure(ifigure); 
-plot(SST,PC,'x')
-hold on
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%%%%%%   REGRESSION LINEAIRE
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-%[b0,b1,s,R,sigb0,sigb1] = linreg(SST,serie);
-
-SSTmean=mean(SST);PCmean=mean(PC);
-n=length(SST);
-
-
-% la fin du programme est a completer (les formules sont dans le cours ou
-% dans les livres)
-
-% calcul de b1 et b0, coefficient directeur et ordonnee a l origine de la
-% droite de regression
-
-
-%trace de la droite de regression
-%---------------------------------
-
-
-%Calcul du coefficient de détermination
-%----------------------------------------
-
-%significativite statistique comparaison de la correlation r et de la
-%valeur critique rc donnee par la table de student en fonction du nombre de
-%degres de liberte ndf
-%------------------------------
-
-
-% significativite statistique de b0 et b1
-%--------------------------------------------
-
-% sauvegarde de la figure
-printer='eps';
-print_printer=['-d', printer];
-fullfilename=['.' filesep file '_' mfilename '_'  num2str(ifigure) '.' printer];
-print(print_printer,fullfilename);
-clear print_printer;
-clear fullfilename;
-clear printer;
-%
-clear ifigure;
-
-return
Index: /trunk/PROGRAMMES/carte_SST.m
===================================================================
--- /trunk/PROGRAMMES/carte_SST.m	(revision 20)
+++ /trunk/PROGRAMMES/carte_SST.m	(revision 30)
@@ -1,3 +1,3 @@
-% CARTE_SST  Ce programme permet de visualiser les cartes mensuelles moyennes de temperature de surface       
+%CARTE_SST Ce programme permet de visualiser les cartes mensuelles moyennes de temperature de surface
 %
 % Date: Avril-Mai 2008
@@ -15,16 +15,18 @@
 % ===========
 %
-% ``carte_SST`` visualise les cartes mensuelles moyennes de température 
+% ``carte_SST`` visualise les cartes mensuelles moyennes de température
 % de surface.
 %
-% ``carte_SST`` lit les fichiers 
+% ``carte_SST`` lit les fichiers
 % ``../DONNEES/lsmask_tropatl_PacandMedblancs.mat``
 % et ``../DONNEES/skt.mon.tropatl.mat``.
 %
-% ``carte_SST`` calcule ++blabla++
-%
-% ``carte_SST`` affiche à l'écran les 12 mois de ++
-%
-% ``carte_SST`` sauve cette image dans ``./carte_SST.eps``.
+% ``carte_SST`` calcule les moyennes mensuelles de SST dans
+% l'Atlantique tropical a partir d'un champ mensuel de 60 ans
+% fourni en entree
+%
+% ``carte_SST`` affiche à l'écran les moyennes mensuelles de SST.
+%
+% ``carte_SST`` sauve cette image dans ``./carte_SST.ps``.
 %
 % EXAMPLES
@@ -47,4 +49,16 @@
 % .. _initfig.m : initfig.m.html
 %
+% whorldmap.m_
+%
+% .. _whorldmap.m : whorldmap.m.html
+%
+% colorbartype.m_
+%
+% .. _colorbartype.m : colorbartype.m.html
+%
+% contlab.m_
+%
+% .. _contlab.m : contlab.m.html
+%
 % TODO
 % ====
@@ -52,4 +66,20 @@
 % improve description
 %
+% KNOWN PROBLEMS
+% ==============
+%
+% - fplod 2009-08-26T08:46:48Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * ne fonctionne pas sous octave 3.2.2
+%     ::
+%     error: gnuplot (as of v4.2) only supports 3D filled triangular patches
+%     error: called from:
+%     error:   /usr/local_linux/octave/64/3.2.2//share/octave/3.2.2/m/plot/__go_draw_axes__.m at line 563, column 6
+%     error:   /usr/local_linux/octave/64/3.2.2//share/octave/3.2.2/m/plot/__go_draw_figure__.m at line 93, column 3
+%     error:   /usr/local_linux/octave/64/3.2.2//share/octave/3.2.2/m/plot/gnuplot_drawnow.m at line 88, column 5
+%     error:   /usr/local_linux/octave/64/3.2.2//share/octave/3.2.2/m/plot/print.m at line 696, column 7
+%     error:   /.autofs/home/fplod/incas/tpoctave/tpacpandreg_ws/PROGRAMMES/carte_SST.m at line 212, column 1
+%
+%
 % EVOLUTIONS
 % ==========
@@ -57,4 +87,40 @@
 % $Id$
 %
+% - fplod 2009-08-26T09:00:11Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * titre: nom de chaque mois
+%   * remplacement de ``eps`` par ``ps`` pour pouvoir imprimer les figures
+%     produites par ``octave``
+%
+%     Dans le fichier .eps produit par ``octave`` on pouvait lire
+%     ``!PS-Adobe-2.0 EPSF-2.0``.
+%
+%     Aucun message d'erreur n'apparaît suite à la demande d'impression
+%     ``zeus$ lpr -o job-sheets='none' carte_SST.ps``, mais aucune feuille ne
+%     sort.
+% 
+%     Pour info dans le fichier .eps produit par matlab, on pouvait lire
+%     ``%!PS-Adobe-2.0 EPSF-1.2``. Aucun souci d'impression avec ce fichier !
+%
+%     En remplaçant ``eps`` par ``ps``, on n'a plus de souci d'impression.
+%
+%
+% - jmignot 2009-08-25
+%
+%   * suppression d'appel à contlab et colorbartype sous octave
+% 
+% - fplod 2009-08-25T10:40:16Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * modif commentaires liée au dernier chgt de Juliette
+%
+% - jmignot 2009-08-25
+%
+%   * remplacement de ``colorbar`` par ``colorbartype``
+%     parce pas beau ni en matlab ni en octave (dégradé, noir et blanc dans
+%     les .jpg)
+%   * contraintes sur les limites x et y de axes car matlab et octave ne
+%     fixent pas les même limites par défaut
+%   * utilisation de ``contlab`` pour ++
+% 
 % - fplod 2009-08-17T12:39:51Z aedon.locean-ipsl.upmc.fr (Darwin)
 %
@@ -63,9 +129,9 @@
 %   * suppression de ``clear all`` présent dans ``tpacpandreg_startup.m``
 %   * contournement pb octave pour ``contourf`` :
-%     Dans le fil de discussion [Pkg-octave-devel] Bug#492223: octave3.0: contourf only accepts vector X and Y if same+++, 
-%     je comprends que j'ai affaire à un bug de octave 3.0 (donc à surveiller 
-%     pour les prochaines versions) qui ne veut pas faire de contourf si les 
-%     variables X et Y n'ont pas la même dimension et c'est le cas ici 
-%     (lat = 42x1, lon= 70x1 vue avec size). 
+%     Dans le fil de discussion [Pkg-octave-devel] Bug#492223: octave3.0: contourf only accepts vector X and Y if same+++,
+%     je comprends que j'ai affaire à un bug de octave 3.0 (donc à surveiller
+%     pour les prochaines versions) qui ne veut pas faire de contourf si les
+%     variables X et Y n'ont pas la même dimension et c'est le cas ici
+%     (lat = 42x1, lon= 70x1 vue avec size).
 %     http://lists.alioth.debian.org/pipermail/pkg-octave-devel/2008-July/004672.html
 %   * sauvegarde de la figure dans ``./carte_SST.eps``
@@ -85,5 +151,5 @@
 %Mask
 %-----
-   file='lsmask_tropatl_PacandMedblancs'
+   file='lsmask_tropatl_PacandMedblancs.mat'
    load([repertoire file]);
 % Data
@@ -101,5 +167,5 @@
 
 %---------------------------------------------------
-%  MANIPULATIONS TEMPORELLES: MOYENNES 
+%  MANIPULATIONS TEMPORELLES: MOYENNES
 %---------------------------------------------------
 tab=reshape(tab,size(tab,1),size(tab,2),12,nannees); %organisation du tableau sous la forme [lat lon mois annees]
@@ -114,9 +180,24 @@
 %%%%               TRACE
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-
+%parametres geographiques du graphe
+  geo=[-100 20 -30 30]; %[lonmin lonmax latmin latmax];
+  yti=[-30:10:30]; xti=[-100:30:20]; %labels de latitude et de longitude
+
+		
 % Variabilité saisonnière
-figure; 
-initfig; 
+%-------------------------
+%niveaux de couleurs
+lev=[0:1:30];%contours a tracer, en unite du champ a dessiner, ici
+             %la SST en degC
+	     %a adapter eventuellement
+clip_lev=[20 30];%definit les limites des couleurs: bleu fonce au
+                 %rouge fonce
+	     %a adapter eventuellement
+
+
+ifigure=1;
+figure(ifigure);
+clf();
+initfig;
 for iplot=1:12
   subplot(4,3,iplot)
@@ -128,15 +209,43 @@
   end
   z=squeeze(tab_mens(:,:,iplot));
-  contourf(x,y,z,[0:1:30])
-  caxis([20 30])
+  contourf(x,y,z,lev)
+  caxis(clip_lev)
   hold on
-  ylim([-30 30])
-  xlim([-100 20])
-  title(['mois ' num2str(iplot)])
-end; 
-colorbar('East')
+  whorldmap
+
+  % limites géographiques et labels de longitude et de latitude
+  ylim([geo(3) geo(4)])
+  xlim([geo(1) geo(2)])
+  set(gca, 'Xtick', xti, 'Ytick', yti);
+  if (run_octave == 0)
+    [longstr, latstr] = contlab(xti, yti);
+    set(gca, 'Xticklabel', longstr, 'Yticklabel', latstr);
+  end;
+  if iplot>=10
+    xlabel('longitude')
+  end;
+  if iplot==1|iplot==4|iplot==7|iplot==10
+    ylabel('latitude')
+  end;
+
+  %titre: nom de chaque mois
+  cmois = datestr(datenum([2000 iplot 1 00 00 00]), 'mmm');
+  title(cmois);
+  clear cmois;
+
+end;
+
+%barre de couleur
+if (run_octave == 0)
+  pos=[.2 .03 .6 .013];
+  hpal=jet(100); %carte des couleurs: du bleu au rouge par defaut.
+  colorbartype(pos,lev,1,clip_lev,hpal,0)
+  title('SST mensuelle moyenne en degC')
+else
+ colorbar('East') %colorbar par defaut
+end;
 
 % sauvegarde de la figure
-printer='eps';
+printer='ps';
 print_printer=['-d', printer];
 fullfilename=['./', mfilename, '.', printer];
Index: /trunk/PROGRAMMES/colorbartype.m
===================================================================
--- /trunk/PROGRAMMES/colorbartype.m	(revision 30)
+++ /trunk/PROGRAMMES/colorbartype.m	(revision 30)
@@ -0,0 +1,51 @@
+function [h0]=colorbartype(pos,lev,lev_pas,clip_lev,hpal,orien,xxx);
+
+% 30/10/97- NS: option "orien" pour positionner la barre verticalement
+%               ou horizontalement--> 1 ou 0
+% 08/09/00- EK: si les valeurs sont petites, ecriture au format exponentiel
+
+h0 = axes('Position',pos,'Box','on');
+ind_lev=find( lev>=(clip_lev(1)-eps) & (lev<=clip_lev(2)+eps) );
+tab_lev=[lev(ind_lev);lev(ind_lev)];
+
+% horizontale par defaut
+if (nargin < 6), orien = 0; end 
+
+% labels en bas ou a droite par defaut
+if nargin < 7, 
+   if orien == 0, xxx=-.5; end
+   if orien == 1, xxx=4; end
+end
+
+% Format ecriture des labels (expo ou non)
+if abs(clip_lev(end)) <= 1e-2, iexp = 1; else iexp = 0; end
+
+if orien == 0,
+
+    pcolor(tab_lev);
+    set(h0,'xticklabel', '','yticklabel','');
+    abs_bar=0:1/(length(ind_lev)-1):1;
+    for j=1:lev_pas:length(ind_lev)
+        if iexp == 0, s=num2str(lev(ind_lev(j))); else
+                      s = sprintf('%0.1e', lev(ind_lev(j)));  end
+        text(abs_bar(j),xxx,s,'Fontsize', 9, 'HorizontalAlignment', ...
+                              'Center', 'Units', 'normalized');
+    end
+
+else
+
+    pcolor(tab_lev');
+    set(h0,'xticklabel', '','yticklabel','');
+    ord_bar=0:1/(length(ind_lev)-1):1;
+    for j=1:lev_pas:length(ind_lev)
+        if iexp == 0, s=num2str(lev(ind_lev(j))); else
+                      s = sprintf('%0.1e', lev(ind_lev(j)));  end
+        tt=text(xxx,ord_bar(j),s,'Fontsize', 9, ...
+            'HorizontalAlignment','right','Units', 'normalized');
+%        tt=text(xxx,ord_bar(j),s,'Fontsize', 6,'HorizontalAlignment','right','Units', 'normalized');
+    end
+
+end
+colormap(hpal);
+caxis(clip_lev);
+
Index: /trunk/PROGRAMMES/contlab.m
===================================================================
--- /trunk/PROGRAMMES/contlab.m	(revision 30)
+++ /trunk/PROGRAMMES/contlab.m	(revision 30)
@@ -0,0 +1,43 @@
+function [ xti_str, yti_str] = contlab(xti, yti)
+%
+%------------------------------------------------------------------------
+% function [ xti_str yti_str] = contlab(xti, yti)
+%
+% fonction pour labels et ticks sur les axes longitude et latitude
+% longitude --> xti tableau defini par lon0:pas:lon1;
+% latitude  --> yti tableau defini par lat0:pas:lat1;
+%
+% 24/10/96 - EK
+%------------------------------------------------------------------------
+ 
+% Tick labels pour longitude
+set(gca, 'Xticklabel', xti);
+xlab =  get(gca, 'Xticklabel');
+[mxlab nxlab] = size(xlab);
+nxlab_str = int2str(nxlab);
+%xti_str = zeros(mxlab, nxlab+1);
+for iti = 1:mxlab
+    xxx = rem(str2num(xlab(iti,:))+540, 360)-180;
+    if (xxx<0), xti_str(iti,:) = sprintf(['%' nxlab_str '.0fW'],abs(xxx));
+        elseif (xxx>0) ...
+                xti_str(iti,:)=sprintf(['%' nxlab_str '.0fE'],(xxx));
+        else xti_str(iti,:)=sprintf(['%' nxlab_str '.0f '],(xxx));
+    end;
+end;
+
+% Tick labels pour latitude
+set(gca, 'Yticklabel', yti);
+ylab = get(gca, 'Yticklabel');
+[mylab nylab] = size(ylab);
+nylab_str = int2str(nylab);
+%yti_str = zeros(mylab, nylab+1);
+for iti = 1:mylab
+    yyy = str2num(ylab(iti,:));
+    if (yyy<0), yti_str(iti,:)=sprintf(['%' nylab_str '.0fS'],abs(yyy));
+        elseif (yyy>0) ...
+                yti_str(iti,:)=sprintf(['%' nylab_str '.0fN'],(yyy)); 
+        else yti_str(iti,:)=sprintf(['%' nylab_str '.0f '],(yyy));
+    end;
+end;
+
+return
Index: /trunk/PROGRAMMES/eof_NCEP_sais.m
===================================================================
--- /trunk/PROGRAMMES/eof_NCEP_sais.m	(revision 20)
+++ /trunk/PROGRAMMES/eof_NCEP_sais.m	(revision 30)
@@ -8,5 +8,5 @@
 %
 % le programme permet egalement de visualiser les resultats. Voir le doc
-% .pdf joint pour l'interpretation de cette visualisation. 
+% .pdf joint pour l'interpretation de cette visualisation.
 %
 % Date: Avril-Mai 2008
@@ -19,10 +19,10 @@
 % ======
 %
-% ``eof_NCEP_sais`` 
+% ``eof_NCEP_sais``
 %
 % DESCRIPTION
 % ===========
 %
-% ``eof_NCEP_sais`` lit les fichiers 
+% ``eof_NCEP_sais`` lit les fichiers
 % ``../DONNEES/lsmask_tropatl_PacandMedblancs.mat``
 % et ``../DONNEES/skt.mon.tropatl.mat``.
@@ -30,15 +30,16 @@
 % ``eof_NCEP_sais`` calcule ++blabla++
 %
-% ``eof_NCEP_sais`` sauvegarde dans ``../DONNEES/EOF_SSTmens_tropAtl30N-20S.mat``.
-%
-% Ces fichiers sont lus par ``Reg_serietemporelles.m_``\ .
+% ``eof_NCEP_sais`` sauvegarde dans
+% ``../DONNEES/eof_SSTmens_tropAtl30N-20S.mat``.
+%
+% Ces fichiers sont lus par ``reg_serietemporelles.m_``\ .
 %
 % ``eof_NCEP_sais`` affiche trois figures :
-%
 %  - ++
 %  - ++
 %  - valeurs propres de la matrice de covariance
 %
-% ``eof_NCEP_sais`` sauve ces images dans ``./EOF_SSTmens_tropAtl30N-20S_[123].eps``.
+% ``eof_NCEP_sais`` sauve ces images dans
+% ``./eof_SSTmens_tropAtl30N-20S_[123].ps``.
 %
 % ++
@@ -63,11 +64,19 @@
 % .. _initfig.m : initfig.m.html
 %
-% Reg_seriestemporelles.m_
-%
-% .. _Reg_seriestemporelles.m : Reg_seriestemporelles.m.html
+% reg_seriestemporelles.m_
+%
+% .. _reg_seriestemporelles.m : reg_seriestemporelles.m.html
 %
 % whorldmap.m_
 %
 % .. _whorldmap.m : whorldmap.m.html
+%
+% colorbartype.m_
+%
+% .. _colorbartype.m : colorbartype.m.html
+%
+% contlab.m_
+%
+% .. _contlab.m : contlab.m.html
 %
 % TODO
@@ -78,4 +87,6 @@
 % tester la branche NetCDF
 %
+% comprendre le pb octave 3.2.2
+%
 % EVOLUTIONS
 % ==========
@@ -83,12 +94,48 @@
 % $Id$
 %
-% - fplod 2009-08-25T09:37:04Z aedon.locean-ipsl.upmc.fr (Darwin)
-%
-%   * typo in ReStructured Text header
+% - fplod 2009-08-26T16:34:59Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * pb contourf sous octave 3.2.2
+%
+% - fplod 2009-08-26T11:05:36Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * remplacement de ``eps`` par ``ps`` pour pouvoir imprimer les figures
+%     produites par ``octave``
+%
+%     Dans le fichier .eps produit par ``octave`` on pouvait lire
+%     ``!PS-Adobe-2.0 EPSF-2.0``.
+%
+%     Aucun message d'erreur n'apparaît suite à la demande d'impression
+%     ``zeus$ lpr -o job-sheets='none' eof_SSTmens_tropAtl30N-20S_1.octave.ps``,
+%     mais aucune feuille ne sort.
+%
+%     Pour info dans le fichier .eps produit par matlab, on pouvait lire
+%     ``%!PS-Adobe-2.0 EPSF-1.2``. Aucun souci d'impression avec ce fichier !
+%
+%     En remplaçant ``eps`` par ``ps``, on n'a plus de souci d'impression.
+%
+% - jmignot 2009-08-25
+%
+%   * suppression d'appel à contlab et colorbartype sous octave
+%
+%
+% - fplod 2009-08-25T10:40:16Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * modif commentaires liée au dernier chgt de Juliette
+%
+% - jmignot 2009-08-25
+%
+%   * chgt de noms de fichiers : ``EOF_`` devient ``eof_``
+%   * remplacement de ``colorbar`` par ``colorbartype``
+%     parce pas beau ni en matlab ni en octave (dégradé, noir et blanc dans
+%     les .jpg)
+%   * contraintes sur les limites x et y de axes car matlab et octave ne
+%     fixent pas les même limites par défaut
+%   * utisation de ``contlab`` pour ++
 %
 % - fplod 2009-08-21T15:04:16Z aedon.locean-ipsl.upmc.fr (Darwin)
 %
-%   * as octave save default format is ASCII and no default extension 
-%     while matlab save default format is MAT v5 mat-file (little endian) and 
+%   * as octave save default format is ASCII and no default extension
+%     while matlab save default format is MAT v5 mat-file (little endian) and
 %     default extension is ``.mat``, write more precise save instruction
 %
@@ -100,10 +147,11 @@
 %   * contournement pb octave pour ``contourf`` :
 %     Dans le fil de discussion [Pkg-octave-devel] Bug#492223: octave3.0: contourf only accepts vector X and Y if same+++,
-%     je comprends que j'ai affaire à un bug de octave 3.0 (donc à surveiller 
-%     pour les prochaines versions) qui ne veut pas faire de contourf si les 
-%     variables X et Y n'ont pas la même dimension et c'est le cas ici 
-%     (lat = 42x1, lon= 70x1 vue avec size). 
+%     je comprends que j'ai affaire à un bug de octave 3.0 (donc à surveiller
+%     pour les prochaines versions) qui ne veut pas faire de contourf si les
+%     variables X et Y n'ont pas la même dimension et c'est le cas ici
+%     (lat = 42x1, lon= 70x1 vue avec size).
 %     http://lists.alioth.debian.org/pipermail/pkg-octave-devel/2008-July/004672.html
 %   * sauvegarde des figures
+
 %
 %-
@@ -111,6 +159,6 @@
 N=4; %nb d'EOF a tracer
 
-nomfic='EOF_SSTmens_tropAtl'
-% mettre 'EOF_SSTan_tropAtl' quand on travaille en moyennes annuelles. Ce
+nomfic='eof_SSTmens_tropAtl'
+% mettre 'eof_SSTan_tropAtl' quand on travaille en moyennes annuelles. Ce
 % fichier sera rappele dans les programme de regression lineaire
 %%
@@ -125,7 +173,7 @@
 % Mask
 %------
-  file='lsmask_tropatl_PacandMedblancs'
+  file='lsmask_tropatl_PacandMedblancs.mat'
    load([repertoire file]);
-  ny=size(mask,1); nx=size(mask,2); % dimensions spatiales des tableaux 
+  ny=size(mask,1); nx=size(mask,2); % dimensions spatiales des tableaux
 
 % Data
@@ -136,5 +184,5 @@
   load([repertoire2 file]); % chargement des donnees
 
-   
+
   %SI NETCDF
   %repertoire2=[repertoire 'NETCDF' filesep]% a adapter a chaque cas: chemin d acces aux donnees
@@ -142,5 +190,5 @@
   % f=netcdf([repertoire file],'nowrite');
   % tab=f{'skt'}(:,:,:);
-  
+
 
    tab=permute(tab,[2 3 1]);% permutation des dimensions pour organiser le tableau sous la forme [lat lon temps]
@@ -169,5 +217,5 @@
 %----------------
 %Si on veut éventuellement attribuer un poids différent aux différents
-%points spatiaux, on le fait ici. 
+%points spatiaux, on le fait ici.
 
 % Préparation du zoom géographique: sur le masque, j'identifie les points
@@ -177,5 +225,5 @@
 limite_latS=-20; %limite du domaine d'etuade en latitude
 limite_latN=30; %limite du domaine d'etuade en latitude
-lat2D=repmat(lat,[1 size(mask,2)]); 
+lat2D=repmat(lat,[1 size(mask,2)]);
 mask(lat2D>limite_latN&mask==0)=4;
 mask(lat2D<limite_latS&mask==0)=4;
@@ -209,24 +257,24 @@
 % matrice de covariance:
   C = Z * Z';
-%  eval(['save ' repertoire 'MatCov_SST_tropAtl' num2str(limite_latN) 'N' num2str(limite_latS) 'S.asc C -ascii']) 
-  
-% diagonalisation  
-  [E,Lambda] = eig(C); 
-   % Lambda est la matrice diagonale des valeurs propres et E la matrice dont les colonnes sont les vecteurs propres correspondants.  
+%  eval(['save ' repertoire 'MatCov_SST_tropAtl' num2str(limite_latN) 'N' num2str(limite_latS) 'S.asc C -ascii'])
+
+% diagonalisation
+  [E,Lambda] = eig(C);
+   % Lambda est la matrice diagonale des valeurs propres et E la matrice dont les colonnes sont les vecteurs propres correspondants.
    % Soit Z*E = E*Lambda
-   % Les vecteurs propres sont appelés les EOF. Chaque EOF peut être considéré comme une carte de l'espace physique.    
-
-   
- %% 
+   % Les vecteurs propres sont appelés les EOF. Chaque EOF peut être considéré comme une carte de l'espace physique.
+
+
+ %%
  %%%% calcul de la matrice donnant les pourcentages de variance
 Lambda = diag(Lambda);
 pctg_var = round(Lambda / sum(Lambda) * 1000) / 10;
 
- 
+
  %On definit les composantes principales de Q comme la decomposition du signal decrit par Q sur les EOF E.
 
    P = Z' * E;
-    
-%%   
+
+%%
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
@@ -271,35 +319,71 @@
 
  %%               TRACE
-figure(1); 
+
+close all
+%parametres geographiques des graphes
+  geo=[-100 20 -30 30]; %[lonmin lonmax latmin latmax];
+  yti=[-30:10:30]; xti=[-100:30:20]; %labels de latitude et de longitude
+%niveaux de couleurs
+  lev=[-2:.02:2];
+  clip_lev=[-0.06 0.06];
+
+figure(1);
 initfig;
-figure(2); 
+figure(2);
 initfig;
 
 for i=1:4
-    
+
    %subplot(4,2,2*i-1)
    ifigure=1;
    figure(ifigure);
    subplot(2,2,i)
+   z=squeeze(TABEOF(:,:,i));
    if (run_octave == 0)
     x=lon;
     y=lat;
+    contourf(x,y,z,lev);
    else
-    [x,y]=meshgrid(lon,lat);
+    switch OCTAVE_VERSION
+     case {'3.0.2'}
+      [x,y]=meshgrid(lon,lat);
+      contourf(x,y,z,lev);
+     case {'3.2.2'}
+      warning('www : do no know about octave 3.2.2 contourf usage');
+%  A(I): Index exceeds matrix dimension ++
+     otherwise
+      warning('www : do no know about octave contourf usage');
+     end
    end
-   z=squeeze(TABEOF(:,:,i));
-   contourf(x,y,z);
-   caxis([-0.06 0.06])
+   caxis(clip_lev);
    hold on
-   xlim([-100 30]); xlabel('longitude')
-   ylim([-40 40]); ylabel('latitude')
-   whorldmap([-100, 30, -40, 40])
-   colorbar
-   title([num2str(p_var(i)) '%  -  (unites: K)'])
-   set(gca,'xtick',[-100;-70;-40;-10;20],'xticklabel',['100W';' 70W';' 40W';' 10W';' 20E'])
-   set(gca,'ytick',[-30;-15;0;15;30],'yticklabel',['30S';'15S';' 0 ';'15N';'30N'])
-   %
+   whorldmap
+   % limites geographies et labels de longitude et de latitude
+  ylim([geo(3) geo(4)])
+  xlim([geo(1) geo(2)])
+  set(gca, 'Xtick', xti, 'Ytick', yti)
+  if (run_octave == 0)
+    [longstr, latstr] = contlab(xti, yti);
+    set(gca, 'Xticklabel', longstr, 'Yticklabel', latstr)
+  end;
+  xlabel('longitude')
+  ylabel('latitude')
+
+   %titre de la figure
+   title([num2str(p_var(i)) '%'])
+
+  %barre de couleur
+  if (run_octave == 0)
+    hpal=jet(100); %carte des couleurs: du bleu au rouge par defaut.
+    pos1=get(gca,'Position');
+    pos=[pos1(1) pos1(2)-.075 pos1(3) .013];
+    colorbartype(pos,lev,1,clip_lev,hpal,0);
+    title('degC')
+  else
+    %++colorbar('East')
+  end;
+
    % sauvegarde de la figure
-   printer='eps';
+   printer='ps';
    print_printer=['-d', printer];
    fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer];
@@ -307,5 +391,5 @@
    clear fullfilename;
    clear printer;
-   
+
    %subplot(4,2,2*i)
    ifigure=2;
@@ -320,5 +404,5 @@
    %
    % sauvegarde de la figure
-   printer='eps';
+   printer='ps';
    print_printer=['-d', printer];
    fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer];
@@ -330,5 +414,6 @@
 % Graphe pourcentages de variance expliqué
 ifigure=3;
-figure(ifigure); 
+figure(ifigure);
+initfig;
 bar(flipud(pctg_var))
 xlim([0 10])
@@ -338,5 +423,5 @@
 
 % sauvegarde de la figure
-printer='eps';
+printer='ps';
 print_printer=['-d', printer];
 fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer];
Index: /trunk/PROGRAMMES/initfig.m
===================================================================
--- /trunk/PROGRAMMES/initfig.m	(revision 20)
+++ /trunk/PROGRAMMES/initfig.m	(revision 30)
@@ -42,4 +42,9 @@
 % $Id$
 %
+% - fplod 2009-08-25T16:41:34Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * ajout pour octave 3.2 (uniquement sous 
+%     /usr/local_linux/octave/64/3.2.2/bin/octave donc dedale)
+%
 % - fplod 2009-08-17T12:39:51Z aedon.locean-ipsl.upmc.fr (Darwin)
 %
@@ -56,5 +61,16 @@
  set(gcf,'PaperPosition',[0,0,1,1]);
 else
- warning('www : octave can not initialise Units, position, PaperType PaperUnits PaperPosition');
+ switch OCTAVE_VERSION
+  case {'3.0.2'}
+   warning('www : octave 3.0.2 can not initialise Units, position, PaperType PaperUnits PaperPosition');
+  case {'3.2.2'}
+   set(gcf,'Units','normalized');
+   set(gcf,'position',[0.1 0.08,0.5,0.8])
+   set(gcf,'PaperType', 'a4');
+   set(gcf,'PaperUnits','normalized');
+   set(gcf,'PaperPosition',[0,0,1,1]);
+  otherwise
+   warning('www : do no know about octave initialise Units, position, PaperType PaperUnits PaperPosition');
+  end
 end
 
Index: /trunk/PROGRAMMES/reg_seriestemporelles.m
===================================================================
--- /trunk/PROGRAMMES/reg_seriestemporelles.m	(revision 30)
+++ /trunk/PROGRAMMES/reg_seriestemporelles.m	(revision 30)
@@ -0,0 +1,227 @@
+%REG_SERIESTEMPORELLES calcule la regression lineaire d'une serie temporelle issue du champ de donnees initial
+% (selectionne en un point d'espace particulier) et une composante principale choisie
+%
+% Date: Avril-Mai 2008
+% Auteur: juliette.mignot@locean-ipsl.upmc.fr
+% Contexte Ecole d'ete ENVI-STAT 2008
+%%
+
+%+
+%
+% module
+% ======
+%
+% ``reg_seriestemporelles``
+%
+% DESCRIPTION
+% ===========
+%
+% ``reg_seriestemporelles`` lit les fichiers fichiers
+% ``../DONNEES/lsmask_tropatl_PacandMedblancs.mat``,
+% ``../DONNEES/skt.mon.tropatl.mat``.
+%
+% ``reg_seriestemporelles`` lit aussi le fichier
+% ``../DONNEES//eof_SSTmens_tropAtl30N-20S.mat``
+% produit par eof_NCEP_sais.m_.
+%
+% ``reg_seriestemporelles`` calcule ++
+%
+% ``reg_seriestemporelles`` affiche ++
+%
+% ``reg_seriestemporelles`` sauve cette image dans
+% ``./eof_SSTmens_tropAtl30N-20S.mat_reg_seriestemporelles_1.ps``.
+%
+% EXAMPLES
+% ========
+%
+% ::
+%
+% >> tpacpandreg_startup
+% >> reg_seriestemporelles
+%
+%
+% SEE ALSO
+% ========
+%
+% tpacpandreg_startup.m_
+%
+% .. _tpacpandreg_startup.m : tpacpandreg_startup.m.html
+%
+% eof_NCEP_sais.m_
+%
+% .. _eof_NCEP_sais.m : eof_NCEP_sais.m.html
+%
+% TODO
+% ====
+%
+% improve description
+%
+% pourquoi ``return`` à la fin
+%
+% ajout xlabel et ylabel
+%
+% EVOLUTIONS
+% ==========
+%
+% $Id$
+%
+% - fplod 2009-08-26T09:49:43Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * remplacement de ``eps`` par ``ps`` pour pouvoir imprimer les figures 
+%     produites par ``octave``
+%
+% - fplod 2009-08-25T10:40:16Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * modif commentaires liée au renommage
+%
+% - jmignot 2009-08-25
+%
+%   * PROGRAMMES/Reg_seriestemporelles.m renommé en
+%     PROGRAMMES/reg_seriestemporelles.m pour éviter confusion
+%     avec RESULTATS/Reg_seriestemporelles.m
+%   * chgt de noms de fichiers : ``EOF_`` devient ``eof_``
+%
+% - fplod 2009-08-24T08:46:31Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * add header
+%   * répertoire de données relatifs au répertoire courant
+%   * suppression de ``clear all`` présent dans ``tpacpandreg_startup.m``
+%   * as octave save default format is ASCII and no default extension
+%     while matlab save default format is MAT v5 mat-file (little endian) and
+%     default extension is ``.mat``, write more precise save instruction
+%   * save the figure
+%
+%-
+
+%% PARAMETRES DU PROGRAMME
+% selection de la composante principale a considerer
+ipc=1; %rang de la composante principale que l'on va considerer ici
+
+% les coordonnees du point d'espace auquel on veut comparer cette
+% composante principale doit etre choisi plus bas, apres chargement des
+% longitudes et latitudes
+%%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%              CHARGEMENT ET PREPARATION DES DONNEES
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+nannees=60;
+
+repertoire=[ '..' filesep 'DONNEES' filesep ]; % a adapter a chaque cas: chemin d acces aux donnees
+
+% Mask
+%-----
+   file='lsmask_tropatl_PacandMedblancs.mat'
+   load([repertoire file]);
+% Data
+%--------
+  repertoire2=[repertoire 'MAT' filesep] % a adapter a chaque cas: chemin d acces aux donnees
+  file=['skt.mon.tropatl.mat']
+  load([repertoire2 file]); % chargement des donnees
+  tab=permute(tab,[2 3 1]); % permutation des dimensions pour organiser le tableau sous la forme [lat lon temps]
+
+% Application du masque
+%----------------------
+mask3D=repmat(mask,[1 1 size(tab,3)]);
+tab(mask3D~=0)=NaN;
+
+
+%  MANIPULATIONS TEMPORELLES: MOYENNES
+
+% Construction de moyennes mensuelles saisonnieres
+%---------------------------------------------------
+tab=reshape(tab,size(tab,1),size(tab,2),12,nannees);
+tab=squeeze(mean(tab,4));
+
+%Calcul des anomalies
+%tab_mean=squeeze(mean(tab,3));
+%tab_mean=repmat(tab_mean,[1 1 size(tab,3)]);
+%anom=tab-tab_mean;
+
+%%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%% selection d'un point d'espace
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% le point d'espace est selectionne en reperant son indice i,j dans le
+% tableau de longitude et de latitude
+
+%i=find(lon>-96 & lon<-94);
+%j=find(lat>19 & lat<21);
+i=find(lon>-96 & lon<-94);
+j=find(lat>19 & lat<21);
+%i=find(lon>-21 & lon<-19);
+%j=find(lat>19 & lat<21);
+%i=find(lon>4 & lon<6)
+%j=find(lat>-11 & lat<-9)
+
+SST=squeeze(tab(j,i,:));
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%% chargement des composantes principales
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+file=['eof_SSTmens_tropAtl30N-20S.mat']; % nom du fichier sous lequel on avait enregistre le resultat de l'ACP en moyenne annuelle
+titledataeof='SST' % je definis ce nom pour pouvoir l'inclure dans le graphique plus bas.
+if (run_octave == 0)
+ load([repertoire file]);
+else
+ load('-v7',[repertoire file]);
+end
+PC=PC(:,ipc);
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%               NUAGE DE POINTS
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+ifigure=1;
+figure(ifigure);
+plot(SST,PC,'x')
+hold on
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%%%   REGRESSION LINEAIRE
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+%[b0,b1,s,R,sigb0,sigb1] = linreg(SST,serie);
+
+SSTmean=mean(SST);PCmean=mean(PC);
+n=length(SST);
+
+
+% la fin du programme est a completer (les formules sont dans le cours ou
+% dans les livres)
+
+% calcul de b1 et b0, coefficient directeur et ordonnee a l origine de la
+% droite de regression
+
+
+%trace de la droite de regression
+%---------------------------------
+
+
+%Calcul du coefficient de détermination
+%----------------------------------------
+
+%significativite statistique comparaison de la correlation r et de la
+%valeur critique rc donnee par la table de student en fonction du nombre de
+%degres de liberte ndf
+%------------------------------
+
+
+% significativite statistique de b0 et b1
+%--------------------------------------------
+
+% sauvegarde de la figure
+printer='ps';
+print_printer=['-d', printer];
+fullfilename=['.' filesep file '_' mfilename '_'  num2str(ifigure) '.' printer];
+print(print_printer,fullfilename);
+clear print_printer;
+clear fullfilename;
+clear printer;
+%
+clear ifigure;
+
+return
Index: /trunk/PROGRAMMES/tpacpandreg_startup.m
===================================================================
--- /trunk/PROGRAMMES/tpacpandreg_startup.m	(revision 20)
+++ /trunk/PROGRAMMES/tpacpandreg_startup.m	(revision 30)
@@ -1,3 +1,3 @@
-% TPACPANDREG_STARTUP initialisation of matlab/octave environment for TPACPANDREG project
+%TPACPANDREG_STARTUP initialisation of matlab/octave environment for TPACPANDREG project
 
 %+
@@ -33,12 +33,12 @@
 % ========
 %
-% tpacpandreg_profile.sh_
-%
-% .. _tpacpandreg_profile.sh : tpacpandreg_profile.sh.html
-%
 % EVOLUTIONS
 % ==========
 %
 % $Id$
+%
+% - fplod 2009-08-25T10:36:26Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * suppression de ``more on`` (parce que pas pratique pour les Ã©tudiants)
 %
 % - fplod 2009-08-21T10:37:54Z aedon.locean-ipsl.upmc.fr (Darwin)
@@ -52,6 +52,4 @@
 %  - path for matlab/octave commands used in TPACPANDREG project
 %  - IO directories in global variables
-%
-% not yet done because of specific meaning of private and toolbox in matlab
 %
 % check if set_param(0,'CharacterEncoding','ISO-8859-1') is a better solution
@@ -74,3 +72,2 @@
 end
 %
-more on
Index: /trunk/PROGRAMMES/whorldmap.m
===================================================================
--- /trunk/PROGRAMMES/whorldmap.m	(revision 20)
+++ /trunk/PROGRAMMES/whorldmap.m	(revision 30)
@@ -40,4 +40,8 @@
 % 
 % $Id$
+%
+% - fplod 2009-08-26T09:02:22Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * add missing octave warning on Erasemode
 %
 % - fplod 2009-08-25T09:32:28Z aedon.locean-ipsl.upmc.fr (Darwin)
@@ -139,5 +143,5 @@
     set(lh,'Erasemode','none');
    else
-    warning('www : octave do not know Erasemode property')
+    warning('www : octave do not know Erasemode property');
    end
    set(gca,'aspect',[1 1]);
@@ -145,6 +149,11 @@
 
 if (drawtwice),
-   line(coastlines(1,:),coastlines(2,:),'Linestyle','-',...
-   'Erasemode','none');
+   if ~ exist('OCTAVE_VERSION','builtin')
+    line(coastlines(1,:),coastlines(2,:),'Linestyle','-',...
+    'Erasemode','none');
+   else
+    warning('www : octave do not know Erasemode property');
+    line(coastlines(1,:),coastlines(2,:),'Linestyle','-')
+   end
 end;
 
Index: unk/RESULTATS/Reg_seriestemporelles.m
===================================================================
--- /trunk/RESULTATS/Reg_seriestemporelles.m	(revision 20)
+++ 	(revision )
@@ -1,234 +1,0 @@
-%REG_SERIESTEMPORELLES calcule la regression lineaire d'une serie temporelle issue du champ de donnees initial
-% (selectionne en un point d'espace particulier) et une composante principale choisie      
-%
-% Date: Avril-Mai 2008
-% Auteur: juliette.mignot@locean-ipsl.upmc.fr
-% Contexte Ecole d'ete ENVI-STAT 2008
-%+
-%
-% module
-% ======
-%
-% ``Reg_seriestemporelles``
-%
-% DESCRIPTION
-% ===========
-%
-% ``Reg_seriestemporelles`` lit les fichiers fichiers 
-% ``../DONNEES/lsmask_tropatl_PacandMedblancs.mat``,
-% ``../DONNEES/skt.mon.tropatl.mat``.
-%
-% ``Reg_seriestemporelles`` lit aussi le fichier 
-% ``../DONNEES//EOF_SSTan_tropAtl30N-20S.mat``
-% produit par eof_NCEP_an.m_.
-%
-% ``Reg_seriestemporelles`` calcule ++
-%
-% ``Reg_seriestemporelles`` affiche ++
-%
-% ``Reg_seriestemporelles`` sauve cette image dans 
-% ``./EOF_SSTan_tropAtl30N-20S.mat_Reg_seriestemporelles_1.eps``.
-%
-% EXAMPLES
-% ========
-%
-% ::
-%
-% >> addpath('../PROGRAMMES/')
-% >> tpacpandreg_startup
-% >> Reg_seriestemporelles
-%
-%
-% SEE ALSO
-% ========
-%
-% tpacpandreg_startup.m_
-%
-% .. _tpacpandreg_startup.m : tpacpandreg_startup.m.html
-%
-% eof_NCEP_an.m_
-%
-% .. _eof_NCEP_an.m : eof_NCEP_an.m.html
-%
-% TODO
-% ====
-%
-% improve description
-%
-% pourquoi ``return`` à la fin
-%
-% ajout xlabel et ylabel
-%
-% EVOLUTIONS
-% ==========
-%
-% $Id$
-%
-% - fplod 2009-08-24T11:44:38Z aedon.locean-ipsl.upmc.fr (Darwin)
-%
-%   * add header
-%   * répertoire de données relatifs au répertoire courant
-%   * suppression de ``clear all`` présent dans ``tpacpandreg_startup.m``
-%   * as octave save default format is ASCII and no default extension
-%     while matlab save default format is MAT v5 mat-file (little endian) and
-%     default extension is ``.mat``, write more precise save instruction
-%   * save the figure
-%   * replace file radical ``EOF_SST_tropAtl`` by ``EOF_SSTan_tropAtl``
-%   * split lines with multiple statements
-%   * handle ``corrcoef`` octave/matlab difference of output
-%
-%-
-
-%%
-%% PARAMETRES DU PROGRAMME
-% selection de la composante principale a considerer
-ipc=1; %rang de la composante principale que l'on va considerer ici
-
-% les coordonnees du point d'espace auquel on veut comparer cette
-% composante principale doit etre choisi plus bas, apres chargement des
-% longitudes et latitudes
-%%
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%%%%              CHARGEMENT ET PREPARATION DES DONNEES
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-nannees=60;
-repertoire=[ '..' filesep 'DONNEES' filesep ]; % a adapter a chaque cas: chemin d acces aux donnees
-
-% Mask
-%------
-file='lsmask_tropatl_PacandMedblancs'
-   load([repertoire file]);
-% Data
-%--------
-  repertoire2=[repertoire 'MAT' filesep] % a adapter a chaque cas: chemin d acces aux donnees
-  file=['skt.mon.tropatl.mat']
-  load([repertoire2 file]); % chargement des donnees
-  tab=permute(tab,[2 3 1]); % permutation des dimensions pour organiser le tableau sous la forme [lat lon temps]
-
-% Application du masque
-%----------------------
-mask3D=repmat(mask,[1 1 size(tab,3)]);
-tab(mask3D~=0)=NaN;
-
-
-%  MANIPULATIONS TEMPORELLES: MOYENNES 
-
-% Construction de moyennes annuelles
-%---------------------------------------
-tab=reshape(tab,size(tab,1),size(tab,2),12,nannees);
-tab_an=squeeze(mean(tab,3));
-%Calcul des anomalies
-%tab_an_mean=squeeze(mean(tab_an,3));
-%tab_an_mean=repmat(tab_an_mean,[1 1 size(tab_an,3)]);
-%tab_an=tab_an-tab_an_mean;
-
-%%
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%%%% selection d'un point d'espace
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-% le point d'espace est selectionne en reperant son indice i,j dans le
-% tableau de longitude et de latitude
-
-%i=find(lon>-96 & lon<-94);
-%j=find(lat>19 & lat<21);
-i=find(lon>-96 & lon<-94);
-j=find(lat>19 & lat<21);
-%i=find(lon>-21 & lon<-19);
-%j=find(lat>19 & lat<21);
-%i=find(lon>4 & lon<6)
-%j=find(lat>-11 & lat<-9)
-
-SST=squeeze(tab_an(j,i,:));
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%%%%% chargement des composantes principales
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-file=['EOF_SSTan_tropAtl30N-20S.mat'];
-titledataeof='SST' % je definis ce nom pour pouvoir l'inclure dans le graphique plus bas. 
-if (run_octave == 0)
- load([repertoire file]);
-else
- load('-v7',[repertoire file]);
-end
-PC=PC(:,ipc);
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%%%%               NUAGE DE POINTS
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-ifigure=1;
-figure(ifigure); 
-plot(SST,PC,'x')
-hold on
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%%%%%%   REGRESSION LINEAIRE
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-%[b0,b1,s,R,sigb0,sigb1] = linreg(SST,serie);
-
-SSTmean=mean(SST);PCmean=mean(PC);
-n=length(SST);
-
-% calcul de b1 et b0, coefficient directeur et ordonnee a l origine de la
-% droite de regression
-b1=sum((PC-PCmean).*(SST-SSTmean))/(sum(SST.*SST)-n*SSTmean^2);
-b0=PCmean-b1*SSTmean;
-
-PCc=b0+b1*SST;
-
-%trace de la droite de regression
-%---------------------------------
-plot(SST,PCc,'r')
-
-%Calcul du coefficient de détermination
-%----------------------------------------
-s=sqrt(sum((PC-PCc).^2)/(n-2));
-R=sum((PCc-PCmean).^2)/sum((PC-PCmean).^2);
-
-%significativite statistique comparaison de la correlation r et de la
-%valeur critique rc donnee par la table de student en fonction du nombre de
-%degres de liberte ndf
-%------------------------------
-
-r=corrcoef(SST,PC);
-if (run_octave == 0)
- r=r(2,1);
-end
-
-% remarque: b1=r*std(PC)/std(SST)
-
-% on estime le nombre de degrés de liberté
-ndf=nannees; %(l'atmosphère est approximativement un bruit blanc à l'échelle de temps interanuellee)
-
-% on calcule le seuil de significativité au percentile P en inversant
-% l'écart réduit T=y*sqrt(nn)./(1-y.^2).^0.5 qui est distribue selon
-% un loi de student avec nn=ndf-2 degres de liberte.
-P=.95; 
-nn=ndf-2;
-t=tinv(P,nn);
-rc = t ./ ( nn + t.^2 ).^0.5;
-% rc est la valeur ccritique de la correlation pour rejeter l'hypothese nulle
-
-title(['r = ' num2str(r) '  valeur critique: rc= ' num2str(rc)])
-
-% significativite statistique de b0 et b1
-%--------------------------------------------
-
-% sauvegarde de la figure
-printer='eps';
-print_printer=['-d', printer];
-fullfilename=['.' filesep file '_' mfilename '_'  num2str(ifigure) '.' printer];
-print(print_printer,fullfilename);
-clear print_printer;
-clear fullfilename;
-clear printer;
-%
-clear ifigure;
-
-return
-
Index: /trunk/RESULTATS/carteetstd_SST.m
===================================================================
--- /trunk/RESULTATS/carteetstd_SST.m	(revision 20)
+++ /trunk/RESULTATS/carteetstd_SST.m	(revision 30)
@@ -1,5 +1,5 @@
 %CARTEETSTD_SST permet de visualiser les cartes mensuelles moyennes
 % de temperature de surface       
-% ainsi que les cartes annuelles moyennes et les ecart type
+% ainsi que les cartes annuelles moyennes et les écart type
 % correspondants
 
@@ -19,5 +19,5 @@
 % ``carteetstd_SST`` permet de visualiser les cartes mensuelles moyennes
 % de temperature de surface
-% ainsi que les cartes annuelles moyennes et les ecart type
+% ainsi que les cartes annuelles moyennes et les écart type
 % correspondants.
 %
@@ -26,9 +26,15 @@
 % et ``../DONNEES/skt.mon.tropatl.mat``.
 %
-% ``carteetstd_SST`` calcule ++blabla++
-%
-% ``carteetstd_SST`` affiche à l'écran les 12 mois de ++
-%
-% ``carteetstd_SST`` sauve ces images dans ``./carteetstd_SST_[123].eps``.
+% ``carteetstd_SST`` calcule les moyennes mensuelles d'une part et
+% les moyennes interannuelles d'autre part de la SST dans
+% l'Atlantique tropical à partir d'un champ mensuel de 60 ans
+% fourni en entrée.
+%
+% ``carteetstd_SST`` affiche à l'écran les 12 valeurs de moyennes
+% mensuelles (figure 1), une sélection de 12 valeurs annuelles
+% moyennes (figure 2) et les écart-type mensuels et interannuels
+% (figure 3).
+%
+% ``carteetstd_SST`` sauve ces images dans ``./carteetstd_SST_[123].ps``.
 %
 % EXAMPLES
@@ -49,5 +55,25 @@
 %
 % initfig.m_
+%
 % .. _initfig.m : initfig.m.html
+%
+% whorldmap.m_
+%
+% .. _whorldmap.m : whorldmap.m.html
+%
+% colorbartype.m_
+%
+% .. _colorbartype.m : colorbartype.m.html
+%
+% contlab.m_
+%
+% .. _contlab.m : contlab.m.html
+%
+%
+% KNOWN PROBLEMS
+% ==============
+%
+% pb avec les diacritiques à l'impression de .ps dans le titres des figures 
+% sous octave : é se transforme en o barré.
 %
 % TODO
@@ -61,4 +87,28 @@
 % $Id$
 %
+% - fplod 2009-08-26T09:56:51Z aedon.locean-ipsl.upmc.fr (Darwin)
+%   
+%   * nom des mois à la place du numéro dans les titres
+%   * remplacement de ``eps`` par ``ps`` pour pouvoir imprimer les figures 
+%     produites par ``octave``
+%
+%
+% - jmignot 2009-08-25
+%
+%   * suppression d'appel à contlab et colorbartype sous octave
+%
+% - fplod 2009-08-25T10:40:16Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * modif commentaires liée au dernier chgt de Juliette
+%
+% - jmignot 2009-08-25
+%
+%   * remplacement de ``colorbar`` par ``colorbartype``
+%     parce pas beau ni en matlab ni en octave (dégradé, noir et blanc dans 
+%     les .jpg) 
+%   * contraintes sur les limites x et y de axes car matlab et octave ne
+%     fixent pas les même limites par défaut
+%   * utisation de ``contlab`` pour ++
+% 
 % - fplod 2009-08-17T12:39:51Z aedon.locean-ipsl.upmc.fr (Darwin)
 %
@@ -92,5 +142,5 @@
 %Mask
 %------
-file='lsmask_tropatl_PacandMedblancs'
+file='lsmask_tropatl_PacandMedblancs.mat'
    load([repertoire file]);
 % Data
@@ -110,5 +160,5 @@
 %  MANIPULATIONS TEMPORELLES: MOYENNES 
 %---------------------------------------------------
-tab=reshape(tab,size(tab,1),size(tab,2),12,nannees); %organisation du tableau sous la forme [lat lon mois annees]
+tab=reshape(tab,size(tab,1),size(tab,2),12,nannees); %organisation du tableau sous la forme [lat lon mois années]
 
 % Construction de moyennes annuelles
@@ -120,4 +170,11 @@
 tab_mens=squeeze(mean(tab,4));
 
+%calcul des écart-types
+%------------------------
+%écart-type mensuel moyen: 
+std_mens=squeeze(std(tab_mens,0,3));
+
+%écart-type interannuel moyen: 
+std_an=squeeze(std(tab_an,0,3));
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
@@ -125,7 +182,18 @@
 %%%%               TRACE
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-
+%parametres geographiques du graphe
+  geo=[-100 20 -30 30]; %[lonmin lonmax latmin latmax];
+  yti=[-30:10:30]; xti=[-100:30:20]; %labels de latitude et de longitude
+
+		 
 % Variabilité saisonnière
+%-------------------------
+%niveaux de couleurs
+lev=[0:1:30];%contours a tracer, en unite du champ a dessiner, ici
+             %la SST en degC
+clip_lev=[20 30];%definit les limites des couleurs: bleu fonce au
+                 %rouge fonce
+
+
 ifigure=1;
 figure(ifigure); 
@@ -141,15 +209,43 @@
   end
   z=squeeze(tab_mens(:,:,iplot));
-  contourf(x,y,z,[0:1:30])
-  caxis([20 30])
+  contourf(x,y,z,lev)
+  caxis(clip_lev)
   hold on
-  ylim([-30 30])
-  xlim([-100 20])
-  title(['mois ' num2str(iplot)])
-end; 
-colorbar('East')
+  whorldmap
+
+  % limites geographies et labels de longitude et de latitude
+  ylim([geo(3) geo(4)])
+  xlim([geo(1) geo(2)])
+  set(gca, 'Xtick', xti, 'Ytick', yti)
+  if (run_octave == 0)
+   [longstr, latstr] = contlab(xti, yti);
+   set(gca, 'Xticklabel', longstr, 'Yticklabel', latstr)
+  end;
+  if iplot>=10
+    xlabel('longitude')
+  end;
+  if iplot==1|iplot==4|iplot==7|iplot==10
+    ylabel('latitude')
+  end;
+  
+  %titre: nom de chaque mois
+  cmois = datestr(datenum([2000 iplot 1 00 00 00]), 'mmm');
+  title(cmois);
+  clear cmois;
+  
+end;
+
+%barre de couleur
+if (run_octave == 0)
+  pos=[.2 .03 .6 .013];
+  hpal=jet(100); %carte des couleurs: du bleu au rouge par defaut.
+  colorbartype(pos,lev,1,clip_lev,hpal,0);
+  title('SST mensuelle moyenne en degC')
+else
+  colorbar('East') %colorbar par defaut
+end;
 
 % sauvegarde de la figure
-printer='eps';
+printer='ps';
 print_printer=['-d', printer];
 fullfilename=['./',mfilename,'_',num2str(ifigure),'.',printer];
@@ -159,5 +255,13 @@
 clear fullfilename;
 
+
 % Variabilité interannuelle
+%-------------------------
+%niveaux de couleurs
+lev=[0:1:30];%contours a tracer, en unite du champ a dessiner, ici
+             %la SST en degC
+clip_lev=[20 30];%definit les limites des couleurs: bleu fonce au
+                 %rouge fonce
+
 % Cartes individuelles
 
@@ -175,15 +279,40 @@
   end
   z=squeeze(tab_an(:,:,iplot));
-  contourf(x,y,z,[0:1:30])
-  caxis([20 30])
+  contourf(x,y,z,lev)
+  caxis(clip_lev)
   hold on
-  ylim([-30 30])
-  xlim([-100 20])
+  whorldmap
+
+  % limites geographies et labels de longitude et de latitude
+  ylim([geo(3) geo(4)])
+  xlim([geo(1) geo(2)])
+  set(gca, 'Xtick', xti, 'Ytick', yti)
+  if (run_octave == 0)
+   [longstr, latstr] = contlab(xti, yti);
+   set(gca, 'Xticklabel', longstr, 'Yticklabel', latstr)
+  end;
+  if iplot>=10
+    xlabel('longitude')
+  end;
+  if iplot==1|iplot==4|iplot==7|iplot==10
+    ylabel('latitude')
+  end;
+
+  %titre: numéro de l'année
   title(['annee ' num2str(iplot+1948-1)])
 end; 
-colorbar('East')
+
+%barre de couleur
+if (run_octave == 0)
+  pos=[.2 .03 .6 .013];
+  hpal=jet(100); %carte des couleurs: du bleu au rouge par defaut.
+  colorbartype(pos,lev,1,clip_lev,hpal,0);
+  title('SST mensuelle moyenne en degC')
+else
+  colorbar('East') %colorbar par defaut
+end;
 
 % sauvegarde de la figure
-printer='eps';
+printer='ps';
 print_printer=['-d', printer];
 fullfilename=['./',mfilename,'_',num2str(ifigure),'.',printer];
@@ -193,9 +322,21 @@
 clear fullfilename;
 
-% ecart type
+
+% écart type
+%-----------
+
+
+%trace
 ifigure=ifigure+1;
 figure(ifigure);
 clf();
 initfig; 
+
+%écart-type mensuel moyen
+
+%niveaux de couleurs
+lev=[0:.4:5];
+clip_lev=[0 3.6];
+
 subplot(2,2,1)
 if (run_octave == 0)
@@ -205,25 +346,72 @@
  [x,y]=meshgrid(lon,lat);
 end
-z=squeeze(std(tab_mens,0,3));
-contourf(x,y,z,[0:.4:5])
-caxis([0 3.5])
+z=std_mens;
+contourf(x,y,z,lev)
+caxis([clip_lev])
 hold on
-ylim([-30 30])
-xlim([-100 20])
+whorldmap
+  % limites geographies et labels de longitude et de latitude
+  ylim([geo(3) geo(4)])
+  xlim([geo(1) geo(2)])
+  set(gca, 'Xtick', xti, 'Ytick', yti)
+  if (run_octave == 0)
+   [longstr, latstr] = contlab(xti, yti);
+   set(gca, 'Xticklabel', longstr, 'Yticklabel', latstr)
+  end;
+  xlabel('longitude')
+  ylabel('latitude')
+
 title('ecart-type mensuel moyen')
-colorbar('SouthOutside')
+%barre de couleur
+if (run_octave == 0)
+  pos=[.1 .5 .35 .013];
+  hpal=jet(100); %carte des couleurs: du bleu au rouge par defaut.
+  colorbartype(pos,lev,1,clip_lev,hpal,0);
+  title('degC')
+else
+  colorbar('East') %colorbar par defaut
+end;
+
+
+
+%écart-type interannuel  moyen
+
+%niveaux de couleurs
+lev=[0:.05:.7];%la variabilite interannuelle est moins forte que la
+               %variabilite saisonniere moyenne. 
+clip_lev=[0.2 .55];
 
 subplot(2,2,2)
-z=squeeze(std(tab_an,0,3));
-contourf(x,y,z,[0:.05:.7])
-caxis([0.2 .55])
+z=std_an;
+contourf(x,y,z,lev)
+caxis(clip_lev)
 hold on
-ylim([-30 30])
-xlim([-100 20])
-title('ecart-type interannuel')
-colorbar('SouthOutside')
+whorldmap
+  % limites geographies et labels de longitude et de latitude
+  ylim([geo(3) geo(4)])
+  xlim([geo(1) geo(2)])
+  set(gca, 'Xtick', xti, 'Ytick', yti)
+  if (run_octave == 0)
+   [longstr, latstr] = contlab(xti, yti);
+   set(gca, 'Xticklabel', longstr, 'Yticklabel', latstr)
+  end;
+  xlabel('longitude')
+  ylabel('latitude')
+
+title('écart-type interannuel')
+%barre de couleur
+if (run_octave == 0)
+  pos=[.55 .5 .35 .013];
+  hpal=jet(100); %carte des couleurs: du bleu au rouge par defaut.
+  colorbartype(pos,lev,1,clip_lev,hpal,0);
+  title('degC')
+else
+  colorbar('East') %colorbar par defaut
+end;
+
+
 
 % sauvegarde de la figure
-printer='eps';
+printer='ps';
 print_printer=['-d', printer];
 fullfilename=['./',mfilename,'_',num2str(ifigure),'.',printer];
Index: /trunk/RESULTATS/eof_NCEP_an.m
===================================================================
--- /trunk/RESULTATS/eof_NCEP_an.m	(revision 20)
+++ /trunk/RESULTATS/eof_NCEP_an.m	(revision 30)
@@ -8,5 +8,5 @@
 %
 % le programme permet egalement de visualiser les resultats. Voir le doc
-% .pdf joint pour l'interpretation de cette visualisation. 
+% .pdf joint pour l'interpretation de cette visualisation.
 %
 % Date: Avril-Mai 2008
@@ -18,10 +18,10 @@
 % ======
 %
-% ``eof_NCEP_an`` 
+% ``eof_NCEP_an``
 %
 % DESCRIPTION
 % ===========
 %
-% ``eof_NCEP_an`` lit les fichiers 
+% ``eof_NCEP_an`` lit les fichiers
 % ``../DONNEES/lsmask_tropatl_PacandMedblancs.mat``
 % et ``../DONNEES/skt.mon.tropatl.mat``.
@@ -29,15 +29,17 @@
 % ``eof_NCEP_an`` calcule ++blabla++
 %
-% ``eof_NCEP_an`` sauvegarde dans ``../DONNEES/EOF_SSTan_tropAtl++.mat`` 
-% ou ``++/DONNEES/EOF_SSTan++.mat`` ++
-%
-% Ces fichiers sont lus par ``Reg_serietemporelles.m_``.
+% ``eof_NCEP_an`` sauvegarde dans ``../DONNEES/eof_SSTan_tropAlt30N-20S.mat``.
+%
+% Ces fichiers sont lus par ``regandsignif_seriestemporelles.m_``.
 %
 % ``eof_NCEP_an`` affiche trois figures :
+%
 %  - ++
+%
 %  - ++
+%
 %  - valeurs propres de la matrice de covariance
 %
-% ``eof_NCEP_an`` sauve ces images dans ``./EOF_SSTan_tropAtl++_[123].eps``.
+% ``eof_NCEP_an`` sauve ces images dans ``./eof_SSTan_tropAlt30N-20S_[123].ps``.
 %
 % ++
@@ -60,9 +62,10 @@
 %
 % initfig.m_
+%
 % .. _initfig.m : initfig.m.html
 %
-% Reg_seriestemporelles.m_
-%
-% .. _Reg_seriestemporelles.m : Reg_seriestemporelles.m.html
+% regandsignif_seriestemporelles.m_
+%
+% .. _regandsignif_seriestemporelles.m : regandsignif_seriestemporelles.m.html
 %
 % whorldmap.m_
@@ -77,5 +80,5 @@
 % tester la branche NetCDF
 %
-% sous octave, triangle bizarre dans la dernière barre du graphique valeurs 
+% sous octave, triangle bizarre dans la dernière barre du graphique valeurs
 % propres.
 %
@@ -85,8 +88,28 @@
 % $Id$
 %
+% - fplod 2009-08-26T10:26:35Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * remplacement de ``eps`` par ``ps`` pour pouvoir imprimer les figures
+%     produites par ``octave``
+%
+% - jmignot 2009-08-25
+%
+%   * suppression d'appel à contlab et colorbartype sous octave
+%
+% - fplod 2009-08-25T10:40:16Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * modif commentaires liée au renommage
+%
+% - jmignot 2009-08-25
+%
+%   * RESULTATS/Reg_seriestemporelles.m renommé en
+%     RESULTATS/regandsignif_seriestemporelles.m pour éviter confusion
+%     avec PROGRAMMES/Reg_seriestemporelles.m
+%   * chgt de noms de fichiers : ``EOF_`` devient ``eof_``
+
 % - fplod 2009-08-21T15:04:16Z aedon.locean-ipsl.upmc.fr (Darwin)
 %
-%   * as octave save default format is ASCII and no default extension 
-%     while matlab save default format is MAT v5 mat-file (little endian) and 
+%   * as octave save default format is ASCII and no default extension
+%     while matlab save default format is MAT v5 mat-file (little endian) and
 %     default extension is ``.mat``, write more precise save instruction
 %   * add header
@@ -109,6 +132,6 @@
 N=4; %nb d'EOF a tracer
 
-nomfic='EOF_SSTan_tropAtl'
-% mettre 'EOF_SSTan_tropAtl' quand on travaille en moyennes annuelles. Ce
+nomfic='eof_SSTan_tropAtl'
+% mettre 'eof_SSTan_tropAtl' quand on travaille en moyennes annuelles. Ce
 % fichier sera rappele dans les programme de regression lineaire
 %%
@@ -119,12 +142,12 @@
 
 nannees=60;
-   
+
 repertoire=[ '..' filesep 'DONNEES' filesep ] % a adapter a chaque cas: chemin d acces aux donnees
 
 % Mask
 %-----
-file='lsmask_tropatl_PacandMedblancs'
+file='lsmask_tropatl_PacandMedblancs.mat'
    load([repertoire file]);
-  ny=size(mask,1); nx=size(mask,2); % dimensions spatiales des tableaux 
+  ny=size(mask,1); nx=size(mask,2); % dimensions spatiales des tableaux
 
 % Data
@@ -135,5 +158,5 @@
   load([repertoire2 file]); % chargement des donnees
 
-   
+
   %SI NETCDF
   %repertoire2=[repertoire 'NETCDF' filesep]% a adapter a chaque cas: chemin d acces aux donnees
@@ -141,5 +164,5 @@
   % f=netcdf([repertoire file],'nowrite');
   % tab=f{'skt'}(:,:,:);
-  
+
 
    tab=permute(tab,[2 3 1]);% permutation des dimensions pour organiser le tableau sous la forme [lat lon temps]
@@ -168,5 +191,5 @@
 %----------------
 %Si on veut éventuellement attribuer un poids différent aux différents
-%points spatiaux, on le fait ici. 
+%points spatiaux, on le fait ici.
 
 % Préparation du zoom géographique: sur le masque, j'identifie les points
@@ -176,5 +199,5 @@
 limite_latS=-20; %limite du domaine d'etuade en latitude
 limite_latN=30; %limite du domaine d'etuade en latitude
-lat2D=repmat(lat,[1 size(mask,2)]); 
+lat2D=repmat(lat,[1 size(mask,2)]);
 mask(lat2D>limite_latN&mask==0)=4;
 mask(lat2D<limite_latS&mask==0)=4;
@@ -208,24 +231,24 @@
 % matrice de covariance:
   C = Z * Z';
-%  eval(['save ' repertoire 'MatCov_SST_tropAtl' num2str(limite_latN) 'N' num2str(limite_latS) 'S.asc C -ascii']) 
-  
-% diagonalisation  
-  [E,Lambda] = eig(C); 
-   % Lambda est la matrice diagonale des valeurs propres et E la matrice dont les colonnes sont les vecteurs propres correspondants.  
+%  eval(['save ' repertoire 'MatCov_SST_tropAtl' num2str(limite_latN) 'N' num2str(limite_latS) 'S.asc C -ascii'])
+
+% diagonalisation
+  [E,Lambda] = eig(C);
+   % Lambda est la matrice diagonale des valeurs propres et E la matrice dont les colonnes sont les vecteurs propres correspondants.
    % Soit Z*E = E*Lambda
-   % Les vecteurs propres sont appelés les EOF. Chaque EOF peut être considéré comme une carte de l'espace physique.    
-
-   
- %% 
+   % Les vecteurs propres sont appelés les EOF. Chaque EOF peut être considéré comme une carte de l'espace physique.
+
+
+ %%
  %%%% calcul de la matrice donnant les pourcentages de variance
 Lambda = diag(Lambda);
 pctg_var = round(Lambda / sum(Lambda) * 1000) / 10;
 
- 
+
  %On definit les composantes principales de Q comme la decomposition du signal decrit par Q sur les EOF E.
 
    P = Z' * E;
-    
-%%   
+
+%%
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
@@ -271,11 +294,19 @@
 
  %%               TRACE
+close all
+%parametres geographiques des graphes
+  geo=[-100 20 -30 30]; %[lonmin lonmax latmin latmax];
+  yti=[-30:10:30]; xti=[-100:30:20]; %labels de latitude et de longitude
+%niveaux de couleurs
+  lev=[-2:.02:2];
+  clip_lev=[-0.06 0.06];
+
 figure(1);
 initfig;
-figure(2); 
+figure(2);
 initfig;
 
 for i=1:4
-    
+
    %subplot(4,2,2*i-1)
    ifigure=1;
@@ -289,17 +320,35 @@
    end
    z=squeeze(TABEOF(:,:,i));
-   contourf(x,y,z);
-   caxis([-0.06 0.06])
+   contourf(x,y,z,lev);
+   caxis(clip_lev);
    hold on
-   xlim([-100 30]); xlabel('longitude')
-   ylim([-40 40]); ylabel('latitude')
-   whorldmap([-100, 30, -40, 40])
-   colorbar
-   title([num2str(p_var(i)) '%  -  (unites: K)'])
-   set(gca,'xtick',[-100;-70;-40;-10;20],'xticklabel',['100W';' 70W';' 40W';' 10W';' 20E'])
-   set(gca,'ytick',[-30;-15;0;15;30],'yticklabel',['30S';'15S';' 0 ';'15N';'30N'])
+   whorldmap
+   % limites geographies et labels de longitude et de latitude
+  ylim([geo(3) geo(4)])
+  xlim([geo(1) geo(2)])
+  set(gca, 'Xtick', xti, 'Ytick', yti)
+   if (run_octave == 0)
+    [longstr, latstr] = contlab(xti, yti);
+    set(gca, 'Xticklabel', longstr, 'Yticklabel', latstr)
+   end;
+  xlabel('longitude')
+  ylabel('latitude')
+
+   %titre de la figure
+   title([num2str(p_var(i)) '%'])
+
+  %barre de couleur
+  if (run_octave == 0)
+    hpal=jet(100); %carte des couleurs: du bleu au rouge par defaut.
+    pos1=get(gca,'Position');
+    pos=[pos1(1) pos1(2)-.075 pos1(3) .013];
+    colorbartype(pos,lev,1,clip_lev,hpal,0);
+    title('degC')
+  else
+    colorbar('East')
+  end;
 
    % sauvegarde de la figure
-   printer='eps';
+   printer='ps';
    print_printer=['-d', printer];
    fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer];
@@ -307,5 +356,5 @@
    clear fullfilename;
    clear printer;
-   
+
    %subplot(4,2,2*i)
    ifigure=2;
@@ -320,5 +369,5 @@
 
    % sauvegarde de la figure
-   printer='eps';
+   printer='ps';
    print_printer=['-d', printer];
    fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer];
@@ -332,5 +381,5 @@
 % Graphe pourcentages de variance expliqué
 ifigure=3;
-figure(ifigure); 
+figure(ifigure);
 bar(flipud(pctg_var))
 xlim([0 10])
@@ -340,5 +389,5 @@
 
 % sauvegarde de la figure
-printer='eps';
+printer='ps';
 print_printer=['-d', printer];
 fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer];
Index: /trunk/RESULTATS/regandsignif_seriestemporelles.m
===================================================================
--- /trunk/RESULTATS/regandsignif_seriestemporelles.m	(revision 30)
+++ /trunk/RESULTATS/regandsignif_seriestemporelles.m	(revision 30)
@@ -0,0 +1,251 @@
+%REGANDSIGNIF_SERIESTEMPORELLES calcule la regression lineaire d'une serie temporelle issue du champ de donnees initial
+% (selectionne en un point d'espace particulier) et une composante principale choisie      
+%
+% Date: Avril-Mai 2008
+% Auteur: juliette.mignot@locean-ipsl.upmc.fr
+% Contexte Ecole d'ete ENVI-STAT 2008
+%+
+%
+% module
+% ======
+%
+% ``regandsignif_seriestemporelles``
+%
+% DESCRIPTION
+% ===========
+%
+% ``regandsignif_seriestemporelles`` lit les fichiers fichiers 
+% ``../DONNEES/lsmask_tropatl_PacandMedblancs.mat``,
+% ``../DONNEES/skt.mon.tropatl.mat``.
+%
+% ``regandsignif_seriestemporelles`` lit aussi le fichier 
+% ``../DONNEES//eof_SSTan_tropAtl30N-20S.mat``
+% produit par eof_NCEP_an.m_.
+%
+% ``regandsignif_seriestemporelles`` calcule ++
+%
+% ``regandsignif_seriestemporelles`` affiche ++
+%
+% ``regandsignif_seriestemporelles`` sauve cette image dans 
+% ``./eof_SSTan_tropAtl30N-20S.mat_regandsignif_seriestemporelles_1.ps``.
+%
+% EXAMPLES
+% ========
+%
+% ::
+%
+% >> addpath('../PROGRAMMES/')
+% >> tpacpandreg_startup
+% >> regandsignif_seriestemporelles
+%
+%
+% SEE ALSO
+% ========
+%
+% tpacpandreg_startup.m_
+%
+% .. _tpacpandreg_startup.m : tpacpandreg_startup.m.html
+%
+% eof_NCEP_an.m_
+%
+% .. _eof_NCEP_an.m : eof_NCEP_an.m.html
+%
+% TODO
+% ====
+%
+% improve description
+%
+% pourquoi ``return`` à la fin
+%
+% ajout xlabel et ylabel
+%
+% EVOLUTIONS
+% ==========
+%
+% $Id$
+%
+% 
+% - fplod 2009-08-26T11:00:35Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * remplacement de ``eps`` par ``ps`` pour pouvoir imprimer les figures 
+%     produites par ``octave``
+% 
+% - fplod 2009-08-25T10:40:16Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * modif commentaires liée au renommage
+%
+% - jmignot 2009-08-25
+%
+%   * RESULTATS/Reg_seriestemporelles.m renommé en 
+%     RESULTATS/regandsignif_seriestemporelles.m pour éviter confusion
+%     avec PROGRAMMES/Reg_seriestemporelles.m
+%   * chgt de noms de fichiers : ``EOF_`` devient ``eof_``
+%
+% - fplod 2009-08-24T11:44:38Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * add header
+%   * répertoire de données relatifs au répertoire courant
+%   * suppression de ``clear all`` présent dans ``tpacpandreg_startup.m``
+%   * as octave save default format is ASCII and no default extension
+%     while matlab save default format is MAT v5 mat-file (little endian) and
+%     default extension is ``.mat``, write more precise save instruction
+%   * save the figure
+%   * replace file radical ``EOF_SST_tropAtl`` by ``EOF_SSTan_tropAtl``
+%   * split lines with multiple statements
+%   * handle ``corrcoef`` octave/matlab difference of output
+%
+%-
+
+%%
+%% PARAMETRES DU PROGRAMME
+% selection de la composante principale a considerer
+ipc=1; %rang de la composante principale que l'on va considerer ici
+
+% les coordonnees du point d'espace auquel on veut comparer cette
+% composante principale doit etre choisi plus bas, apres chargement des
+% longitudes et latitudes
+%%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%              CHARGEMENT ET PREPARATION DES DONNEES
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+nannees=60;
+repertoire=[ '..' filesep 'DONNEES' filesep ]; % a adapter a chaque cas: chemin d acces aux donnees
+
+% Mask
+%------
+file='lsmask_tropatl_PacandMedblancs.mat'
+   load([repertoire file]);
+% Data
+%--------
+  repertoire2=[repertoire 'MAT' filesep] % a adapter a chaque cas: chemin d acces aux donnees
+  file=['skt.mon.tropatl.mat']
+  load([repertoire2 file]); % chargement des donnees
+  tab=permute(tab,[2 3 1]); % permutation des dimensions pour organiser le tableau sous la forme [lat lon temps]
+
+% Application du masque
+%----------------------
+mask3D=repmat(mask,[1 1 size(tab,3)]);
+tab(mask3D~=0)=NaN;
+
+
+%  MANIPULATIONS TEMPORELLES: MOYENNES 
+
+% Construction de moyennes annuelles
+%---------------------------------------
+tab=reshape(tab,size(tab,1),size(tab,2),12,nannees);
+tab_an=squeeze(mean(tab,3));
+%Calcul des anomalies
+%tab_an_mean=squeeze(mean(tab_an,3));
+%tab_an_mean=repmat(tab_an_mean,[1 1 size(tab_an,3)]);
+%tab_an=tab_an-tab_an_mean;
+
+%%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%% selection d'un point d'espace
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% le point d'espace est selectionne en reperant son indice i,j dans le
+% tableau de longitude et de latitude
+
+%i=find(lon>-96 & lon<-94);
+%j=find(lat>19 & lat<21);
+i=find(lon>-96 & lon<-94);
+j=find(lat>19 & lat<21);
+%i=find(lon>-21 & lon<-19);
+%j=find(lat>19 & lat<21);
+%i=find(lon>4 & lon<6)
+%j=find(lat>-11 & lat<-9)
+
+SST=squeeze(tab_an(j,i,:));
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%% chargement des composantes principales
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+file=['eof_SSTan_tropAtl30N-20S.mat'];
+titledataeof='SST' % je definis ce nom pour pouvoir l'inclure dans le graphique plus bas. 
+if (run_octave == 0)
+ load([repertoire file]);
+else
+ load('-v7',[repertoire file]);
+end
+PC=PC(:,ipc);
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%               NUAGE DE POINTS
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+ifigure=1;
+figure(ifigure); 
+plot(SST,PC,'x')
+hold on
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%%%   REGRESSION LINEAIRE
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+%[b0,b1,s,R,sigb0,sigb1] = linreg(SST,serie);
+
+SSTmean=mean(SST);PCmean=mean(PC);
+n=length(SST);
+
+% calcul de b1 et b0, coefficient directeur et ordonnee a l origine de la
+% droite de regression
+b1=sum((PC-PCmean).*(SST-SSTmean))/(sum(SST.*SST)-n*SSTmean^2);
+b0=PCmean-b1*SSTmean;
+
+PCc=b0+b1*SST;
+
+%trace de la droite de regression
+%---------------------------------
+plot(SST,PCc,'r')
+
+%Calcul du coefficient de détermination
+%----------------------------------------
+s=sqrt(sum((PC-PCc).^2)/(n-2));
+R=sum((PCc-PCmean).^2)/sum((PC-PCmean).^2);
+
+%significativite statistique comparaison de la correlation r et de la
+%valeur critique rc donnee par la table de student en fonction du nombre de
+%degres de liberte ndf
+%------------------------------
+
+r=corrcoef(SST,PC);
+if (run_octave == 0)
+ r=r(2,1);
+end
+
+% remarque: b1=r*std(PC)/std(SST)
+
+% on estime le nombre de degrés de liberté
+ndf=nannees; %(l'atmosphère est approximativement un bruit blanc à l'échelle de temps interanuellee)
+
+% on calcule le seuil de significativité au percentile P en inversant
+% l'écart réduit T=y*sqrt(nn)./(1-y.^2).^0.5 qui est distribue selon
+% un loi de student avec nn=ndf-2 degres de liberte.
+P=.95; 
+nn=ndf-2;
+t=tinv(P,nn);
+rc = t ./ ( nn + t.^2 ).^0.5;
+% rc est la valeur ccritique de la correlation pour rejeter l'hypothese nulle
+
+title(['r = ' num2str(r) '  valeur critique: rc= ' num2str(rc)])
+
+% significativite statistique de b0 et b1
+%--------------------------------------------
+
+% sauvegarde de la figure
+printer='ps';
+print_printer=['-d', printer];
+fullfilename=['.' filesep file '_' mfilename '_'  num2str(ifigure) '.' printer];
+print(print_printer,fullfilename);
+clear print_printer;
+clear fullfilename;
+clear printer;
+%
+clear ifigure;
+
+return
+
Index: /trunk/adm/makefile
===================================================================
--- /trunk/adm/makefile	(revision 20)
+++ /trunk/adm/makefile	(revision 30)
@@ -76,8 +76,8 @@
 $(DIRSRC)/PROGRAMMES/eof_NCEP_sais.m \
 $(DIRSRC)/PROGRAMMES/whorldmap.m \
-$(DIRSRC)/PROGRAMMES/Reg_seriestemporelles.m \
+$(DIRSRC)/PROGRAMMES/reg_seriestemporelles.m \
 $(DIRSRC)/RESULTATS/carteetstd_SST.m \
 $(DIRSRC)/RESULTATS/eof_NCEP_an.m \
-$(DIRSRC)/RESULTATS/Reg_seriestemporelles.m
+$(DIRSRC)/RESULTATS/regandsignif_seriestemporelles.m
 
 LIST_SRCMATLAB_RST = \
