function h=plotell(x2,y2,z,fac,isub);
%function h=plotell(x2,y2,z,fac,isub);
[m,n,p]=size(z);
zs=z(1:isub:m,1:isub:n,:);
xs=x2(1:isub:m,1:isub:n);
ys=y2(1:isub:m,1:isub:n);
[m,n,p]=size(zs);
x3=repmat(xs,[1 1 p])+real(zs)*fac;
y3=repmat(ys,[1 1 p])+imag(zs)*fac;
x3=reshape(x3,m*n,p).';
y3=reshape(y3,m*n,p).';
h=line(x3,y3);
