Ignore:
Timestamp:
01/06/09 10:53:18 (15 years ago)
Author:
pinsard
Message:

remove trailing blanks, split lines with multiple statements

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/mode_sahelien/reconstitution_1030_std.m

    r15 r16  
     1% 
    12 
    2 clear;  
     3clear; 
    34close all; 
    45 
     
    1011 
    1112cosa=zeros(4,3416,425); 
    12 for a=1:4; % boucle sur les VP 
     13% boucle sur les VP 
     14for a=1:4; 
    1315      for b=1:3416 % boucle sur le temps 
    1416          for c=1:425 % boucle sur l espace 
     
    1820end; 
    1921 
    20 cosa1(:,:)=cosa(1,:,:);   
    21 cosa1=reshape(cosa1,3416,25,17);  
     22cosa1(:,:)=cosa(1,:,:); 
     23cosa1=reshape(cosa1,3416,25,17); 
    2224 
    23 cosabis=cosa(2:4,:,:);  
    24 cosa24(:,:)=sum(cosabis); clear cosabis; 
    25 cosa24=reshape(cosa24,3416,25,17);  
     25cosabis=cosa(2:4,:,:); 
     26cosa24(:,:)=sum(cosabis); 
     27clear cosabis; 
     28cosa24=reshape(cosa24,3416,25,17); 
    2629 
    27  
    28  
    29 figure(1); orient('landscape') 
    30 figure(2); orient('landscape') 
    31 figure(3); orient('landscape') 
     30figure(1); 
     31orient('landscape') 
     32figure(2); 
     33orient('landscape') 
     34figure(3); 
     35orient('landscape') 
    3236 
    3337% On selectionne la demaine geographique 
    3438 
    3539% la latitude : 7.5N (8) 12.5N (10) 
    36 % la longitude : 10W (9) 10E(17)  
    37 lat1=8; lat2=10; 
    38 lon1=9; lon2=17; 
     40% la longitude : 10W (9) 10E(17) 
     41lat1=8; 
     42lat2=10; 
     43lon1=9; 
     44lon2=17; 
    3945taille=(lat2-lat1+1)*(lon2-lon1+1); 
    4046 
    41 cosa1=cosa1(:,lon1:lon2,lat1:lat2); clear cosi2; 
    42 cosa1=reshape(cosa1,3416,taille);  
     47cosa1=cosa1(:,lon1:lon2,lat1:lat2); 
     48clear cosi2; 
     49cosa1=reshape(cosa1,3416,taille); 
    4350cosa1=cosa1'; 
    44 cosa1=mean(cosa1);  
     51cosa1=mean(cosa1); 
    4552cosa1=cosa1'; 
    4653std1=std(cosa1); 
     
    4855cosa1=reshape(cosa1,122,28); 
    4956 
    50 cosa24=cosa24(:,lon1:lon2,lat1:lat2); clear cosi2; 
    51 cosa24=reshape(cosa24,3416,taille);  
     57cosa24=cosa24(:,lon1:lon2,lat1:lat2); 
     58clear cosi2; 
     59cosa24=reshape(cosa24,3416,taille); 
    5260cosa24=cosa24'; 
    53 cosa24=mean(cosa24);  
     61cosa24=mean(cosa24); 
    5462cosa24=cosa24'; 
    5563std24=std(cosa24); 
     
    6068for a=1:28; 
    6169     if a<12; 
    62         figure(1); subplot(4,3,a); 
     70        figure(1); 
     71        subplot(4,3,a); 
    6372     elseif a<23 
    64         figure(2); subplot(4,3,a-11); 
     73        figure(2); 
     74        subplot(4,3,a-11); 
    6575    else 
    66         figure(3); subplot(4,3,a-22); 
     76        figure(3); 
     77        subplot(4,3,a-22); 
    6778     end 
    6879     y1=cosa1(:,a); 
    6980     y2=cosa24(:,a); 
    7081 
    71      x=[1:122]';  
    72      plot(x,y1,'k','Linewidth',1); hold on; 
     82     x=[1:122]'; 
     83     plot(x,y1,'k','Linewidth',1); 
     84     hold on; 
    7385     plot(x,y2,'r','Linewidth',1); 
    7486     y5=zeros(122,1); 
    7587     plot(x,y5,'k:'); 
    76       
     88 
    7789     title(a+1978,'Fontsize',8); 
    7890     if a==2; 
     
    96108   set (gca,'XTick',[10:10:122],'XTickLabel',[10:10:122],'fontname','Arial','fontsize',6); 
    97109 
    98 end;  
     110end; 
    99111 
    100 figure(1); print -depsc2 reconstitution1_eof23_filtreER_1030.eps; 
    101 figure(2); print -depsc2 reconstitution2_eof23_filtreER_1030.eps; 
    102 figure(3); print -depsc2 reconstitution3_eof23_filtreER_1030.eps; 
    103  
     112figure(1); 
     113print -depsc2 reconstitution1_eof23_filtreER_1030.eps; 
     114figure(2); 
     115print -depsc2 reconstitution2_eof23_filtreER_1030.eps; 
     116figure(3); 
     117print -depsc2 reconstitution3_eof23_filtreER_1030.eps; 
Note: See TracChangeset for help on using the changeset viewer.