function [UU,WW]=wing(S,theta,gR,gL,XX,ZZ)
[x,y,u,v]=vtheta();
xx=x+(S*cosd(theta)/2)*ones(size(x));
zz=y+(S*sind(theta)/2)*ones(size(y));
UR=interp2(xx,zz,u,XX,ZZ,'linear',0.0);
WR=interp2(xx,zz,v,XX,ZZ,'linear',0.0);
xx=x-(S*cosd(theta)/2)*ones(size(x));
zz=y+(S*sind(theta)/2)*ones(size(y));
UL=interp2(xx,zz,-u,XX,ZZ,'linear',0.0);
WL=interp2(xx,zz,-v,XX,ZZ,'linear',0.0);
UU=gR*UL+gL*UR;
WW=gR*WL+gL*WR;
VV=zeros(size(XX));
return