Ignore:
Timestamp:
01/06/09 10:53:18 (15 years ago)
Author:
pinsard
Message:

remove trailing blanks, split lines with multiple statements

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/mode_sahelien/composite_wind925_EOF12_ER30100.m

    r15 r16  
     1% 
     2 
    13% initialisation 
    24close all; 
    … …  
    810uwnd925=zeros(122,49,144,28); 
    911 
    10 siz=31+28+31+30+31+1; sizbis=31+28+31+30+31+30+31+31+30; 
     12siz=31+28+31+30+31+1; 
     13sizbis=31+28+31+30+31+30+31+31+30; 
    1114ncload('uwnd925.1979.nc'); uwnd925(:,:,:,1)=uwnd(siz:sizbis,1,:,:); clear uwnd; 
    1215ncload('uwnd925.1980.nc'); uwnd925(:,:,:,2)=uwnd(siz+1:sizbis+1,1,:,:); clear uwnd; 
    … …  
    4649 
    4750%%%%%%%%%%%%%%% FICHIER OLR 
    48 % ouverture / lecture du fichier  
    49 load eof12_rossby_30100.txt; ind_olr=reshape(eof12_rossby_30100,122,28); 
     51% ouverture / lecture du fichier 
     52load eof12_rossby_30100.txt; 
     53ind_olr=reshape(eof12_rossby_30100,122,28); 
    5054clear eof12_rossby_30100; 
    5155 
    … …  
    5458% on calcule l ecart type sur la serie entiere de l OLR 
    5559ind_olr1=reshape(ind_olr, k,1); 
    56 ind_olr2=std(ind_olr1); clear ind_olr1; 
     60ind_olr2=std(ind_olr1); 
     61clear ind_olr1; 
    5762 
    5863 
    5964% intialisation des tableaux 
    60 uwnd925_filtre_compomax=NaN*ones(60,31,49,144);  
     65uwnd925_filtre_compomax=NaN*ones(60,31,49,144); 
    6166uwnd925_filtre_compomin=NaN*ones(60,31,49,144); 
    6267 
    63 compteurmax=0;  
     68compteurmax=0; 
    6469compteurmin=0; 
    6570 
    66 for a=1:28; % boucle sur les annees 
    67    clear delta;  
     71% boucle sur les annees 
     72for a=1:28; 
     73   clear delta; 
    6874   delta=zeros(122,1); 
    69    for b=1:122; % boucle sur les jours : on prepare avant    
     75   % boucle sur les jours : on prepare avant 
     76   for b=1:122; 
    7077        if b>1; 
    7178            delta(b)=ind_olr(b,a)-ind_olr(b-1,a); 
    7279        end; 
    7380   end; 
    74    for b=1:122; % boucle sur les jours 
    75       if (b>20 && b<112); % 1ere condition (laisser une marge de 10)  
    76 % cas 1 : max deltas apres negetifs et deltas avant postifs    
    77          if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2);  
     81   % boucle sur les jours 
     82   for b=1:122; 
     83      % 1ere condition (laisser une marge de 10) 
     84      if (b>20 && b<112); 
     85% cas 1 : max deltas apres negetifs et deltas avant postifs 
     86         if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2); 
    7887            compteurmax=compteurmax+1; 
    7988            uwnd925_filtre_compomax(compteurmax,1,:,:)=uwnd(b-20,a,:,:); 
    … …  
    110119            uwnd925_filtre_compomax(compteurmax,30,:,:)=uwnd(b+9,a,:,:); 
    111120            uwnd925_filtre_compomax(compteurmax,31,:,:)=uwnd(b+10,a,:,:); 
    112 % cas 2 : min deltas apres positifs et deltas avant negatifs    
     121% cas 2 : min deltas apres positifs et deltas avant negatifs 
    113122         elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && ind_olr(b,a)<(-ind_olr2)); 
    114123            compteurmin=compteurmin+1; 
    … …  
    149158         end; 
    150159      end; 
    151    end;          
     160   end; 
    152161end; 
    153162 
    … …  
    155164% on fait la somme des valeurs de uwnd925_filtreentiel sur le nbre d occurences 
    156165% on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
    157 uwnd925_filtre_compomin1(:,:,:)=nanmean(uwnd925_filtre_compomin); clear uwnd925_filtre_compomin; 
    158 uwnd925_filtre_compomax1(:,:,:)=nanmean(uwnd925_filtre_compomax); clear uwnd925_filtre_compomax; 
    159  
    160  
    161  
    162  
    163 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
     166uwnd925_filtre_compomin1(:,:,:)=nanmean(uwnd925_filtre_compomin); 
     167clear uwnd925_filtre_compomin; 
     168uwnd925_filtre_compomax1(:,:,:)=nanmean(uwnd925_filtre_compomax); 
     169clear uwnd925_filtre_compomax; 
    164170 
    165171% etape 1 :on selectionne les mois d ete (attention aux annees bixestiles!) 
    … …  
    168174vwnd925=zeros(122,49,144,28); 
    169175 
    170 siz=31+28+31+30+31+1; sizbis=31+28+31+30+31+30+31+31+30; 
     176siz=31+28+31+30+31+1; 
     177sizbis=31+28+31+30+31+30+31+31+30; 
    171178ncload('vwnd925.1979.nc'); vwnd925(:,:,:,1)=vwnd(siz:sizbis,1,:,:); clear vwnd; 
    172179ncload('vwnd925.1980.nc'); vwnd925(:,:,:,2)=vwnd(siz+1:sizbis+1,1,:,:); clear vwnd; 
    … …  
    204211 
    205212%%%%%%%%%%%%%% FICHIER OLR 
    206 % ouverture / lecture du fichier  
    207 load eof12_rossby_30100.txt; ind_olr=reshape(eof12_rossby_30100,122,28); 
     213% ouverture / lecture du fichier 
     214load eof12_rossby_30100.txt; 
     215ind_olr=reshape(eof12_rossby_30100,122,28); 
    208216clear eof12_rossby_30100 
    209217k=122*28; 
    … …  
    211219% on calcule l ecart type sur la serie entiere de l OLR 
    212220ind_olr1=reshape(ind_olr, k,1); 
    213 ind_olr2=std(ind_olr1); clear ind_olr1; 
     221ind_olr2=std(ind_olr1); 
     222clear ind_olr1; 
    214223 
    215224 
    216225% intialisation des tableaux 
    217 vwnd925_filtre_compomax=NaN*ones(60,31,49,144);  
     226vwnd925_filtre_compomax=NaN*ones(60,31,49,144); 
    218227vwnd925_filtre_compomin=NaN*ones(60,31,49,144); 
    219228 
    220 compteurmax=0;  
     229compteurmax=0; 
    221230compteurmin=0; 
    222231 
    223 for a=1:28; % boucle sur les annees 
    224    clear delta;  
     232% boucle sur les annees 
     233for a=1:28; 
     234   clear delta; 
    225235   delta=zeros(122,1); 
    226    for b=1:122; % boucle sur les jours : on prepare avant    
     236   % boucle sur les jours : on prepare avant 
     237   for b=1:122; 
    227238        if b>1; 
    228239            delta(b)=ind_olr(b,a)-ind_olr(b-1,a); 
    229240        end; 
    230241   end; 
    231    for b=1:122; % boucle sur les jours 
    232       if (b>20 && b<112); % 1ere condition (laisser une marge de 10)  
    233 % cas 1 : max deltas apres negetifs et deltas avant postifs    
    234          if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2);  
     242   % boucle sur les jours 
     243   for b=1:122; 
     244      % 1ere condition (laisser une marge de 10) 
     245      if (b>20 && b<112); 
     246% cas 1 : max deltas apres negetifs et deltas avant postifs 
     247         if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2); 
    235248            compteurmax=compteurmax+1; 
    236249            vwnd925_filtre_compomax(compteurmax,1,:,:)=vwnd(b-20,a,:,:); 
    … …  
    267280            vwnd925_filtre_compomax(compteurmax,30,:,:)=vwnd(b+9,a,:,:); 
    268281            vwnd925_filtre_compomax(compteurmax,31,:,:)=vwnd(b+10,a,:,:); 
    269 % cas 2 : min deltas apres positifs et deltas avant negatifs    
     282% cas 2 : min deltas apres positifs et deltas avant negatifs 
    270283         elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && ind_olr(b,a)<(-ind_olr2)); 
    271284            compteurmin=compteurmin+1; 
    … …  
    306319         end; 
    307320      end; 
    308    end;          
     321   end; 
    309322end; 
    310323 
    … …  
    312325% on fait la somme des valeurs de vwnd925_filtreentiel sur le nbre d occurences 
    313326% on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
    314 vwnd925_filtre_compomin1(:,:,:)=nanmean(vwnd925_filtre_compomin); clear vwnd925_filtre_compomin; 
    315 vwnd925_filtre_compomax1(:,:,:)=nanmean(vwnd925_filtre_compomax); clear vwnd925_filtre_compomax; 
    316  
    317 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
     327vwnd925_filtre_compomin1(:,:,:)=nanmean(vwnd925_filtre_compomin); 
     328clear vwnd925_filtre_compomin; 
     329vwnd925_filtre_compomax1(:,:,:)=nanmean(vwnd925_filtre_compomax); 
     330clear vwnd925_filtre_compomax; 
    318331 
    319332figure(1); 
    … …  
    333346             delta_uwnd925_filtre_min(b,c)=uwnd925_filtre_compomin1(f,b,c)-uwnd925_filtre_compomax1(f,b,c); 
    334347             delta_vwnd925_filtre_min(b,c)=vwnd925_filtre_compomin1(f,b,c)-vwnd925_filtre_compomax1(f,b,c); 
    335   
     348 
    336349         end; 
    337350       end; 
    … …  
    344357       delta_vwnd925_filtre(:,1:72)=delta_vwnd925_filtre_min(:,73:144); 
    345358       delta_vwnd925_filtre(:,73:144)=delta_vwnd925_filtre_min(:,1:72); 
    346        xi=[-180:3.75:176.25]; xi=xi'; yi=[-87.1591:3.708895:87.1591];  
     359       xi=[-180:3.75:176.25]; 
     360       xi=xi'; 
     361       yi=[-87.1591:3.708895:87.1591]; 
    347362       clear carteu cartev; 
    348        carteu=interp2(lon-180,lat,delta_uwnd925_filtre,xi,yi);  
     363       carteu=interp2(lon-180,lat,delta_uwnd925_filtre,xi,yi); 
    349364       cartev=interp2(lon-180,lat,delta_vwnd925_filtre,xi,yi); 
    350365       clear coco ; 
    351366       coco(:,:)=pression(5-e,:,:); 
    352       hgt=interp2(lon-180,lat,coco,xi,yi);  
     367       hgt=interp2(lon-180,lat,coco,xi,yi); 
    353368       palette 
    354  
    355        aa=contourf(xi,yi,hgt,[-20:2:20],'LineStyle','none'); hold on; 
     369       aa=contourf(xi,yi,hgt,[-20:2:20],'LineStyle','none'); 
     370       hold on; 
    356371       caxis([-20 20]); 
    357372       aa=quiver(xi,yi,carteu,cartev,'k'); 
    … …  
    360375       else 
    361376set (gca,'XTick',[-60:20:120],'XTickLabel',[' ';' ';' ';' ';' ';' ';' ';' ';' ';' '],'fontsize',6); 
    362        end    
     377       end 
    363378       set (gca,'YTick',[-10:10:30],'YTickLabel',['10S';' 0 ';'10N';'20N';'30N'],'fontname','Arial','fontsize',6); 
    364379       cartemonde1; 
    365380       if e==2; 
    366           co=colorbar;             
     381          co=colorbar; 
    367382          set (co,'ytick',[-20:10:20],'yticklabel',[-20:10:20],'fontname','arial','fontsize',6); 
    368383       end; 
    … …  
    376391       text(-150,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 
    377392       pos_vert=e/6.1; 
    378 %       set (gca,'position',[0.2 pos_vert 0.3 1/6.1])  
    379        set (gca,'position',[0.2 pos_vert 0.4  1/6.1])  
     393%       set (gca,'position',[0.2 pos_vert 0.3 1/6.1]) 
     394       set (gca,'position',[0.2 pos_vert 0.4  1/6.1]) 
    380395end; 
    381396 
    382  
    383 print -depsc2 composite_wind925_EOF12_ER30100_part1.eps;  
    384  
    385  
    386  
    387 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
     397print -depsc2 composite_wind925_EOF12_ER30100_part1.eps; 
    388398 
    389399figure(2); 
    … …  
    403413             delta_uwnd925_filtre_min(b,c)=uwnd925_filtre_compomin1(f,b,c)-uwnd925_filtre_compomax1(f,b,c); 
    404414             delta_vwnd925_filtre_min(b,c)=vwnd925_filtre_compomin1(f,b,c)-vwnd925_filtre_compomax1(f,b,c); 
    405   
     415 
    406416         end; 
    407417       end; 
    … …  
    414424       delta_vwnd925_filtre(:,1:72)=delta_vwnd925_filtre_min(:,73:144); 
    415425       delta_vwnd925_filtre(:,73:144)=delta_vwnd925_filtre_min(:,1:72); 
    416        xi=[-180:3.75:176.25]; xi=xi'; yi=[-87.1591:3.708895:87.1591];  
     426       xi=[-180:3.75:176.25]; 
     427       xi=xi'; 
     428       yi=[-87.1591:3.708895:87.1591]; 
    417429       clear carteu cartev; 
    418        carteu=interp2(lon-180,lat,delta_uwnd925_filtre,xi,yi);  
     430       carteu=interp2(lon-180,lat,delta_uwnd925_filtre,xi,yi); 
    419431       cartev=interp2(lon-180,lat,delta_vwnd925_filtre,xi,yi); 
    420432       clear coco ; 
    421433       coco(:,:)=pression(5-e+4,:,:); 
    422       hgt=interp2(lon-180,lat,coco,xi,yi);  
     434       hgt=interp2(lon-180,lat,coco,xi,yi); 
    423435       palette 
    424        aa=contourf(xi,yi,hgt,[-20:2:20],'LineStyle','none');  
     436       aa=contourf(xi,yi,hgt,[-20:2:20],'LineStyle','none'); 
    425437       caxis([-20 20]); 
    426438       hold on; 
    … …  
    430442       else 
    431443set (gca,'XTick',[-60:20:120],'XTickLabel',[' ';' ';' ';' ';' ';' ';' ';' ';' ';' '],'fontsize',6); 
    432        end   
     444       end 
    433445       set (gca,'YTick',[-10:10:30],'YTickLabel',['10S';' 0 ';'10N';'20N';'30N'],'fontname','Arial','fontsize',6); 
    434446 
    435447       cartemonde1; 
    436448       if e==2; 
    437           co=colorbar;             
     449          co=colorbar; 
    438450          set (co,'ytick',[-20:10:20],'yticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 
    439451       end; 
    … …  
    447459       text(-150,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 
    448460       pos_vert=e/6.1; 
    449        set (gca,'position',[0.2 pos_vert 0.4 1/6.1])  
     461       set (gca,'position',[0.2 pos_vert 0.4 1/6.1]) 
    450462end; 
    451463 
    452464 
    453 print -depsc2 composite_wind925_EOF12_ER30100_part2.eps;  
    454  
    455  
    456 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
     465print -depsc2 composite_wind925_EOF12_ER30100_part2.eps; 
    457466 
    458467figure(3); 
    … …  
    471480             delta_uwnd925_filtre_min(b,c)=uwnd925_filtre_compomin1(f,b,c)-uwnd925_filtre_compomax1(f,b,c); 
    472481             delta_vwnd925_filtre_min(b,c)=vwnd925_filtre_compomin1(f,b,c)-vwnd925_filtre_compomax1(f,b,c); 
    473   
     482 
    474483         end; 
    475484       end; 
    … …  
    482491       delta_vwnd925_filtre(:,1:72)=delta_vwnd925_filtre_min(:,73:144); 
    483492       delta_vwnd925_filtre(:,73:144)=delta_vwnd925_filtre_min(:,1:72); 
    484        xi=[-180:3.75:176.25]; xi=xi'; yi=[-87.1591:3.708895:87.1591];  
     493       xi=[-180:3.75:176.25]; 
     494       xi=xi'; 
     495       yi=[-87.1591:3.708895:87.1591]; 
    485496       clear carteu cartev; 
    486        carteu=interp2(lon-180,lat,delta_uwnd925_filtre,xi,yi);  
     497       carteu=interp2(lon-180,lat,delta_uwnd925_filtre,xi,yi); 
    487498       cartev=interp2(lon-180,lat,delta_vwnd925_filtre,xi,yi); 
    488499       clear coco ; 
    489500       coco(:,:)=pression(5-e+8,:,:); 
    490       hgt=interp2(lon-180,lat,coco,xi,yi);  
     501      hgt=interp2(lon-180,lat,coco,xi,yi); 
    491502       palette 
    492        aa=contourf(xi,yi,hgt,[-20:2:20],'LineStyle','none');  
     503       aa=contourf(xi,yi,hgt,[-20:2:20],'LineStyle','none'); 
    493504       caxis([-20 20]); 
    494505       hold on; 
    … …  
    498509       else 
    499510set (gca,'XTick',[-60:20:120],'XTickLabel',[' ';' ';' ';' ';' ';' ';' ';' ';' ';' '],'fontsize',6); 
    500        end   
     511       end 
    501512       set (gca,'YTick',[-10:10:30],'YTickLabel',['10S';' 0 ';'10N';'20N';'30N'],'fontname','Arial','fontsize',6); 
    502513 
    503514       cartemonde1; 
    504515       if e==2; 
    505           co=colorbar;             
     516          co=colorbar; 
    506517          set (co,'ytick',[-20:10:20],'yticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 
    507518       end; 
    … …  
    515526       text(-150,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 
    516527       pos_vert=e/6.1; 
    517        set (gca,'position',[0.2 pos_vert 0.4 1/6.1])  
     528       set (gca,'position',[0.2 pos_vert 0.4 1/6.1]) 
    518529end; 
    519530 
    520531 
    521 print -depsc2 composite_wind925_EOF12_ER30100_part3.eps;  
    522  
    523  
    524 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
     532print -depsc2 composite_wind925_EOF12_ER30100_part3.eps; 
    525533 
    526534figure(4); 
    … …  
    539547             delta_uwnd925_filtre_min(b,c)=uwnd925_filtre_compomin1(f,b,c)-uwnd925_filtre_compomax1(f,b,c); 
    540548             delta_vwnd925_filtre_min(b,c)=vwnd925_filtre_compomin1(f,b,c)-vwnd925_filtre_compomax1(f,b,c); 
    541   
     549 
    542550         end; 
    543551       end; 
    … …  
    550558       delta_vwnd925_filtre(:,1:72)=delta_vwnd925_filtre_min(:,73:144); 
    551559       delta_vwnd925_filtre(:,73:144)=delta_vwnd925_filtre_min(:,1:72); 
    552        xi=[-180:3.75:176.25]; xi=xi'; yi=[-87.1591:3.708895:87.1591];  
     560       xi=[-180:3.75:176.25]; 
     561       xi=xi'; 
     562       yi=[-87.1591:3.708895:87.1591]; 
    553563       clear carteu cartev; 
    554        carteu=interp2(lon-180,lat,delta_uwnd925_filtre,xi,yi);  
     564       carteu=interp2(lon-180,lat,delta_uwnd925_filtre,xi,yi); 
    555565       cartev=interp2(lon-180,lat,delta_vwnd925_filtre,xi,yi); 
    556566       clear coco ; 
    557567       coco(:,:)=pression(5-e+10,:,:); 
    558       hgt=interp2(lon-180,lat,coco,xi,yi);  
     568      hgt=interp2(lon-180,lat,coco,xi,yi); 
    559569       palette 
    560        aa=contourf(xi,yi,hgt,[-20:2:20],'LineStyle','none');  
     570       aa=contourf(xi,yi,hgt,[-20:2:20],'LineStyle','none'); 
    561571        caxis([-20 20]); 
    562572       hold on; 
    … …  
    566576       else 
    567577set (gca,'XTick',[-60:20:120],'XTickLabel',[' ';' ';' ';' ';' ';' ';' ';' ';' ';' '],'fontsize',6); 
    568        end     
     578       end 
    569579       set (gca,'YTick',[-10:10:30],'YTickLabel',['10S';' 0 ';'10N';'20N';'30N'],'fontname','Arial','fontsize',6); 
    570580 
    571581       cartemonde1; 
    572582       if e==2; 
    573           co=colorbar;             
     583          co=colorbar; 
    574584          set (co,'ytick',[-20:10:20],'yticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 
    575585       end; 
    … …  
    583593       text(-150,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 
    584594       pos_vert=e/6.1; 
    585        set (gca,'position',[0.2 pos_vert 0.4 1/6.1])  
     595       set (gca,'position',[0.2 pos_vert 0.4 1/6.1]) 
    586596end; 
    587597 
    588 print -depsc2 composite_wind925_EOF12_ER30100_part4.eps;  
    589  
    590  
    591  
     598print -depsc2 composite_wind925_EOF12_ER30100_part4.eps; 
Note: See TracChangeset for help on using the changeset viewer.