--- MITgcm_contrib/osse/utils/cyl2cartuv.m	2004/07/12 21:43:16	1.2
+++ MITgcm_contrib/osse/utils/cyl2cartuv.m	2005/01/13 15:50:36	1.3
@@ -12,7 +12,7 @@
 % >> xi=-179:2:180;yi=-89:2:90;
 % >> ti=cyl2cart(x,y,t,xi,yi);
 %
-% $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/osse/utils/cyl2cartuv.m,v 1.2 2004/07/12 21:43:16 afe Exp $
+% $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/osse/utils/cyl2cartuv.m,v 1.3 2005/01/13 15:50:36 afe Exp $
 
 if ~isequal(size(thetav),size(rhov))
   error('Theta and rho vector arrays must be same size');
@@ -23,19 +23,11 @@
 [ntheta nrho nz]=size(thetav);
 [RHO,THETA,NZ] = meshgrid(1:nrho,-pi+2*pi/ntheta:2*pi/ntheta:pi,1:nz);
 [x,y] = pol2cart(THETA(:,:,1),RHO(:,:,1));
-%[nx ny nz]=size(c);
 nx=ntheta;ny=nrho;
 
 % break out components
-%vv=-(thetav.*-cos(THETA)+rhov.*sin(THETA));
-%uv=-(thetav.*-sin(THETA)+rhov.*cos(THETA));
 uv=thetav.*cos(THETA)+rhov.*sin(THETA);
 vv=thetav.*-sin(THETA)+rhov.*cos(THETA);
-%uv=thetav.*cos(THETA); %+rhov.*sin(THETA);
-%vv=thetav.*-sin(THETA); %+rhov.*cos(THETA);
-%uv=rhov.*sin(THETA);
-%vv=rhov.*cos(THETA);
-
 
 X=reshape(x,[1 nx*ny]);
 Y=reshape(y,[1 nx*ny]);
@@ -46,7 +38,6 @@
  VV=reshape(vv(:,:,k),[1 nx*ny]);
 u(:,:,k)=griddata(Y,X,UV,yi,xi',varargin{:});
 v(:,:,k)=griddata(Y,X,VV,yi,xi',varargin{:});
-% z(:,:,k)=griddata_fast(del,[C C(il) C(ig)],varargin{:});
 end % k
 
 % Split vertical and time dimensions