profilesMatlabProcessing/ecco_v4/insitu_misfit.m

function []=insitu_misfit(dirData,doComp,varargin);
%[]=insitu_misfit(dirData,listData,listVar,varargin);
%object:        computes or displays cost function statistics
%input:         dirData is the data main directory name
%               doComp states whether to compute & save (1) or load & display (0)
%optional must take the following form {'name',param1,param2,...} Active ones are
%               {'subdirData',subdirData} is the data subdirectory name ('profiles/output/' by default)
%               {'listData',data1,data2,...} where data1 e.g. is 'argo_indian.nc'
%               {'listVar',var1,var2} where var1 and var2 are 'T' and 'S'
%               {'years',year0,year1} to limit computation to [year0 year1] interval
%               {'subdirMat',subdirMat} where subdirMat is the output subdirectory ('mat/' by default)
%               {'suffMat',suffMat} is the mat file suffix ('all' by default)
%               {'subdirTex',subdirTex} where subdirTex is the tex directory ('tex/' by default)
%               {'addToTex',addToTex} states whether (1) or not (0 default) to 
%                       augment an ongoing tex file (see write2tex.m, subdirTex)
%
%example:        insitu_misfit('./',1,{'listData','argo_in*'});

vecDepth=[[75 125];[250 350];[900 1100]];
vecCaxis=[[3 3 1];[0.5 0.5 0.2]];
%if nargin>6; vecCaxis=varargin{1}; else; vecCaxis=[]; end;

%load paths and grid:
%--------------------
global mygrid; if isempty(mygrid); grid_load('GRID/',5,'compact'); end;

%set more optional paramaters to default values
global myparms; year0=myparms.yearInAve(1); year1=myparms.yearInAve(2);
%year0=1992; year1=2008; 
listData={'WOD09_XBT*','WOD09_CTD*','argo_in*','argo_at*','argo_pa*','seals*'}; 
listVar={'T','S'}; 
%shorter test case: listData={'argo_in*'}; listVar={'T'};
subdirData='/profiles/output/';
subdirMat='mat/'; suffMat='_all';
subdirTex='tex/'; addToTex=0;
%set more optional paramaters to user defined values
for ii=1:nargin-2;
   if strcmp(varargin{ii}{1},'listData'); listData={varargin{ii}{2:end}};
   elseif strcmp(varargin{ii}{1},'listVar'); listVar={varargin{ii}{2:end}};
   elseif strcmp(varargin{ii}{1},'years'); year0=varargin{ii}{2}; year1=varargin{ii}{3};
   elseif strcmp(varargin{ii}{1},'subdirMat'); subdirMat=varargin{ii}{2};
   elseif strcmp(varargin{ii}{1},'suffMat'); suffMat=['_' varargin{ii}{2}];
   elseif strcmp(varargin{ii}{1},'subdirTex'); subdirTex=varargin{ii}{2};
   elseif strcmp(varargin{ii}{1},'addToTex'); addToTex=varargin{ii}{2};
   elseif strcmp(varargin{ii}{1},'subdirData'); subdirData=varargin{ii}{2};
   else;
       warning('inputCheck:insitu_misfit',...
           ['unknown option ''' varargin{ii}{1} ''' was ignored']);
   end;
end;

if doComp;

%initialize delaunay triangulation:
gcmfaces_bindata;
%time limits:
date0=datenum(year0,1,1); date1=datenum(year1,12,31);
%number of levels:
nk=size(vecDepth,1);

for vv=1:length(listVar);
    varCur=listVar{vv};
    %get the misfits (here not normalized):
    [MITprof]=MITprof_stats_load([dirData '/' subdirData],listData,varCur,1);
    %mask out values that are not in year range:
    ii=find(MITprof.prof_date<date0|MITprof.prof_date>date1);
    MITprof.prof(ii,:)=NaN;
    %misfit maps:
    for kkk=1:size(vecDepth,1);
        kk=find(MITprof.prof_depth>=vecDepth(kkk,1)&MITprof.prof_depth<=vecDepth(kkk,2));
        lon=MITprof.prof_lon; lat=MITprof.prof_lat; obs=nanmean(MITprof.prof(:,kk),2);
        ii=find(~isnan(obs)); lon=lon(ii); lat=lat(ii); obs=obs(ii);
        misfit_map=gcmfaces_bindata(lon,lat,obs);
        eval(['misfit_map_' varCur num2str(kkk) '=misfit_map;']);
    end;
end;

eval(['save ' dirData '/' subdirMat 'insitu_misfit' suffMat '.mat misfit_map_* vecDepth  year* listData listVar;']);

else;

eval(['load ' dirData '/' subdirMat 'insitu_misfit' suffMat '.mat;']);

nk=size(vecDepth,1);
figureL; 
for vv=1:length(listVar);
    varCur=listVar{vv};
    for kkk=1:size(vecDepth,1);
        eval(['misfit_map=misfit_map_' varCur num2str(kkk) ';']); 
        dep=mean(vecDepth(kkk,:));

        if isempty(vecCaxis); cc=3*sqrt(nanmean(misfit_map.^2)); else; cc=vecCaxis(vv,kkk); end;
        subplot(nk,length(listVar),vv+(kkk-1)*length(listVar));
        m_map_gcmfaces(misfit_map,1,{'myCaxis',[-1 1]*cc});
        title(sprintf('%s misfits at %dm   (%d to %d)',varCur,dep,year0,year1));
    end;
end;

myCaption={'Time mean misfit (model-data) for in situ profiles, ',...
         'at various depths (rows), for T (left; in K) and S (right; in psu).'};
if addToTex; write2tex([dirData '/' subdirTex '/myPlots.tex'],2,myCaption,gcf); end;

end;


1	gforget	1.1	function []=insitu_misfit(dirData,doComp,varargin);
2			%[]=insitu_misfit(dirData,listData,listVar,varargin);
3			%object: computes or displays cost function statistics
4			%input: dirData is the data main directory name
5			% doComp states whether to compute & save (1) or load & display (0)
6			%optional must take the following form {'name',param1,param2,...} Active ones are
7			% {'subdirData',subdirData} is the data subdirectory name ('profiles/output/' by default)
8			% {'listData',data1,data2,...} where data1 e.g. is 'argo_indian.nc'
9			% {'listVar',var1,var2} where var1 and var2 are 'T' and 'S'
10			% {'years',year0,year1} to limit computation to [year0 year1] interval
11			% {'subdirMat',subdirMat} where subdirMat is the output subdirectory ('mat/' by default)
12			% {'suffMat',suffMat} is the mat file suffix ('all' by default)
13			% {'subdirTex',subdirTex} where subdirTex is the tex directory ('tex/' by default)
14			% {'addToTex',addToTex} states whether (1) or not (0 default) to
15			% augment an ongoing tex file (see write2tex.m, subdirTex)
16			%
17			%example: insitu_misfit('./',1,{'listData','argo_in*'});
18
19			vecDepth=[[75 125];[250 350];[900 1100]];
20			vecCaxis=[[3 3 1];[0.5 0.5 0.2]];
21			%if nargin>6; vecCaxis=varargin{1}; else; vecCaxis=[]; end;
22
23			%load paths and grid:
24			%--------------------
25			global mygrid; if isempty(mygrid); grid_load('GRID/',5,'compact'); end;
26
27			%set more optional paramaters to default values
28	gforget	1.5	global myparms; year0=myparms.yearInAve(1); year1=myparms.yearInAve(2);
29			%year0=1992; year1=2008;
30	gforget	1.2	listData={'WOD09_XBT','WOD09_CTD','argo_in','argo_at','argo_pa','seals'};
31			listVar={'T','S'};
32			%shorter test case: listData={'argo_in*'}; listVar={'T'};
33	gforget	1.1	subdirData='/profiles/output/';
34			subdirMat='mat/'; suffMat='_all';
35			subdirTex='tex/'; addToTex=0;
36			%set more optional paramaters to user defined values
37			for ii=1:nargin-2;
38			if strcmp(varargin{ii}{1},'listData'); listData={varargin{ii}{2:end}};
39			elseif strcmp(varargin{ii}{1},'listVar'); listVar={varargin{ii}{2:end}};
40			elseif strcmp(varargin{ii}{1},'years'); year0=varargin{ii}{2}; year1=varargin{ii}{3};
41			elseif strcmp(varargin{ii}{1},'subdirMat'); subdirMat=varargin{ii}{2};
42			elseif strcmp(varargin{ii}{1},'suffMat'); suffMat=['_' varargin{ii}{2}];
43			elseif strcmp(varargin{ii}{1},'subdirTex'); subdirTex=varargin{ii}{2};
44			elseif strcmp(varargin{ii}{1},'addToTex'); addToTex=varargin{ii}{2};
45			elseif strcmp(varargin{ii}{1},'subdirData'); subdirData=varargin{ii}{2};
46			else;
47			warning('inputCheck:insitu_misfit',...
48			['unknown option ''' varargin{ii}{1} ''' was ignored']);
49			end;
50			end;
51
52			if doComp;
53
54			%initialize delaunay triangulation:
55			gcmfaces_bindata;
56			%time limits:
57			date0=datenum(year0,1,1); date1=datenum(year1,12,31);
58			%number of levels:
59			nk=size(vecDepth,1);
60
61			for vv=1:length(listVar);
62			varCur=listVar{vv};
63			%get the misfits (here not normalized):
64			[MITprof]=MITprof_stats_load([dirData '/' subdirData],listData,varCur,1);
65			%mask out values that are not in year range:
66			ii=find(MITprof.prof_date<date0\|MITprof.prof_date>date1);
67			MITprof.prof(ii,:)=NaN;
68			%misfit maps:
69			for kkk=1:size(vecDepth,1);
70			kk=find(MITprof.prof_depth>=vecDepth(kkk,1)&MITprof.prof_depth<=vecDepth(kkk,2));
71			lon=MITprof.prof_lon; lat=MITprof.prof_lat; obs=nanmean(MITprof.prof(:,kk),2);
72			ii=find(~isnan(obs)); lon=lon(ii); lat=lat(ii); obs=obs(ii);
73			misfit_map=gcmfaces_bindata(lon,lat,obs);
74			eval(['misfit_map_' varCur num2str(kkk) '=misfit_map;']);
75			end;
76			end;
77
78			eval(['save ' dirData '/' subdirMat 'insitu_misfit' suffMat '.mat misfit_map_* vecDepth year* listData listVar;']);
79
80			else;
81
82			eval(['load ' dirData '/' subdirMat 'insitu_misfit' suffMat '.mat;']);
83
84			nk=size(vecDepth,1);
85	gforget	1.3	figureL;
86	gforget	1.1	for vv=1:length(listVar);
87			varCur=listVar{vv};
88			for kkk=1:size(vecDepth,1);
89			eval(['misfit_map=misfit_map_' varCur num2str(kkk) ';']);
90			dep=mean(vecDepth(kkk,:));
91
92			if isempty(vecCaxis); cc=3*sqrt(nanmean(misfit_map.^2)); else; cc=vecCaxis(vv,kkk); end;
93			subplot(nk,length(listVar),vv+(kkk-1)*length(listVar));
94			m_map_gcmfaces(misfit_map,1,{'myCaxis',[-1 1]*cc});
95			title(sprintf('%s misfits at %dm (%d to %d)',varCur,dep,year0,year1));
96			end;
97			end;
98
99	gforget	1.4	myCaption={'Time mean misfit (model-data) for in situ profiles, ',...
100	gforget	1.5	'at various depths (rows), for T (left; in K) and S (right; in psu).'};
101	gforget	1.1	if addToTex; write2tex([dirData '/' subdirTex '/myPlots.tex'],2,myCaption,gcf); end;
102
103			end;
104
105
106