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/composite_pluie_ird_rossby1030.m

    r15 r16  
     1% 
     2 
    13%%%%%%%%%%%%%%% FICHIER pluie IRD 
    24% ouverture / lecture du fichier 
    … …  
    1921       for c=1:365; 
    2022          if v(a,b,c)>990; 
    21               v(a,b,c)=NaN;         
     23              v(a,b,c)=NaN; 
    2224          elseif v(a,b,c)<0; 
    2325              v(a,b,c)=NaN; 
    … …  
    3638       for c=1:366; 
    3739          if v(a,b,c)>990; 
    38               v(a,b,c)=NaN;         
     40              v(a,b,c)=NaN; 
    3941          elseif v(a,b,c)<0; 
    4042              v(a,b,c)=NaN; 
    … …  
    5355       for c=1:365; 
    5456          if v(a,b,c)>990; 
    55               v(a,b,c)=NaN;         
     57              v(a,b,c)=NaN; 
    5658          elseif v(a,b,c)<0; 
    5759              v(a,b,c)=NaN; 
    … …  
    7072       for c=1:365; 
    7173          if v(a,b,c)>990; 
    72               v(a,b,c)=NaN;         
     74              v(a,b,c)=NaN; 
    7375          elseif v(a,b,c)<0; 
    7476              v(a,b,c)=NaN; 
    … …  
    8789       for c=1:365; 
    8890          if v(a,b,c)>990; 
    89               v(a,b,c)=NaN;         
     91              v(a,b,c)=NaN; 
    9092          elseif v(a,b,c)<0; 
    9193              v(a,b,c)=NaN; 
    … …  
    106108       for c=1:366; 
    107109          if v(a,b,c)>990; 
    108               v(a,b,c)=NaN;         
     110              v(a,b,c)=NaN; 
    109111          elseif v(a,b,c)<0; 
    110112              v(a,b,c)=NaN; 
    … …  
    123125       for c=1:365; 
    124126          if v(a,b,c)>990; 
    125               v(a,b,c)=NaN;         
     127              v(a,b,c)=NaN; 
    126128          elseif v(a,b,c)<0; 
    127129              v(a,b,c)=NaN; 
    … …  
    140142       for c=1:365; 
    141143          if v(a,b,c)>990; 
    142               v(a,b,c)=NaN;         
     144              v(a,b,c)=NaN; 
    143145          elseif v(a,b,c)<0; 
    144146              v(a,b,c)=NaN; 
    … …  
    157159       for c=1:365; 
    158160          if v(a,b,c)>990; 
    159               v(a,b,c)=NaN;         
     161              v(a,b,c)=NaN; 
    160162          elseif v(a,b,c)<0; 
    161163              v(a,b,c)=NaN; 
    … …  
    174176       for c=1:366; 
    175177          if v(a,b,c)>990; 
    176               v(a,b,c)=NaN;         
     178              v(a,b,c)=NaN; 
    177179          elseif v(a,b,c)<0; 
    178180              v(a,b,c)=NaN; 
    … …  
    191193       for c=1:365; 
    192194          if v(a,b,c)>990; 
    193               v(a,b,c)=NaN;         
     195              v(a,b,c)=NaN; 
    194196          elseif v(a,b,c)<0; 
    195197              v(a,b,c)=NaN; 
    … …  
    208210       for c=1:365; 
    209211          if v(a,b,c)>990; 
    210               v(a,b,c)=NaN;         
     212              v(a,b,c)=NaN; 
    211213          elseif v(a,b,c)<0; 
    212214              v(a,b,c)=NaN; 
    … …  
    223225clear 
    224226 
    225 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
    226227load pluie.mat; 
    227228% pluie dimension : lon lat jours annee (41 17 122 12) => jour annee lat lon  (17 41 122 12) 
    228229 
    229 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
    230230load eof12_rossby_1030.txt; 
    231231 
    232232% on met sous la forme annee mois 
    233233eof12_rossby_1030=reshape(eof12_rossby_1030,122,28); 
    234 olr_eof=eof12_rossby_1030(:,1:12); % selection des annees avant 1990  
     234% selection des annees avant 1990 
     235olr_eof=eof12_rossby_1030(:,1:12); 
    235236clear eof12_rossby_1030; 
    236237 
    237238% on calcule l ecart type sur la serie temporelle de l'olr 
    238 k=122*12;  
     239k=122*12; 
    239240olrbis=reshape(olr_eof,k,1); 
    240 olrbis=std(olrbis);  
     241olrbis=std(olrbis); 
    241242 
    242243% intialisation des tableaux 
    243 pluie_filtre_compomax=NaN*ones(30,31,17,41);  
     244pluie_filtre_compomax=NaN*ones(30,31,17,41); 
    244245pluie_filtre_compomin=NaN*ones(30,31,17,41); 
    245246pluie=permute(pluie,[3 4 2 1]); 
    246 compteurmax=0;  
     247compteurmax=0; 
    247248compteurmin=0; 
    248249 
    249 for a=1:12; % boucle sur les annees 
    250    clear delta;  
     250% boucle sur les annees 
     251for a=1:12; 
     252   clear delta; 
    251253   delta=zeros(122,1); 
    252    for b=1:122; % boucle sur les jours : on prepare avant    
     254   % boucle sur les jours : on prepare avant 
     255   for b=1:122; 
    253256        if b>1; 
    254257            delta(b)=olr_eof(b,a)-olr_eof(b-1,a); 
    255258        end; 
    256259   end; 
    257    for b=1:122; % boucle sur les jours 
    258       if (b>20 && b<112); % 1ere condition (laisser une marge de 10)  
    259 % cas 1 : max deltas apres negetifs et deltas avant postifs    
    260          if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && olr_eof(b,a)>olrbis);  
     260   % boucle sur les jours 
     261   for b=1:122; 
     262      % 1ere condition (laisser une marge de 10) 
     263      if (b>20 && b<112); 
     264% cas 1 : max deltas apres negetifs et deltas avant postifs 
     265         if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && olr_eof(b,a)>olrbis); 
    261266            compteurmax=compteurmax+1; 
    262267            pluie_filtre_compomax(compteurmax,1,:,:)=pluie(b-20,a,:,:); 
    … …  
    291296            pluie_filtre_compomax(compteurmax,30,:,:)=pluie(b+9,a,:,:); 
    292297            pluie_filtre_compomax(compteurmax,31,:,:)=pluie(b+10,a,:,:); 
    293 % cas 2 : min deltas apres positifs et deltas avant negatifs    
     298% cas 2 : min deltas apres positifs et deltas avant negatifs 
    294299         elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && olr_eof(b,a)<(-olrbis)); 
    295300            compteurmin=compteurmin+1; 
    … …  
    328333         end; 
    329334      end; 
    330    end;          
     335   end; 
    331336end; 
    332337 
    … …  
    334339% on fait la somme des valeurs de pluie_filtreentiel sur le nbre d occurences 
    335340% on obtent un tableau de la forme (jours-avant-apres,lon,lat) 
    336 pluie_filtre_compomin1(:,:,:)=nanmean(pluie_filtre_compomin); clear pluie_filtre_compomin; 
    337 pluie_filtre_compomax1(:,:,:)=nanmean(pluie_filtre_compomax); clear pluie_filtre_compomax; 
    338  
    339 figure(1); orient('landscape'); 
    340 figure(2); orient('landscape'); 
    341 figure(3); orient('landscape'); 
     341pluie_filtre_compomin1(:,:,:)=nanmean(pluie_filtre_compomin); 
     342clear pluie_filtre_compomin; 
     343pluie_filtre_compomax1(:,:,:)=nanmean(pluie_filtre_compomax); 
     344clear pluie_filtre_compomax; 
     345 
     346figure(1); 
     347orient('landscape'); 
     348figure(2); 
     349orient('landscape'); 
     350figure(3); 
     351orient('landscape'); 
    342352 
    343353for e=1:7; 
    … …  
    368378       pluie_max(1,1)=-0.0001; 
    369379       pluie_max(1,2)=10.0001; 
    370        lat=[-10:2.5:30]; lon=[-70:2.5:30]'; 
    371        figure(1); subplot(7,1,8-e); palette 
     380       lat=[-10:2.5:30]; 
     381       lon=[-70:2.5:30]'; 
     382       figure(1); 
     383       subplot(7,1,8-e); 
     384       palette 
    372385       aa=contourf(lon,lat,pluie_min,[0:0.5:10],'LineStyle','none'); 
    373386       if e==1; 
    … …  
    375388       else 
    376389         set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 
    377        end;        
     390       end; 
    378391       set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6); 
    379392       if e==4; 
    380           co=colorbar;             
     393          co=colorbar; 
    381394          set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 
    382395       end; 
    383        hold on; cartemonde1; 
     396       hold on; 
     397       cartemonde1; 
    384398       set (gca,'tickDir','out'); 
    385399       axis equal; 
    … …  
    390404       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 
    391405       pos_vert=e/10; 
    392        set (gca,'position',[0.2 pos_vert 0.5 0.1]);  
     406       set (gca,'position',[0.2 pos_vert 0.5 0.1]); 
    393407%         *************************** 
    394        figure(2); subplot(7,1,8-e); palette 
     408       figure(2); 
     409       subplot(7,1,8-e); 
     410       palette 
    395411       aa=contourf(lon,lat,pluie_max,[0:0.5:10],'LineStyle','none'); 
    396412       if e==1; 
    … …  
    398414       else 
    399415         set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 
    400        end;        
     416       end; 
    401417       set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6); 
    402418       if e==4; 
    403           co=colorbar;             
     419          co=colorbar; 
    404420          set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 
    405421       end; 
    406        hold on; cartemonde1; 
     422       hold on; 
     423       cartemonde1; 
    407424       set (gca,'tickDir','out'); 
    408425       axis equal; 
    … …  
    413430       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 
    414431       pos_vert=e/10; 
    415        set (gca,'position',[0.2 pos_vert 0.5 0.1]);  
     432       set (gca,'position',[0.2 pos_vert 0.5 0.1]); 
    416433%******************************************** 
    417434 
    418        figure(3); subplot(7,1,8-e); palette 
     435       figure(3); 
     436       subplot(7,1,8-e); 
     437       palette 
    419438       diff=pluie_min-pluie_max; 
    420439       diff(1,1)=-5.01; 
    … …  
    425444       else 
    426445         set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 
    427        end;        
     446       end; 
    428447       set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6); 
    429448       if e==4; 
    430           co=colorbar;             
     449          co=colorbar; 
    431450          set (co,'xtick',[-20:10:20],'xticklabel',[-20:10:20],'fontname','Arial','fontsize',6); 
    432451       end; 
    433        hold on; cartemonde1; 
     452       hold on; 
     453       cartemonde1; 
    434454       set (gca,'tickDir','out'); 
    435455       axis equal; 
    … …  
    440460       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 
    441461       pos_vert=e/10; 
    442        set (gca,'position',[0.2 pos_vert 0.5 0.1]);  
    443  
    444  
    445  
    446 end; 
    447 figure(1); print -depsc2 composite_pluie_min_eof2a.eps; 
    448 figure(2); print -depsc2 composite_pluie_max_eof2a.eps; 
    449 figure(3); print -depsc2 composite_pluie_diff_eof2a.eps; 
    450 close all 
    451  
    452 %******************************* 
    453  
    454  
    455 figure(1); orient('landscape'); 
    456 figure(2); orient('landscape'); 
    457 figure(3); orient('landscape'); 
     462       set (gca,'position',[0.2 pos_vert 0.5 0.1]); 
     463 
     464 
     465 
     466end; 
     467figure(1); 
     468print -depsc2 composite_pluie_min_eof2a.eps; 
     469figure(2); 
     470print -depsc2 composite_pluie_max_eof2a.eps; 
     471figure(3); 
     472print -depsc2 composite_pluie_diff_eof2a.eps; 
     473close all; 
     474 
     475figure(1); 
     476orient('landscape'); 
     477figure(2); 
     478orient('landscape'); 
     479figure(3); 
     480orient('landscape'); 
    458481 
    459482for e=1:7; 
    … …  
    484507       pluie_max(1,1)=-0.0001; 
    485508       pluie_max(1,2)=10.0001; 
    486        lat=[-10:2.5:30]; lon=[-70:2.5:30]'; 
    487        figure(1); subplot(7,1,8-e); palette 
     509       lat=[-10:2.5:30]; 
     510       lon=[-70:2.5:30]'; 
     511       figure(1); 
     512       subplot(7,1,8-e); 
     513       palette 
    488514       aa=contourf(lon,lat,pluie_min,[0:0.5:10],'LineStyle','none'); 
    489515       if e==1; 
    … …  
    491517       else 
    492518         set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 
    493        end;        
     519       end; 
    494520       set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6); 
    495521       if e==4; 
    496           co=colorbar;             
     522          co=colorbar; 
    497523          set (co,'xtick',[0:2:10],'xticklabel',[0:2:10],'fontname','Arial','fontsize',6); 
    498524       end; 
    499        hold on; cartemonde1; 
     525       hold on; 
     526       cartemonde1; 
    500527       set (gca,'tickDir','out'); 
    501528       axis equal; 
    … …  
    506533       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 
    507534       pos_vert=e/10; 
    508        set (gca,'position',[0.2 pos_vert 0.5 0.1]);  
     535       set (gca,'position',[0.2 pos_vert 0.5 0.1]); 
    509536%         *************************** 
    510        figure(2); subplot(7,1,8-e); palette 
     537       figure(2); 
     538       subplot(7,1,8-e); 
     539       palette 
    511540       aa=contourf(lon,lat,pluie_max,[0:0.5:10],'LineStyle','none'); 
    512541       if e==1; 
    … …  
    514543       else 
    515544         set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 
    516        end;        
     545       end; 
    517546       set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6); 
    518547       if e==4; 
    519           co=colorbar;             
     548          co=colorbar; 
    520549          set (co,'xtick',[0:2:10],'xticklabel',[0:2:10],'fontname','Arial','fontsize',6); 
    521550       end; 
    522        hold on; cartemonde1; 
     551       hold on; 
     552       cartemonde1; 
    523553       set (gca,'tickDir','out'); 
    524554       axis equal; 
    … …  
    529559       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 
    530560       pos_vert=e/10; 
    531        set (gca,'position',[0.2 pos_vert 0.5 0.1]);  
     561       set (gca,'position',[0.2 pos_vert 0.5 0.1]); 
    532562%******************************************** 
    533563 
    534        figure(3); subplot(7,1,8-e); palette 
     564       figure(3); 
     565       subplot(7,1,8-e); 
     566       palette 
    535567       diff=pluie_min-pluie_max; 
    536568       diff(1,1)=-5.01; 
    … …  
    541573       else 
    542574         set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 
    543        end;        
     575       end; 
    544576       set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6); 
    545577       if e==4; 
    546           co=colorbar;             
     578          co=colorbar; 
    547579          set (co,'xtick',[-5:1:5],'xticklabel',[-5:1:5],'fontname','Arial','fontsize',6); 
    548580       end; 
    549        hold on; cartemonde1; 
     581       hold on; 
     582       cartemonde1; 
    550583       set (gca,'tickDir','out'); 
    551584       axis equal; 
    … …  
    556589       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 
    557590       pos_vert=e/10; 
    558        set (gca,'position',[0.2 pos_vert 0.5 0.1]);  
    559  
    560  
    561  
    562 end; 
    563 figure(1); print -depsc2 composite_pluie_min_eof2b.eps; 
    564 figure(2); print -depsc2 composite_pluie_max_eof2b.eps; 
    565 figure(3); print -depsc2 composite_pluie_diff_eof2b.eps; 
    566  
    567  
    568 %******************************* 
    569  
     591       set (gca,'position',[0.2 pos_vert 0.5 0.1]); 
     592end; 
     593figure(1); 
     594print -depsc2 composite_pluie_min_eof2b.eps; 
     595figure(2); 
     596print -depsc2 composite_pluie_max_eof2b.eps; 
     597figure(3); 
     598print -depsc2 composite_pluie_diff_eof2b.eps; 
Note: See TracChangeset for help on using the changeset viewer.