Ignore:
File:
1 edited

Legend:

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

    r26 r18  
    55close all; 
    66 
     7cd('..'); 
    78ncload('olr.day.mean.nc'); 
    8  
     9cd('mode_sahelien'); 
    910% taille 12327x73x144 1er juin 1974-31 decembre2007 
    1011% on enleve les annees 1974 a 1978 nombre de jours a enlever au debut 214+365+366+365+365 
     
    4142      for c=1:siz; 
    4243           vecteur(c)=olr1(a,b,c); 
    43       end; 
     44      end 
    4445      std_c=nanmean(vecteur); 
    4546      for c=1:siz; 
    4647          olr1(a,b,c)=olr1(a,b,c)-std_c; 
    47       end; 
    48    end; 
     48      end 
     49   end 
    4950end; 
    5051 
     
    145146clear olr; 
    146147 
    147 status=load('eof24_1025.txt'); 
    148 eof24_1025=status; 
    149 clear status; 
     148%load ind_olr_10_25_eof234.txt 
     149load eof24_1025.txt 
    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.