gmaze_pv/subfct/diagCatH.m

% Ch = diagCatH(C,depth,h)
%
% Get field C(depth,lat,lon) at depth h(lat,lon)
%
% depth < 0
% h     < 0 
%
% G. Maze, MIT, June 2007
%y

function varargout = diagCatH(C,Z,h)

% 0 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% PREPROC
[nz,ny,nx] = size(C);
Ch = zeros(ny,nx);

% 1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% COMPUTING
warning off
 for ix = 1 : nx
   for iy = 1 : ny
     Ch(iy,ix) = interp1( Z, squeeze(C(:,iy,ix)) , h(iy,ix) , 'linear');
   end
 end
warning on
 
% 2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% OUTPUTS
switch nargout
 case 1
  varargout(1) = {Ch};
end
1	gmaze	1.1	% Ch = diagCatH(C,depth,h)
2			%
3			% Get field C(depth,lat,lon) at depth h(lat,lon)
4			%
5			% depth < 0
6			% h < 0
7			%
8			% G. Maze, MIT, June 2007
9			%y
10
11			function varargout = diagCatH(C,Z,h)
12
13			% 0 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% PREPROC
14			[nz,ny,nx] = size(C);
15			Ch = zeros(ny,nx);
16
17			% 1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% COMPUTING
18			warning off
19			for ix = 1 : nx
20			for iy = 1 : ny
21			Ch(iy,ix) = interp1( Z, squeeze(C(:,iy,ix)) , h(iy,ix) , 'linear');
22			end
23			end
24			warning on
25
26			% 2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% OUTPUTS
27			switch nargout
28			case 1
29			varargout(1) = {Ch};
30			end