function z=getdepth(cdf,sta);
%function z=getdepth(cdf,sta);
d=mcvgt(cdf,'depth',sta-1,1);
z=mcvgt(cdf,'sigma');
n=length(z);
z(n+1)=-1;
z=.5*(z(1:n)+z(2:n+1));
z=z*d;
