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

trunk/src/mode_sahelien/composite_hgt925_sahel.m

@@ +1 @@
+%
+
 % initialisation
-clear; close all;
-
+clear;
+close all;
 
 % etape 1 :on selectionne les mois d ete (attention aux annees bixestiles!)
@@ -8 +10 @@
 hgt700=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('hgt700.1979.nc'); hgt700(:,:,:,1)=hgt(siz:sizbis,1,:,:); clear hgt;
 ncload('hgt700.1980.nc'); hgt700(:,:,:,2)=hgt(siz+1:sizbis+1,1,:,:); clear hgt;
@@ -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;
+
 % GEOPOTENTIEL
 % intialisation des tableaux
-hgt700_filtre_compomax=NaN*ones(60,31,49,144); 
+hgt700_filtre_compomax=NaN*ones(60,31,49,144);
 hgt700_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;
             hgt700_filtre_compomax(compteurmax,1,:,:)=hgt(b-20,a,:,:);
@@ -109 +120 @@
             hgt700_filtre_compomax(compteurmax,30,:,:)=hgt(b+9,a,:,:);
             hgt700_filtre_compomax(compteurmax,31,:,:)=hgt(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;
@@ -148 +159 @@
          end;
       end;
-   end;         
+   end;
 end;
 
@@ -154 +165 @@
 % on fait la somme des valeurs de hgt700_filtreentiel sur le nbre d occurences
 % on obtent un tableau de la forme (jours-avant-apres,lon,lat)
-hgt700_filtre_compomin1(:,:,:)=nanmean(hgt700_filtre_compomin); clear hgt700_filtre_compomin;
-hgt700_filtre_compomax1(:,:,:)=nanmean(hgt700_filtre_compomax); clear hgt700_filtre_compomax;
-
-x=[-10:1:10]'; %'
-
-figure(1); 
+hgt700_filtre_compomin1(:,:,:)=nanmean(hgt700_filtre_compomin);
+clear hgt700_filtre_compomin;
+hgt700_filtre_compomax1(:,:,:)=nanmean(hgt700_filtre_compomax);
+clear hgt700_filtre_compomax;
+
+x=[-10:1:10]';
+
+figure(1);
 orient('landscape');
 lat=lat';
@@ -183 +196 @@
           end;
        end
-       subplot(7,1,8-e); hold on;
+       subplot(7,1,8-e);
+       hold on;
        palette
        clear delta_hgt700;
@@ -197 +211 @@
        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);
        if e==4;
-          co=colorbar;            
+          co=colorbar;
           set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6);
        end;
-       hold on; cartemonde1;
+       hold on;
+       cartemonde1;
        set (gca,'tickDir','out');
        axis equal;
@@ -212 +227 @@
        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 compo_hgt700_sahel.eps
 
-figure(2); 
+figure(2);
 orient('landscape');
 lat=lat';
@@ -240 +255 @@
           end;
        end
-       subplot(7,1,8-e); hold on;
+       subplot(7,1,8-e);
+       hold on;
        palette
        clear delta_hgt700;
@@ -254 +270 @@
        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);
        if e==4;
-          co=colorbar;            
+          co=colorbar;
           set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6);
        end;
-       hold on; cartemonde1;
+       hold on;
+       cartemonde1;
        set (gca,'tickDir','out');
        axis equal;
@@ -269 +286 @@
        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;