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_olr_phase5.m

    r15 r16  
     1% 
     2 
    13% initialisation 
    2 clear;  
     4clear; 
    35 
    46cd('..'); 
     
    2527olr1=zeros(siz,21,144); 
    2628olr1(:,:,1:72)=olr(:,:,73:144); 
    27 olr1(:,:,73:144)=olr(:,:,1:72); clear olr; 
     29olr1(:,:,73:144)=olr(:,:,1:72); 
     30clear olr; 
    2831 
    2932% on permture le temps et l'espce 
     
    3740      vecteur=zeros(siz,1); 
    3841      for c=1:siz; 
    39            vecteur(c)=olr1(a,b,c);     
     42           vecteur(c)=olr1(a,b,c); 
    4043      end 
    4144      std_c=nanmean(vecteur); 
    4245      for c=1:siz; 
    43           olr1(a,b,c)=olr1(a,b,c)-std_c;     
    44       end  
     46          olr1(a,b,c)=olr1(a,b,c)-std_c; 
     47      end 
    4548   end 
    46 end;    
    47  
    48 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
     49end; 
     50 
    4951% on met sous la forme longitude latitude jour 
    50  
    5152 
    5253% on selectionne les mois d ete (attention aux annees bixestiles!) 
     
    151152clear olr; 
    152153olr2=reshape(olr1,3416,21,144); 
    153 olr3(:,:)=mean(olr2); clear olr2; 
    154 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
     154olr3(:,:)=mean(olr2); 
     155clear olr2; 
    155156 
    156157load phase5.txt; 
    157158phase5=reshape(phase5,122,28,3); 
    158159 olr_filtre_compomax=zeros(90,21,144); 
    159 compteurmax=0;  
    160  
    161 for a=1:28; % boucle sur les annees 
    162    for b=1:122; % boucle sur les jours 
    163       if (b>11 && b<112); % 1ere condition (laisser une marge de 10)  
    164          if phase5(b,a,1)==1;  
     160compteurmax=0; 
     161 
     162% boucle sur les annees 
     163for a=1:28; 
     164   % boucle sur les jours 
     165   for b=1:122; 
     166      % 1ere condition (laisser une marge de 10) 
     167      if (b>11 && b<112); 
     168         if phase5(b,a,1)==1; 
    165169            compteurmax=compteurmax+1; 
    166170            olr_filtre_compomax(compteurmax,:,:)=olr1(b,a,:,:); 
    167171         end; 
    168172      end; 
    169    end;          
     173   end; 
    170174end; 
    171175 
     
    175179% on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 
    176180% on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
    177 olr_filtre_compomax1(:,:)=mean(olr_filtre_compomax); clear olr_filtre_compomax; 
    178  
    179  
    180 x=[-10:1:10]'; %' 
    181  
    182 lat=lat(21:41)';  
     181olr_filtre_compomax1(:,:)=mean(olr_filtre_compomax); 
     182clear olr_filtre_compomax; 
     183 
     184x=[-10:1:10]'; 
     185 
     186lat=lat(21:41)'; 
    183187lon=lon-180; 
    184 delta=olr_filtre_compomax1-olr3; clear olr3; 
     188delta=olr_filtre_compomax1-olr3; 
     189clear olr3; 
    185190for b=1:21; 
    186191    for c=1:144; 
     
    198203set (gca,'XTick',[-180:60:180],'XTickLabel',[' 180';'120W';' 60W';'   0';' 60E';'120E';' 180'],'fontsize',6); 
    199204set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6); 
    200 co=colorbar;             
     205co=colorbar; 
    201206set(co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 
    202 hold on; cartemonde1; 
     207hold on; 
     208cartemonde1; 
    203209set (gca,'tickDir','out'); 
    204210axis equal; 
    205211axis([-180 180 -10 40]); 
    206  
    207  
    208  
    209  
    210  
    211  
    212  
Note: See TracChangeset for help on using the changeset viewer.