profilesMatlabProcessing/profiles_misc/MITprof_subset.m

function [MITprofSub]=MITprof_subset(MITprof,varargin);
%   [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;
    eval(['prof_sub=MITprof.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 strcmp('depth',varargin{(iSub-1)*2+1});
        KK=find(prof_sub>=range_sub(1)&prof_sub<range_sub(2));
    else;
        KK=find(prof_sub>=range_sub(1)&prof_sub<range_sub(2));
    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(MITprof.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;
    
end;

%add a couple things:
%--------------------
MITprofSub.np=length(MITprofSub.prof_lon);
MITprofSub.nr=length(MITprofSub.prof_depth);

1	gforget	1.1	function [MITprofSub]=MITprof_subset(MITprof,varargin);
2			% [MITprofSub]=MITprof_subset(MITprof,'depth',[50 150]);
3			% returns the subset of MITprof that has prof_depth>=50 and <150
4			% [MITprofSub]=MITprof_subset(MITprof,'depth',[50 150],'descr',platformName);
5			% ... and prof_descr=platformName;
6
7			MITprofSub=MITprof;
8			fldNames=fieldnames(MITprof);
9
10			nSub=(nargin-1)/2;
11			for iSub=1:nSub;
12			eval(['prof_sub=MITprof.prof_' varargin{(iSub-1)*2+1} ';']);
13			range_sub=varargin{(iSub-1)*2+2};
14			if strcmp('descr',varargin{(iSub-1)*2+1});
15			KK=find(strcmp(prof_sub,range_sub));
16			elseif strcmp('depth',varargin{(iSub-1)*2+1});
17			KK=find(prof_sub>=range_sub(1)&prof_sub<range_sub(2));
18			else;
19			KK=find(prof_sub>=range_sub(1)&prof_sub<range_sub(2));
20			end;
21			% [varargin{(iSub-1)*2+1} ' -- ' num2str(length(KK))]
22			for iFld=1:length(fldNames);
23			eval(['tmp1=MITprofSub.' fldNames{iFld} ';']);
24			if strcmp('depth',varargin{(iSub-1)*2+1}) & strcmp(fldNames{iFld},'prof_depth');
25			tmp1=tmp1(KK);
26			elseif strcmp('depth',varargin{(iSub-1)*2+1});
27			if size(tmp1,2)==length(MITprof.prof_depth); tmp1=tmp1(:,KK); end;
28			else;
29			if size(tmp1,1)==size(prof_sub,1); tmp1=tmp1(KK,:); end;
30			end;
31			eval(['MITprofSub.' fldNames{iFld} '=tmp1;']);
32			end;
33
34			end;
35
36			%add a couple things:
37			%--------------------
38			MITprofSub.np=length(MITprofSub.prof_lon);
39			MITprofSub.nr=length(MITprofSub.prof_depth);
40