--- MITgcm_contrib/osse/utils/cyl2cart.m 2004/07/12 21:43:16 1.2 +++ MITgcm_contrib/osse/utils/cyl2cart.m 2005/01/13 15:50:36 1.3 @@ -17,7 +17,6 @@ [theta rho nz]=size(c); [RHO,THETA] = meshgrid(1:rho,-pi+2*pi/theta:2*pi/theta:pi); [x,y] = pol2cart(THETA,RHO); -%[nx ny nz]=size(c); nx=theta;ny=rho; X=reshape(x,[1 nx*ny]); @@ -27,7 +26,6 @@ for k=1:nz; C=reshape(c(:,:,k),[1 nx*ny]); z(:,:,k)=griddata(Y,X,C,yi,xi',varargin{:}); -% z(:,:,k)=griddata_fast(del,[C C(il) C(ig)],varargin{:}); end % k % Split vertical and time dimensions