profilesMatlabProcessing/profiles_misc/MITprof_subset.m

function [MITprofSub]=MITprof_subset(MITprof,varargin);
% [MITprofSub]=MITprof_subset(MITprof,'PropertyName',PropertyValue,...)
%   load a subset of profiles from MITprof into MITprofSub
%
%   if PropertyName = 'list' : PropertyValue contains the list of profiles
%   otherwise PropertyName = 'field' : selection on MITprof.prof_field,
%       range or unique value given in PropertyValue.
%   if PropertyName = 'list' is used, then it needs to be the first one.
%
%   [MITprofSub]=MITprof_subset(MITprof,'list',1:50);
%       returns the 50th first profiles in MITprof
%   [MITprofSub]=MITprof_subset(MITprof,'depth',[50 150]);
%       returns the subset of MITprof that has prof_depth>=50 and <150
%   [MITprofSub]=MITprof_subset(MITprof,'depth',[50 150],'descr',platformName);
%       ... and prof_descr=platformName;

MITprofSub=MITprof;
fldNames=fieldnames(MITprof);

nSub=(nargin-1)/2;
for iSub=1:nSub;
    if ~strcmp('list',varargin{(iSub-1)*2+1});
        eval(['prof_sub=MITprofSub.prof_' varargin{(iSub-1)*2+1} ';']);
        range_sub=varargin{(iSub-1)*2+2};
        if strcmp('descr',varargin{(iSub-1)*2+1});
            KK=find(strcmp(prof_sub,range_sub));
        elseif length(range_sub)==2 & range_sub(1)<range_sub(2);
            KK=find(prof_sub>=range_sub(1)&prof_sub<range_sub(2));
        elseif length(range_sub)==1;
            KK=find(prof_sub==range_sub);
        else
            error('MITprof_subset: wrong arguments');
        end;
    else;
        KK=varargin{(iSub-1)*2+2};
        prof_sub=[1:MITprofSub.np]';
        if iSub~=1; error('''list'' subset should come first'); end;
    end;
    
    %  [varargin{(iSub-1)*2+1} ' -- ' num2str(length(KK))]
    for iFld=1:length(fldNames);
        eval(['tmp1=MITprofSub.' fldNames{iFld} ';']);
        if strcmp('depth',varargin{(iSub-1)*2+1}) & strcmp(fldNames{iFld},'prof_depth');
            tmp1=tmp1(KK);
        elseif strcmp('depth',varargin{(iSub-1)*2+1});
            if size(tmp1,2)==length(MITprofSub.prof_depth); tmp1=tmp1(:,KK); end;
        else;
            if size(tmp1,1)==size(prof_sub,1); tmp1=tmp1(KK,:); end;
        end;
        eval(['MITprofSub.' fldNames{iFld} '=tmp1;']);
    end;
        
    %add a couple things:
    %--------------------
    MITprofSub.np=length(MITprofSub.prof_lon);
    MITprofSub.nr=length(MITprofSub.prof_depth);
    MITprofSub.list_descr=unique(MITprofSub.prof_descr);
    MITprofSub.nd=length(MITprofSub.list_descr);
    
end;
1	function [MITprofSub]=MITprof_subset(MITprof,varargin);
2	% [MITprofSub]=MITprof_subset(MITprof,'PropertyName',PropertyValue,...)
3	% load a subset of profiles from MITprof into MITprofSub
4	%
5	% if PropertyName = 'list' : PropertyValue contains the list of profiles
6	% otherwise PropertyName = 'field' : selection on MITprof.prof_field,
7	% range or unique value given in PropertyValue.
8	% if PropertyName = 'list' is used, then it needs to be the first one.
9	%
10	% [MITprofSub]=MITprof_subset(MITprof,'list',1:50);
11	% returns the 50th first profiles in MITprof
12	% [MITprofSub]=MITprof_subset(MITprof,'depth',[50 150]);
13	% returns the subset of MITprof that has prof_depth>=50 and <150
14	% [MITprofSub]=MITprof_subset(MITprof,'depth',[50 150],'descr',platformName);
15	% ... and prof_descr=platformName;
16
17	MITprofSub=MITprof;
18	fldNames=fieldnames(MITprof);
19
20	nSub=(nargin-1)/2;
21	for iSub=1:nSub;
22	if ~strcmp('list',varargin{(iSub-1)*2+1});
23	eval(['prof_sub=MITprofSub.prof_' varargin{(iSub-1)*2+1} ';']);
24	range_sub=varargin{(iSub-1)*2+2};
25	if strcmp('descr',varargin{(iSub-1)*2+1});
26	KK=find(strcmp(prof_sub,range_sub));
27	elseif length(range_sub)==2 & range_sub(1)<range_sub(2);
28	KK=find(prof_sub>=range_sub(1)&prof_sub<range_sub(2));
29	elseif length(range_sub)==1;
30	KK=find(prof_sub==range_sub);
31	else
32	error('MITprof_subset: wrong arguments');
33	end;
34	else;
35	KK=varargin{(iSub-1)*2+2};
36	prof_sub=[1:MITprofSub.np]';
37	if iSub~=1; error('''list'' subset should come first'); end;
38	end;
39
40	% [varargin{(iSub-1)*2+1} ' -- ' num2str(length(KK))]
41	for iFld=1:length(fldNames);
42	eval(['tmp1=MITprofSub.' fldNames{iFld} ';']);
43	if strcmp('depth',varargin{(iSub-1)*2+1}) & strcmp(fldNames{iFld},'prof_depth');
44	tmp1=tmp1(KK);
45	elseif strcmp('depth',varargin{(iSub-1)*2+1});
46	if size(tmp1,2)==length(MITprofSub.prof_depth); tmp1=tmp1(:,KK); end;
47	else;
48	if size(tmp1,1)==size(prof_sub,1); tmp1=tmp1(KK,:); end;
49	end;
50	eval(['MITprofSub.' fldNames{iFld} '=tmp1;']);
51	end;
52
53	%add a couple things:
54	%--------------------
55	MITprofSub.np=length(MITprofSub.prof_lon);
56	MITprofSub.nr=length(MITprofSub.prof_depth);
57	MITprofSub.list_descr=unique(MITprofSub.prof_descr);
58	MITprofSub.nd=length(MITprofSub.list_descr);
59
60	end;