function s = meyer_peter_muller(Acal,D50,delt,g,mu,shld,xi,shld_cr)
s = ( Acal *8 .* D50 .* sqrt(delt .* g .* D50).*(max(mu.*shld - xi.*shld_cr,0)) .^ (3/2) );