function [x,z,u,w]=vthetaR(R)

[x,z]=meshgrid([-750.5:1:750.5],[-750.5:1:750.5]);
Xp=complex(x,z);
        rhop=abs(Xp);
        ang=angle(Xp);
        Vtheta=(1-exp(-1.26*((rhop/R).^2)))./(2*pi*rhop);
        u=-Vtheta.*sin(ang);
        w=Vtheta.*cos(ang);
return
end
        