gmaze_pv/subfct/subfct_getisoS.m

%
% [Iiso MASK] = SUBFCT_GETISOS(TRACER,ISO)
%
% This function determines the iso-surface ISO of the 
% 3D field TRACER(Z,Y,X)
% Iiso(Y,X) contains indices of axis Z of the iso-surface.
% MASK(Y,X) contains 1 where the iso-surface is defined
% and 0 elsewhere. Useful to plot other fields.
%
% 06/19/06
% gmaze@mit.edu

function varargout = subfct_getisoS(CHP,ISO)


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% PRE-PROCESS and ERROR CHECK %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
pv_checkpath

if nargin ~= 2
  help subfct_getisoS.m
  error('subfct_getisoS.m : Wrong number of parameters')
  return
end

[nz ny nx] = size(CHP);
if nz ==1 | nx == 1 | ny == 1
  help subfct_getisoS.m
  error('subfct_getisoS.m : 3DFIELD must be a 3D field')
  return
end


%%%%%%%%%%%%%%%%%%%%%%%%%
% DETERMINE ISO-SURFACE %
%%%%%%%%%%%%%%%%%%%%%%%%%

mask = ones(ny,nx).*NaN;
%isoS = ones(ny,nx).*NaN;
Iiso = ones(ny,nx).*NaN;

for ix = 1 : nx
  for iy = 1 : ny    
    p = squeeze(CHP(:,iy,ix));
    p = find(p<=ISO);
    if isempty(p)
%      isoS(iy,ix) = NaN;
      Iiso(iy,ix) = NaN;
    else
%      isoS(iy,ix) = CHP(max(p),iy,ix);
      Iiso(iy,ix) = max(p);
      mask(iy,ix) = 1;
    end %if
  end %for ix
end %for iy


%%%%%%%%%%%
% OUTPUTS %
%%%%%%%%%%%
switch nargout
 case {0,1}
  varargout(1) = {Iiso};
 case 2
  varargout(1) = {Iiso};
  varargout(2) = {mask};
end %switch nargout


1	gmaze	1.1	%
2			% [Iiso MASK] = SUBFCT_GETISOS(TRACER,ISO)
3			%
4			% This function determines the iso-surface ISO of the
5			% 3D field TRACER(Z,Y,X)
6			% Iiso(Y,X) contains indices of axis Z of the iso-surface.
7			% MASK(Y,X) contains 1 where the iso-surface is defined
8			% and 0 elsewhere. Useful to plot other fields.
9			%
10			% 06/19/06
11			% gmaze@mit.edu
12
13			function varargout = subfct_getisoS(CHP,ISO)
14
15
16			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
17			% PRE-PROCESS and ERROR CHECK %
18			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
19			pv_checkpath
20
21			if nargin ~= 2
22			help subfct_getisoS.m
23			error('subfct_getisoS.m : Wrong number of parameters')
24			return
25			end
26
27			[nz ny nx] = size(CHP);
28			if nz ==1 \| nx == 1 \| ny == 1
29			help subfct_getisoS.m
30			error('subfct_getisoS.m : 3DFIELD must be a 3D field')
31			return
32			end
33
34
35
36			%%%%%%%%%%%%%%%%%%%%%%%%%
37			% DETERMINE ISO-SURFACE %
38			%%%%%%%%%%%%%%%%%%%%%%%%%
39
40			mask = ones(ny,nx).*NaN;
41			%isoS = ones(ny,nx).*NaN;
42			Iiso = ones(ny,nx).*NaN;
43
44			for ix = 1 : nx
45			for iy = 1 : ny
46			p = squeeze(CHP(:,iy,ix));
47			p = find(p<=ISO);
48			if isempty(p)
49			% isoS(iy,ix) = NaN;
50			Iiso(iy,ix) = NaN;
51			else
52			% isoS(iy,ix) = CHP(max(p),iy,ix);
53			Iiso(iy,ix) = max(p);
54			mask(iy,ix) = 1;
55			end %if
56			end %for ix
57			end %for iy
58
59
60
61			%%%%%%%%%%%
62			% OUTPUTS %
63			%%%%%%%%%%%
64			switch nargout
65			case {0,1}
66			varargout(1) = {Iiso};
67			case 2
68			varargout(1) = {Iiso};
69			varargout(2) = {mask};
70			end %switch nargout
71
72