--- MITgcm_contrib/osse/utils/cart2cyl.m	2004/07/12 23:24:35	1.1
+++ MITgcm_contrib/osse/utils/cart2cyl.m	2004/07/13 16:26:51	1.2
@@ -20,6 +20,7 @@
 %[x,y] = pol2cart(THETA,RHO);
 [theta,rho] = cart2pol(X,Y);
 theta=(theta./(pi*2)+0.5).*size(thetai,2);
+rho=max(rhoi)*rho/((x-1)/2);
 %theta=rot90(theta.*180./pi,-1);
 %[nx ny nz]=size(c);
 %nx=theta;ny=rho;