function [Vrm,scale,ee]=vpairR(x,y,u,w,pd1,zd1,pd2,zd2,pd,zd,az,el,Vreff,Ratio)

xx=x+pd2*ones(size(x));
yy=y+zd2*ones(size(y));
UUR=interp2(xx,yy,u,pd,zd,'linear',0.0);
WWR=interp2(xx,yy,w,pd,zd,'linear',0.0);
xx=x+pd1*ones(size(x));
yy=y+zd1*ones(size(y));
UUL=interp2(xx,yy,-u,pd,zd,'linear',0.0);
WWL=interp2(xx,yy,-w,pd,zd,'linear',0.0);
UU=UUL+Ratio*UUR;
WW=WWL+Ratio*WWR;
Vrm=sind(az).*cosd(el).*UU+sind(el).*WW;
in=find(~isnan(Vreff)&abs(Vrm)>0.001);
%scale=mean(abs(Vreff(in))./abs(Vrm(in)));
ee=pcorr(Vreff(in),Vrm(in));
scale=std(Vreff(in))/std(Vrm(in));
%ee=sqrt(sum((Vreff(in)-Vrm(in)*scale).^2)/length(Vrm(in)));
return
end