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 |
|
|
25 |
|
|
26 |
%set more optional paramaters to default values |
%set more optional paramaters to default values |
27 |
global myparms; year0=myparms.yearInAve(1); year1=myparms.yearInAve(2); |
global myparms; year0=myparms.yearInAve(1); year1=myparms.yearInAve(2); |
28 |
%year0=1992; year1=2008; |
listData={'argo_june2012_1992_to_2007*','argo_june2012_2008_to_2010*',... |
29 |
listData={'WOD09_XBT*','WOD09_CTD*','argo_in*','argo_at*','argo_pa*','seals*'}; |
'argo_june2012_2011_to_2012*','climode_june2012*',... |
30 |
|
'ctd_june2012*','itp_june2012*','seals_june2012*','xbt_june2012*'}; |
31 |
listVar={'T','S'}; |
listVar={'T','S'}; |
32 |
%shorter test case: listData={'argo_in*'}; listVar={'T'}; |
%shorter test case: listData={'argo_in*'}; listVar={'T'}; |
33 |
subdirData='/profiles/output/'; |
dirData='mat/profiles/output/'; |
34 |
subdirMat='mat/'; suffMat='_all'; |
dirMat='mat/'; suffMat='_all'; |
35 |
subdirTex='tex/'; addToTex=0; |
dirTex='tex/'; addToTex=0; nameTex='myPlots'; |
36 |
%set more optional paramaters to user defined values |
%set more optional paramaters to user defined values |
37 |
for ii=1:nargin-2; |
for ii=1:nargin-1; |
38 |
if strcmp(varargin{ii}{1},'listData'); listData={varargin{ii}{2:end}}; |
if strcmp(varargin{ii}{1},'listData'); listData={varargin{ii}{2:end}}; |
39 |
elseif strcmp(varargin{ii}{1},'listVar'); listVar={varargin{ii}{2:end}}; |
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}; |
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}; |
elseif strcmp(varargin{ii}{1},'dirMat'); dirMat=varargin{ii}{2}; |
42 |
elseif strcmp(varargin{ii}{1},'suffMat'); suffMat=['_' varargin{ii}{2}]; |
elseif strcmp(varargin{ii}{1},'suffMat'); suffMat=['_' varargin{ii}{2}]; |
43 |
elseif strcmp(varargin{ii}{1},'subdirTex'); subdirTex=varargin{ii}{2}; |
elseif strcmp(varargin{ii}{1},'dirTex'); dirTex=varargin{ii}{2}; |
44 |
elseif strcmp(varargin{ii}{1},'addToTex'); addToTex=varargin{ii}{2}; |
elseif strcmp(varargin{ii}{1},'addToTex'); addToTex=varargin{ii}{2}; |
45 |
elseif strcmp(varargin{ii}{1},'subdirData'); subdirData=varargin{ii}{2}; |
elseif strcmp(varargin{ii}{1},'nameTex'); nameTex=varargin{ii}{2}; |
46 |
|
elseif strcmp(varargin{ii}{1},'dirData'); dirData=varargin{ii}{2}; |
47 |
else; |
else; |
48 |
warning('inputCheck:insitu_misfit',... |
warning('inputCheck:insitu_misfit',... |
49 |
['unknown option ''' varargin{ii}{1} ''' was ignored']); |
['unknown option ''' varargin{ii}{1} ''' was ignored']); |
62 |
for vv=1:length(listVar); |
for vv=1:length(listVar); |
63 |
varCur=listVar{vv}; |
varCur=listVar{vv}; |
64 |
%get the misfits (here not normalized): |
%get the misfits (here not normalized): |
65 |
[MITprof]=MITprof_stats_load([dirData '/' subdirData],listData,varCur,1); |
[MITprof]=MITprof_stats_load(dirData,listData,varCur,1); |
66 |
%mask out values that are not in year range: |
%mask out values that are not in year range: |
67 |
ii=find(MITprof.prof_date<date0|MITprof.prof_date>date1); |
ii=find(MITprof.prof_date<date0|MITprof.prof_date>date1); |
68 |
MITprof.prof(ii,:)=NaN; |
MITprof.prof(ii,:)=NaN; |
76 |
end; |
end; |
77 |
end; |
end; |
78 |
|
|
79 |
eval(['save ' dirData '/' subdirMat 'insitu_misfit' suffMat '.mat misfit_map_* vecDepth year* listData listVar;']); |
eval(['save ' dirMat 'insitu_misfit' suffMat '.mat misfit_map_* vecDepth year* listData listVar;']); |
80 |
|
|
81 |
else; |
else; |
82 |
|
|
83 |
eval(['load ' dirData '/' subdirMat 'insitu_misfit' suffMat '.mat;']); |
eval(['load ' dirMat 'insitu_misfit' suffMat '.mat;']); |
84 |
|
|
85 |
nk=size(vecDepth,1); |
nk=size(vecDepth,1); |
86 |
figureL; |
figureL; |
99 |
|
|
100 |
myCaption={'Time mean misfit (model-data) for in situ profiles, ',... |
myCaption={'Time mean misfit (model-data) for in situ profiles, ',... |
101 |
'at various depths (rows), for T (left; in K) and S (right; in psu).'}; |
'at various depths (rows), for T (left; in K) and S (right; in psu).'}; |
102 |
if addToTex; write2tex([dirData '/' subdirTex '/myPlots.tex'],2,myCaption,gcf); end; |
if addToTex; write2tex([dirTex '/' nameTex '.tex'],2,myCaption,gcf); end; |
103 |
|
|
104 |
end; |
end; |
105 |
|
|