function s = slowfunctionsmooth(depth)
v=(3.0+sqrt(depth)/sqrt(1000))*1000;
%true 1-d slowness structure
s=1./v;


