Changeset 16 for trunk/data_out/plot_results.m
- Timestamp:
- 03/22/11 10:35:04 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/data_out/plot_results.m
r15 r16 1 f=netcdf('exp_T.nc'); 2 3 addpath('../scripts/matlab_toolbox'); 1 4 sparam=cell(0,5); 2 5 sparam(end+1,:)={'t',[12:1:24],[0:2:24],[-1:0.2:1],[-1:0.2:1]}; … … 4 7 sparam(end+1,:)={'u',[-0.2:0.02:0.2],[-0.2:0.015:0.2],[-0.14:0.02:0.14],[-0.14:0.02:0.14]}; 5 8 sparam(end+1,:)={'v',[-0.2:0.02:0.2],[-0.2:0.015:0.2],[-0.14:0.02:0.14],[-0.14:0.02:0.14]}; 6 f=netcdf('exp_T.nc');9 sparam(end+1,:)={'ssh',[-0.8:0.05:0.2],[-0.2:0.015:0.2],[-0.14:0.02:0.14],[-0.14:0.02:0.14]}; 7 10 8 11 suff=cell(0,2); … … 13 16 suff(end+1,:)={'a_c_true','parametre vrai (TP)'}; 14 17 suff(end+1,:)={'b_obs_48','forward du TP'}; 18 19 suff2=cell(0,2); 20 suff2(end+1,:)={'n_c_init','first guess (FG)'}; 21 suff2(end+1,:)={'b_forw0_48','forward du FG'}; 22 suff2(end+1,:)={'n_c_fin','param control final (AP)'}; 23 suff2(end+1,:)={'b_forwfin_48','forward du PA'}; 24 suff2(end+1,:)={'n_c_true','parametre vrai (TP)'}; 25 suff2(end+1,:)={'b_obs_48','forward du TP'}; 26 15 27 16 28 delta=cell(0,2); … … 38 50 for k=1:length(suff) 39 51 52 if strcmp(sparam{j,1},'ssh') 53 param=f{[sparam{j,1} suff2{k,1}]}(:,:,:); 54 55 param=squeeze(param); 56 paramy{k}=mat2yao(param); 57 paramy{k}=mat2yao(param); 58 paramy{k}(:,1:2)=fliplr(paramy{k}(:,1:2)); 59 paramy{k}=sortrows(paramy{k},[1 2]); 60 figure(j); 61 subplot(3,2,k) 62 % contourf(longi(2:end-1,2:end-1),lati(2:end-1,2:end-1),param(2:end-1,2:end-1),sparam{j,2}); 63 % caxis([sparam{j,2}(1) sparam{j,2}(end)]) 64 %xlabel('longitude','Fontsize',9) 65 %ylabel('latitude','Fontsize',9) 66 plot_nemo_3D(paramy{k},longi,lati,zlev,1,sparam{j,2},1); 67 68 69 title(suff2{k,2}); 70 if (k==1) 71 colorbar 72 end 40 73 41 param=f{[sparam{j,1} suff{k,1}]}(:,:,:,:);42 param=squeeze(param);43 paramy{k}=mat2yao(param);44 paramy{k}(:,1:3)=fliplr(paramy{k}(:,1:3));45 paramy{k}=sortrows(paramy{k},[1 2 3]);46 74 else 75 param=f{[sparam{j,1} suff{k,1}]}(:,:,:,:); 76 param=squeeze(param); 77 paramy{k}=mat2yao(param); 78 paramy{k}(:,1:3)=fliplr(paramy{k}(:,1:3)); 79 paramy{k}=sortrows(paramy{k},[1 2 3]); 47 80 figure(j); 48 81 subplot(3,2,k) 49 82 plot_nemo_3D(paramy{k},longi,lati,zlev,1,sparam{j,2},1); 50 title(suff{k,2}) 83 title(suff{k,2}) 84 85 51 86 if (k==1) 52 87 colorbar … … 59 94 colorbar 60 95 end 61 96 end %if ~ssh 62 97 63 end 98 end %for k 99 64 100 %%Graphe des différences 65 101 for k=1:size(delta,1) 66 102 67 param=paramy{delta{k,1}(1)}; 68 param(:,end)=param(:,end)-paramy{delta{k,1}(2)}(:,end); 69 70 figure(j+2*size(sparam,1)) 71 subplot(2,2,k) 72 plot_nemo_3D(param,longi,lati,zlev,1,sparam{j,4},1); 73 colormap(cmapdif); 74 title(delta{k,2}); 75 if (k==2) 76 colorbar 77 end 78 79 figure(j+3*size(sparam,1)) 80 subplot(2,2,k) 81 plot_nemo_3D(param,longi,lati,zlev,lat0,sparam{j,5},2); 82 colormap(cmapdif); 83 title(delta{k,2}); 84 if (k==2) 85 colorbar 86 end 87 fprintf(1,[sparam{j,1} ' ' delta{k,2} '\n']); 88 [N,RPD,URPD,MAD,RMS,slope,intercept,r2]=make_stats(paramy{delta{k,1}(1)}(:,end),paramy{delta{k,1}(2)}(:,end)) 103 104 param=paramy{delta{k,1}(1)}; 105 param(:,end)=param(:,end)-paramy{delta{k,1}(2)}(:,end); 106 107 figure(j+2*size(sparam,1)) 108 subplot(2,2,k) 109 plot_nemo_3D(param,longi,lati,zlev,1,sparam{j,4},1); 110 colormap(cmapdif); 111 title(delta{k,2}); 112 if (k==2) 113 colorbar 114 end 115 if ~strcmp(sparam{j,1},'ssh') 116 figure(j+3*size(sparam,1)) 117 subplot(2,2,k) 118 plot_nemo_3D(param,longi,lati,zlev,lat0,sparam{j,5},2); 119 colormap(cmapdif); 120 title(delta{k,2}); 121 if (k==2) 122 colorbar 123 end 124 end % if ~ssh 125 126 if strcmp(sparam{j,1},'ssh') 127 n1=max(paramy{delta{k}(1)}(:,1)); 128 n2=max(paramy{delta{k}(1)}(:,2)); 129 iok=find(paramy{delta{k}(1)}(:,1)>1 & paramy{delta{k}(1)}(:,1)<n1 & paramy{delta{k}(1)}(:,2)>1 & ... 130 paramy{delta{k}(1)}(:,2)<n2); 131 else 132 n1=max(paramy{delta{k}(1)}(:,1)); 133 n2=max(paramy{delta{k}(1)}(:,2)); 134 n3=max(paramy{delta{k}(1)}(:,3)); 135 iok=find(paramy{delta{k}(1)}(:,1)>1 & paramy{delta{k}(1)}(:,1)<n1 & paramy{delta{k}(1)}(:,2)>1 & ... 136 paramy{delta{k}(1)}(:,2)<n2 & paramy{delta{k}(1)}(:,3)<n3); 137 end 138 139 140 % fprintf(1,[int2str(j) sparam{j,1} ' ' delta{k,2} '\n']); 141 [N(j,k),RPD(j,k),URPD(j,k),MAD(j,k),RMS(j,k),slope(j,k),intercept(j,k),r2(j,k) ]=make_stats(paramy{delta{k,1}(1)}(iok,end),paramy{delta{k,1}(2)} (iok,end)); 89 142 90 91 end 143 end %for k 92 144 93 145 … … 102 154 end 103 155 close(f) 156 157 %sauvegarde du tableau 158 159 fid=fopen('tab.txt','w'); 160 for j=1:size(sparam,1) 161 162 fprintf(fid,'\n * ParamÚtre %s \n',sparam{j,1}); 163 fprintf(fid,'||Name ||bias|| RMS || rel. err || r2 ||\n'); 164 for k=1:size(delta,1) 165 166 167 fprintf(fid,'||%s %s || %4.3e || %4.3e || %4.3f || %4.3f ||\n',sparam{j,1},delta{k,2}, MAD(j,k),RMS(j,k), ... 168 URPD(j,k),r2(j,k)); 169 170 end 171 172 fprintf(fid,'\n'); 173 174 end 175 fclose(fid);
Note: See TracChangeset
for help on using the changeset viewer.