Changeset 75 for trunk/data_out/plot_results.m

Timestamp:

03/02/12 09:38:08 (12 years ago)

Author:

jbrlod

Message:

new option for xchangmode

File:

• trunk/data_out/plot_results.m (modified) (3 diffs)

trunk/data_out/plot_results.m

@@ -1 @@
-f=netcdf('exp_T.nc');
+f=netcdf('exp_T_3D_all_ebxhz.nc');
 
 addpath('../scripts/matlab_toolbox');
@@ -11 +11 @@
 suff=cell(0,2);
 suff(end+1,:)={'a_c_init','first guess (FG)'};
-suff(end+1,:)={'b_forw0_48','forward du FG'};
+suff(end+1,:)={'b_forw0_49','forward du FG'};
 suff(end+1,:)={'a_c_fin','param control final (AP)'};
-suff(end+1,:)={'b_forwfin_48','forward du PA'};
+suff(end+1,:)={'b_forwfin_99','forward du PA'};
 suff(end+1,:)={'a_c_true','parametre vrai (TP)'};
-suff(end+1,:)={'b_obs_48','forward du TP'};
+suff(end+1,:)={'b_obs_99','forward du TP'};
 
 suff2=cell(0,2);
 suff2(end+1,:)={'n_c_init','first guess (FG)'};
-suff2(end+1,:)={'b_forw0_48','forward du FG'};
+suff2(end+1,:)={'b_forw0_49','forward du FG'};
 suff2(end+1,:)={'n_c_fin','param control final (AP)'};
-suff2(end+1,:)={'b_forwfin_48','forward du PA'};
+suff2(end+1,:)={'b_forwfin_99','forward du PA'};
 suff2(end+1,:)={'n_c_true','parametre vrai (TP)'};
-suff2(end+1,:)={'b_obs_48','forward du TP'};
+suff2(end+1,:)={'b_obs_99','forward du TP'};
 
 
@@ -133 +133 @@
         n2=max(paramy{delta{k}(1)}(:,2));
         n3=max(paramy{delta{k}(1)}(:,3));
+       
+      
+      for iz=1:length(unique(paramy{delta{k}(1)}(:,3)))
         iok=find(paramy{delta{k}(1)}(:,1)>1 & paramy{delta{k}(1)}(:,1)<n1 & paramy{delta{k}(1)}(:,2)>1 & ...
+                 paramy{delta{k}(1)}(:,2)<n2 & paramy{delta{k}(1)}(:,3)==iz);
+        
+ [Nz(j,k,iz),RPDz(j,k,iz),URPDz(j,k,iz),MADz(j,k,iz),RMSz(j,k,iz),slopez(j,k,iz),interceptz(j,k,iz),r2z(j,k,iz) ]=make_stats(paramy{delta{k,1}(1)}(iok,end),paramy{delta{k,1}(2)} (iok,end));
+      end 
+       iok=find(paramy{delta{k}(1)}(:,1)>1 & paramy{delta{k}(1)}(:,1)<n1 & paramy{delta{k}(1)}(:,2)>1 & ...
                  paramy{delta{k}(1)}(:,2)<n2 & paramy{delta{k}(1)}(:,3)<n3);
       end