summaryrefslogtreecommitdiff
path: root/2087/CH2/EX2.3/example2_3.sce
blob: fcb14146a0c6958db52e434480d71edb12c77396 (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
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
//example 2.3            

//calculate 
//time of water application
//optimum length of each border strip
//dischrge for each border strip

clc;
//Given
d=0.05;//depth of root zone
I=1.25D-5;//average infiltration rate
s=0.0035//slope of border strip
t=d/(I*3600);
t=round(t*1000)/1000;
mprintf("Time of water application=%f hours.",t);

//Part (a)
q=2D-3;//discharge entering water source
qdash=q*100^2*60;
n=0.55425-(0.0001386*qdash);
yo=(n*q/(s^0.5))^0.6;
y=0.665*yo;
L=(q/I)*(1-%e^(-d/y));
L=round(10*L)/10;
mprintf("\nPart (a):");
mprintf("\nOptimum length of each border strip=%f m.",L);

//Part (b)
Lgiven=150//given value of length
//First Trial
q=3D-3;
qdash=q*100^2*60;
n=0.55425-(0.0001386*qdash);
yo=(n*q/(s^0.5))^0.6;
y=0.665*yo;
L=(q/I)*(1-%e^(-d/y));
//second trial
q=3.15D-3;
qdash=q*100^2*60;
n=0.55425-(0.0001386*qdash);
yo=(n*q/(s^0.5))^0.6;
y=0.665*yo;
L=(q/I)*(1-%e^(-d/y));
q=9*Lgiven*q*1000/L;
q=round(q*10)/10;
mprintf("\nPart (b):");
mprintf("\nDischarge for each border strip=%f lps.",q);