/[MITgcm]/MITgcm_contrib/osse/utils/cyl2cart.m
ViewVC logotype

Diff of /MITgcm_contrib/osse/utils/cyl2cart.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by afe, Mon Jul 12 21:43:16 2004 UTC revision 1.3 by afe, Thu Jan 13 15:50:36 2005 UTC
# Line 17  NN=size(c); Line 17  NN=size(c);
17  [theta rho nz]=size(c);  [theta rho nz]=size(c);
18  [RHO,THETA] = meshgrid(1:rho,-pi+2*pi/theta:2*pi/theta:pi);  [RHO,THETA] = meshgrid(1:rho,-pi+2*pi/theta:2*pi/theta:pi);
19  [x,y] = pol2cart(THETA,RHO);  [x,y] = pol2cart(THETA,RHO);
 %[nx ny nz]=size(c);  
20  nx=theta;ny=rho;  nx=theta;ny=rho;
21    
22  X=reshape(x,[1 nx*ny]);  X=reshape(x,[1 nx*ny]);
# Line 27  del=griddata_preprocess(Y,X,yi,xi',varar Line 26  del=griddata_preprocess(Y,X,yi,xi',varar
26  for k=1:nz;  for k=1:nz;
27   C=reshape(c(:,:,k),[1 nx*ny]);   C=reshape(c(:,:,k),[1 nx*ny]);
28  z(:,:,k)=griddata(Y,X,C,yi,xi',varargin{:});  z(:,:,k)=griddata(Y,X,C,yi,xi',varargin{:});
 % z(:,:,k)=griddata_fast(del,[C C(il) C(ig)],varargin{:});  
29  end % k  end % k
30    
31  % Split vertical and time dimensions  % Split vertical and time dimensions
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22