Changeset 13

Timestamp:

08/24/09 13:04:41 (15 years ago)

Author:

pinsard

Message:

eof_NCEP_an : : header and octave compatibility and save figures

File:

• trunk/RESULTATS/eof_NCEP_an.m (modified) (6 diffs)

trunk/RESULTATS/eof_NCEP_an.m

@@ -1 @@
-% Ce programme calcule les EOFs et composantes principales du champs de temperature de surface en moyennes mensuelles climatologiques
+%EOF_NCEP_AN calcule les EOFs et composantes principales du champs de temperature de surface en moyennes mensuelles climatologiques
 % (12 points temporels)
-
+%
 % le programme enregistre automatiquement les resultats du calcul
 % (attention: penser a changer le nom du fichier de sortie en cas de
 % modification de parametres du programme: domaine spatial, moyenne
 % temporelle effectuee: annee / mois, etc. parametre "nomfic" ci dessous)
-
+%
 % le programme permet egalement de visualiser les resultats. Voir le doc
 % .pdf joint pour l'interpretation de cette visualisation. 
-
+%
 % Date: Avril-Mai 2008
 % Auteur: juliette.mignot@locean-ipsl.upmc.fr
 % Contexte Ecole d'ete ENVI-STAT 2008
+%+
+%
+% module
+% ======
+%
+% ``eof_NCEP_an`` 
+%
+% DESCRIPTION
+% ===========
+%
+% ``eof_NCEP_an`` lit les fichiers 
+% ``../DONNEES/lsmask_tropatl_PacandMedblancs.mat``
+% et ``../DONNEES/skt.mon.tropatl.mat``.
+%
+% ``eof_NCEP_an`` calcule ++blabla++
+%
+% ``eof_NCEP_an`` sauvegarde dans ``../DONNEES/EOF_SSTmens_tropAtl++.mat`` 
+% ou ``++/DONNEES/EOF_SSTan++.mat`` ++
+%
+% Ces fichiers sont lus par ``Reg_serietemporelles.m_``.
+%
+% ``eof_NCEP_an`` affiche trois figures :
+%  - ++
+%  - ++
+%  - valeurs propres de la matrice de covariance
+%
+% ``eof_NCEP_an`` sauve ces images dans ``./EOF_SSTmens_tropAtl++_[123].eps``.
+%
+% ++
+%
+% EXAMPLES
+% ========
+%
+% ::
+%
+% >> addpath('../PROGRAMMES/')
+% >> tpacpandreg_startup
+% >> eof_NCEP_an
+%
+% SEE ALSO
+% ========
+%
+% tpacpandreg_startup.m_
+%
+% .. _tpacpandreg_startup.m : tpacpandreg_startup.m.html
+%
+% initfig.m_
+% .. _initfig.m : initfig.m.html
+%
+% Reg_seriestemporelles.m_
+%
+% .. _Reg_seriestemporelles.m : Reg_seriestemporelles.m.html
+%
+% whorldmap.m_
+%
+% .. _whorldmap.m : whorldmap.m.html
+%
+% TODO
+% ====
+%
+% improve description
+%
+% tester la branche NetCDF
+%
+% sous octave, triangle bizarre dans la dernière barre du graphique valeurs 
+% propres.
+%
+% EVOLUTIONS
+% ==========
+%
+% $Id$
+%
+% - 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 
+%     default extension is ``.mat``, write more precise save instruction
+%
+%
+%   * add header
+%   * répertoire de données relatifs au répertoire courant
+%   * 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'on 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
+%
+%-
+
 
 %%
-clear all;
 N=4; %nb d'EOF a tracer
 
@@ -29 +121 @@
 nannees=60;
    
-%repertoire=['C:' filesep 'TPACPandReg' filesep 'DONNEES' filesep] % a adapter a chaque cas: chemin d acces aux donnees
-   repertoire=['/usr/temp/jmignot/TPOCTAVE/TPACPandReg/DONNEES/'] % a adapter a chaque cas: chemin d acces aux donnees
+repertoire=[ '..' filesep 'DONNEES' filesep ] % a adapter a chaque cas: chemin d acces aux donnees
 
 % Mask
@@ -171 +262 @@
 end;
 
-eval(['save ' repertoire nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S TABEOF PC p_var lon lat mask'])
+fullfilename=[repertoire nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S.mat'];
+if (run_octave == 0)
+ save(fullfilename,'TABEOF','PC','p_var','lon','lat','mask','-v7');
+else
+ save('-v7',fullfilename,'TABEOF','PC','p_var','lon','lat','mask');
+end
+clear fullfilename;
 
 
  %%               TRACE
-figure(1); initfig;
-figure(2); initfig;
+figure(1);
+initfig;
+figure(2); 
+initfig;
 
 for i=1:4
     
    %subplot(4,2,2*i-1)
-   figure(1)
+   ifigure=1;
+   figure(ifigure);
    subplot(2,2,i)
-   contourf(lon,lat,squeeze(TABEOF(:,:,i)))
+   if (run_octave == 0)
+    x=lon;
+    y=lat;
+   else
+    [x,y]=meshgrid(lon,lat);
+   end
+   z=squeeze(TABEOF(:,:,i));
+   contourf(x,y,z);
    caxis([-0.06 0.06])
    hold on
@@ -193 +300 @@
    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'])
+
+   % sauvegarde de la figure
+   printer='eps';
+   print_printer=['-d', printer];
+   fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer];
+   print(print_printer,fullfilename);
+   clear fullfilename;
+   clear printer;
    
    %subplot(4,2,2*i)
-   figure(2)
+   ifigure=2;
+   figure(ifigure);
    subplot(2,2,i)
    plot(PC(:,i),'x-')
@@ -203 +319 @@
    xlim([1 ntemp])
    xlabel('mois')
+
+   % sauvegarde de la figure
+   printer='eps';
+   print_printer=['-d', printer];
+   fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer];
+   print(print_printer,fullfilename);
+   clear fullfilename;
+   clear printer;
+
 end;
 
 
 % Graphe pourcentages de variance expliqué
-figure; 
+ifigure=3;
+figure(ifigure); 
 bar(flipud(pctg_var))
 xlim([0 10])
@@ -214 +340 @@
 xlabel('indice du vecteur propre')
 
+% sauvegarde de la figure
+printer='eps';
+print_printer=['-d', printer];
+fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer];
+print(print_printer,fullfilename);
+clear fullfilename;
+clear printer;
+