Changeset 21 for trunk/src/mode_sahelien


Ignore:
Timestamp:
01/06/09 16:22:54 (15 years ago)
Author:
pinsard
Message:

small homogenization

Location:
trunk/src/mode_sahelien
Files:
23 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/mode_sahelien/composite_olr_phase2.m

    r18 r21  
    149149clear olr; 
    150150olr2=reshape(olr1,3416,21,144); 
    151 olr3(:,:)=mean(olr2); clear olr2; 
     151olr3(:,:)=mean(olr2); 
     152clear olr2; 
    152153 
    153154load phase2.txt; 
    154155phase2=reshape(phase2,122,28,3); 
    155  olr_filtre_compomax=zeros(90,21,144); 
     156olr_filtre_compomax=zeros(90,21,144); 
    156157compteurmax=0; 
    157158 

  • trunk/src/mode_sahelien/composite_olr_phase3.m

    r18 r21  
    154154load phase3.txt; 
    155155phase3=reshape(phase3,122,28,3); 
    156  olr_filtre_compomax=zeros(90,21,144); 
     156olr_filtre_compomax=zeros(90,21,144); 
    157157compteurmax=0; 
    158158 

  • trunk/src/mode_sahelien/composite_olr_phase4.m

    r18 r21  
    149149clear olr; 
    150150olr2=reshape(olr1,3416,21,144); 
    151 olr3(:,:)=mean(olr2); clear olr2; 
     151olr3(:,:)=mean(olr2); 
     152clear olr2; 
    152153 
    153154load phase4.txt; 
    154155phase4=reshape(phase4,122,28,3); 
    155  olr_filtre_compomax=zeros(90,21,144); 
     156olr_filtre_compomax=zeros(90,21,144); 
    156157compteurmax=0; 
    157158 

  • trunk/src/mode_sahelien/composite_olr_phase5.m

    r18 r21  
    154154load phase5.txt; 
    155155phase5=reshape(phase5,122,28,3); 
    156  olr_filtre_compomax=zeros(90,21,144); 
     156olr_filtre_compomax=zeros(90,21,144); 
    157157compteurmax=0; 
    158158 

  • trunk/src/mode_sahelien/composite_olr_phase6.m

    r18 r21  
    154154load phase6.txt; 
    155155phase6=reshape(phase6,122,28,3); 
    156  olr_filtre_compomax=zeros(90,21,144); 
     156olr_filtre_compomax=zeros(90,21,144); 
    157157compteurmax=0; 
    158158 

  • trunk/src/mode_sahelien/composite_olr_phase7.m

    r18 r21  
    149149clear olr; 
    150150olr2=reshape(olr1,3416,21,144); 
    151 olr3(:,:)=mean(olr2); clear olr2; 
     151olr3(:,:)=mean(olr2); 
     152clear olr2; 
    152153 
    153154load phase7.txt; 
    154155phase7=reshape(phase7,122,28,3); 
    155  olr_filtre_compomax=zeros(90,21,144); 
     156olr_filtre_compomax=zeros(90,21,144); 
    156157compteurmax=0; 
    157158 

  • trunk/src/mode_sahelien/composite_olr_phase8.m

    r18 r21  
    154154load phase8.txt; 
    155155phase8=reshape(phase8,122,28,3); 
    156  olr_filtre_compomax=zeros(90,21,144); 
     156olr_filtre_compomax=zeros(90,21,144); 
    157157compteurmax=0; 
    158158 

  • trunk/src/mode_sahelien/time_serie_ER_rossby_1030.m

    r18 r21  
    1313    for b=1:28; 
    1414       year(a,b)=1979+(b-1); 
    15     end 
    16 end 
     15    end; 
     16end; 
    1717year=reshape(year,3416,1); 
    1818 
    … …  
    2121    for b=1:28; 
    2222       day(a,b)=31+28+31+30+31+a; 
    23     end 
    24 end 
     23    end; 
     24end; 
    2525day=reshape(day,3416,1); 
    2626 
    … …  
    5858set (gca,'YTick',[-100:20:100],'YTickLabel',[-100:20:100],'fontname','Arial','fontsize',6); 
    5959set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 
    60  

  • trunk/src/mode_sahelien/time_serie_ER_rossby_1030_complet.m

    r18 r21  
    2828      figure(3); 
    2929      subplot(4,3,a-22); 
    30    end 
     30   end; 
    3131   clear koro1 
    3232   debut=(a-1)*122+1; 
    … …  
    4949   elseif a==24; 
    5050       text(0,183,'time series EOF1 and EOF2 OLR ER 10-30d'); 
    51    end 
     51   end; 
    5252   title(a+1978,'FontSize',8); 
    5353   set (gca,'tickDir','out'); 
    … …  
    5656set (gca,'YTick',[-100:20:100],'YTickLabel',[-100:20:100],'fontname','Arial','fontsize',6); 
    5757set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 
    58 end 
     58end; 
    5959 
    6060figure(1); 

  • trunk/src/mode_sahelien/time_serie_ER_rossby_30100.m

    r18 r21  
    1313    for b=1:28; 
    1414       year(a,b)=1979+(b-1); 
    15     end 
    16 end 
     15    end; 
     16end; 
    1717year=reshape(year,3416,1); 
    1818 
    … …  
    2121    for b=1:28; 
    2222       day(a,b)=31+28+31+30+31+a; 
    23     end 
    24 end 
     23    end; 
     24end; 
    2525day=reshape(day,3416,1); 
    2626 

  • trunk/src/mode_sahelien/time_serie_ER_rossby_30100_complet.m

    r18 r21  
    2828      figure(3); 
    2929      subplot(4,3,a-22); 
    30    end 
     30   end; 
    3131   clear koro1 
    3232   debut=(a-1)*122+1; 
    … …  
    4949   elseif a==24; 
    5050       text(0,183,'time series EOF1 and EOF2 OLR ER 30-100d'); 
    51    end 
     51   end; 
    5252   title(a+1978,'FontSize',8); 
    5353   set (gca,'tickDir','out'); 
    … …  
    5656set (gca,'YTick',[-100:20:100],'YTickLabel',[-100:20:100],'fontname','Arial','fontsize',6); 
    5757set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 
    58 end 
     58end; 
    5959 
    6060figure(1); 

  • trunk/src/mode_sahelien/time_serie_histog_phase1.m

    r18 r21  
    8888    for b=1:4 
    8989       kara_diff(a,b)=kara_max(a,b)-kara_moy(b); 
    90     end 
    91 end 
     90    end; 
     91end; 
    9292 
    9393for a=1:4; 
    … …  
    105105    else 
    106106       bar(jour,diff,'r'); 
    107     end 
     107    end; 
    108108    set (gca,'XTick',[-10:2:10],'XTickLabel',[-10:2:10],'fontname','Arial','fontsize',6); 
    109109    set (gca,'YTick',[-20:10:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); 
    … …  
    120120       title('Blocking','Fontsize',8,'Fontname','Arial'); 
    121121    end; 
    122 end 
     122end; 

  • trunk/src/mode_sahelien/time_serie_histog_phase2.m

    r18 r21  
    8888    for b=1:4 
    8989       kara_diff(a,b)=kara_max(a,b)-kara_moy(b); 
    90     end 
    91 end 
     90    end; 
     91end; 
    9292 
    9393for a=1:4; 
    … …  
    105105    else 
    106106       bar(jour,diff,'r'); 
    107     end 
     107    end; 
    108108    set (gca,'XTick',[-10:2:10],'XTickLabel',[-10:2:10],'fontname','Arial','fontsize',6); 
    109109    set (gca,'YTick',[-20:10:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); 
    … …  
    120120       title('Blocking','Fontsize',8,'Fontname','Arial'); 
    121121    end; 
    122 end 
     122end; 

  • trunk/src/mode_sahelien/time_serie_histog_phase3.m

    r17 r21  
    2323% boucle sur les annees 
    2424for a=1:28; 
    25    % 1ere condition (laisser une marge de 10) 
    2625   for b=1:122; 
     26      % 1ere condition (laisser une marge de 10) 
    2727      if (b>11 && b<112); 
    2828%        condition 1 on est a un maximum deltas apres sont negetifs et 
    … …  
    8888    for b=1:4 
    8989       kara_diff(a,b)=kara_max(a,b)-kara_moy(b); 
    90     end 
    91 end 
     90    end; 
     91end; 
    9292 
    9393for a=1:4; 
    … …  
    105105    else 
    106106       bar(jour,diff,'r'); 
    107     end 
     107    end; 
    108108    set (gca,'XTick',[-10:2:10],'XTickLabel',[-10:2:10],'fontname','Arial','fontsize',6); 
    109109    set (gca,'YTick',[-20:10:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); 

  • trunk/src/mode_sahelien/time_serie_histog_phase4.m

    r18 r21  
    8888    for b=1:4 
    8989       kara_diff(a,b)=kara_max(a,b)-kara_moy(b); 
    90     end 
    91 end 
     90    end; 
     91end; 
    9292 
    9393for a=1:4; 
    … …  
    105105    else 
    106106       bar(jour,diff,'r'); 
    107     end 
     107    end; 
    108108    set (gca,'XTick',[-10:2:10],'XTickLabel',[-10:2:10],'fontname','Arial','fontsize',6); 
    109109    set (gca,'YTick',[-20:10:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); 
    … …  
    120120       title('Blocking','Fontsize',8,'Fontname','Arial'); 
    121121    end; 
    122 end 
     122end; 

  • trunk/src/mode_sahelien/time_serie_histog_phase5.m

    r18 r21  
    8888    for b=1:4 
    8989       kara_diff(a,b)=kara_max(a,b)-kara_moy(b); 
    90     end 
    91 end 
     90    end; 
     91end; 
    9292 
    9393for a=1:4; 
    … …  
    105105    else 
    106106       bar(jour,diff,'r'); 
    107     end 
     107    end; 
    108108    set (gca,'XTick',[-10:2:10],'XTickLabel',[-10:2:10],'fontname','Arial','fontsize',6); 
    109109    set (gca,'YTick',[-20:10:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); 
    … …  
    120120       title('Blocking','Fontsize',8,'Fontname','Arial'); 
    121121    end; 
    122 end 
     122end; 

  • trunk/src/mode_sahelien/time_serie_histog_phase6.m

    r18 r21  
    8888    for b=1:4 
    8989       kara_diff(a,b)=kara_max(a,b)-kara_moy(b); 
    90     end 
    91 end 
     90    end; 
     91end; 
    9292 
    9393for a=1:4; 
    … …  
    105105    else 
    106106       bar(jour,diff,'r'); 
    107     end 
     107    end; 
    108108    set (gca,'XTick',[-10:2:10],'XTickLabel',[-10:2:10],'fontname','Arial','fontsize',6); 
    109109    set (gca,'YTick',[-20:10:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); 
    … …  
    120120       title('Blocking','Fontsize',8,'Fontname','Arial'); 
    121121    end; 
    122 end 
     122end; 

  • trunk/src/mode_sahelien/time_serie_histog_phase7.m

    r17 r21  
    8888    for b=1:4 
    8989       kara_diff(a,b)=kara_max(a,b)-kara_moy(b); 
    90     end 
    91 end 
     90    end; 
     91end; 
    9292 
    9393for a=1:4; 
    … …  
    105105    else 
    106106       bar(jour,diff,'r'); 
    107     end 
     107    end; 
    108108    set (gca,'XTick',[-10:2:10],'XTickLabel',[-10:2:10],'fontname','Arial','fontsize',6); 
    109109    set (gca,'YTick',[-20:10:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); 
    … …  
    120120       title('Blocking','Fontsize',8,'Fontname','Arial'); 
    121121    end; 
    122 end 
     122end; 

  • trunk/src/mode_sahelien/time_serie_histog_phase8.m

    r18 r21  
    8888    for b=1:4 
    8989       kara_diff(a,b)=kara_max(a,b)-kara_moy(b); 
    90     end 
    91 end 
     90    end; 
     91end; 
    9292 
    9393for a=1:4; 
    … …  
    105105    else 
    106106       bar(jour,diff,'r'); 
    107     end 
     107    end; 
    108108    set (gca,'XTick',[-10:2:10],'XTickLabel',[-10:2:10],'fontname','Arial','fontsize',6); 
    109109    set (gca,'YTick',[-20:10:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); 
    … …  
    117117    elseif a==3; 
    118118       title('NAO-','Fontsize',8,'Fontname','Arial'); 
    119        else 
     119    else 
    120120       title('Blocking','Fontsize',8,'Fontname','Arial'); 
    121121    end; 
    122 end 
     122end; 

  • trunk/src/mode_sahelien/time_serie_jjas.m

    r18 r21  
    2020   time_serie_olr2(a)=time_serie_olr2(a)/dt_olr2; 
    2121   time_serie_olr3(a)=time_serie_olr3(a)/dt_olr3; 
    22 end 
     22end; 
    2323 
    2424save eof1_1030.txt time_serie_olr1 -ascii; 
    … …  
    3737      figure(3); 
    3838      subplot(4,3,a-22); 
    39    end 
    40    clear koro1 
     39   end; 
     40   clear koro1; 
    4141   debut=(a-1)*122+1; 
    4242   fin=a*122; 
    … …  
    4949   y=zeros(122,1); 
    5050   plot(x,y,'k:'); 
    51    title(a+1978,'FontSize',8) 
     51   title(a+1978,'FontSize',8); 
    5252   if a==2; 
    5353       text(0,5.5,'serie temporelle VP1'); 
    54     elseif a==13; 
     54   elseif a==13; 
    5555       text(0,5.5,'serie temporelle VP1'); 
    5656   elseif a==24; 
    5757       text(0,5.5,'serie temporelle VP1'); 
    58    end 
     58   end; 
    5959   set (gca,'tickDir','out'); 
    6060   axis([1 122 -3 3]); 
    6161 
    62 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 
    63 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 
     62   set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 
     63   set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 
    6464end; 

  • trunk/src/mode_sahelien/time_serie_jjas_10_25.m

    r18 r21  
    2020   time_serie_olr2(a)=time_serie_olr2(a)/dt_olr2; 
    2121   time_serie_olr3(a)=time_serie_olr3(a)/dt_olr3; 
    22 end 
     22end; 
    2323 
    2424save eof1_olr_rossby_jjas_1025.txt time_serie_olr1 -ascii; 
    … …  
    3737      figure(3); 
    3838      subplot(4,3,a-22); 
    39    end 
    40    clear koro1 
     39   end; 
     40   clear koro1; 
    4141   debut=(a-1)*122+1; 
    4242   fin=a*122; 
    … …  
    5151   if a==2; 
    5252       text(0,5.5,'serie temporelle VP1 filtrage rossby 10-25j'); 
    53     elseif a==13; 
     53   elseif a==13; 
    5454       text(0,5.5,'serie temporelle VP1 filtrage rossby 10-25j'); 
    5555   elseif a==24; 
    … …  
    6060 
    6161   set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 
    62 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 
     62   set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 
    6363end; 
    6464 

  • trunk/src/mode_sahelien/time_serie_jjas_25_90.m

    r18 r21  
    2020   time_serie_olr2(a)=time_serie_olr2(a)/dt_olr2; 
    2121   time_serie_olr3(a)=time_serie_olr3(a)/dt_olr3; 
    22 end 
     22end; 
    2323 
    2424save eof1_olr_rossby_jjas_2590.txt time_serie_olr1 -ascii; 
    … …  
    3737      figure(3); 
    3838      subplot(4,3,a-22); 
    39    end 
    40    clear koro1 
     39   end; 
     40   clear koro1; 
    4141   debut=(a-1)*122+1; 
    4242   fin=a*122; 
    … …  
    5252   if a==2; 
    5353       text(0,5.5,'serie temporelle VP1 filtrage rossby 25-90j'); 
    54     elseif a==13; 
     54   elseif a==13; 
    5555       text(0,5.5,'serie temporelle VP1 filtrage rossby 25-90j'); 
    5656   elseif a==24; 
    … …  
    6161 
    6262   set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 
    63 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 
     63   set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 
    6464end; 
    6565 

  • trunk/src/mode_sahelien/time_serie_jjas_30100.m

    r18 r21  
    2020   time_serie_olr2(a)=time_serie_olr2(a)/dt_olr2; 
    2121   time_serie_olr3(a)=time_serie_olr3(a)/dt_olr3; 
    22 end 
     22end; 
    2323 
    2424save eof1_30100.txt time_serie_olr1 -ascii; 
    … …  
    3737      figure(3); 
    3838      subplot(4,3,a-22); 
    39    end 
    40    clear koro1 
     39   end; 
     40   clear koro1; 
    4141   debut=(a-1)*122+1; 
    4242   fin=a*122; 
    … …  
    4949   y=zeros(122,1); 
    5050   plot(x,y,'k:'); 
    51    title(a+1978,'FontSize',8) 
     51   title(a+1978,'FontSize',8); 
    5252   if a==2; 
    5353       text(0,5.5,'serie temporelle VP1 filtrage 30100'); 
    54     elseif a==13; 
     54   elseif a==13; 
    5555       text(0,5.5,'serie temporelle VP1 filtrage 30100'); 
    5656   elseif a==24; 
    5757       text(0,5.5,'serie temporelle VP1 filtrage 30100'); 
    58    end 
     58   end; 
    5959   set (gca,'tickDir','out'); 
    6060   axis([1 122 -3 3]); 
    6161 
    62 set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 
    63 set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 
     62   set (gca,'YTick',[-3:1:3],'YTickLabel',[-3:1:3],'fontname','Arial','fontsize',6); 
     63   set (gca,'XTick',[0:20:120],'XTickLabel',[0:20:120],'fontname','Arial','fontsize',6); 
    6464end; 
    6565 
Note: See TracChangeset for help on using the changeset viewer.