Changeset 16 for trunk/src/mode_sahelien/composite_olr_eof234_1030_filtre.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_olr_eof234_1030_filtre.m (modified) (7 diffs)

trunk/src/mode_sahelien/composite_olr_eof234_1030_filtre.m

@@ +1 @@
+%
+
 
 %%%%%%%%%%%%%%% FICHIER OLR
 % ouverture / lecture du fichier
-clear; 
+clear;
 
 fid=fopen('olr_filtre1030_bis.dat','r');
@@ -8 +10 @@
 v=fread(fid,siz,'float');
 fclose(fid);
-olr=reshape(v,144,25,122,28); clear v;
-olr1=permute(olr,[3 4 2 1]); clear olr
+olr=reshape(v,144,25,122,28);
+clear v;
+olr1=permute(olr,[3 4 2 1]);
+clear olr
 
 cd('..');
@@ -17 +21 @@
 cd('mode_sahelien');
 
-
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 load eof24_1030.txt
 
 
 % on met sous la forme annee mois
-olr_eof=reshape(eof24_1030,122,28); clear eof24_1030;
+olr_eof=reshape(eof24_1030,122,28);
+clear eof24_1030;
 
 % on calcule l ecart type sur la serie temporelle de l OLR
-k=122*28; 
+k=122*28;
 olrbis=reshape(olr_eof,k,1);
-olrbis=std(olrbis); 
+olrbis=std(olrbis);
 
 % intialisation des tableaux
-olr_filtre_compomax=NaN*ones(50,25,144); 
+olr_filtre_compomax=NaN*ones(50,25,144);
 olr_filtre_compomin=NaN*ones(50,25,144);
 
-compteurmax=0; 
+compteurmax=0;
 compteurmin=0;
 
-for a=1:28; % boucle sur les annees
-   clear delta; 
+% boucle sur les annees
+for a=1:28;
+   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;
             olr_filtre_compomax(compteurmax,:,:)=olr1(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;
@@ -58 +64 @@
          end;
       end;
-   end;         
+   end;
 end;
 compteurmax_eof234_1030=compteurmax;
@@ -66 +72 @@
 % on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences
 % on obtent un tableau de la forme (jours-avant-apres,lon,lat)
-olr_filtre_compomin1(:,:)=nanmean(olr_filtre_compomin); clear olr_filtre_compomin;
-olr_filtre_compomax1(:,:)=nanmean(olr_filtre_compomax); clear olr_filtre_compomax;
+olr_filtre_compomin1(:,:)=nanmean(olr_filtre_compomin);
+clear olr_filtre_compomin;
+olr_filtre_compomax1(:,:)=nanmean(olr_filtre_compomax);
+clear olr_filtre_compomax;
 
-x=[-10:1:10]'; %'
+x=[-10:1:10]';
 
 delta_olr_filtre_min=zeros(25,144);
@@ -90 +98 @@
 
 
-clear fin; fin=-10+16*2.5;
+clear fin;
+fin=-10+16*2.5;
 [cs,h]=contour(lon,lat,delta_olr_filtre_min,[5:5:25],'k');
 clabel(cs,h,[10 20 30],'color','k','rotation',0,'fontsize',5);
 hold on;
 
-clear cs h; hold on
+clear cs h;
+hold on;
 [cs,h]=contour(lon,lat,delta_olr_filtre_min,[-25:5:-5],'k:','Linewidth',0.3);
 clabel(cs,h,[-30 -20 -10],'color','k','rotation',0,'fontsize',5);
@@ -101 +111 @@
 cartemonde1;
 
-axis equal; axis([-80 100 -10 30]);
+axis equal;
+axis([-80 100 -10 30]);
 
 set (gca,'XTick',[-180:60:180],'XTickLabel',[' 180';'120W';' 60W';'  0 ';' 60E';'120E';'180 '],'fontname','Arial','fontsize',6);