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

trunk/src/mode_sahelien/reconstitution_1025.m

@@ +1 @@
+%
 
-clear; close all;
+clear;
+close all;
 
 %lon=[-30:2.5:30], lat=[-10:2.5:30]
@@ -9 +11 @@
 
 cosa=zeros(4,3416,425);
-for a=1:4; % boucle sur les VP
+% boucle sur les VP
+for a=1:4;
       for b=1:3416 % boucle sur le temps
           for c=1:425 % boucle sur l espace
@@ -17 +20 @@
 end;
 
-cosa1(:,:)=cosa(1,:,:);  
-cosa1=reshape(cosa1,3416,25,17); 
+cosa1(:,:)=cosa(1,:,:);
+cosa1=reshape(cosa1,3416,25,17);
 
-cosa2(:,:)=cosa(2,:,:);  
-cosa2=reshape(cosa2,3416,25,17); 
+cosa2(:,:)=cosa(2,:,:);
+cosa2=reshape(cosa2,3416,25,17);
 
-cosabis=cosa(2:3,:,:); 
-cosa23(:,:)=sum(cosabis); clear cosabis;
-cosa23=reshape(cosa23,3416,25,17); 
+cosabis=cosa(2:3,:,:);
+cosa23(:,:)=sum(cosabis);
+clear cosabis;
+cosa23=reshape(cosa23,3416,25,17);
 
-cosabis=cosa(2:4,:,:); 
-cosa24(:,:)=sum(cosabis); clear cosabis;
-cosa24=reshape(cosa24,3416,25,17); 
+cosabis=cosa(2:4,:,:);
+cosa24(:,:)=sum(cosabis);
+clear cosabis;
+cosa24=reshape(cosa24,3416,25,17);
 
-figure(1); orient('landscape')
-figure(2); orient('landscape')
-figure(3); orient('landscape')
+figure(1);
+orient('landscape')
+figure(2);
+orient('landscape')
+figure(3);
+orient('landscape')
 
 % On selectionne la demaine geographique
 
 % la latitude : 5N (7) 12.5N (10)
-% la longitude : 10W (9) 10E(17) 
-lat1=7; lat2=10;
-lon1=9; lon2=17;
+% la longitude : 10W (9) 10E(17)
+lat1=7;
+lat2=10;
+lon1=9;
+lon2=17;
 taille=(lat2-lat1+1)*(lon2-lon1+1);
 
-
-cosa23=cosa23(:,lon1:lon2,lat1:lat2); clear cosi2;
-cosa23=reshape(cosa23,3416,taille); 
+cosa23=cosa23(:,lon1:lon2,lat1:lat2);
+clear cosi2;
+cosa23=reshape(cosa23,3416,taille);
 cosa23=cosa23';
-cosa23=mean(cosa23); 
+cosa23=mean(cosa23);
 cosa23=cosa23';
 save eof23_1025.txt cosa23 -ascii;
 cosa23=reshape(cosa23,122,28);
 
-cosa24=cosa24(:,lon1:lon2,lat1:lat2); clear cosi2;
-cosa24=reshape(cosa24,3416,taille); 
+cosa24=cosa24(:,lon1:lon2,lat1:lat2);
+clear cosi2;
+cosa24=reshape(cosa24,3416,taille);
 cosa24=cosa24';
-cosa24=mean(cosa24); 
+cosa24=mean(cosa24);
 cosa24=cosa24';
 save eof24_1025.txt cosa24 -ascii;
@@ -66 +77 @@
 
 % 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 2006
-ind_olr=ind_olr(:,1:28); 
-
-
+ind_olr=ind_olr(:,1:28);
 
 for a=1:28;
      if a<12;
-        figure(1); subplot(4,3,a);
+        figure(1);
+        subplot(4,3,a);
      elseif a<24
-        figure(2); subplot(4,3,a-11);
+        figure(2);
+        subplot(4,3,a-11);
     else
-        figure(3); subplot(4,3,a-23);
+        figure(3);
+        subplot(4,3,a-23);
      end
      y2=ind_olr(:,a);
      y4=cosa24(:,a);
-     x=[1:122]'; 
-     plot(x,y2); hold on;
+     x=[1:122]';
+     plot(x,y2);
+     hold on;
      plot(x,y4,'k');
      y5=zeros(122,1);
      plot(x,y5,'k:');
-     
+
      title(a+1978,'Fontsize',8);
      if a==2;
@@ -110 +124 @@
    set (gca,'XTick',[10:10:122],'XTickLabel',[10:10:122],'fontname','Arial','fontsize',6);
 
-end; 
+end;
 
-figure(1); print -depsc2 reconstitution1_eof234_1025.eps;
-figure(2); print -depsc2 reconstitution2_eof234_1025.eps;
-figure(3); print -depsc2 reconstitution3_eof234_1025.eps;
-
+figure(1);
+print -depsc2 reconstitution1_eof234_1025.eps;
+figure(2);
+print -depsc2 reconstitution2_eof234_1025.eps;
+figure(3);
+print -depsc2 reconstitution3_eof234_1025.eps;