Changeset 75 for trunk/data_out/plot_results.m
- Timestamp:
- 03/02/12 09:38:08 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/data_out/plot_results.m
r38 r75 1 f=netcdf('exp_T .nc');1 f=netcdf('exp_T_3D_all_ebxhz.nc'); 2 2 3 3 addpath('../scripts/matlab_toolbox'); … … 11 11 suff=cell(0,2); 12 12 suff(end+1,:)={'a_c_init','first guess (FG)'}; 13 suff(end+1,:)={'b_forw0_4 8','forward du FG'};13 suff(end+1,:)={'b_forw0_49','forward du FG'}; 14 14 suff(end+1,:)={'a_c_fin','param control final (AP)'}; 15 suff(end+1,:)={'b_forwfin_ 48','forward du PA'};15 suff(end+1,:)={'b_forwfin_99','forward du PA'}; 16 16 suff(end+1,:)={'a_c_true','parametre vrai (TP)'}; 17 suff(end+1,:)={'b_obs_ 48','forward du TP'};17 suff(end+1,:)={'b_obs_99','forward du TP'}; 18 18 19 19 suff2=cell(0,2); 20 20 suff2(end+1,:)={'n_c_init','first guess (FG)'}; 21 suff2(end+1,:)={'b_forw0_4 8','forward du FG'};21 suff2(end+1,:)={'b_forw0_49','forward du FG'}; 22 22 suff2(end+1,:)={'n_c_fin','param control final (AP)'}; 23 suff2(end+1,:)={'b_forwfin_ 48','forward du PA'};23 suff2(end+1,:)={'b_forwfin_99','forward du PA'}; 24 24 suff2(end+1,:)={'n_c_true','parametre vrai (TP)'}; 25 suff2(end+1,:)={'b_obs_ 48','forward du TP'};25 suff2(end+1,:)={'b_obs_99','forward du TP'}; 26 26 27 27 … … 133 133 n2=max(paramy{delta{k}(1)}(:,2)); 134 134 n3=max(paramy{delta{k}(1)}(:,3)); 135 136 137 for iz=1:length(unique(paramy{delta{k}(1)}(:,3))) 135 138 iok=find(paramy{delta{k}(1)}(:,1)>1 & paramy{delta{k}(1)}(:,1)<n1 & paramy{delta{k}(1)}(:,2)>1 & ... 139 paramy{delta{k}(1)}(:,2)<n2 & paramy{delta{k}(1)}(:,3)==iz); 140 141 [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)); 142 end 143 iok=find(paramy{delta{k}(1)}(:,1)>1 & paramy{delta{k}(1)}(:,1)<n1 & paramy{delta{k}(1)}(:,2)>1 & ... 136 144 paramy{delta{k}(1)}(:,2)<n2 & paramy{delta{k}(1)}(:,3)<n3); 137 145 end
Note: See TracChangeset
for help on using the changeset viewer.