function [l]=hydrau(m,x,wh,wl,dcol)
[mwa,densa]=mwdens(x);
const=183.2*m*mwa/(densa*dcol*dcol)-wh/12;
if(const<=0)
    l=0;
else
    l=densa*wl*999*((183.2*m*mwa/(densa*dcol*dcol)-wh/12)^1.5)/mwa;
end