utils/matlab/cube2latlon.m

function [z] = cube2latlon(x,y,c,xi,yi)
% 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/models/MITgcmUV/utils/matlab/Attic/cube2latlon.m,v 1.1.2.1 2001/04/13 18:28:23 adcroft Exp $

[nx ny nz]=size(c);

for k=1:nz;
X=x;
Y=y;
C=c(:,:,k);

i=3*ny+(1:ny);j=floor(ny/2);
X(end+1,:)=X(i,j)'-360; Y(end+1,:)=Y(i,j)'; C(end+1,:)=C(i,j)';
i=3*ny+(1:ny);j=floor(ny/2)+1;
X(end+1,:)=X(i,j)'+360; Y(end+1,:)=Y(i,j)'; C(end+1,:)=C(i,j)';
i=5*ny+round(ny/2);j=1:floor(ny/2);
X(end+1,j)=X(i,j)+360;
Y(end+1,j)=Y(i,j);
C(end+1,j)=C(i,j);
i=5*ny+round(ny/2)+1;j=1:floor(ny/2);
X(end,j+ny/2)=X(i,j)-360;
Y(end,j+ny/2)=Y(i,j);
C(end,j+ny/2)=C(i,j);
i=2*32+(ny/2+1:ny);j=floor(ny/2);
X(end+1,1:ny/2)=X(i,j)'-360;
Y(end+1,1:ny/2)=Y(i,j)';
C(end+1,1:ny/2)=C(i,j)';
i=2*32+(ny/2+1:ny);j=floor(ny/2)+1;
X(end,ny/2+1:ny)=X(i,j)'+360;
Y(end,ny/2+1:ny)=Y(i,j)';
C(end,ny/2+1:ny)=C(i,j)';

z(:,:,k)=griddata(Y,X,C,yi,xi');
end % k
1	function [z] = cube2latlon(x,y,c,xi,yi)
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/models/MITgcmUV/utils/matlab/Attic/cube2latlon.m,v 1.1.2.1 2001/04/13 18:28:23 adcroft Exp $
18
19	[nx ny nz]=size(c);
20
21	for k=1:nz;
22	X=x;
23	Y=y;
24	C=c(:,:,k);
25
26	i=3*ny+(1:ny);j=floor(ny/2);
27	X(end+1,:)=X(i,j)'-360; Y(end+1,:)=Y(i,j)'; C(end+1,:)=C(i,j)';
28	i=3*ny+(1:ny);j=floor(ny/2)+1;
29	X(end+1,:)=X(i,j)'+360; Y(end+1,:)=Y(i,j)'; C(end+1,:)=C(i,j)';
30	i=5*ny+round(ny/2);j=1:floor(ny/2);
31	X(end+1,j)=X(i,j)+360;
32	Y(end+1,j)=Y(i,j);
33	C(end+1,j)=C(i,j);
34	i=5*ny+round(ny/2)+1;j=1:floor(ny/2);
35	X(end,j+ny/2)=X(i,j)-360;
36	Y(end,j+ny/2)=Y(i,j);
37	C(end,j+ny/2)=C(i,j);
38	i=2*32+(ny/2+1:ny);j=floor(ny/2);
39	X(end+1,1:ny/2)=X(i,j)'-360;
40	Y(end+1,1:ny/2)=Y(i,j)';
41	C(end+1,1:ny/2)=C(i,j)';
42	i=2*32+(ny/2+1:ny);j=floor(ny/2)+1;
43	X(end,ny/2+1:ny)=X(i,j)'+360;
44	Y(end,ny/2+1:ny)=Y(i,j)';
45	C(end,ny/2+1:ny)=C(i,j)';
46
47	z(:,:,k)=griddata(Y,X,C,yi,xi');
48	end % k