Changeset 24 for trunk/src/indice_flore

Timestamp:

01/07/09 16:49:58 (15 years ago)

Author:

pinsard

Message:

correction of usage of print function, fix side effect of load function

Location:

trunk/src/indice_flore

Files:

• time_serie_histog.m (modified) (4 diffs)
• time_serie_histog3.m (modified) (8 diffs)
• time_serie_histog3_max_moy.m (modified) (8 diffs)
• time_serie_histog3_min_moy.m (modified) (6 diffs)

trunk/src/indice_flore/time_serie_histog.m

@@ -27 +27 @@
 % on ouvre le fichier
 status=load('occu.txt');
+occu=status;
+clear status;
 
 % on reorganise le fichier de type de temps (jours,annees,cluster)
@@ -147 +149 @@
     if a==1;
        bar(jour,diff,'r');
-c=[1 0.55 0];
-bar(jour,diff),colormap(c);
+       c=[1 0.55 0];
+       bar(jour,diff),colormap(c);
     elseif a==2;
        bar(jour,diff,'g');
@@ -155 +157 @@
     else
        bar(jour,diff,'r');
-    end
+    end;
     set (gca,'XTick',[-10:10],'XTickLabel',[-10:10],'fontname','Arial','fontsize',6);
     set (gca,'YTick',[-40:20:40],'YTickLabel',[-40:20:40],'fontname','Arial','fontsize',6);
@@ -169 +171 @@
        title('Blocking','Fontsize',8,'Fontname','Arial');
     end;
-end
+end;
 
 figure(1);
-status=print('-depsc2','time_serie_type_tps_pb.eps');
+print('-depsc2','time_serie_type_tps_pb.eps');

trunk/src/indice_flore/time_serie_histog3.m

@@ -21 +21 @@
 
 for a=1:3538;
-   ind_olr(a)=(ind_olr(a)-olr_moy)/olr_std;
+    ind_olr(a)=(ind_olr(a)-olr_moy)/olr_std;
 end;
 ind_olr=reshape(ind_olr,122,29);
@@ -28 +28 @@
 % on ouvre le fichier
 status=load('occu.txt');
+occu=status;
+clear status;
 
 % on reorganise le fichier de type de temps (jours,annees,cluster)
@@ -44 +46 @@
 
 % boucle sur les annees
-for a=1:29
+for a=1:29;
    clear delta;
    delta=zeros(122,1);
    % boucle sur les jours
-   for b=1:122
+   for b=1:122;
         if b>1;
             delta(b)=ind_olr(b,a)-ind_olr(b-1,a);
@@ -55 +57 @@
    for b=1:122;
       % 1ere condition (laisser une marge de 10)
-      if (b>11 && b<112);
+      if (b>11 && b<112)
 %        condition 1 on est a un maximum deltas apres sont negetifs et
 %             delta avant postifs
@@ -218 +220 @@
        y=type_temps_compomax(1:sizemax,a,b);
        student1(a,b)=ttest(x,y);
-    end
-end
-
-
+    end;
+end;
 
 % on fait la somme des types de temps sur toute la periode de temps
@@ -256 +256 @@
            table_ttest(a,b)=NaN;
         end;
-    end
-end
+    end;
+end;
 for a=1:4;
     subplot(2,2,a);
@@ -280 +280 @@
        hold on;
        plot(jour,diff2,'k*');
-    end
+    end;
     set (gca,'XTick',[-10:10],'XTickLabel',[-10:10],'fontname','Arial','fontsize',6);
     set (gca,'YTick',[-40:20:40],'YTickLabel',[-40:20:40],'fontname','Arial','fontsize',6);
@@ -294 +294 @@
        title('Frequence Blocking (filtre PB) Differences min - max  Seuil 1.5','Fontsize',8,'Fontname','Arial');
     end;
-end
+end;
 
 figure(1);
-status=print('-depsc2','time_serie_type_tps_pb.eps');
+print('-depsc2','time_serie_type_tps_pb.eps');

trunk/src/indice_flore/time_serie_histog3_max_moy.m

@@ -24 +24 @@
 
 for a=1:3538;
-   ind_olr(a)=(ind_olr(a)-olr_moy)/olr_std;
+    ind_olr(a)=(ind_olr(a)-olr_moy)/olr_std;
 end;
 ind_olr=reshape(ind_olr,122,29);
+
 %%%%%%%%%%%% FICHIER TYPE DE TEMPS
 % on ouvre le fichier
 status=load('occu.txt');
+occu=status;
+clear status;
 
 % on reorganise le fichier de type de temps (jours,annees,cluster)
@@ -46 +49 @@
 
 % boucle sur les annees
-for a=1:29
+for a=1:29;
    clear delta;
    delta=zeros(122,1);
    % boucle sur les jours
-   for b=1:122
+   for b=1:122;
         if b>1;
             delta(b)=ind_olr(b,a)-ind_olr(b-1,a);
@@ -223 +226 @@
        y=type_temps_compomax(1:sizemax,a,b);
        student1(a,b)=ttest(x,y);
-    end
-end
-
+    end;
+end;
 
 % on fait la somme des types de temps sur toute la periode de temps
@@ -258 +260 @@
     for b=1:4
        kara_diff(a,b)=-kara_max(a,b)+kara_moy(b);
-    end
-end
+    end;
+end;
 
 table_ttest=zeros(21,4);
@@ -269 +271 @@
            table_ttest(a,b)=NaN;
         end;
-    end
-end
+    end;
+end;
 
 for a=1:4;
@@ -282 +284 @@
        hold on;
        plot(jour,diff2,'k*');
-     elseif a==2;
+    elseif a==2;
        bar(jour,diff,'g');
-        hold on;
-        plot(jour,diff2,'k*');
+       hold on;
+       plot(jour,diff2,'k*');
     elseif a==3;
        bar(jour,diff,'b');
@@ -294 +296 @@
        hold on;
        plot(jour,diff2,'k*');
-    end
+    end;
     set (gca,'XTick',[-10:10],'XTickLabel',[-10:10],'fontname','Arial','fontsize',6);
     set (gca,'YTick',[-20:5:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6);
@@ -308 +310 @@
        title('Blocking','Fontsize',8,'Fontname','Arial');
     end;
-end
+end;
 
 figure(1);
-status=print('-depsc2','time_serie_type_tps_max_moy_pb.eps');
+print('-depsc2','time_serie_type_tps_max_moy_pb.eps');

trunk/src/indice_flore/time_serie_histog3_min_moy.m

@@ -25 +25 @@
 for a=1:3538;
     ind_olr(a)=(ind_olr(a)-olr_moy)/olr_std;
-end
+end;
 ind_olr=reshape(ind_olr,122,29);
 
@@ -31 +31 @@
 % on ouvre le fichier
 status=load('occu.txt');
+occu=status;
+clear status;
 
 % on reorganise le fichier de type de temps (jours,annees,cluster)
@@ -58 +60 @@
    for b=1:122;
       % 1ere condition (laisser une marge de 10)
-      if (b>11 && b<112);
+      if (b>11 && b<112)
 %        condition 1 on est a un maximum deltas apres sont negetifs et
 %             delta avant postifs
@@ -239 +241 @@
     for b=1:4
        kara_diff(a,b)=kara_min(a,b)-kara_moy(b);
-    end
-end
+    end;
+end;
 
 for a=1:4;
@@ -256 +258 @@
     else
        bar(jour,diff,'r');
-    end
+    end;
     set (gca,'XTick',[-10:10],'XTickLabel',[-10:10],'fontname','Arial','fontsize',6);
     set (gca,'YTick',[-20:5:20],'YTickLabel',[-20:5:20],'fontname','Arial','fontsize',6);
@@ -270 +272 @@
        title('Blocking','Fontsize',8,'Fontname','Arial');
     end;
-end
+end;
 
 figure(1);
-status=print('-depsc2','time_serie_type_tps_min_moy_pb.eps');
+print('-depsc2','time_serie_type_tps_min_moy_pb.eps');