function f0=hslice_1(z,f,z0);
%function f0=hslice(z,f,z0);

[nn,nv]=size(f);
f0=zeros(nn,1);
for k=1:nn
    f0(k)=interp1(z(k,:),f(k,:),z0);
end
