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