function [xd,zd,UU,WW]=vpair_plane_reference(wtv)


SS1=wtv.vortex.port.core;
SS2=wtv.vortex.starboard.core;
sc1=wtv.vortex.port.circulation;
sc2=wtv.vortex.starboard.circulation;
xd1=wtv.vortex.port.position(1);
zd1=wtv.vortex.port.position(2);
xd2=wtv.vortex.starboard.position(1);
zd2=wtv.vortex.starboard.position(2);
downwash=wtv.vortex.downwash;

xl=min(xd1,xd2);
xh=max(xd1,xd2);
zl=min(zd1,zd2);
zh=max(zd1,zd2);

dd=max(SS1,SS2);
cenx=(xl+xh)/2;
cenz=(zl+zh)/2;

[xd,zd]=meshgrid([xl-4*dd:3:xh+4*dd],[zl-4*dd:3:zh+4*dd]);

xdp=xd-xd1*ones(size(xd));
zdp=zd-zd1*ones(size(zd));
%Left vortex (clockwise rotation)
Xp=complex(xdp,zdp);
        rhop=abs(Xp);
        ang=angle(Xp);
        Vtheta=(1-exp(-1.26*((rhop/SS1).^2)))./(2*pi*rhop);
        up=Vtheta.*sin(ang);
        wp=-Vtheta.*cos(ang);
        
UUL=up;
WWL=wp;

xdp=xd-xd2*ones(size(xd));
zdp=zd-zd2*ones(size(zd));

Xp=complex(xdp,zdp);
        rhop=abs(Xp);
        ang=angle(Xp);
        Vtheta=(1-exp(-1.26*((rhop/SS2).^2)))./(2*pi*rhop);
        up=-Vtheta.*sin(ang);
        wp=Vtheta.*cos(ang);
        
UUR=up;
WWR=wp;

UU=sc1*UUL+sc2*UUR;

WW=sc1*WWL+sc2*WWR;

return