[thetav,rhov] = cart2cyluv(u(:,:,15),v(:,:,15),thetai,rhoi);
figure(1);
imagesc(rhov);
%imagesc(v(:,:,15));
colorbar;
title('thetav (u)');

figure(2);
imagesc(V(:,:,15));
colorbar;
title('rhov (v)');
title('V');

figure(3)
imagesc(rhov-V(:,:,15))
colorbar;
title('err');