Changeset 15 for trunk/src/indice_flore


Ignore:
Timestamp:
12/23/08 10:36:27 (15 years ago)
Author:
pinsard
Message:

modification according to mlint diagnostic

Location:
trunk/src/indice_flore
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/indice_flore/time_serie_histog.m

    r10 r15  
    3939   end; 
    4040   for b=1:122; 
    41       if (b>10 & b<112); % 1ere condition (laisser une marge de 10)  
     41      if (b>10 && b<112); % 1ere condition (laisser une marge de 10)  
    4242%        condition 1 on est a un maximum deltas apres sont negetifs et 
    4343%             delta avant postifs    
    44          if (delta(b+1)<0 & delta(b)<0 & delta(b-1)>0 & delta(b-2)>0 & ind_olr(b,a)>ind_olr2);  
     44         if (delta(b+1)<0 && delta(b)<0 && delta(b-1)>0 && delta(b-2)>0 && ind_olr(b,a)>ind_olr2);  
    4545             
    4646            compteurmax=compteurmax+1; 
    … …  
    6666            type_temps_compomax(compteurmax,20,:)=occu_type_tps(b+9,a,:); 
    6767            type_temps_compomax(compteurmax,21,:)=occu_type_tps(b+10,a,:); 
    68          elseif (delta(b+1)>0 & delta(b)>0 & delta(b-1)<0 & delta(b-2)<0 & ind_olr(b,a)<(-ind_olr2)); 
     68         elseif (delta(b+1)>0 && delta(b)>0 && delta(b-1)<0 && delta(b-2)<0 && ind_olr(b,a)<(-ind_olr2)); 
    6969            compteurmin=compteurmin+1; 
    7070            type_temps_compomin(compteurmin,1,:)=occu_type_tps(b-10,a,:); 
    … …  
    112112 
    113113% on traduit le nbre d occurence en frequences relatives 
     114kara_min=zeros(21,4); 
     115kara_max=zeros(21,4); 
    114116for a=1:21; 
    115117   for b=1:4; 
    … …  
    139141    set (gca,'YTick',[-40:20:40],'YTickLabel',[-40:20:40],'fontname','Arial','fontsize',6); 
    140142    axis([-10 10 -40 40]); 
    141     set(gca,'tickDir','out'); 
     143    set (gca,'tickDir','out'); 
    142144    if a==1; 
    143145       title('Atl low','Fontsize',8,'Fontname','Arial'); 
    … …  
    146148    elseif a==3; 
    147149       title('NAO-','Fontsize',8,'Fontname','Arial'); 
    148     else;  
     150    else 
    149151       title('Blocking','Fontsize',8,'Fontname','Arial'); 
    150152    end;  

  • trunk/src/indice_flore/time_serie_histog3.m

    r10 r15  
    4444   end; 
    4545   for b=1:122; 
    46       if (b>11 & b<112); % 1ere condition (laisser une marge de 10)  
     46      if (b>11 && b<112); % 1ere condition (laisser une marge de 10)  
    4747%        condition 1 on est a un maximum deltas apres sont negetifs et 
    4848%             delta avant postifs    
    49          if (delta(b+1)<0 & delta(b)>0 & ind_olr(b,a)>1.5);  
     49         if (delta(b+1)<0 && delta(b)>0 && ind_olr(b,a)>1.5);  
    5050             
    5151            compteurmax=compteurmax+1; 
    … …  
    117117 
    118118 
    119          elseif (delta(b+1)>0 & delta(b)<0 & ind_olr(b,a)<-1.5); 
     119         elseif (delta(b+1)>0 && delta(b)<0 && ind_olr(b,a)<-1.5); 
    120120            compteurmin=compteurmin+1; 
    121121            type_temps_compomin(1,compteurmin,1,:)=occu_type_tps(b-10,a,:); 
    … …  
    197197type_temps_compomin=reshape(type_temps_compomin,sizemin,21,4); 
    198198type_temps_compomax=reshape(type_temps_compomax,sizemax,21,4); 
    199 compteurmin 
    200 compteurmax 
     199disp(['iii : compteurmin =', int2str(compteurmin)]); 
     200disp(['iii : compteurmax =', int2str(compteurmax)]); 
    201201 
    202202% on fait le test de Student 
    … …  
    226226 
    227227% on traduit le nbre d occurence en frequences relatives 
     228kara_min=zeros(21,4); 
     229kara_max=zeros(21,4); 
    228230for a=1:21; 
    229231   for b=1:4; 
    … …  
    240242        if student1(a,b)==1; 
    241243           table_ttest(a,b)=kara_diff(a,b); 
    242         else; 
     244        else 
    243245           table_ttest(a,b)=NaN; 
    244246        end; 
    … …  
    266268    set (gca,'YTick',[-40:20:40],'YTickLabel',[-40:20:40],'fontname','Arial','fontsize',6); 
    267269    axis([-10 10 -40 40]); 
    268     set(gca,'tickDir','out'); 
     270    set (gca,'tickDir','out'); 
    269271    if a==1; 
    270272       title('Frequence Atl Low (filtre PB) Differences min - max  Seuil 1.5','Fontsize',8,'Fontname','Arial'); 
    … …  
    273275    elseif a==3; 
    274276       title('Frequence NAO- (filtre PB) Differences min - max  Seuil 1.5','Fontsize',8,'Fontname','Arial'); 
    275     else;  
     277    else 
    276278       title('Frequence Blocking (filtre PB) Differences min - max  Seuil 1.5','Fontsize',8,'Fontname','Arial'); 
    277279    end;  

  • trunk/src/indice_flore/time_serie_histog3_max_moy.m

    r10 r15  
    4545   end; 
    4646   for b=1:122; 
    47       if (b>11 & b<112); % 1ere condition (laisser une marge de 10)  
     47      if (b>11 && b<112); % 1ere condition (laisser une marge de 10)  
    4848%        condition 1 on est a un maximum deltas apres sont negetifs et 
    4949%             delta avant postifs    
    50          if (delta(b+1)<0 & delta(b)>0 & ind_olr(b,a)>1.5);  
     50         if (delta(b+1)<0 && delta(b)>0 && ind_olr(b,a)>1.5);  
    5151             
    5252            compteurmax=compteurmax+1; 
    … …  
    118118 
    119119 
    120          elseif (delta(b+1)>0 & delta(b)<0 & ind_olr(b,a)<-1.5); 
     120         elseif (delta(b+1)>0 && delta(b)<0 && ind_olr(b,a)<-1.5); 
    121121            compteurmin=compteurmin+1; 
    122122            type_temps_compomin(1,compteurmin,1,:)=occu_type_tps(b-10,a,:); 
    … …  
    202202type_temps_moy=reshape(occu_type_tps,3538,4); 
    203203 
    204 compteurmin 
    205 compteurmax 
     204disp(['iii : compteurmin =', int2str(compteurmin)]); 
     205disp(['iii : compteurmax =', int2str(compteurmax)]); 
    206206 
    207207student1=zeros(21,4); 
    … …  
    229229 
    230230% on traduit le nbre d occurence en frequences relatives 
     231kara_min=zeros(21,4); 
     232kara_max=zeros(21,4); 
     233kara_moy=zeros(4); 
     234 
    231235for a=1:21; 
    232236   for b=1:4; 
    … …  
    250254        if student1(a,b)==1; 
    251255           table_ttest(a,b)=kara_diff(a,b); 
    252         else; 
     256        else 
    253257           table_ttest(a,b)=NaN; 
    254258        end; 
    … …  
    277281    set (gca,'YTick',[-20:5:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); 
    278282    axis([-10 10 -20 20]); 
    279     set(gca,'tickDir','out'); 
     283    set (gca,'tickDir','out'); 
    280284    if a==1; 
    281285       title('Atl low','Fontsize',8,'Fontname','Arial'); 
    … …  
    284288    elseif a==3; 
    285289       title('NAO-','Fontsize',8,'Fontname','Arial'); 
    286     else;  
     290    else 
    287291       title('Blocking','Fontsize',8,'Fontname','Arial'); 
    288292    end;  

  • trunk/src/indice_flore/time_serie_histog3_min_moy.m

    r10 r15  
    4646   end; 
    4747   for b=1:122; 
    48       if (b>11 & b<112); % 1ere condition (laisser une marge de 10)  
     48      if (b>11 && b<112); % 1ere condition (laisser une marge de 10)  
    4949%        condition 1 on est a un maximum deltas apres sont negetifs et 
    5050%             delta avant postifs    
    51          if (delta(b+1)<0 & delta(b)>0 & ind_olr(b,a)>1.5);  
     51         if (delta(b+1)<0 && delta(b)>0 && ind_olr(b,a)>1.5);  
    5252             
    5353            compteurmax=compteurmax+1; 
    … …  
    119119 
    120120 
    121          elseif (delta(b+1)>0 & delta(b)<0 & ind_olr(b,a)<-1.5); 
     121         elseif (delta(b+1)>0 && delta(b)<0 && ind_olr(b,a)<-1.5); 
    122122            compteurmin=compteurmin+1; 
    123123            type_temps_compomin(1,compteurmin,1,:)=occu_type_tps(b-10,a,:); 
    … …  
    195195type_temps_compomax=reshape(type_temps_compomax,240,21,4); 
    196196type_temps_moy=reshape(occu_type_tps,3538,4); 
    197 compteurmin 
    198 compteurmax 
     197disp(['iii : compteurmin =', int2str(compteurmin)]); 
     198disp(['iii : compteurmax =', int2str(compteurmax)]); 
    199199% on fait la somme des types de temps sur toute la periode de temps 
    200200% on obtent un tableau de la forme (jours-avant-apres,cluster) 
    … …  
    211211 
    212212% on traduit le nbre d occurence en frequences relatives 
     213kara_min=zeros(21,4); 
     214kara_max=zeros(21,4); 
     215kara_moy=zeros(4); 
    213216for a=1:21; 
    214217   for b=1:4; 
    … …  
    244247    set (gca,'YTick',[-20:5:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6); 
    245248    axis([-10 10 -20 20]); 
    246     set(gca,'tickDir','out'); 
     249    set (gca,'tickDir','out'); 
    247250    if a==1; 
    248251       title('Atl low','Fontsize',8,'Fontname','Arial'); 
    … …  
    251254    elseif a==3; 
    252255       title('NAO-','Fontsize',8,'Fontname','Arial'); 
    253     else;  
     256    else 
    254257       title('Blocking','Fontsize',8,'Fontname','Arial'); 
    255258    end;  
Note: See TracChangeset for help on using the changeset viewer.