Changeset 16 for trunk/src/mode_sahelien/composite_pluie_eof23_1030.m

Timestamp:

01/06/09 10:53:18 (15 years ago)

Author:

pinsard

Message:

remove trailing blanks, split lines with multiple statements

File:

• trunk/src/mode_sahelien/composite_pluie_eof23_1030.m (modified) (17 diffs)

trunk/src/mode_sahelien/composite_pluie_eof23_1030.m

@@ +1 @@
+%
+
 %%%%%%%%%%%%%%% FICHIER pluie IRD
 % ouverture / lecture du fichier
@@ -10 +12 @@
 pluie=zeros(41,17,122,12);
 
-clear v fid; fid=fopen('pluvio.79.haut','r','b'); siz=41*17*365;
-v=fread(fid,siz,'float'); fclose(fid);
-v=reshape(v,41,17,365);
-clear a b c;
-for a=1:41;
-   for b=1:17;
-       for c=1:365;
-          if v(a,b,c)>990;
-              v(a,b,c)=NaN;        
+clear v fid;
+fid=fopen('pluvio.79.haut','r','b');
+siz=41*17*365;
+v=fread(fid,siz,'float');
+fclose(fid);
+v=reshape(v,41,17,365);
+clear a b c;
+for a=1:41;
+   for b=1:17;
+       for c=1:365;
+          if v(a,b,c)>990;
+              v(a,b,c)=NaN;
           elseif v(a,b,c)<0;
               v(a,b,c)=NaN;
@@ -35 +40 @@
        for c=1:366;
           if v(a,b,c)>990;
-              v(a,b,c)=NaN;        
+              v(a,b,c)=NaN;
           elseif v(a,b,c)<0;
               v(a,b,c)=NaN;
@@ -52 +57 @@
        for c=1:365;
           if v(a,b,c)>990;
-              v(a,b,c)=NaN;        
+              v(a,b,c)=NaN;
           elseif v(a,b,c)<0;
               v(a,b,c)=NaN;
@@ -69 +74 @@
        for c=1:365;
           if v(a,b,c)>990;
-              v(a,b,c)=NaN;        
+              v(a,b,c)=NaN;
           elseif v(a,b,c)<0;
               v(a,b,c)=NaN;
@@ -86 +91 @@
        for c=1:365;
           if v(a,b,c)>990;
-              v(a,b,c)=NaN;        
+              v(a,b,c)=NaN;
           elseif v(a,b,c)<0;
               v(a,b,c)=NaN;
@@ -105 +110 @@
        for c=1:366;
           if v(a,b,c)>990;
-              v(a,b,c)=NaN;        
+              v(a,b,c)=NaN;
           elseif v(a,b,c)<0;
               v(a,b,c)=NaN;
@@ -122 +127 @@
        for c=1:365;
           if v(a,b,c)>990;
-              v(a,b,c)=NaN;        
+              v(a,b,c)=NaN;
           elseif v(a,b,c)<0;
               v(a,b,c)=NaN;
@@ -139 +144 @@
        for c=1:365;
           if v(a,b,c)>990;
-              v(a,b,c)=NaN;        
+              v(a,b,c)=NaN;
           elseif v(a,b,c)<0;
               v(a,b,c)=NaN;
@@ -156 +161 @@
        for c=1:365;
           if v(a,b,c)>990;
-              v(a,b,c)=NaN;        
+              v(a,b,c)=NaN;
           elseif v(a,b,c)<0;
               v(a,b,c)=NaN;
@@ -173 +178 @@
        for c=1:366;
           if v(a,b,c)>990;
-              v(a,b,c)=NaN;        
+              v(a,b,c)=NaN;
           elseif v(a,b,c)<0;
               v(a,b,c)=NaN;
@@ -190 +195 @@
        for c=1:365;
           if v(a,b,c)>990;
-              v(a,b,c)=NaN;        
+              v(a,b,c)=NaN;
           elseif v(a,b,c)<0;
               v(a,b,c)=NaN;
@@ -207 +212 @@
        for c=1:365;
           if v(a,b,c)>990;
-              v(a,b,c)=NaN;        
+              v(a,b,c)=NaN;
           elseif v(a,b,c)<0;
               v(a,b,c)=NaN;
@@ -222 +227 @@
 clear
 
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 load pluie.mat;
 % pluie dimension : lon lat jours annee (41 17 122 12) => jour annee lat lon  (17 41 122 12)
 
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 load eof23_1030.txt
-
 
 % on met sous la forme annee mois
 eof23_1030=reshape(eof23_1030,122,28);
-olr_eof=eof23_1030(:,1:12); % selection des annees avant 1990 
+olr_eof=eof23_1030(:,1:12); % selection des annees avant 1990
 clear eof23_1030
 
 % on calcule l ecart type sur la serie temporelle de l'olr
-k=122*12; 
+k=122*12;
 olrbis=reshape(olr_eof,k,1);
-olrbis=std(olrbis); 
+olrbis=std(olrbis);
 
 % intialisation des tableaux
-pluie_filtre_compomax=NaN*ones(30,17,41); 
+pluie_filtre_compomax=NaN*ones(30,17,41);
 pluie_filtre_compomin=NaN*ones(30,17,41);
 pluie=permute(pluie,[3 4 2 1]);
-compteurmax=0; 
+compteurmax=0;
 compteurmin=0;
 
-for a=1:12; % boucle sur les annees
-   clear delta; 
+% boucle sur les annees
+for a=1:12;
+   clear delta;
    delta=zeros(122,1);
-   for b=1:122; % boucle sur les jours : on prepare avant   
+   % boucle sur les jours : on prepare avant
+   for b=1:122;
         if b>1;
             delta(b)=olr_eof(b,a)-olr_eof(b-1,a);
         end;
    end;
-   for b=1:122; % boucle sur les jours
-      if (b>20 && b<112); % 1ere condition (laisser une marge de 10) 
-% cas 1 : max deltas apres negetifs et deltas avant postifs   
-         if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && olr_eof(b,a)>olrbis); 
+   % boucle sur les jours
+   for b=1:122;
+      % 1ere condition (laisser une marge de 10)
+      if (b>20 && b<112);
+% cas 1 : max deltas apres negetifs et deltas avant postifs
+         if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && olr_eof(b,a)>olrbis);
             compteurmax=compteurmax+1;
             pluie_filtre_compomax(compteurmax,:,:)=pluie(b,a,:,:);
-% cas 2 : min deltas apres positifs et deltas avant negatifs   
+% cas 2 : min deltas apres positifs et deltas avant negatifs
          elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && olr_eof(b,a)<(-olrbis));
             compteurmin=compteurmin+1;
@@ -267 +273 @@
          end;
       end;
-   end;         
+   end;
 end;
 
@@ -273 +279 @@
 % on fait la somme des valeurs de pluie_filtreentiel sur le nbre d occurences
 % on obtent un tableau de la forme (jours-avant-apres,lon,lat)
-pluie_filtre_compomin1(:,:)=nanmean(pluie_filtre_compomin); clear pluie_filtre_compomin;
-pluie_filtre_compomax1(:,:)=nanmean(pluie_filtre_compomax); clear pluie_filtre_compomax;
+pluie_filtre_compomin1(:,:)=nanmean(pluie_filtre_compomin);
+clear pluie_filtre_compomin;
+pluie_filtre_compomax1(:,:)=nanmean(pluie_filtre_compomax);
+clear pluie_filtre_compomax;
 
 clear pluie_min pluie_max;
@@ -298 +306 @@
 
 palette
-lat=[-10:2.5:30]; lon=[-70:2.5:30]';
+lat=[-10:2.5:30];
+lon=[-70:2.5:30]';
 
 contourf(lon,lat,diff,[-5:0.5:5],'LineStyle','none');
 caxis([-5 5])
-co=colorbar('horiz');            
+co=colorbar('horiz');
 set(co,'position',[0.15 0.26 0.3 0.02]);
 set(co,'xtick',[-4:2:4],'xticklabel',[-4:2:4],'fontname','Arial','fontsize',4);
 
-hold on; 
-
-
+hold on;