function [zi,di,fi]=zslice1(g,z,f,p,N);
%function [zi,di,fi]=zslice1(g,z,f,p,N);

x0=p(1);
y0=p(2);
x1=p(3);
y1=p(4);
xi=x0:(x1-x0)/N:x1;
yi=y0:(y1-y0)/N:y1;
N=length(xi);
[nn,nnv]=size(z);
fi=zeros(N,nnv);
zi=zeros(N,nnv);

fi=scattered_interpm(g.x,g.y,f,g.e,xi,yi);
zi=scattered_interpm(g.x,g.y,z,g.e,xi,yi);

d(1)=0;
for n=2:N,
    d(n)=d(n-1)+abs( xi(n)+i*yi(n)-xi(n-1)-i*yi(n-1) );
end
d=d(:);
di=d*ones(1,nnv);
