Diagnostics/DiagUtility/cube2latlon.m

function [z] = cube2latlon(x,y,c,xi,yi,varargin)
% z=cube2latlon(x,y,c,xi,yi);
%
% Re-grids model output on expanded spherical cube to lat-lon grid.
%  x,y   are 2-D arrays of the cell-centered coordinates 
%  c     is a 2-D or 3-D scalar field
%  xi,yi are vectors of the new regular lat-lon grid to interpolate to.
%  z     is the interpolated data with dimensions of size(xi) by size(yi).
%
% e.g.
% >> x=rdmds('XC');
% >> y=rdmds('YC');
% >> t=rdmds('Ttave.0000513360');
% >> xi=-179:2:180;yi=-89:2:90;
% >> ti=cube2latlon(x,y,t,xi,yi);
%
% $Header: /u/gcmpack/MITgcm/utils/matlab/cube2latlon.m,v 1.4 2004/06/04 17:03:50 adcroft Exp $

NN=size(c);
[nx ny nz]=size(c);

X=reshape(x,[1 nx*ny]);
Y=reshape(y,[1 nx*ny]);
ig=find(X>90);
il=find(X<-90);
del=griddata_preprocess([Y Y(il) Y(ig)],[X X(il)+360 X(ig)-360],yi,xi',varargin{:});

for k=1:nz;
 C=reshape(c(:,:,k),[1 nx*ny]);
%z(:,:,k)=griddata([Y Y(il) Y(ig)],[X X(il)+360 X(ig)-360],[C C(il) C(ig)],yi,xi',varargin{:});
 z(:,:,k)=griddata_fast(del,[C C(il) C(ig)],varargin{:});
end % k

% Split vertical and time dimensions
if size(NN,2)>2
z=reshape(z,[size(z,1) size(z,2) NN(3:end)]);
end
1	enderton	1.1	function [z] = cube2latlon(x,y,c,xi,yi,varargin)
2			% z=cube2latlon(x,y,c,xi,yi);
3			%
4			% Re-grids model output on expanded spherical cube to lat-lon grid.
5			% x,y are 2-D arrays of the cell-centered coordinates
6			% c is a 2-D or 3-D scalar field
7			% xi,yi are vectors of the new regular lat-lon grid to interpolate to.
8			% z is the interpolated data with dimensions of size(xi) by size(yi).
9			%
10			% e.g.
11			% >> x=rdmds('XC');
12			% >> y=rdmds('YC');
13			% >> t=rdmds('Ttave.0000513360');
14			% >> xi=-179:2:180;yi=-89:2:90;
15			% >> ti=cube2latlon(x,y,t,xi,yi);
16			%
17			% $Header: /u/gcmpack/MITgcm/utils/matlab/cube2latlon.m,v 1.4 2004/06/04 17:03:50 adcroft Exp $
18
19			NN=size(c);
20			[nx ny nz]=size(c);
21
22			X=reshape(x,[1 nx*ny]);
23			Y=reshape(y,[1 nx*ny]);
24			ig=find(X>90);
25			il=find(X<-90);
26			del=griddata_preprocess([Y Y(il) Y(ig)],[X X(il)+360 X(ig)-360],yi,xi',varargin{:});
27
28			for k=1:nz;
29			C=reshape(c(:,:,k),[1 nx*ny]);
30			%z(:,:,k)=griddata([Y Y(il) Y(ig)],[X X(il)+360 X(ig)-360],[C C(il) C(ig)],yi,xi',varargin{:});
31			z(:,:,k)=griddata_fast(del,[C C(il) C(ig)],varargin{:});
32			end % k
33
34			% Split vertical and time dimensions
35			if size(NN,2)>2
36			z=reshape(z,[size(z,1) size(z,2) NN(3:end)]);
37			end