%function f = sismo_Stoneley_disp_left(k,sigma,g,alphaw,h)
function f = sismo_Stoneley_disp_left(k,sigma,alphaw,h)
% left-hand side of the Stoneley dispersion relation
Q=0; %g/(2*alphaw);


l=sqrt(-1)*Q+sqrt(sigma.^2./alphaw.^2-k.^2);
f=tan(l.*h);
end

