blob: 4f8e9641f64b14618cfd14b22d114b014a87947e (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
//example 4.28
//calculate ordinates of storm hydrograph
clc;funcprot(0);
//given
fi=2.5; //infiltration index
B=10; //Base flow
O=[0 110 365 500 390 310 250 235 175 130 95 65 40 22 10 0 0 0]; //ordinates of unit hydrograph
R1=2;R2=6.75;R3=3.75;
r1=(R1*10-(fi*3)-5)/10; //rainfall excess in first three hour
r2=(R2*10-(fi*3))/10; //rainfall excess in second three hour
r3=(R3*10-(fi*3))/10; //rainfall excess in third three hour
for i=1:18
s1(i)=r1*O(i);
end
for i=2:18
s2(i)=r2*O(i-1);
end
for i=3:18
s3(i)=r3*O(i-2);
end //surface run-off from rainfall excess during succesive unit periods
mprintf("ordinates of storm hydrograph");
for i=1:18
T(i)=s1(i)+s2(i)+s3(i);
t(i)=T(i)+B;
t(i)=round(t(i)*10)/10;
mprintf("\n%f",t(i));
end
|
