/[MITgcm]/MITgcm_contrib/gael/profilesMatlabProcessing/ecco_v4/insitu_misfit.m
ViewVC logotype

Diff of /MITgcm_contrib/gael/profilesMatlabProcessing/ecco_v4/insitu_misfit.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by gforget, Fri Jul 1 18:03:43 2011 UTC revision 1.12 by gforget, Mon Feb 15 18:07:32 2016 UTC
# Line 1  Line 1 
1  function []=insitu_misfit(dirData,doComp,varargin);  function []=insitu_misfit(doComp,varargin);
 %[]=insitu_misfit(dirData,listData,listVar,varargin);  
2  %object:        computes or displays cost function statistics  %object:        computes or displays cost function statistics
3  %input:         dirData is the data main directory name  %input:         doComp states whether to compute & save (1) or load & display (0)
 %               doComp states whether to compute & save (1) or load & display (0)  
4  %optional must take the following form {'name',param1,param2,...} Active ones are  %optional must take the following form {'name',param1,param2,...} Active ones are
5  %               {'subdirData',subdirData} is the data subdirectory name ('profiles/output/' by default)  %               {'dirData',dirData} is the data directory name ('mat/profiles/output/' by default)
6  %               {'listData',data1,data2,...} where data1 e.g. is 'argo_indian.nc'  %               {'listData',data1,data2,...} where data1 e.g. is 'argo_indian.nc'
7  %               {'listVar',var1,var2} where var1 and var2 are 'T' and 'S'  %               {'listVar',var1,var2} where var1 and var2 are 'T' and 'S'
8  %               {'years',year0,year1} to limit computation to [year0 year1] interval  %               {'years',year0,year1} to limit computation to [year0 year1] interval
9  %               {'subdirMat',subdirMat} where subdirMat is the output subdirectory ('mat/' by default)  %               {'dirMat',dirMat} where dirMat is the output directory ('mat/' by default)
10  %               {'suffMat',suffMat} is the mat file suffix ('all' by default)  %               {'suffMat',suffMat} is the mat file suffix ('all' by default)
11  %               {'subdirTex',subdirTex} where subdirTex is the tex directory ('tex/' by default)  %               {'dirTex',dirTex} where dirTex is the tex directory ('tex/' by default)
12    %               {'nameTex',nameTex} is the tex file name ('myPlots' by default)
13  %               {'addToTex',addToTex} states whether (1) or not (0 default) to  %               {'addToTex',addToTex} states whether (1) or not (0 default) to
14  %                       augment an ongoing tex file (see write2tex.m, subdirTex)  %                       augment an ongoing tex file (see write2tex.m, dirTex)
15  %  %
16  %example:        insitu_misfit('./',1,{'listData','argo_in*'});  %example:        insitu_misfit('./',1,{'listData','argo_in*'});
17    
# Line 20  vecDepth=[[75 125];[250 350];[900 1100]] Line 19  vecDepth=[[75 125];[250 350];[900 1100]]
19  vecCaxis=[[3 3 1];[0.5 0.5 0.2]];  vecCaxis=[[3 3 1];[0.5 0.5 0.2]];
20  %if nargin>6; vecCaxis=varargin{1}; else; vecCaxis=[]; end;  %if nargin>6; vecCaxis=varargin{1}; else; vecCaxis=[]; end;
21    
22  %load paths and grid:  gcmfaces_global;
23  %--------------------  global myparms;
24  global mygrid; if isempty(mygrid); grid_load('GRID/',5,'compact'); end;  
25    listVar={'T','S'};
26  %set more optional paramaters to default values  %shorter test case: listData={'argo_in*'}; listVar={'T'};
27  year0=1992; year1=2008;  dirData='mat/profiles/output/';
28  %listData={'WOD09_XBT*','WOD09_CTD*','argo_in*','argo_at*','argo_pa*','seals*'};  dirMat='mat/'; suffMat='_all';
29  %listVar={'T','S'};  dirTex='tex/'; addToTex=0; nameTex='myPlots';
 listData={'argo_in*'}; listVar={'T'};  
 subdirData='/profiles/output/';  
 subdirMat='mat/'; suffMat='_all';  
 subdirTex='tex/'; addToTex=0;  
30  %set more optional paramaters to user defined values  %set more optional paramaters to user defined values
31  for ii=1:nargin-2;  for ii=1:nargin-1;
32     if strcmp(varargin{ii}{1},'listData'); listData={varargin{ii}{2:end}};     if strcmp(varargin{ii}{1},'listData'); listData={varargin{ii}{2:end}};
33     elseif strcmp(varargin{ii}{1},'listVar'); listVar={varargin{ii}{2:end}};     elseif strcmp(varargin{ii}{1},'listVar'); listVar={varargin{ii}{2:end}};
34     elseif strcmp(varargin{ii}{1},'years'); year0=varargin{ii}{2}; year1=varargin{ii}{3};     elseif strcmp(varargin{ii}{1},'years'); year0=varargin{ii}{2}; year1=varargin{ii}{3};
35     elseif strcmp(varargin{ii}{1},'subdirMat'); subdirMat=varargin{ii}{2};     elseif strcmp(varargin{ii}{1},'dirMat'); dirMat=varargin{ii}{2};
36     elseif strcmp(varargin{ii}{1},'suffMat'); suffMat=['_' varargin{ii}{2}];     elseif strcmp(varargin{ii}{1},'suffMat'); suffMat=['_' varargin{ii}{2}];
37     elseif strcmp(varargin{ii}{1},'subdirTex'); subdirTex=varargin{ii}{2};     elseif strcmp(varargin{ii}{1},'dirTex'); dirTex=varargin{ii}{2};
38     elseif strcmp(varargin{ii}{1},'addToTex'); addToTex=varargin{ii}{2};     elseif strcmp(varargin{ii}{1},'addToTex'); addToTex=varargin{ii}{2};
39     elseif strcmp(varargin{ii}{1},'subdirData'); subdirData=varargin{ii}{2};     elseif strcmp(varargin{ii}{1},'nameTex'); nameTex=varargin{ii}{2};
40       elseif strcmp(varargin{ii}{1},'dirData'); dirData=varargin{ii}{2};
41     else;     else;
42         warning('inputCheck:insitu_misfit',...         warning('inputCheck:insitu_misfit',...
43             ['unknown option ''' varargin{ii}{1} ''' was ignored']);             ['unknown option ''' varargin{ii}{1} ''' was ignored']);
44     end;     end;
45  end;  end;
46    
47    if isempty(myparms)|isempty(mygrid);
48      load([dirMat 'diags_grid_parms.mat']);
49    end;
50    
51    if isempty(whos('year0'));
52      year0=myparms.yearInAve(1); year1=myparms.yearInAve(2);
53    end;
54    
55    if isempty(whos('listData'));
56      listData=dir([dirMat 'profiles/output/*.nc']);
57      listData={listData(:).name};
58      for ff=1:length(listData); listData{ff}=[listData{ff}(1:end-3) '*']; end;
59    end;
60    
61  if doComp;  if doComp;
62    
63  %initialize delaunay triangulation:  %initialize delaunay triangulation:
# Line 60  nk=size(vecDepth,1); Line 70  nk=size(vecDepth,1);
70  for vv=1:length(listVar);  for vv=1:length(listVar);
71      varCur=listVar{vv};      varCur=listVar{vv};
72      %get the misfits (here not normalized):      %get the misfits (here not normalized):
73      [MITprof]=MITprof_stats_load([dirData '/' subdirData],listData,varCur,1);      [MITprof]=MITprof_stats_load(dirData,listData,varCur,1);
74      %mask out values that are not in year range:      %mask out values that are not in year range:
75      ii=find(MITprof.prof_date<date0|MITprof.prof_date>date1);      ii=find(MITprof.prof_date<date0|MITprof.prof_date>date1);
76      MITprof.prof(ii,:)=NaN;      MITprof.prof(ii,:)=NaN;
# Line 74  for vv=1:length(listVar); Line 84  for vv=1:length(listVar);
84      end;      end;
85  end;  end;
86    
87  eval(['save ' dirData '/' subdirMat 'insitu_misfit' suffMat '.mat misfit_map_* vecDepth  year* listData listVar;']);  if ~isdir([dirMat 'cost/']); mkdir([dirMat 'cost/']); end;
88    eval(['save ' dirMat 'cost/insitu_misfit' suffMat '.mat misfit_map_* vecDepth  year* listData listVar;']);
89    
90  else;  else;
91    
92  eval(['load ' dirData '/' subdirMat 'insitu_misfit' suffMat '.mat;']);  if isdir([dirMat 'cost/']); dirMat=[dirMat 'cost/']; end;
93    
94    eval(['load ' dirMat 'insitu_misfit' suffMat '.mat;']);
95    
96  nk=size(vecDepth,1);  nk=size(vecDepth,1);
97  figureL; orient tall;  figureL;
98  for vv=1:length(listVar);  for vv=1:length(listVar);
99      varCur=listVar{vv};      varCur=listVar{vv};
100      for kkk=1:size(vecDepth,1);      for kkk=1:size(vecDepth,1);
# Line 90  for vv=1:length(listVar); Line 103  for vv=1:length(listVar);
103    
104          if isempty(vecCaxis); cc=3*sqrt(nanmean(misfit_map.^2)); else; cc=vecCaxis(vv,kkk); end;          if isempty(vecCaxis); cc=3*sqrt(nanmean(misfit_map.^2)); else; cc=vecCaxis(vv,kkk); end;
105          subplot(nk,length(listVar),vv+(kkk-1)*length(listVar));          subplot(nk,length(listVar),vv+(kkk-1)*length(listVar));
106          m_map_gcmfaces(misfit_map,1,{'myCaxis',[-1 1]*cc});          m_map_gcmfaces(misfit_map,1,{'myCaxis',[-1 1]*cc},{'myShading','flat'});
107          title(sprintf('%s misfits at %dm   (%d to %d)',varCur,dep,year0,year1));          title(sprintf('%s misfits at %dm   (%d to %d)',varCur,dep,year0,year1));
108      end;      end;
109  end;  end;
110    
111  myCaption={'model minus in situ profiles'};  myCaption={'Time mean misfit (model-data) for in situ profiles, ',...
112  if addToTex; write2tex([dirData '/' subdirTex '/myPlots.tex'],2,myCaption,gcf); end;           'at various depths (rows), for T (left; in K) and S (right; in psu).'};
113    if addToTex; write2tex([dirTex '/' nameTex '.tex'],2,myCaption,gcf); end;
114    
115  end;  end;
116    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22