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

    r15 r16  
     1% 
     2 
    13% initialisation 
    2 clear; close all; 
    3  
     4clear; 
     5close all; 
    46 
    57% etape 1 :on selectionne les mois d ete (attention aux annees bixestiles!) 
     
    810vwnd700=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('vwnd700.1979.nc'); vwnd700(:,:,:,1)=vwnd(siz:sizbis,1,:,:); clear vwnd; 
    1215ncload('vwnd700.1980.nc'); vwnd700(:,:,:,2)=vwnd(siz+1:sizbis+1,1,:,:); clear vwnd; 
     
    4851 
    4952%%%%%%%%%%%%%%% FICHIER OLR 
    50 % ouverture / lecture du fichier  
     53% ouverture / lecture du fichier 
    5154load eof12_rossby_30100.txt 
    5255ind_olr=reshape(eof12_rossby_30100,122,28); 
     
    5558% on calcule l ecart type sur la serie entiere de l OLR 
    5659ind_olr1=reshape(ind_olr, k,1); 
    57 ind_olr2=std(ind_olr1); clear ind_olr1; 
    58  
    59 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
     60ind_olr2=std(ind_olr1); 
     61clear ind_olr1; 
     62 
    6063% GEOPOTENTIEL 
    6164% intialisation des tableaux 
    62 vwnd700_filtre_compomax=NaN*ones(60,31,49,144);  
     65vwnd700_filtre_compomax=NaN*ones(60,31,49,144); 
    6366vwnd700_filtre_compomin=NaN*ones(60,31,49,144); 
    6467 
    65 compteurmax=0;  
     68compteurmax=0; 
    6669compteurmin=0; 
    6770 
    68 for a=1:28; % boucle sur les annees 
    69    clear delta;  
     71% boucle sur les annees 
     72for a=1:28; 
     73   clear delta; 
    7074   delta=zeros(122,1); 
    71    for b=1:122; % boucle sur les jours : on prepare avant    
     75   % boucle sur les jours : on prepare avant 
     76   for b=1:122; 
    7277        if b>1; 
    7378            delta(b)=ind_olr(b,a)-ind_olr(b-1,a); 
    7479        end; 
    7580   end; 
    76    for b=1:122; % boucle sur les jours 
    77       if (b>22 && b<112); % 1ere condition (laisser une marge de 10)  
    78 % cas 1 : max deltas apres negetifs et deltas avant postifs    
    79          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>22 && 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); 
    8087            compteurmax=compteurmax+1; 
    8188            vwnd700_filtre_compomax(compteurmax,1,:,:)=vwnd(b-21,a,:,:); 
     
    113120            vwnd700_filtre_compomax(compteurmax,31,:,:)=vwnd(b+9,a,:,:); 
    114121 
    115 % cas 2 : min deltas apres positifs et deltas avant negatifs    
     122% cas 2 : min deltas apres positifs et deltas avant negatifs 
    116123         elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && ind_olr(b,a)<(-ind_olr2)); 
    117124            compteurmin=compteurmin+1; 
     
    149156            vwnd700_filtre_compomin(compteurmin,30,:,:)=vwnd(b+8,a,:,:); 
    150157            vwnd700_filtre_compomin(compteurmin,31,:,:)=vwnd(b+9,a,:,:); 
    151             
     158 
    152159         end; 
    153160      end; 
    154    end;          
     161   end; 
    155162end; 
    156163 
     
    158165% on fait la somme des valeurs de vwnd700_filtreentiel sur le nbre d occurences 
    159166% on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
    160 vwnd700_filtre_compomin1(:,:,:)=nanmean(vwnd700_filtre_compomin); clear vwnd700_filtre_compomin; 
    161 vwnd700_filtre_compomax1(:,:,:)=nanmean(vwnd700_filtre_compomax); clear vwnd700_filtre_compomax; 
    162  
    163 x=[-10:1:10]'; %' 
    164  
    165 figure(1);  
     167vwnd700_filtre_compomin1(:,:,:)=nanmean(vwnd700_filtre_compomin); 
     168clear vwnd700_filtre_compomin; 
     169vwnd700_filtre_compomax1(:,:,:)=nanmean(vwnd700_filtre_compomax); 
     170clear vwnd700_filtre_compomax; 
     171 
     172x=[-10:1:10]'; 
     173 
     174figure(1); 
    166175orient('landscape'); 
    167176lat=lat'; 
Note: See TracChangeset for help on using the changeset viewer.