main_macrorugo.sce 1.65 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
17
ks = 0.01 // Rugosité de fond (m)
Cd0=1.5 // Forme
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)
18
19

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

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

// *****************************************************************************
printf("\n*** Submerged conditions ***\n")
// *****************************************************************************
38
k = 0.7
39
h = 0.8
40
Cd0 = 1.5
41
C = 0.13
42
macrorugo_searchQ(ks, D, k, Cd0, S, B, h, C, z_amont, long, bDbg)
43

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


macrorugo_searchQ(ks, D, 0.5, 2, S, 2, 0.85, C, z_amont, long, bDbg)