top of page

MatLab Code Showcase

per=extremewave(:,3);
dpt=extremewave(:,7);
length0=extremewave(:,6);
angle=90-extremewave(:,4);
h=extremewave(:,5);
d=1;l=431;i=1;length=[];alpha=[];
for i=1:l
    for l=1:100
    if abs(d-l)<1
       
        d=1.56*per(i)^2*tanh(2*3.14/l*dpt(i));
        length(i)=d;
    end
    end
end
length=reshape(length,431,1);
alpha=asin(sin(angle/180*3.14)).*length./length0;
Q=237600.*h.^2.5.*sin(2*alpha);
sedimenttransport=Q;

© 2023 by Zhoufeng Dai and Boyuan Lu. Proudly created with Wix.com

bottom of page