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

trunk/src/mode_sahelien/composite_wind925_sahel.m

@@ +1 @@
+%
+
 % initialisation
 close all
@@ -8 +10 @@
 uwnd700=zeros(122,49,144,22);
 
-siz=31+28+31+30+31+1; sizbis=31+28+31+30+31+30+31+31+30;
+siz=31+28+31+30+31+1;
+sizbis=31+28+31+30+31+30+31+31+30;
 ncload('uwnd700.1979.nc'); uwnd700(:,:,:,1)=uwnd(siz:sizbis,1,:,:); clear uwnd;
 ncload('uwnd700.1980.nc'); uwnd700(:,:,:,2)=uwnd(siz+1:sizbis+1,1,:,:); clear uwnd;
@@ -38 +41 @@
 
 %%%%%%%%%%%%%%% FICHIER OLR
-% ouverture / lecture du fichier 
+% ouverture / lecture du fichier
 fid=fopen('sahel_pb.dat','r');
 v=fread(fid,3538,'float');
@@ -44 +47 @@
 
 % on met sous la forme annee mois
-ind_olr=reshape(v,122,29); clear v; 
+ind_olr=reshape(v,122,29);
+clear v;
 
 % on selectionne les annees de 1979 a 2000
-ind_olr1=ind_olr(:,1:22); clear ind_olr;
-k=122*22; ind_olr=ind_olr1; clear ind_olr1;
+ind_olr1=ind_olr(:,1:22);
+clear ind_olr;
+k=122*22;
+ind_olr=ind_olr1;
+clear ind_olr1;
 
 % on calcule l ecart type sur la serie entiere de l OLR
 ind_olr1=reshape(ind_olr, k,1);
-ind_olr2=std(ind_olr1); clear ind_olr1;
+ind_olr2=std(ind_olr1);
+clear ind_olr1;
 
 
 % intialisation des tableaux
-uwnd700_filtre_compomax=NaN*ones(60,31,49,144); 
+uwnd700_filtre_compomax=NaN*ones(60,31,49,144);
 uwnd700_filtre_compomin=NaN*ones(60,31,49,144);
 
-compteurmax=0; 
+compteurmax=0;
 compteurmin=0;
 
-for a=1:22; % boucle sur les annees
-   clear delta; 
+% boucle sur les annees
+for a=1:22;
+   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)=ind_olr(b,a)-ind_olr(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 && ind_olr(b,a)>ind_olr2); 
+   % 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 && ind_olr(b,a)>ind_olr2);
             compteurmax=compteurmax+1;
             uwnd700_filtre_compomax(compteurmax,1,:,:)=uwnd(b-20,a,:,:);
@@ -108 +120 @@
             uwnd700_filtre_compomax(compteurmax,30,:,:)=uwnd(b+9,a,:,:);
             uwnd700_filtre_compomax(compteurmax,31,:,:)=uwnd(b+10,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 && ind_olr(b,a)<(-ind_olr2));
             compteurmin=compteurmin+1;
@@ -147 +159 @@
          end;
       end;
-   end;         
+   end;
 end;
 
@@ -153 +165 @@
 % on fait la somme des valeurs de uwnd700_filtreentiel sur le nbre d occurences
 % on obtent un tableau de la forme (jours-avant-apres,lon,lat)
-uwnd700_filtre_compomin1(:,:,:)=nanmean(uwnd700_filtre_compomin); clear uwnd700_filtre_compomin;
-uwnd700_filtre_compomax1(:,:,:)=nanmean(uwnd700_filtre_compomax); clear uwnd700_filtre_compomax;
-
-
-
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+uwnd700_filtre_compomin1(:,:,:)=nanmean(uwnd700_filtre_compomin);
+clear uwnd700_filtre_compomin;
+uwnd700_filtre_compomax1(:,:,:)=nanmean(uwnd700_filtre_compomax);
+clear uwnd700_filtre_compomax;
 
 % etape 1 :on selectionne les mois d ete (attention aux annees bixestiles!)
@@ -166 +175 @@
 vwnd700=zeros(122,49,144,22);
 
-siz=31+28+31+30+31+1; sizbis=31+28+31+30+31+30+31+31+30;
+siz=31+28+31+30+31+1;
+sizbis=31+28+31+30+31+30+31+31+30;
 ncload('vwnd700.1979.nc'); vwnd700(:,:,:,1)=vwnd(siz:sizbis,1,:,:); clear vwnd;
 ncload('vwnd700.1980.nc'); vwnd700(:,:,:,2)=vwnd(siz+1:sizbis+1,1,:,:); clear vwnd;
@@ -195 +205 @@
 
 %%%%%%%%%%%%%% FICHIER OLR
-% ouverture / lecture du fichier 
+% ouverture / lecture du fichier
 fid=fopen('sahel_pb.dat','r');
 v=fread(fid,3538,'float');
@@ -201 +211 @@
 
 % on met sous la forme annee mois
-ind_olr=reshape(v,122,29); clear v; 
+ind_olr=reshape(v,122,29);
+clear v;
 
 % on selectionne les annees de 1979 a 2000
-ind_olr1=ind_olr(:,1:22); clear ind_olr;
-k=122*22; ind_olr=ind_olr1; clear ind_olr1;
+ind_olr1=ind_olr(:,1:22);
+clear ind_olr;
+k=122*22;
+ind_olr=ind_olr1;
+clear ind_olr1;
 
 % on calcule l ecart type sur la serie entiere de l OLR
 ind_olr1=reshape(ind_olr, k,1);
-ind_olr2=std(ind_olr1); clear ind_olr1;
-
+ind_olr2=std(ind_olr1);
+clear ind_olr1;
 
 % intialisation des tableaux
-vwnd700_filtre_compomax=NaN*ones(60,31,49,144); 
+vwnd700_filtre_compomax=NaN*ones(60,31,49,144);
 vwnd700_filtre_compomin=NaN*ones(60,31,49,144);
 
-compteurmax=0; 
+compteurmax=0;
 compteurmin=0;
 
-for a=1:22; % boucle sur les annees
-   clear delta; 
+% boucle sur les annees
+for a=1:22;
+   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)=ind_olr(b,a)-ind_olr(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 && ind_olr(b,a)>ind_olr2); 
+   % 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 && ind_olr(b,a)>ind_olr2);
             compteurmax=compteurmax+1;
             vwnd700_filtre_compomax(compteurmax,1,:,:)=vwnd(b-20,a,:,:);
@@ -265 +283 @@
             vwnd700_filtre_compomax(compteurmax,30,:,:)=vwnd(b+9,a,:,:);
             vwnd700_filtre_compomax(compteurmax,31,:,:)=vwnd(b+10,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 && ind_olr(b,a)<(-ind_olr2));
             compteurmin=compteurmin+1;
@@ -304 +322 @@
          end;
       end;
-   end;         
+   end;
 end;
 
@@ -310 +328 @@
 % on fait la somme des valeurs de vwnd700_filtreentiel sur le nbre d occurences
 % on obtent un tableau de la forme (jours-avant-apres,lon,lat)
-vwnd700_filtre_compomin1(:,:,:)=nanmean(vwnd700_filtre_compomin); clear vwnd700_filtre_compomin;
-vwnd700_filtre_compomax1(:,:,:)=nanmean(vwnd700_filtre_compomax); clear vwnd700_filtre_compomax;
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+vwnd700_filtre_compomin1(:,:,:)=nanmean(vwnd700_filtre_compomin);
+clear vwnd700_filtre_compomin;
+vwnd700_filtre_compomax1(:,:,:)=nanmean(vwnd700_filtre_compomax);
+clear vwnd700_filtre_compomax;
 
 figure(1)
@@ -331 +349 @@
              delta_uwnd700_filtre_min(b,c)=uwnd700_filtre_compomin1(f,b,c)-uwnd700_filtre_compomax1(f,b,c);
              delta_vwnd700_filtre_min(b,c)=vwnd700_filtre_compomin1(f,b,c)-vwnd700_filtre_compomax1(f,b,c);
- 
+
          end;
        end;
@@ -342 +360 @@
        delta_uwnd700_filtre(:,1:72)=delta_uwnd700_filtre_min(:,73:144);
        delta_uwnd700_filtre(:,73:144)=delta_uwnd700_filtre_min(:,1:72);
-       xi=[-180:3.75:176.25]; xi=xi'; yi=[-87.1591:3.708895:87.1591]; 
+       xi=[-180:3.75:176.25];
+       xi=xi';
+       yi=[-87.1591:3.708895:87.1591];
        clear carteu cartev;
-       carteu=interp2(lon-180,lat,delta_uwnd700_filtre,xi,yi); 
+       carteu=interp2(lon-180,lat,delta_uwnd700_filtre,xi,yi);
        cartev=interp2(lon-180,lat,delta_vwnd700_filtre,xi,yi);
        clear coco ;
        coco(:,:)=pression(e,:,:);
-      hgt=interp2(lon-180,lat,coco,xi,yi); 
+       hgt=interp2(lon-180,lat,coco,xi,yi);
        palette
-       aa=contourf(xi,yi,hgt,[-20:2:20],'LineStyle','none'); hold on;
+       aa=contourf(xi,yi,hgt,[-20:2:20],'LineStyle','none');
+       hold on;
        aa=quiver(xi,yi,carteu,cartev,'k');
        if e==1;
@@ -356 +377 @@
        else
       set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6);
-       end       
+       end
        set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6);
        cartemonde1;
        if e==4;
-          co=colorbar;            
+          co=colorbar;
           set (co,'xtick',[-40:20:40],'xticklabel',[-40:20:40],'fontname','Arial','fontsize',6);
        end;
@@ -372 +393 @@
        text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10);
        pos_vert=e/10;
-       set (gca,'position',[0.2 pos_vert 0.5 0.1]) 
+       set (gca,'position',[0.2 pos_vert 0.5 0.1])
 end;
 
-
 print -depsc2 composite_sahel_1_u700_hgt.eps;
 
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%ù
 figure(2)
 orient('landscape')
@@ -395 +414 @@
              delta_uwnd700_filtre_min(b,c)=uwnd700_filtre_compomin1(f,b,c)-uwnd700_filtre_compomax1(f,b,c);
              delta_vwnd700_filtre_min(b,c)=vwnd700_filtre_compomin1(f,b,c)-vwnd700_filtre_compomax1(f,b,c);
- 
+
          end;
        end;
@@ -406 +425 @@
        delta_uwnd700_filtre(:,1:72)=delta_uwnd700_filtre_min(:,73:144);
        delta_uwnd700_filtre(:,73:144)=delta_uwnd700_filtre_min(:,1:72);
-       xi=[-180:3.75:176.25]; xi=xi'; yi=[-87.1591:3.708895:87.1591]; 
+       xi=[-180:3.75:176.25];
+       xi=xi';
+       yi=[-87.1591:3.708895:87.1591];
        clear carteu cartev;
-       carteu=interp2(lon-180,lat,delta_uwnd700_filtre,xi,yi); 
+       carteu=interp2(lon-180,lat,delta_uwnd700_filtre,xi,yi);
        cartev=interp2(lon-180,lat,delta_vwnd700_filtre,xi,yi);
        clear coco ;
        coco(:,:)=pression(e+7,:,:);
-      hgt=interp2(lon-180,lat,coco,xi,yi); 
+       hgt=interp2(lon-180,lat,coco,xi,yi);
        palette
-       aa=contourf(xi,yi,hgt,[-20:2:20],'LineStyle','none'); hold on;
+       aa=contourf(xi,yi,hgt,[-20:2:20],'LineStyle','none');
+       hold on;
        aa=quiver(xi,yi,carteu,cartev,'k');
        if e==1;
@@ -420 +442 @@
        else
       set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6);
-       end       
+       end
        set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6);
        cartemonde1;
        if e==4;
-          co=colorbar;            
+          co=colorbar;
           set (co,'xtick',[-40:20:40],'xticklabel',[-40:20:40],'fontname','Arial','fontsize',6);
        end;
@@ -436 +458 @@
        text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10);
        pos_vert=e/10;
-       set (gca,'position',[0.2 pos_vert 0.5 0.1]) 
+       set (gca,'position',[0.2 pos_vert 0.5 0.1])
 end;
 
-
 print -depsc2 composite_sahel_2_u700_hgt.eps;
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%ù