main_macrorugo.sce 2.51 KB
Newer Older
1
2
3
4
5
// Test des passes à macro rugosité pour JaLHyd
clear
sCurrentPath = get_absolute_file_path("main_macrorugo.sce");
getd(sCurrentPath);

6
bDbg = %f;
7
bDbg_RA = %f;
8

9
// Tests parameters
10
11
12
13
14
15
16
ks = 0.01 // Rugosité de fond (m)
S = 0.05 // Pente
C = 0.13; // Concentration
D = 0.5; // Diamètre
B = 1 // Largeur
z_amont = 12.5; // Cote amont (m)
long = 6; // Longueur rampe (m)
17
18

// *****************************************************************************
19
printf("\n*** Emergent conditions Cd=1.1***\n")
20
21
// *****************************************************************************
h = 0.6
22
k = 0.7
23
Cd0 = 1.1 // Forme ronde
24
macrorugo_searchQ(ks, D, k, Cd0, S, B, h, C, z_amont, long, bDbg)
25

26
// *****************************************************************************
27
printf("\n*** Emergent conditions Cd=2.6***\n")
28
29
30
// *****************************************************************************
h = 0.6
k = 0.7
31
Cd0 = 2.6 // Forme plane
32
macrorugo_searchQ(ks, D, k, Cd0, S, B, h, C, z_amont, long, bDbg)
33
34

// *****************************************************************************
35
printf("\n*** Submerged conditions  Cd=1.1***\n")
36
// *****************************************************************************
37
k = 0.7
38
h = 0.8
39
Cd0 = 1.1 // Forme ronde
40
C = 0.13
41
macrorugo_searchQ(ks, D, k, Cd0, S, B, h, C, z_amont, long, bDbg)
42

43
44
45
// *****************************************************************************
printf("\n*** JalHyd #85 ***\n")
// *****************************************************************************
46
Cd0 = 1.1
47
48
k = 0.8
C = 0.2; // Concentration
49
Q = []
50
for h = 0.7:0.1:1.2
51
    Q = [Q, macrorugo_searchQ(ks, D, k, Cd0, S, B, h, C, z_amont, long, bDbg)]
52
end
53
disp Q
54

55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
macrorugo_searchQ(ks, D, 0.5, 2.6, S, 2, 0.85, C, z_amont, long, bDbg)

// *****************************************************************************
printf("\n*** JalHyd #144 ***\n")
// rZF1, rL, rB,  rIf, rQ,   rY,  rRF,  rCB, rPBD, rPBH, rCd0
//  100, 25,  5, 0.05,  1, 0.05, 0.15, 0.13,  0.5,  0.4,  1.1
// *****************************************************************************
// macrorugo_searchQ(ks, D, k, Cd0, S, B, h, C, z_amont, long, bDbg)
//macrorugo_searchQ(0.15, 0.5, 0.4, 1.1, 0.05, 5, 0.05, 0.13, z_amont, long, bDbg)

// *****************************************************************************
printf("\n*** JalHyd #154 ***\n")
// *****************************************************************************
//Q = []
//for h = 0.34:0.01:0.41
//    Q = [Q, macrorugo_searchQ(0.1, 0.4, 0.4, 1, S, 1, h, 0.13, z_amont, long, bDbg)]
//end
//print_r("Q")