Changeset 16 for trunk/data_out

Timestamp:

03/22/11 10:35:04 (13 years ago)

Author:

jbrlod

Message:

make v1.1

Location:

trunk/data_out

Files:

• plot_expref.m (deleted)
• plot_results.m (modified) (6 diffs)
• sa_c_fin.dat (deleted)
• sa_c_init.dat (deleted)
• sa_c_true.dat (deleted)
• sb_obs_48.dat (deleted)
• sb_obs_8.dat (deleted)
• sb_sbd100.dat (deleted)
• sshb_obs20.dat (deleted)
• sshb_obs20_0_1.dat (deleted)
• sshb_obs_48.dat (deleted)
• sshb_obs_8.dat (deleted)
• sshb_obs_std_424.dat (deleted)
• sshb_sbd100.dat (deleted)
• sshn_c_fin.dat (deleted)
• sshn_c_init.dat (deleted)
• sshn_c_true.dat (deleted)
• ta_c_fin.dat (deleted)
• ta_c_init.dat (deleted)
• ta_c_true.dat (deleted)
• ta_c_yao.dat (deleted)
• tb100.dat (deleted)
• tb_50.dat (deleted)
• tb_50_0.dat (deleted)
• tb_forw_init_424.dat (deleted)
• tb_obs20.dat (deleted)
• tb_obs20_0_1.dat (deleted)
• tb_obs20_0_2.dat (deleted)
• tb_obs20_2.dat (deleted)
• tb_obs_48.dat (deleted)
• tb_obs_8.dat (deleted)
• tb_obs_std_16_12_1_424.dat (deleted)
• tb_obs_std_424.dat (deleted)
• tb_ret_424.dat (deleted)
• tb_sbd100.dat (deleted)
• tb_st10.dat (deleted)
• ua_c_fin.dat (deleted)
• ua_c_init.dat (deleted)
• ua_c_true.dat (deleted)
• ub_obs_48.dat (deleted)
• ub_obs_8.dat (deleted)
• ub_sbd100.dat (deleted)
• va_c_fin.dat (deleted)
• va_c_init.dat (deleted)
• va_c_true.dat (deleted)
• vb_obs_48.dat (deleted)
• vb_obs_8.dat (deleted)
• vb_sbd100.dat (deleted)

trunk/data_out/plot_results.m

@@ +1 @@
+f=netcdf('exp_T.nc');
+
+addpath('../scripts/matlab_toolbox');
 sparam=cell(0,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]};
 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]};
-f=netcdf('exp_T.nc');
+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]};
 
 suff=cell(0,2);
@@ -13 +16 @@
 suff(end+1,:)={'a_c_true','parametre vrai (TP)'};
 suff(end+1,:)={'b_obs_48','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,:)={'n_c_fin','param control final (AP)'};
+suff2(end+1,:)={'b_forwfin_48','forward du PA'};
+suff2(end+1,:)={'n_c_true','parametre vrai (TP)'};
+suff2(end+1,:)={'b_obs_48','forward du TP'};
+
 
 delta=cell(0,2);
@@ -38 +50 @@
   for k=1:length(suff)
     
+   if strcmp(sparam{j,1},'ssh')
+     param=f{[sparam{j,1} suff2{k,1}]}(:,:,:);
+     
+      param=squeeze(param);
+     paramy{k}=mat2yao(param);
+     paramy{k}=mat2yao(param);
+     paramy{k}(:,1:2)=fliplr(paramy{k}(:,1:2));
+     paramy{k}=sortrows(paramy{k},[1 2]);
+     figure(j);
+    subplot(3,2,k)
+    % 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});
+    % caxis([sparam{j,2}(1) sparam{j,2}(end)])
+     %xlabel('longitude','Fontsize',9)
+     %ylabel('latitude','Fontsize',9)
+      plot_nemo_3D(paramy{k},longi,lati,zlev,1,sparam{j,2},1);
+
+    
+    title(suff2{k,2});
+     if (k==1) 
+      colorbar
+     end
    
-    param=f{[sparam{j,1} suff{k,1}]}(:,:,:,:);
-    param=squeeze(param);
-    paramy{k}=mat2yao(param);
-    paramy{k}(:,1:3)=fliplr(paramy{k}(:,1:3));
-    paramy{k}=sortrows(paramy{k},[1 2 3]);
-   
+   else
+     param=f{[sparam{j,1} suff{k,1}]}(:,:,:,:);
+     param=squeeze(param);
+     paramy{k}=mat2yao(param);
+     paramy{k}(:,1:3)=fliplr(paramy{k}(:,1:3));
+     paramy{k}=sortrows(paramy{k},[1 2 3]);
     figure(j);
     subplot(3,2,k)
     plot_nemo_3D(paramy{k},longi,lati,zlev,1,sparam{j,2},1);
-    title(suff{k,2})
+    title(suff{k,2}) 
+  
+   
     if (k==1) 
       colorbar
@@ -59 +94 @@
       colorbar
     end
-   
+   end %if ~ssh
     
-    end    
+  end %for k   
+  
      %%Graphe des différences
     for k=1:size(delta,1)
       
-    param=paramy{delta{k,1}(1)};
-    param(:,end)=param(:,end)-paramy{delta{k,1}(2)}(:,end);
-    
-    figure(j+2*size(sparam,1))
-    subplot(2,2,k)
-    plot_nemo_3D(param,longi,lati,zlev,1,sparam{j,4},1);
-    colormap(cmapdif);
-    title(delta{k,2});
-    if (k==2)
-      colorbar
-    end
-    
-      figure(j+3*size(sparam,1))
-    subplot(2,2,k)
-     plot_nemo_3D(param,longi,lati,zlev,lat0,sparam{j,5},2);
-    colormap(cmapdif);
-    title(delta{k,2});
-    if (k==2)
-      colorbar
-    end
-    fprintf(1,[sparam{j,1} ' ' delta{k,2} '\n']);
-    [N,RPD,URPD,MAD,RMS,slope,intercept,r2]=make_stats(paramy{delta{k,1}(1)}(:,end),paramy{delta{k,1}(2)}(:,end))
+      
+      param=paramy{delta{k,1}(1)};
+      param(:,end)=param(:,end)-paramy{delta{k,1}(2)}(:,end);
+      
+      figure(j+2*size(sparam,1))
+      subplot(2,2,k)
+      plot_nemo_3D(param,longi,lati,zlev,1,sparam{j,4},1);
+      colormap(cmapdif);
+      title(delta{k,2});
+      if (k==2)
+        colorbar
+      end
+      if ~strcmp(sparam{j,1},'ssh')
+        figure(j+3*size(sparam,1))
+        subplot(2,2,k)
+        plot_nemo_3D(param,longi,lati,zlev,lat0,sparam{j,5},2);
+        colormap(cmapdif);
+        title(delta{k,2});
+        if (k==2)
+          colorbar
+        end
+      end % if ~ssh
+      
+      if strcmp(sparam{j,1},'ssh')
+        n1=max(paramy{delta{k}(1)}(:,1));
+        n2=max(paramy{delta{k}(1)}(:,2));
+        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);
+      else
+        n1=max(paramy{delta{k}(1)}(:,1));
+        n2=max(paramy{delta{k}(1)}(:,2));
+        n3=max(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)<n3);
+      end
+      
+      
+         %  fprintf(1,[int2str(j) sparam{j,1} ' ' delta{k,2} '\n']);
+      [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));
 
-    
-  end
+    end %for k
   
   
@@ -102 +154 @@
 end
 close(f)
+
+%sauvegarde du tableau
+
+fid=fopen('tab.txt','w');
+for j=1:size(sparam,1)
+
+  fprintf(fid,'\n * ParamÚtre %s \n',sparam{j,1});
+  fprintf(fid,'||Name  ||bias|| RMS || rel. err || r2    ||\n');
+    for k=1:size(delta,1)
+      
+
+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), ...
+        URPD(j,k),r2(j,k));
+
+    end
+
+    fprintf(fid,'\n');
+
+end
+fclose(fid);