profilesMatlabProcessing/profiles_stats/MITprof_stats_misfit.m

function []=MITprof_stats_misfit(dirData,listData,listVar,year0,year1,vecLevel,varargin);
%[]=MITprof_stats_misfit(dirData,listData,listVar,year0,year1,vecLevel);
%object:        maps and displays estim-obs at various levels
%input:         dirData is the data directory name
%               listData is the data file list (e.g. {'argo_in*'} or {'argo_in*','argo_at*'} )
%               listVar is the variable list varCur (e.g. {'T'} or {'T','S'} )
%               year0,year1 is the first,last year pair
%               vecLevel is the vector of depth levels
%               vecCaxis (optional) is the matrix of colorbar range
%output:        (none -- figures only)
%
%for example:
% dirData='./';
% listData={'argo_in*'}; listVar={'T','S'};
% MITprof_stats_misfit(dirData,listData,listVar,2004,2008,[10 25 45],[[3 3 1];[0.5 0.5 0.2]]);

%load paths and grid:
%--------------------
gcmfaces_path; global mygrid; mygrid=[]; grid_load([gcmfaces_dir '/sample_input/GRIDv4/'],5,'compact');

gcmfaces_bindata;%initialize delaunay triangulation

%time limits:
date0=datenum(year0,1,1); date1=datenum(year1,12,31);

%number of levels:
nk=length(vecLevel);

%caxis limits:
if nargin>6; vecCaxis=varargin{1}; else; vecCaxis=[]; end;

figureL; orient tall;

for vv=1:length(listVar);
    varCur=listVar{vv};
    
    %get the cost square root:
    [MITprof]=MITprof_stats_load(dirData,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:length(vecLevel);
        kk=vecLevel(kkk);
        
        lon=MITprof.prof_lon; lat=MITprof.prof_lat; obs=MITprof.prof(:,kk);
        ii=find(~isnan(obs)); lon=lon(ii); lat=lat(ii); obs=obs(ii);
        misfit_map=gcmfaces_bindata(lon,lat,obs);
        
        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,[-1 1]*cc);
        title(sprintf('%s misfits at %dm   (%d to %d)',varCur,MITprof.prof_depth(kk),year0,year1));
        
    end;
    
    
end;

1	gforget	1.1	function []=MITprof_stats_misfit(dirData,listData,listVar,year0,year1,vecLevel,varargin);
2			%[]=MITprof_stats_misfit(dirData,listData,listVar,year0,year1,vecLevel);
3			%object: maps and displays estim-obs at various levels
4			%input: dirData is the data directory name
5			% listData is the data file list (e.g. {'argo_in'} or {'argo_in','argo_at*'} )
6			% listVar is the variable list varCur (e.g. {'T'} or {'T','S'} )
7			% year0,year1 is the first,last year pair
8			% vecLevel is the vector of depth levels
9			% vecCaxis (optional) is the matrix of colorbar range
10			%output: (none -- figures only)
11			%
12			%for example:
13			% dirData='./';
14			% listData={'argo_in*'}; listVar={'T','S'};
15			% MITprof_stats_misfit(dirData,listData,listVar,2004,2008,[10 25 45],[[3 3 1];[0.5 0.5 0.2]]);
16
17			%load paths and grid:
18			%--------------------
19	gforget	1.3	gcmfaces_path; global mygrid; mygrid=[]; grid_load([gcmfaces_dir '/sample_input/GRIDv4/'],5,'compact');
20	gforget	1.1
21			gcmfaces_bindata;%initialize delaunay triangulation
22
23			%time limits:
24			date0=datenum(year0,1,1); date1=datenum(year1,12,31);
25
26			%number of levels:
27			nk=length(vecLevel);
28
29			%caxis limits:
30			if nargin>6; vecCaxis=varargin{1}; else; vecCaxis=[]; end;
31
32			figureL; orient tall;
33
34			for vv=1:length(listVar);
35			varCur=listVar{vv};
36
37			%get the cost square root:
38			[MITprof]=MITprof_stats_load(dirData,listData,varCur,1);
39
40			%mask out values that are not in year range:
41			ii=find(MITprof.prof_date<date0\|MITprof.prof_date>date1);
42			MITprof.prof(ii,:)=NaN;
43
44
45			%misfit maps:
46			for kkk=1:length(vecLevel);
47			kk=vecLevel(kkk);
48
49			lon=MITprof.prof_lon; lat=MITprof.prof_lat; obs=MITprof.prof(:,kk);
50			ii=find(~isnan(obs)); lon=lon(ii); lat=lat(ii); obs=obs(ii);
51			misfit_map=gcmfaces_bindata(lon,lat,obs);
52
53			if isempty(vecCaxis); cc=3*sqrt(nanmean(misfit_map.^2)); else; cc=vecCaxis(vv,kkk); end;
54			subplot(nk,length(listVar),vv+(kkk-1)*length(listVar));
55			m_map_gcmfaces(misfit_map,-1,[-1 1]*cc);
56			title(sprintf('%s misfits at %dm (%d to %d)',varCur,MITprof.prof_depth(kk),year0,year1));
57
58			end;
59
60
61			end;
62