Ignore:
File:
1 edited

Legend:

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

    r18 r26  
    55close all; 
    66 
    7 cd('..'); 
    87ncload('olr.day.mean.nc'); 
    9 cd('mode_sahelien'); 
     8 
    109% taille 12327x73x144 1er juin 1974-31 decembre2007 
    1110% on enleve les annees 1974 a 1978 nombre de jours a enlever au debut 214+365+366+365+365 
     
    4241      for c=1:siz; 
    4342           vecteur(c)=olr1(a,b,c); 
    44       end 
     43      end; 
    4544      std_c=nanmean(vecteur); 
    4645      for c=1:siz; 
    4746          olr1(a,b,c)=olr1(a,b,c)-std_c; 
    48       end 
    49    end 
     47      end; 
     48   end; 
    5049end; 
    5150 
     
    146145clear olr; 
    147146 
    148 %load ind_olr_10_25_eof234.txt 
    149 load eof24_1025.txt 
     147status=load('eof24_1025.txt'); 
     148eof24_1025=status; 
     149clear status; 
    150150% on met sous la forme annee mois 
    151151olr_eof=reshape(eof24_1025,122,28); 
     
    264264           y=olr_filtre_compomax(1:nb_value,a,b,c); 
    265265           student1(a,b,c)=ttest(x,y,.2); 
    266         end 
    267     end 
    268 end 
     266        end; 
     267    end; 
     268end; 
    269269 
    270270% on fait la somme des valeurs de olr_filtreentiel sur le nbre d occurences 
     
    305305             end; 
    306306          end; 
    307        end 
     307       end; 
    308308       subplot(7,1,8-e); 
    309309       palette; 
     
    319319       else 
    320320      set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 
    321        end 
     321       end; 
    322322       set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6); 
    323323       if e==4; 
     
    332332       if e==7; 
    333333           title('composites OLR filtre 10-25 jours EOF234','fontsize',8) ; 
    334        end 
     334       end; 
    335335       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 
    336336       pos_vert=e/10; 
     
    339339 
    340340figure(1); 
    341 print -depsc2 composite_olr_jjas_sahel_10_25_eof234_part1.eps; 
     341print('-depsc2','composite_olr_jjas_sahel_10_25_eof234_part1.eps'); 
    342342 
    343343figure(2); 
     
    367367             end; 
    368368          end; 
    369        end 
     369       end; 
    370370       subplot(7,1,8-e); 
    371371       palette 
     
    381381       else 
    382382      set (gca,'XTick',[-180:60:180],'XTickLabel',[' ';' ';' ';' ';' ';' '],'fontname','Arial','fontsize',6); 
    383        end 
     383       end; 
    384384       set (gca,'YTick',[-10:10:40],'YTickLabel',['   ';'  0';'10N';'20N';'30N';'   '],'fontname','Arial','fontsize',6); 
    385385       if e==4; 
     
    394394       if e==7; 
    395395           title('composites OLR filtre 10-25 jours EOF234','fontsize',8) ; 
    396        end 
     396       end; 
    397397       text(-270,10,['day = ',num2str(f-21)],'fontname','Arial','fontsize',10); 
    398398       pos_vert=e/10; 
     
    401401 
    402402figure(2); 
    403 print -depsc2 composite_olr_jjas_sahel_10_25_eof234_part2.eps; 
     403print('-depsc2','composite_olr_jjas_sahel_10_25_eof234_part2.eps'); 
Note: See TracChangeset for help on using the changeset viewer.