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_ER_rossby_1030.m

    r15 r16  
     1% 
    12 
    23clear; 
     
    1011 
    1112cosa=zeros(3,3416,425); 
    12 for a=1:3; % boucle sur les VP 
     13% boucle sur les VP 
     14for a=1:3; 
    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 cosa2(:,:)=cosa(2,:,:);   
    24 cosa2=reshape(cosa2,3416,25,17);  
     25cosa2(:,:)=cosa(2,:,:); 
     26cosa2=reshape(cosa2,3416,25,17); 
    2527 
    26 cosabis=cosa(1:2,:,:);  
    27 cosa12(:,:)=sum(cosabis); clear cosabis; 
    28 cosa12=reshape(cosa12,3416,25,17);  
     28cosabis=cosa(1:2,:,:); 
     29cosa12(:,:)=sum(cosabis); 
     30clear cosabis; 
     31cosa12=reshape(cosa12,3416,25,17); 
    2932 
    30  
    31  
    32 figure(1); orient('landscape') 
    33 figure(2); orient('landscape') 
    34 figure(3); orient('landscape') 
     33figure(1); 
     34orient('landscape') 
     35figure(2); 
     36orient('landscape') 
     37figure(3); 
     38orient('landscape') 
    3539 
    3640% On selectionne la demaine geographique 
    3741 
    3842% la latitude : 7.5N (8) 12.5N (10) 
    39 % la longitude : 10W (9) 10E(17)  
    40 lat1=8; lat2=10; 
    41 lon1=9; lon2=17; 
     43% la longitude : 10W (9) 10E(17) 
     44lat1=8; 
     45lat2=10; 
     46lon1=9; 
     47lon2=17; 
    4248taille=(lat2-lat1+1)*(lon2-lon1+1); 
    4349 
    44  
    45 cosa12=cosa12(:,lon1:lon2,lat1:lat2); clear cosi2; 
    46 cosa12=reshape(cosa12,3416,taille);  
     50cosa12=cosa12(:,lon1:lon2,lat1:lat2); 
     51clear cosi2; 
     52cosa12=reshape(cosa12,3416,taille); 
    4753cosa12=cosa12'; 
    48 cosa12=mean(cosa12);  
     54cosa12=mean(cosa12); 
    4955cosa12=cosa12'; 
    5056save eof12_rossby_1030.txt cosa12 -ascii; 
     
    5258 
    5359% On selectionne la demaine geographique 
    54 olr=reshape(olr,3416,25,17);  
    55 cosa4=olr(:,lon1:lon2,lat1:lat2); clear olr; 
    56 cosa4=reshape(cosa4,3416,taille); clear taille; 
     60olr=reshape(olr,3416,25,17); 
     61cosa4=olr(:,lon1:lon2,lat1:lat2); 
     62clear olr; 
     63cosa4=reshape(cosa4,3416,taille); 
     64clear taille; 
    5765cosa4=cosa4'; 
    5866cosa4=mean(cosa4); 
     
    6270for a=1:28; 
    6371     if a<12; 
    64         figure(1); subplot(4,3,a); 
     72        figure(1); 
     73        subplot(4,3,a); 
    6574     elseif a<24 
    66         figure(2); subplot(4,3,a-11); 
     75        figure(2); 
     76        subplot(4,3,a-11); 
    6777    else 
    68         figure(3); subplot(4,3,a-23); 
     78        figure(3); 
     79        subplot(4,3,a-23); 
    6980     end 
    7081     y2=cosa12(:,a); 
    7182     y4=cosa4(:,a); 
    72      x=[1:122]';  
    73      plot(x,y2); hold on; 
     83     x=[1:122]'; 
     84     plot(x,y2); 
     85     hold on; 
    7486     plot(x,y4,'k'); 
    7587     y5=zeros(122,1); 
    7688     plot(x,y5,'k:'); 
    77       
     89 
    7890     title(a+1978,'Fontsize',8); 
    7991     if a==2; 
     
    97109   set (gca,'XTick',[10:10:122],'XTickLabel',[10:10:122],'fontname','Arial','fontsize',6); 
    98110 
    99 end;  
     111end; 
    100112 
    101 figure(1); print -depsc2 reconstitution1_eof12_filtreER_1030.eps; 
    102 figure(2); print -depsc2 reconstitution2_eof12_filtreER_1030.eps; 
    103 figure(3); print -depsc2 reconstitution3_eof12_filtreER_1030.eps; 
    104  
     113figure(1); 
     114print -depsc2 reconstitution1_eof12_filtreER_1030.eps; 
     115figure(2); 
     116print -depsc2 reconstitution2_eof12_filtreER_1030.eps; 
     117figure(3); 
     118print -depsc2 reconstitution3_eof12_filtreER_1030.eps; 
Note: See TracChangeset for help on using the changeset viewer.