Changes in / [20:30]


Ignore:
Location:
/trunk
Files:
4 added
2 deleted
8 edited

Legend:

Unmodified
Added
Removed

  • /trunk/PROGRAMMES/carte_SST.m

    • Property svn:keywords set to Id
    r20 r30  
    1 % CARTE_SST  Ce programme permet de visualiser les cartes mensuelles moyennes de temperature de surface        
     1%CARTE_SST Ce programme permet de visualiser les cartes mensuelles moyennes de temperature de surface 
    22% 
    33% Date: Avril-Mai 2008 
    … …  
    1515% =========== 
    1616% 
    17 % ``carte_SST`` visualise les cartes mensuelles moyennes de température  
     17% ``carte_SST`` visualise les cartes mensuelles moyennes de température 
    1818% de surface. 
    1919% 
    20 % ``carte_SST`` lit les fichiers  
     20% ``carte_SST`` lit les fichiers 
    2121% ``../DONNEES/lsmask_tropatl_PacandMedblancs.mat`` 
    2222% et ``../DONNEES/skt.mon.tropatl.mat``. 
    2323% 
    24 % ``carte_SST`` calcule ++blabla++ 
    25 % 
    26 % ``carte_SST`` affiche à l'écran les 12 mois de ++ 
    27 % 
    28 % ``carte_SST`` sauve cette image dans ``./carte_SST.eps``. 
     24% ``carte_SST`` calcule les moyennes mensuelles de SST dans 
     25% l'Atlantique tropical a partir d'un champ mensuel de 60 ans 
     26% fourni en entree 
     27% 
     28% ``carte_SST`` affiche à l'écran les moyennes mensuelles de SST. 
     29% 
     30% ``carte_SST`` sauve cette image dans ``./carte_SST.ps``. 
    2931% 
    3032% EXAMPLES 
    … …  
    4749% .. _initfig.m : initfig.m.html 
    4850% 
     51% whorldmap.m_ 
     52% 
     53% .. _whorldmap.m : whorldmap.m.html 
     54% 
     55% colorbartype.m_ 
     56% 
     57% .. _colorbartype.m : colorbartype.m.html 
     58% 
     59% contlab.m_ 
     60% 
     61% .. _contlab.m : contlab.m.html 
     62% 
    4963% TODO 
    5064% ==== 
    … …  
    5266% improve description 
    5367% 
     68% KNOWN PROBLEMS 
     69% ============== 
     70% 
     71% - fplod 2009-08-26T08:46:48Z aedon.locean-ipsl.upmc.fr (Darwin) 
     72% 
     73%   * ne fonctionne pas sous octave 3.2.2 
     74%     :: 
     75%     error: gnuplot (as of v4.2) only supports 3D filled triangular patches 
     76%     error: called from: 
     77%     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 
     78%     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 
     79%     error:   /usr/local_linux/octave/64/3.2.2//share/octave/3.2.2/m/plot/gnuplot_drawnow.m at line 88, column 5 
     80%     error:   /usr/local_linux/octave/64/3.2.2//share/octave/3.2.2/m/plot/print.m at line 696, column 7 
     81%     error:   /.autofs/home/fplod/incas/tpoctave/tpacpandreg_ws/PROGRAMMES/carte_SST.m at line 212, column 1 
     82% 
     83% 
    5484% EVOLUTIONS 
    5585% ========== 
    … …  
    5787% $Id$ 
    5888% 
     89% - fplod 2009-08-26T09:00:11Z aedon.locean-ipsl.upmc.fr (Darwin) 
     90% 
     91%   * titre: nom de chaque mois 
     92%   * remplacement de ``eps`` par ``ps`` pour pouvoir imprimer les figures 
     93%     produites par ``octave`` 
     94% 
     95%     Dans le fichier .eps produit par ``octave`` on pouvait lire 
     96%     ``!PS-Adobe-2.0 EPSF-2.0``. 
     97% 
     98%     Aucun message d'erreur n'apparaît suite à la demande d'impression 
     99%     ``zeus$ lpr -o job-sheets='none' carte_SST.ps``, mais aucune feuille ne 
     100%     sort. 
     101%  
     102%     Pour info dans le fichier .eps produit par matlab, on pouvait lire 
     103%     ``%!PS-Adobe-2.0 EPSF-1.2``. Aucun souci d'impression avec ce fichier ! 
     104% 
     105%     En remplaçant ``eps`` par ``ps``, on n'a plus de souci d'impression. 
     106% 
     107% 
     108% - jmignot 2009-08-25 
     109% 
     110%   * suppression d'appel à contlab et colorbartype sous octave 
     111%  
     112% - fplod 2009-08-25T10:40:16Z aedon.locean-ipsl.upmc.fr (Darwin) 
     113% 
     114%   * modif commentaires liée au dernier chgt de Juliette 
     115% 
     116% - jmignot 2009-08-25 
     117% 
     118%   * remplacement de ``colorbar`` par ``colorbartype`` 
     119%     parce pas beau ni en matlab ni en octave (dégradé, noir et blanc dans 
     120%     les .jpg) 
     121%   * contraintes sur les limites x et y de axes car matlab et octave ne 
     122%     fixent pas les même limites par défaut 
     123%   * utilisation de ``contlab`` pour ++ 
     124%  
    59125% - fplod 2009-08-17T12:39:51Z aedon.locean-ipsl.upmc.fr (Darwin) 
    60126% 
    … …  
    63129%   * suppression de ``clear all`` présent dans ``tpacpandreg_startup.m`` 
    64130%   * contournement pb octave pour ``contourf`` : 
    65 %     Dans le fil de discussion [Pkg-octave-devel] Bug#492223: octave3.0: contourf only accepts vector X and Y if same+++,  
    66 %     je comprends que j'ai affaire à un bug de octave 3.0 (donc à surveiller  
    67 %     pour les prochaines versions) qui ne veut pas faire de contourf si les  
    68 %     variables X et Y n'ont pas la même dimension et c'est le cas ici  
    69 %     (lat = 42x1, lon= 70x1 vue avec size).  
     131%     Dans le fil de discussion [Pkg-octave-devel] Bug#492223: octave3.0: contourf only accepts vector X and Y if same+++, 
     132%     je comprends que j'ai affaire à un bug de octave 3.0 (donc à surveiller 
     133%     pour les prochaines versions) qui ne veut pas faire de contourf si les 
     134%     variables X et Y n'ont pas la même dimension et c'est le cas ici 
     135%     (lat = 42x1, lon= 70x1 vue avec size). 
    70136%     http://lists.alioth.debian.org/pipermail/pkg-octave-devel/2008-July/004672.html 
    71137%   * sauvegarde de la figure dans ``./carte_SST.eps`` 
    … …  
    85151%Mask 
    86152%----- 
    87    file='lsmask_tropatl_PacandMedblancs' 
     153   file='lsmask_tropatl_PacandMedblancs.mat' 
    88154   load([repertoire file]); 
    89155% Data 
    … …  
    101167 
    102168%--------------------------------------------------- 
    103 %  MANIPULATIONS TEMPORELLES: MOYENNES  
     169%  MANIPULATIONS TEMPORELLES: MOYENNES 
    104170%--------------------------------------------------- 
    105171tab=reshape(tab,size(tab,1),size(tab,2),12,nannees); %organisation du tableau sous la forme [lat lon mois annees] 
    … …  
    114180%%%%               TRACE 
    115181%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
    116  
    117  
     182%parametres geographiques du graphe 
     183  geo=[-100 20 -30 30]; %[lonmin lonmax latmin latmax]; 
     184  yti=[-30:10:30]; xti=[-100:30:20]; %labels de latitude et de longitude 
     185 
     186                 
    118187% Variabilité saisonnière 
    119 figure;  
    120 initfig;  
     188%------------------------- 
     189%niveaux de couleurs 
     190lev=[0:1:30];%contours a tracer, en unite du champ a dessiner, ici 
     191             %la SST en degC 
     192             %a adapter eventuellement 
     193clip_lev=[20 30];%definit les limites des couleurs: bleu fonce au 
     194                 %rouge fonce 
     195             %a adapter eventuellement 
     196 
     197 
     198ifigure=1; 
     199figure(ifigure); 
     200clf(); 
     201initfig; 
    121202for iplot=1:12 
    122203  subplot(4,3,iplot) 
    … …  
    128209  end 
    129210  z=squeeze(tab_mens(:,:,iplot)); 
    130   contourf(x,y,z,[0:1:30]) 
    131   caxis([20 30]) 
     211  contourf(x,y,z,lev) 
     212  caxis(clip_lev) 
    132213  hold on 
    133   ylim([-30 30]) 
    134   xlim([-100 20]) 
    135   title(['mois ' num2str(iplot)]) 
    136 end;  
    137 colorbar('East') 
     214  whorldmap 
     215 
     216  % limites géographiques et labels de longitude et de latitude 
     217  ylim([geo(3) geo(4)]) 
     218  xlim([geo(1) geo(2)]) 
     219  set(gca, 'Xtick', xti, 'Ytick', yti); 
     220  if (run_octave == 0) 
     221    [longstr, latstr] = contlab(xti, yti); 
     222    set(gca, 'Xticklabel', longstr, 'Yticklabel', latstr); 
     223  end; 
     224  if iplot>=10 
     225    xlabel('longitude') 
     226  end; 
     227  if iplot==1|iplot==4|iplot==7|iplot==10 
     228    ylabel('latitude') 
     229  end; 
     230 
     231  %titre: nom de chaque mois 
     232  cmois = datestr(datenum([2000 iplot 1 00 00 00]), 'mmm'); 
     233  title(cmois); 
     234  clear cmois; 
     235 
     236end; 
     237 
     238%barre de couleur 
     239if (run_octave == 0) 
     240  pos=[.2 .03 .6 .013]; 
     241  hpal=jet(100); %carte des couleurs: du bleu au rouge par defaut. 
     242  colorbartype(pos,lev,1,clip_lev,hpal,0) 
     243  title('SST mensuelle moyenne en degC') 
     244else 
     245 colorbar('East') %colorbar par defaut 
     246end; 
    138247 
    139248% sauvegarde de la figure 
    140 printer='eps'; 
     249printer='ps'; 
    141250print_printer=['-d', printer]; 
    142251fullfilename=['./', mfilename, '.', printer]; 

  • /trunk/PROGRAMMES/eof_NCEP_sais.m

    • Property svn:keywords set to Id
    r20 r30  
    88% 
    99% le programme permet egalement de visualiser les resultats. Voir le doc 
    10 % .pdf joint pour l'interpretation de cette visualisation.  
     10% .pdf joint pour l'interpretation de cette visualisation. 
    1111% 
    1212% Date: Avril-Mai 2008 
    … …  
    1919% ====== 
    2020% 
    21 % ``eof_NCEP_sais``  
     21% ``eof_NCEP_sais`` 
    2222% 
    2323% DESCRIPTION 
    2424% =========== 
    2525% 
    26 % ``eof_NCEP_sais`` lit les fichiers  
     26% ``eof_NCEP_sais`` lit les fichiers 
    2727% ``../DONNEES/lsmask_tropatl_PacandMedblancs.mat`` 
    2828% et ``../DONNEES/skt.mon.tropatl.mat``. 
    … …  
    3030% ``eof_NCEP_sais`` calcule ++blabla++ 
    3131% 
    32 % ``eof_NCEP_sais`` sauvegarde dans ``../DONNEES/EOF_SSTmens_tropAtl30N-20S.mat``. 
    33 % 
    34 % Ces fichiers sont lus par ``Reg_serietemporelles.m_``\ . 
     32% ``eof_NCEP_sais`` sauvegarde dans 
     33% ``../DONNEES/eof_SSTmens_tropAtl30N-20S.mat``. 
     34% 
     35% Ces fichiers sont lus par ``reg_serietemporelles.m_``\ . 
    3536% 
    3637% ``eof_NCEP_sais`` affiche trois figures : 
    37 % 
    3838%  - ++ 
    3939%  - ++ 
    4040%  - valeurs propres de la matrice de covariance 
    4141% 
    42 % ``eof_NCEP_sais`` sauve ces images dans ``./EOF_SSTmens_tropAtl30N-20S_[123].eps``. 
     42% ``eof_NCEP_sais`` sauve ces images dans 
     43% ``./eof_SSTmens_tropAtl30N-20S_[123].ps``. 
    4344% 
    4445% ++ 
    … …  
    6364% .. _initfig.m : initfig.m.html 
    6465% 
    65 % Reg_seriestemporelles.m_ 
    66 % 
    67 % .. _Reg_seriestemporelles.m : Reg_seriestemporelles.m.html 
     66% reg_seriestemporelles.m_ 
     67% 
     68% .. _reg_seriestemporelles.m : reg_seriestemporelles.m.html 
    6869% 
    6970% whorldmap.m_ 
    7071% 
    7172% .. _whorldmap.m : whorldmap.m.html 
     73% 
     74% colorbartype.m_ 
     75% 
     76% .. _colorbartype.m : colorbartype.m.html 
     77% 
     78% contlab.m_ 
     79% 
     80% .. _contlab.m : contlab.m.html 
    7281% 
    7382% TODO 
    … …  
    7887% tester la branche NetCDF 
    7988% 
     89% comprendre le pb octave 3.2.2 
     90% 
    8091% EVOLUTIONS 
    8192% ========== 
    … …  
    8394% $Id$ 
    8495% 
    85 % - fplod 2009-08-25T09:37:04Z aedon.locean-ipsl.upmc.fr (Darwin) 
    86 % 
    87 %   * typo in ReStructured Text header 
     96% - fplod 2009-08-26T16:34:59Z aedon.locean-ipsl.upmc.fr (Darwin) 
     97% 
     98%   * pb contourf sous octave 3.2.2 
     99% 
     100% - fplod 2009-08-26T11:05:36Z aedon.locean-ipsl.upmc.fr (Darwin) 
     101% 
     102%   * remplacement de ``eps`` par ``ps`` pour pouvoir imprimer les figures 
     103%     produites par ``octave`` 
     104% 
     105%     Dans le fichier .eps produit par ``octave`` on pouvait lire 
     106%     ``!PS-Adobe-2.0 EPSF-2.0``. 
     107% 
     108%     Aucun message d'erreur n'apparaît suite à la demande d'impression 
     109%     ``zeus$ lpr -o job-sheets='none' eof_SSTmens_tropAtl30N-20S_1.octave.ps``, 
     110%     mais aucune feuille ne sort. 
     111% 
     112%     Pour info dans le fichier .eps produit par matlab, on pouvait lire 
     113%     ``%!PS-Adobe-2.0 EPSF-1.2``. Aucun souci d'impression avec ce fichier ! 
     114% 
     115%     En remplaçant ``eps`` par ``ps``, on n'a plus de souci d'impression. 
     116% 
     117% - jmignot 2009-08-25 
     118% 
     119%   * suppression d'appel à contlab et colorbartype sous octave 
     120% 
     121% 
     122% - fplod 2009-08-25T10:40:16Z aedon.locean-ipsl.upmc.fr (Darwin) 
     123% 
     124%   * modif commentaires liée au dernier chgt de Juliette 
     125% 
     126% - jmignot 2009-08-25 
     127% 
     128%   * chgt de noms de fichiers : ``EOF_`` devient ``eof_`` 
     129%   * remplacement de ``colorbar`` par ``colorbartype`` 
     130%     parce pas beau ni en matlab ni en octave (dégradé, noir et blanc dans 
     131%     les .jpg) 
     132%   * contraintes sur les limites x et y de axes car matlab et octave ne 
     133%     fixent pas les même limites par défaut 
     134%   * utisation de ``contlab`` pour ++ 
    88135% 
    89136% - fplod 2009-08-21T15:04:16Z aedon.locean-ipsl.upmc.fr (Darwin) 
    90137% 
    91 %   * as octave save default format is ASCII and no default extension  
    92 %     while matlab save default format is MAT v5 mat-file (little endian) and  
     138%   * as octave save default format is ASCII and no default extension 
     139%     while matlab save default format is MAT v5 mat-file (little endian) and 
    93140%     default extension is ``.mat``, write more precise save instruction 
    94141% 
    … …  
    100147%   * contournement pb octave pour ``contourf`` : 
    101148%     Dans le fil de discussion [Pkg-octave-devel] Bug#492223: octave3.0: contourf only accepts vector X and Y if same+++, 
    102 %     je comprends que j'ai affaire à un bug de octave 3.0 (donc à surveiller  
    103 %     pour les prochaines versions) qui ne veut pas faire de contourf si les  
    104 %     variables X et Y n'ont pas la même dimension et c'est le cas ici  
    105 %     (lat = 42x1, lon= 70x1 vue avec size).  
     149%     je comprends que j'ai affaire à un bug de octave 3.0 (donc à surveiller 
     150%     pour les prochaines versions) qui ne veut pas faire de contourf si les 
     151%     variables X et Y n'ont pas la même dimension et c'est le cas ici 
     152%     (lat = 42x1, lon= 70x1 vue avec size). 
    106153%     http://lists.alioth.debian.org/pipermail/pkg-octave-devel/2008-July/004672.html 
    107154%   * sauvegarde des figures 
     155 
    108156% 
    109157%- 
    … …  
    111159N=4; %nb d'EOF a tracer 
    112160 
    113 nomfic='EOF_SSTmens_tropAtl' 
    114 % mettre 'EOF_SSTan_tropAtl' quand on travaille en moyennes annuelles. Ce 
     161nomfic='eof_SSTmens_tropAtl' 
     162% mettre 'eof_SSTan_tropAtl' quand on travaille en moyennes annuelles. Ce 
    115163% fichier sera rappele dans les programme de regression lineaire 
    116164%% 
    … …  
    125173% Mask 
    126174%------ 
    127   file='lsmask_tropatl_PacandMedblancs' 
     175  file='lsmask_tropatl_PacandMedblancs.mat' 
    128176   load([repertoire file]); 
    129   ny=size(mask,1); nx=size(mask,2); % dimensions spatiales des tableaux  
     177  ny=size(mask,1); nx=size(mask,2); % dimensions spatiales des tableaux 
    130178 
    131179% Data 
    … …  
    136184  load([repertoire2 file]); % chargement des donnees 
    137185 
    138     
     186 
    139187  %SI NETCDF 
    140188  %repertoire2=[repertoire 'NETCDF' filesep]% a adapter a chaque cas: chemin d acces aux donnees 
    … …  
    142190  % f=netcdf([repertoire file],'nowrite'); 
    143191  % tab=f{'skt'}(:,:,:); 
    144    
     192 
    145193 
    146194   tab=permute(tab,[2 3 1]);% permutation des dimensions pour organiser le tableau sous la forme [lat lon temps] 
    … …  
    169217%---------------- 
    170218%Si on veut éventuellement attribuer un poids différent aux différents 
    171 %points spatiaux, on le fait ici.  
     219%points spatiaux, on le fait ici. 
    172220 
    173221% Préparation du zoom géographique: sur le masque, j'identifie les points 
    … …  
    177225limite_latS=-20; %limite du domaine d'etuade en latitude 
    178226limite_latN=30; %limite du domaine d'etuade en latitude 
    179 lat2D=repmat(lat,[1 size(mask,2)]);  
     227lat2D=repmat(lat,[1 size(mask,2)]); 
    180228mask(lat2D>limite_latN&mask==0)=4; 
    181229mask(lat2D<limite_latS&mask==0)=4; 
    … …  
    209257% matrice de covariance: 
    210258  C = Z * Z'; 
    211 %  eval(['save ' repertoire 'MatCov_SST_tropAtl' num2str(limite_latN) 'N' num2str(limite_latS) 'S.asc C -ascii'])  
    212    
    213 % diagonalisation   
    214   [E,Lambda] = eig(C);  
    215    % Lambda est la matrice diagonale des valeurs propres et E la matrice dont les colonnes sont les vecteurs propres correspondants.   
     259%  eval(['save ' repertoire 'MatCov_SST_tropAtl' num2str(limite_latN) 'N' num2str(limite_latS) 'S.asc C -ascii']) 
     260 
     261% diagonalisation 
     262  [E,Lambda] = eig(C); 
     263   % Lambda est la matrice diagonale des valeurs propres et E la matrice dont les colonnes sont les vecteurs propres correspondants. 
    216264   % Soit Z*E = E*Lambda 
    217    % Les vecteurs propres sont appelés les EOF. Chaque EOF peut être considéré comme une carte de l'espace physique.     
    218  
    219     
    220  %%  
     265   % Les vecteurs propres sont appelés les EOF. Chaque EOF peut être considéré comme une carte de l'espace physique. 
     266 
     267 
     268 %% 
    221269 %%%% calcul de la matrice donnant les pourcentages de variance 
    222270Lambda = diag(Lambda); 
    223271pctg_var = round(Lambda / sum(Lambda) * 1000) / 10; 
    224272 
    225   
     273 
    226274 %On definit les composantes principales de Q comme la decomposition du signal decrit par Q sur les EOF E. 
    227275 
    228276   P = Z' * E; 
    229      
    230 %%    
     277 
     278%% 
    231279%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
    232280%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
    … …  
    271319 
    272320 %%               TRACE 
    273 figure(1);  
     321 
     322close all 
     323%parametres geographiques des graphes 
     324  geo=[-100 20 -30 30]; %[lonmin lonmax latmin latmax]; 
     325  yti=[-30:10:30]; xti=[-100:30:20]; %labels de latitude et de longitude 
     326%niveaux de couleurs 
     327  lev=[-2:.02:2]; 
     328  clip_lev=[-0.06 0.06]; 
     329 
     330figure(1); 
    274331initfig; 
    275 figure(2);  
     332figure(2); 
    276333initfig; 
    277334 
    278335for i=1:4 
    279      
     336 
    280337   %subplot(4,2,2*i-1) 
    281338   ifigure=1; 
    282339   figure(ifigure); 
    283340   subplot(2,2,i) 
     341   z=squeeze(TABEOF(:,:,i)); 
    284342   if (run_octave == 0) 
    285343    x=lon; 
    286344    y=lat; 
     345    contourf(x,y,z,lev); 
    287346   else 
    288     [x,y]=meshgrid(lon,lat); 
     347    switch OCTAVE_VERSION 
     348     case {'3.0.2'} 
     349      [x,y]=meshgrid(lon,lat); 
     350      contourf(x,y,z,lev); 
     351     case {'3.2.2'} 
     352      warning('www : do no know about octave 3.2.2 contourf usage'); 
     353%  A(I): Index exceeds matrix dimension ++ 
     354     otherwise 
     355      warning('www : do no know about octave contourf usage'); 
     356     end 
    289357   end 
    290    z=squeeze(TABEOF(:,:,i)); 
    291    contourf(x,y,z); 
    292    caxis([-0.06 0.06]) 
     358   caxis(clip_lev); 
    293359   hold on 
    294    xlim([-100 30]); xlabel('longitude') 
    295    ylim([-40 40]); ylabel('latitude') 
    296    whorldmap([-100, 30, -40, 40]) 
    297    colorbar 
    298    title([num2str(p_var(i)) '%  -  (unites: K)']) 
    299    set(gca,'xtick',[-100;-70;-40;-10;20],'xticklabel',['100W';' 70W';' 40W';' 10W';' 20E']) 
    300    set(gca,'ytick',[-30;-15;0;15;30],'yticklabel',['30S';'15S';' 0 ';'15N';'30N']) 
    301    % 
     360   whorldmap 
     361   % limites geographies et labels de longitude et de latitude 
     362  ylim([geo(3) geo(4)]) 
     363  xlim([geo(1) geo(2)]) 
     364  set(gca, 'Xtick', xti, 'Ytick', yti) 
     365  if (run_octave == 0) 
     366    [longstr, latstr] = contlab(xti, yti); 
     367    set(gca, 'Xticklabel', longstr, 'Yticklabel', latstr) 
     368  end; 
     369  xlabel('longitude') 
     370  ylabel('latitude') 
     371 
     372   %titre de la figure 
     373   title([num2str(p_var(i)) '%']) 
     374 
     375  %barre de couleur 
     376  if (run_octave == 0) 
     377    hpal=jet(100); %carte des couleurs: du bleu au rouge par defaut. 
     378    pos1=get(gca,'Position'); 
     379    pos=[pos1(1) pos1(2)-.075 pos1(3) .013]; 
     380    colorbartype(pos,lev,1,clip_lev,hpal,0); 
     381    title('degC') 
     382  else 
     383    %++colorbar('East') 
     384  end; 
     385 
    302386   % sauvegarde de la figure 
    303    printer='eps'; 
     387   printer='ps'; 
    304388   print_printer=['-d', printer]; 
    305389   fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer]; 
    … …  
    307391   clear fullfilename; 
    308392   clear printer; 
    309     
     393 
    310394   %subplot(4,2,2*i) 
    311395   ifigure=2; 
    … …  
    320404   % 
    321405   % sauvegarde de la figure 
    322    printer='eps'; 
     406   printer='ps'; 
    323407   print_printer=['-d', printer]; 
    324408   fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer]; 
    … …  
    330414% Graphe pourcentages de variance expliqué 
    331415ifigure=3; 
    332 figure(ifigure);  
     416figure(ifigure); 
     417initfig; 
    333418bar(flipud(pctg_var)) 
    334419xlim([0 10]) 
    … …  
    338423 
    339424% sauvegarde de la figure 
    340 printer='eps'; 
     425printer='ps'; 
    341426print_printer=['-d', printer]; 
    342427fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer]; 

  • /trunk/PROGRAMMES/initfig.m

    r20 r30  
    4242% $Id$ 
    4343% 
     44% - fplod 2009-08-25T16:41:34Z aedon.locean-ipsl.upmc.fr (Darwin) 
     45% 
     46%   * ajout pour octave 3.2 (uniquement sous  
     47%     /usr/local_linux/octave/64/3.2.2/bin/octave donc dedale) 
     48% 
    4449% - fplod 2009-08-17T12:39:51Z aedon.locean-ipsl.upmc.fr (Darwin) 
    4550% 
    … …  
    5661 set(gcf,'PaperPosition',[0,0,1,1]); 
    5762else 
    58  warning('www : octave can not initialise Units, position, PaperType PaperUnits PaperPosition'); 
     63 switch OCTAVE_VERSION 
     64  case {'3.0.2'} 
     65   warning('www : octave 3.0.2 can not initialise Units, position, PaperType PaperUnits PaperPosition'); 
     66  case {'3.2.2'} 
     67   set(gcf,'Units','normalized'); 
     68   set(gcf,'position',[0.1 0.08,0.5,0.8]) 
     69   set(gcf,'PaperType', 'a4'); 
     70   set(gcf,'PaperUnits','normalized'); 
     71   set(gcf,'PaperPosition',[0,0,1,1]); 
     72  otherwise 
     73   warning('www : do no know about octave initialise Units, position, PaperType PaperUnits PaperPosition'); 
     74  end 
    5975end 
    6076 

  • /trunk/PROGRAMMES/tpacpandreg_startup.m

    r20 r30  
    1 % TPACPANDREG_STARTUP initialisation of matlab/octave environment for TPACPANDREG project 
     1%TPACPANDREG_STARTUP initialisation of matlab/octave environment for TPACPANDREG project 
    22 
    33%+ 
    … …  
    3333% ======== 
    3434% 
    35 % tpacpandreg_profile.sh_ 
    36 % 
    37 % .. _tpacpandreg_profile.sh : tpacpandreg_profile.sh.html 
    38 % 
    3935% EVOLUTIONS 
    4036% ========== 
    4137% 
    4238% $Id$ 
     39% 
     40% - fplod 2009-08-25T10:36:26Z aedon.locean-ipsl.upmc.fr (Darwin) 
     41% 
     42%   * suppression de ``more on`` (parce que pas pratique pour les étudiants) 
    4343% 
    4444% - fplod 2009-08-21T10:37:54Z aedon.locean-ipsl.upmc.fr (Darwin) 
    … …  
    5252%  - path for matlab/octave commands used in TPACPANDREG project 
    5353%  - IO directories in global variables 
    54 % 
    55 % not yet done because of specific meaning of private and toolbox in matlab 
    5654% 
    5755% check if set_param(0,'CharacterEncoding','ISO-8859-1') is a better solution 
    … …  
    7472end 
    7573% 
    76 more on 

  • /trunk/PROGRAMMES/whorldmap.m

    • Property svn:keywords set to Id
    r20 r30  
    4040%  
    4141% $Id$ 
     42% 
     43% - fplod 2009-08-26T09:02:22Z aedon.locean-ipsl.upmc.fr (Darwin) 
     44% 
     45%   * add missing octave warning on Erasemode 
    4246% 
    4347% - fplod 2009-08-25T09:32:28Z aedon.locean-ipsl.upmc.fr (Darwin) 
    … …  
    139143    set(lh,'Erasemode','none'); 
    140144   else 
    141     warning('www : octave do not know Erasemode property') 
     145    warning('www : octave do not know Erasemode property'); 
    142146   end 
    143147   set(gca,'aspect',[1 1]); 
    … …  
    145149 
    146150if (drawtwice), 
    147    line(coastlines(1,:),coastlines(2,:),'Linestyle','-',... 
    148    'Erasemode','none'); 
     151   if ~ exist('OCTAVE_VERSION','builtin') 
     152    line(coastlines(1,:),coastlines(2,:),'Linestyle','-',... 
     153    'Erasemode','none'); 
     154   else 
     155    warning('www : octave do not know Erasemode property'); 
     156    line(coastlines(1,:),coastlines(2,:),'Linestyle','-') 
     157   end 
    149158end; 
    150159 

  • /trunk/RESULTATS/carteetstd_SST.m

    • Property svn:keywords set to Id
    r20 r30  
    11%CARTEETSTD_SST permet de visualiser les cartes mensuelles moyennes 
    22% de temperature de surface        
    3 % ainsi que les cartes annuelles moyennes et les ecart type 
     3% ainsi que les cartes annuelles moyennes et les écart type 
    44% correspondants 
    55 
    … …  
    1919% ``carteetstd_SST`` permet de visualiser les cartes mensuelles moyennes 
    2020% de temperature de surface 
    21 % ainsi que les cartes annuelles moyennes et les ecart type 
     21% ainsi que les cartes annuelles moyennes et les écart type 
    2222% correspondants. 
    2323% 
    … …  
    2626% et ``../DONNEES/skt.mon.tropatl.mat``. 
    2727% 
    28 % ``carteetstd_SST`` calcule ++blabla++ 
    29 % 
    30 % ``carteetstd_SST`` affiche à l'écran les 12 mois de ++ 
    31 % 
    32 % ``carteetstd_SST`` sauve ces images dans ``./carteetstd_SST_[123].eps``. 
     28% ``carteetstd_SST`` calcule les moyennes mensuelles d'une part et 
     29% les moyennes interannuelles d'autre part de la SST dans 
     30% l'Atlantique tropical à partir d'un champ mensuel de 60 ans 
     31% fourni en entrée. 
     32% 
     33% ``carteetstd_SST`` affiche à l'écran les 12 valeurs de moyennes 
     34% mensuelles (figure 1), une sélection de 12 valeurs annuelles 
     35% moyennes (figure 2) et les écart-type mensuels et interannuels 
     36% (figure 3). 
     37% 
     38% ``carteetstd_SST`` sauve ces images dans ``./carteetstd_SST_[123].ps``. 
    3339% 
    3440% EXAMPLES 
    … …  
    4955% 
    5056% initfig.m_ 
     57% 
    5158% .. _initfig.m : initfig.m.html 
     59% 
     60% whorldmap.m_ 
     61% 
     62% .. _whorldmap.m : whorldmap.m.html 
     63% 
     64% colorbartype.m_ 
     65% 
     66% .. _colorbartype.m : colorbartype.m.html 
     67% 
     68% contlab.m_ 
     69% 
     70% .. _contlab.m : contlab.m.html 
     71% 
     72% 
     73% KNOWN PROBLEMS 
     74% ============== 
     75% 
     76% pb avec les diacritiques à l'impression de .ps dans le titres des figures  
     77% sous octave : é se transforme en o barré. 
    5278% 
    5379% TODO 
    … …  
    6187% $Id$ 
    6288% 
     89% - fplod 2009-08-26T09:56:51Z aedon.locean-ipsl.upmc.fr (Darwin) 
     90%    
     91%   * nom des mois à la place du numéro dans les titres 
     92%   * remplacement de ``eps`` par ``ps`` pour pouvoir imprimer les figures  
     93%     produites par ``octave`` 
     94% 
     95% 
     96% - jmignot 2009-08-25 
     97% 
     98%   * suppression d'appel à contlab et colorbartype sous octave 
     99% 
     100% - fplod 2009-08-25T10:40:16Z aedon.locean-ipsl.upmc.fr (Darwin) 
     101% 
     102%   * modif commentaires liée au dernier chgt de Juliette 
     103% 
     104% - jmignot 2009-08-25 
     105% 
     106%   * remplacement de ``colorbar`` par ``colorbartype`` 
     107%     parce pas beau ni en matlab ni en octave (dégradé, noir et blanc dans  
     108%     les .jpg)  
     109%   * contraintes sur les limites x et y de axes car matlab et octave ne 
     110%     fixent pas les même limites par défaut 
     111%   * utisation de ``contlab`` pour ++ 
     112%  
    63113% - fplod 2009-08-17T12:39:51Z aedon.locean-ipsl.upmc.fr (Darwin) 
    64114% 
    … …  
    92142%Mask 
    93143%------ 
    94 file='lsmask_tropatl_PacandMedblancs' 
     144file='lsmask_tropatl_PacandMedblancs.mat' 
    95145   load([repertoire file]); 
    96146% Data 
    … …  
    110160%  MANIPULATIONS TEMPORELLES: MOYENNES  
    111161%--------------------------------------------------- 
    112 tab=reshape(tab,size(tab,1),size(tab,2),12,nannees); %organisation du tableau sous la forme [lat lon mois annees] 
     162tab=reshape(tab,size(tab,1),size(tab,2),12,nannees); %organisation du tableau sous la forme [lat lon mois années] 
    113163 
    114164% Construction de moyennes annuelles 
    … …  
    120170tab_mens=squeeze(mean(tab,4)); 
    121171 
     172%calcul des écart-types 
     173%------------------------ 
     174%écart-type mensuel moyen:  
     175std_mens=squeeze(std(tab_mens,0,3)); 
     176 
     177%écart-type interannuel moyen:  
     178std_an=squeeze(std(tab_an,0,3)); 
    122179 
    123180%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
    … …  
    125182%%%%               TRACE 
    126183%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
    127  
    128  
     184%parametres geographiques du graphe 
     185  geo=[-100 20 -30 30]; %[lonmin lonmax latmin latmax]; 
     186  yti=[-30:10:30]; xti=[-100:30:20]; %labels de latitude et de longitude 
     187 
     188                  
    129189% Variabilité saisonnière 
     190%------------------------- 
     191%niveaux de couleurs 
     192lev=[0:1:30];%contours a tracer, en unite du champ a dessiner, ici 
     193             %la SST en degC 
     194clip_lev=[20 30];%definit les limites des couleurs: bleu fonce au 
     195                 %rouge fonce 
     196 
     197 
    130198ifigure=1; 
    131199figure(ifigure);  
    … …  
    141209  end 
    142210  z=squeeze(tab_mens(:,:,iplot)); 
    143   contourf(x,y,z,[0:1:30]) 
    144   caxis([20 30]) 
     211  contourf(x,y,z,lev) 
     212  caxis(clip_lev) 
    145213  hold on 
    146   ylim([-30 30]) 
    147   xlim([-100 20]) 
    148   title(['mois ' num2str(iplot)]) 
    149 end;  
    150 colorbar('East') 
     214  whorldmap 
     215 
     216  % limites geographies et labels de longitude et de latitude 
     217  ylim([geo(3) geo(4)]) 
     218  xlim([geo(1) geo(2)]) 
     219  set(gca, 'Xtick', xti, 'Ytick', yti) 
     220  if (run_octave == 0) 
     221   [longstr, latstr] = contlab(xti, yti); 
     222   set(gca, 'Xticklabel', longstr, 'Yticklabel', latstr) 
     223  end; 
     224  if iplot>=10 
     225    xlabel('longitude') 
     226  end; 
     227  if iplot==1|iplot==4|iplot==7|iplot==10 
     228    ylabel('latitude') 
     229  end; 
     230   
     231  %titre: nom de chaque mois 
     232  cmois = datestr(datenum([2000 iplot 1 00 00 00]), 'mmm'); 
     233  title(cmois); 
     234  clear cmois; 
     235   
     236end; 
     237 
     238%barre de couleur 
     239if (run_octave == 0) 
     240  pos=[.2 .03 .6 .013]; 
     241  hpal=jet(100); %carte des couleurs: du bleu au rouge par defaut. 
     242  colorbartype(pos,lev,1,clip_lev,hpal,0); 
     243  title('SST mensuelle moyenne en degC') 
     244else 
     245  colorbar('East') %colorbar par defaut 
     246end; 
    151247 
    152248% sauvegarde de la figure 
    153 printer='eps'; 
     249printer='ps'; 
    154250print_printer=['-d', printer]; 
    155251fullfilename=['./',mfilename,'_',num2str(ifigure),'.',printer]; 
    … …  
    159255clear fullfilename; 
    160256 
     257 
    161258% Variabilité interannuelle 
     259%------------------------- 
     260%niveaux de couleurs 
     261lev=[0:1:30];%contours a tracer, en unite du champ a dessiner, ici 
     262             %la SST en degC 
     263clip_lev=[20 30];%definit les limites des couleurs: bleu fonce au 
     264                 %rouge fonce 
     265 
    162266% Cartes individuelles 
    163267 
    … …  
    175279  end 
    176280  z=squeeze(tab_an(:,:,iplot)); 
    177   contourf(x,y,z,[0:1:30]) 
    178   caxis([20 30]) 
     281  contourf(x,y,z,lev) 
     282  caxis(clip_lev) 
    179283  hold on 
    180   ylim([-30 30]) 
    181   xlim([-100 20]) 
     284  whorldmap 
     285 
     286  % limites geographies et labels de longitude et de latitude 
     287  ylim([geo(3) geo(4)]) 
     288  xlim([geo(1) geo(2)]) 
     289  set(gca, 'Xtick', xti, 'Ytick', yti) 
     290  if (run_octave == 0) 
     291   [longstr, latstr] = contlab(xti, yti); 
     292   set(gca, 'Xticklabel', longstr, 'Yticklabel', latstr) 
     293  end; 
     294  if iplot>=10 
     295    xlabel('longitude') 
     296  end; 
     297  if iplot==1|iplot==4|iplot==7|iplot==10 
     298    ylabel('latitude') 
     299  end; 
     300 
     301  %titre: numéro de l'année 
    182302  title(['annee ' num2str(iplot+1948-1)]) 
    183303end;  
    184 colorbar('East') 
     304 
     305%barre de couleur 
     306if (run_octave == 0) 
     307  pos=[.2 .03 .6 .013]; 
     308  hpal=jet(100); %carte des couleurs: du bleu au rouge par defaut. 
     309  colorbartype(pos,lev,1,clip_lev,hpal,0); 
     310  title('SST mensuelle moyenne en degC') 
     311else 
     312  colorbar('East') %colorbar par defaut 
     313end; 
    185314 
    186315% sauvegarde de la figure 
    187 printer='eps'; 
     316printer='ps'; 
    188317print_printer=['-d', printer]; 
    189318fullfilename=['./',mfilename,'_',num2str(ifigure),'.',printer]; 
    … …  
    193322clear fullfilename; 
    194323 
    195 % ecart type 
     324 
     325% écart type 
     326%----------- 
     327 
     328 
     329%trace 
    196330ifigure=ifigure+1; 
    197331figure(ifigure); 
    198332clf(); 
    199333initfig;  
     334 
     335%écart-type mensuel moyen 
     336 
     337%niveaux de couleurs 
     338lev=[0:.4:5]; 
     339clip_lev=[0 3.6]; 
     340 
    200341subplot(2,2,1) 
    201342if (run_octave == 0) 
    … …  
    205346 [x,y]=meshgrid(lon,lat); 
    206347end 
    207 z=squeeze(std(tab_mens,0,3)); 
    208 contourf(x,y,z,[0:.4:5]) 
    209 caxis([0 3.5]) 
     348z=std_mens; 
     349contourf(x,y,z,lev) 
     350caxis([clip_lev]) 
    210351hold on 
    211 ylim([-30 30]) 
    212 xlim([-100 20]) 
     352whorldmap 
     353  % limites geographies et labels de longitude et de latitude 
     354  ylim([geo(3) geo(4)]) 
     355  xlim([geo(1) geo(2)]) 
     356  set(gca, 'Xtick', xti, 'Ytick', yti) 
     357  if (run_octave == 0) 
     358   [longstr, latstr] = contlab(xti, yti); 
     359   set(gca, 'Xticklabel', longstr, 'Yticklabel', latstr) 
     360  end; 
     361  xlabel('longitude') 
     362  ylabel('latitude') 
     363 
    213364title('ecart-type mensuel moyen') 
    214 colorbar('SouthOutside') 
     365%barre de couleur 
     366if (run_octave == 0) 
     367  pos=[.1 .5 .35 .013]; 
     368  hpal=jet(100); %carte des couleurs: du bleu au rouge par defaut. 
     369  colorbartype(pos,lev,1,clip_lev,hpal,0); 
     370  title('degC') 
     371else 
     372  colorbar('East') %colorbar par defaut 
     373end; 
     374 
     375 
     376 
     377%écart-type interannuel  moyen 
     378 
     379%niveaux de couleurs 
     380lev=[0:.05:.7];%la variabilite interannuelle est moins forte que la 
     381               %variabilite saisonniere moyenne.  
     382clip_lev=[0.2 .55]; 
    215383 
    216384subplot(2,2,2) 
    217 z=squeeze(std(tab_an,0,3)); 
    218 contourf(x,y,z,[0:.05:.7]) 
    219 caxis([0.2 .55]) 
     385z=std_an; 
     386contourf(x,y,z,lev) 
     387caxis(clip_lev) 
    220388hold on 
    221 ylim([-30 30]) 
    222 xlim([-100 20]) 
    223 title('ecart-type interannuel') 
    224 colorbar('SouthOutside') 
     389whorldmap 
     390  % limites geographies et labels de longitude et de latitude 
     391  ylim([geo(3) geo(4)]) 
     392  xlim([geo(1) geo(2)]) 
     393  set(gca, 'Xtick', xti, 'Ytick', yti) 
     394  if (run_octave == 0) 
     395   [longstr, latstr] = contlab(xti, yti); 
     396   set(gca, 'Xticklabel', longstr, 'Yticklabel', latstr) 
     397  end; 
     398  xlabel('longitude') 
     399  ylabel('latitude') 
     400 
     401title('écart-type interannuel') 
     402%barre de couleur 
     403if (run_octave == 0) 
     404  pos=[.55 .5 .35 .013]; 
     405  hpal=jet(100); %carte des couleurs: du bleu au rouge par defaut. 
     406  colorbartype(pos,lev,1,clip_lev,hpal,0); 
     407  title('degC') 
     408else 
     409  colorbar('East') %colorbar par defaut 
     410end; 
     411 
     412 
    225413 
    226414% sauvegarde de la figure 
    227 printer='eps'; 
     415printer='ps'; 
    228416print_printer=['-d', printer]; 
    229417fullfilename=['./',mfilename,'_',num2str(ifigure),'.',printer]; 

  • /trunk/RESULTATS/eof_NCEP_an.m

    • Property svn:keywords set to Id
    r20 r30  
    88% 
    99% le programme permet egalement de visualiser les resultats. Voir le doc 
    10 % .pdf joint pour l'interpretation de cette visualisation.  
     10% .pdf joint pour l'interpretation de cette visualisation. 
    1111% 
    1212% Date: Avril-Mai 2008 
    … …  
    1818% ====== 
    1919% 
    20 % ``eof_NCEP_an``  
     20% ``eof_NCEP_an`` 
    2121% 
    2222% DESCRIPTION 
    2323% =========== 
    2424% 
    25 % ``eof_NCEP_an`` lit les fichiers  
     25% ``eof_NCEP_an`` lit les fichiers 
    2626% ``../DONNEES/lsmask_tropatl_PacandMedblancs.mat`` 
    2727% et ``../DONNEES/skt.mon.tropatl.mat``. 
    … …  
    2929% ``eof_NCEP_an`` calcule ++blabla++ 
    3030% 
    31 % ``eof_NCEP_an`` sauvegarde dans ``../DONNEES/EOF_SSTan_tropAtl++.mat``  
    32 % ou ``++/DONNEES/EOF_SSTan++.mat`` ++ 
    33 % 
    34 % Ces fichiers sont lus par ``Reg_serietemporelles.m_``. 
     31% ``eof_NCEP_an`` sauvegarde dans ``../DONNEES/eof_SSTan_tropAlt30N-20S.mat``. 
     32% 
     33% Ces fichiers sont lus par ``regandsignif_seriestemporelles.m_``. 
    3534% 
    3635% ``eof_NCEP_an`` affiche trois figures : 
     36% 
    3737%  - ++ 
     38% 
    3839%  - ++ 
     40% 
    3941%  - valeurs propres de la matrice de covariance 
    4042% 
    41 % ``eof_NCEP_an`` sauve ces images dans ``./EOF_SSTan_tropAtl++_[123].eps``. 
     43% ``eof_NCEP_an`` sauve ces images dans ``./eof_SSTan_tropAlt30N-20S_[123].ps``. 
    4244% 
    4345% ++ 
    … …  
    6062% 
    6163% initfig.m_ 
     64% 
    6265% .. _initfig.m : initfig.m.html 
    6366% 
    64 % Reg_seriestemporelles.m_ 
    65 % 
    66 % .. _Reg_seriestemporelles.m : Reg_seriestemporelles.m.html 
     67% regandsignif_seriestemporelles.m_ 
     68% 
     69% .. _regandsignif_seriestemporelles.m : regandsignif_seriestemporelles.m.html 
    6770% 
    6871% whorldmap.m_ 
    … …  
    7780% tester la branche NetCDF 
    7881% 
    79 % sous octave, triangle bizarre dans la dernière barre du graphique valeurs  
     82% sous octave, triangle bizarre dans la dernière barre du graphique valeurs 
    8083% propres. 
    8184% 
    … …  
    8588% $Id$ 
    8689% 
     90% - fplod 2009-08-26T10:26:35Z aedon.locean-ipsl.upmc.fr (Darwin) 
     91% 
     92%   * remplacement de ``eps`` par ``ps`` pour pouvoir imprimer les figures 
     93%     produites par ``octave`` 
     94% 
     95% - jmignot 2009-08-25 
     96% 
     97%   * suppression d'appel à contlab et colorbartype sous octave 
     98% 
     99% - fplod 2009-08-25T10:40:16Z aedon.locean-ipsl.upmc.fr (Darwin) 
     100% 
     101%   * modif commentaires liée au renommage 
     102% 
     103% - jmignot 2009-08-25 
     104% 
     105%   * RESULTATS/Reg_seriestemporelles.m renommé en 
     106%     RESULTATS/regandsignif_seriestemporelles.m pour éviter confusion 
     107%     avec PROGRAMMES/Reg_seriestemporelles.m 
     108%   * chgt de noms de fichiers : ``EOF_`` devient ``eof_`` 
     109 
    87110% - fplod 2009-08-21T15:04:16Z aedon.locean-ipsl.upmc.fr (Darwin) 
    88111% 
    89 %   * as octave save default format is ASCII and no default extension  
    90 %     while matlab save default format is MAT v5 mat-file (little endian) and  
     112%   * as octave save default format is ASCII and no default extension 
     113%     while matlab save default format is MAT v5 mat-file (little endian) and 
    91114%     default extension is ``.mat``, write more precise save instruction 
    92115%   * add header 
    … …  
    109132N=4; %nb d'EOF a tracer 
    110133 
    111 nomfic='EOF_SSTan_tropAtl' 
    112 % mettre 'EOF_SSTan_tropAtl' quand on travaille en moyennes annuelles. Ce 
     134nomfic='eof_SSTan_tropAtl' 
     135% mettre 'eof_SSTan_tropAtl' quand on travaille en moyennes annuelles. Ce 
    113136% fichier sera rappele dans les programme de regression lineaire 
    114137%% 
    … …  
    119142 
    120143nannees=60; 
    121     
     144 
    122145repertoire=[ '..' filesep 'DONNEES' filesep ] % a adapter a chaque cas: chemin d acces aux donnees 
    123146 
    124147% Mask 
    125148%----- 
    126 file='lsmask_tropatl_PacandMedblancs' 
     149file='lsmask_tropatl_PacandMedblancs.mat' 
    127150   load([repertoire file]); 
    128   ny=size(mask,1); nx=size(mask,2); % dimensions spatiales des tableaux  
     151  ny=size(mask,1); nx=size(mask,2); % dimensions spatiales des tableaux 
    129152 
    130153% Data 
    … …  
    135158  load([repertoire2 file]); % chargement des donnees 
    136159 
    137     
     160 
    138161  %SI NETCDF 
    139162  %repertoire2=[repertoire 'NETCDF' filesep]% a adapter a chaque cas: chemin d acces aux donnees 
    … …  
    141164  % f=netcdf([repertoire file],'nowrite'); 
    142165  % tab=f{'skt'}(:,:,:); 
    143    
     166 
    144167 
    145168   tab=permute(tab,[2 3 1]);% permutation des dimensions pour organiser le tableau sous la forme [lat lon temps] 
    … …  
    168191%---------------- 
    169192%Si on veut éventuellement attribuer un poids différent aux différents 
    170 %points spatiaux, on le fait ici.  
     193%points spatiaux, on le fait ici. 
    171194 
    172195% Préparation du zoom géographique: sur le masque, j'identifie les points 
    … …  
    176199limite_latS=-20; %limite du domaine d'etuade en latitude 
    177200limite_latN=30; %limite du domaine d'etuade en latitude 
    178 lat2D=repmat(lat,[1 size(mask,2)]);  
     201lat2D=repmat(lat,[1 size(mask,2)]); 
    179202mask(lat2D>limite_latN&mask==0)=4; 
    180203mask(lat2D<limite_latS&mask==0)=4; 
    … …  
    208231% matrice de covariance: 
    209232  C = Z * Z'; 
    210 %  eval(['save ' repertoire 'MatCov_SST_tropAtl' num2str(limite_latN) 'N' num2str(limite_latS) 'S.asc C -ascii'])  
    211    
    212 % diagonalisation   
    213   [E,Lambda] = eig(C);  
    214    % Lambda est la matrice diagonale des valeurs propres et E la matrice dont les colonnes sont les vecteurs propres correspondants.   
     233%  eval(['save ' repertoire 'MatCov_SST_tropAtl' num2str(limite_latN) 'N' num2str(limite_latS) 'S.asc C -ascii']) 
     234 
     235% diagonalisation 
     236  [E,Lambda] = eig(C); 
     237   % Lambda est la matrice diagonale des valeurs propres et E la matrice dont les colonnes sont les vecteurs propres correspondants. 
    215238   % Soit Z*E = E*Lambda 
    216    % Les vecteurs propres sont appelés les EOF. Chaque EOF peut être considéré comme une carte de l'espace physique.     
    217  
    218     
    219  %%  
     239   % Les vecteurs propres sont appelés les EOF. Chaque EOF peut être considéré comme une carte de l'espace physique. 
     240 
     241 
     242 %% 
    220243 %%%% calcul de la matrice donnant les pourcentages de variance 
    221244Lambda = diag(Lambda); 
    222245pctg_var = round(Lambda / sum(Lambda) * 1000) / 10; 
    223246 
    224   
     247 
    225248 %On definit les composantes principales de Q comme la decomposition du signal decrit par Q sur les EOF E. 
    226249 
    227250   P = Z' * E; 
    228      
    229 %%    
     251 
     252%% 
    230253%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
    231254%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
    … …  
    271294 
    272295 %%               TRACE 
     296close all 
     297%parametres geographiques des graphes 
     298  geo=[-100 20 -30 30]; %[lonmin lonmax latmin latmax]; 
     299  yti=[-30:10:30]; xti=[-100:30:20]; %labels de latitude et de longitude 
     300%niveaux de couleurs 
     301  lev=[-2:.02:2]; 
     302  clip_lev=[-0.06 0.06]; 
     303 
    273304figure(1); 
    274305initfig; 
    275 figure(2);  
     306figure(2); 
    276307initfig; 
    277308 
    278309for i=1:4 
    279      
     310 
    280311   %subplot(4,2,2*i-1) 
    281312   ifigure=1; 
    … …  
    289320   end 
    290321   z=squeeze(TABEOF(:,:,i)); 
    291    contourf(x,y,z); 
    292    caxis([-0.06 0.06]) 
     322   contourf(x,y,z,lev); 
     323   caxis(clip_lev); 
    293324   hold on 
    294    xlim([-100 30]); xlabel('longitude') 
    295    ylim([-40 40]); ylabel('latitude') 
    296    whorldmap([-100, 30, -40, 40]) 
    297    colorbar 
    298    title([num2str(p_var(i)) '%  -  (unites: K)']) 
    299    set(gca,'xtick',[-100;-70;-40;-10;20],'xticklabel',['100W';' 70W';' 40W';' 10W';' 20E']) 
    300    set(gca,'ytick',[-30;-15;0;15;30],'yticklabel',['30S';'15S';' 0 ';'15N';'30N']) 
     325   whorldmap 
     326   % limites geographies et labels de longitude et de latitude 
     327  ylim([geo(3) geo(4)]) 
     328  xlim([geo(1) geo(2)]) 
     329  set(gca, 'Xtick', xti, 'Ytick', yti) 
     330   if (run_octave == 0) 
     331    [longstr, latstr] = contlab(xti, yti); 
     332    set(gca, 'Xticklabel', longstr, 'Yticklabel', latstr) 
     333   end; 
     334  xlabel('longitude') 
     335  ylabel('latitude') 
     336 
     337   %titre de la figure 
     338   title([num2str(p_var(i)) '%']) 
     339 
     340  %barre de couleur 
     341  if (run_octave == 0) 
     342    hpal=jet(100); %carte des couleurs: du bleu au rouge par defaut. 
     343    pos1=get(gca,'Position'); 
     344    pos=[pos1(1) pos1(2)-.075 pos1(3) .013]; 
     345    colorbartype(pos,lev,1,clip_lev,hpal,0); 
     346    title('degC') 
     347  else 
     348    colorbar('East') 
     349  end; 
    301350 
    302351   % sauvegarde de la figure 
    303    printer='eps'; 
     352   printer='ps'; 
    304353   print_printer=['-d', printer]; 
    305354   fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer]; 
    … …  
    307356   clear fullfilename; 
    308357   clear printer; 
    309     
     358 
    310359   %subplot(4,2,2*i) 
    311360   ifigure=2; 
    … …  
    320369 
    321370   % sauvegarde de la figure 
    322    printer='eps'; 
     371   printer='ps'; 
    323372   print_printer=['-d', printer]; 
    324373   fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer]; 
    … …  
    332381% Graphe pourcentages de variance expliqué 
    333382ifigure=3; 
    334 figure(ifigure);  
     383figure(ifigure); 
    335384bar(flipud(pctg_var)) 
    336385xlim([0 10]) 
    … …  
    340389 
    341390% sauvegarde de la figure 
    342 printer='eps'; 
     391printer='ps'; 
    343392print_printer=['-d', printer]; 
    344393fullfilename=['.' filesep nomfic num2str(limite_latN) 'N' num2str(limite_latS) 'S_' , num2str(ifigure) '.' printer]; 

  • /trunk/adm/makefile

    r20 r30  
    7676$(DIRSRC)/PROGRAMMES/eof_NCEP_sais.m \ 
    7777$(DIRSRC)/PROGRAMMES/whorldmap.m \ 
    78 $(DIRSRC)/PROGRAMMES/Reg_seriestemporelles.m \ 
     78$(DIRSRC)/PROGRAMMES/reg_seriestemporelles.m \ 
    7979$(DIRSRC)/RESULTATS/carteetstd_SST.m \ 
    8080$(DIRSRC)/RESULTATS/eof_NCEP_an.m \ 
    81 $(DIRSRC)/RESULTATS/Reg_seriestemporelles.m 
     81$(DIRSRC)/RESULTATS/regandsignif_seriestemporelles.m 
    8282 
    8383LIST_SRCMATLAB_RST = \ 
Note: See TracChangeset for help on using the changeset viewer.