macrorugo_resultats_complementaires.sci 1.22 KB
Newer Older
1
function macrorugo_resultComp(z_amont, S, long, Q, L, pf, C, Cd0, h, D)
2
    g = 9.81
3

4
5
    ZF2=z_amont-S*long;
    print_r("ZF2");
6
7
8
9
10
11

    Vdeb=Q/L/pf;
    print_r("Vdeb");


    Vg=Q/L/pf/(1-C^0.5);
12
    print_r("Vg");
13
14
15
    Fr=Vg./(g.*pf).^0.5;
    print_r("Fr");

16
17
    PV=1000*g*Vdeb*S;
    print_r("PV");
18
19
20

    if pf/h<1
        flowcond = 'emergent'
21
    elseif pf/h<1.1 & pf/h>=1
22
        flowcond = 'quasi emergent'
23
    else
24
25
26
        flowcond = 'immerge'
    end
    print_r("flowcond");
27

28
    if pf/h>1.1
29
        q_technique=0.955*(pf/h)^2.282*S^0.466*C^(-0.23)*(9.81*h)^0.5.*h*L;
30
    else
31
        if Cd0==2
32
33
34
35
36
37
            q_technique= 0.648*(pf/D)^1.084*S^0.56*C^(-0.456)*(9.81*D)^0.5.*D*L;
            V_technique=3.35*(pf/D)^0.27*S^0.53*(9.81*D)^0.5;
        else
            q_technique=0.815*(pf/D)^1.45*S^0.557*C^(-0.456)*(9.81*D)^0.5.*D*L;
            V_technique=4.54*(pf/D)^0.32*S^0.56*(9.81*D)^0.5;
        end
38
39
40
41
42
        if Cd0 > 2
            coeff_contraction=1.5;
        else
            coeff_contraction=1.1;
        end
43
        Vmax=Vg.*coeff_contraction*calcfFr(Fr);
44
        print_r("Vmax");
45
46
47
48
        print_r("V_technique");
    end
    print_r("q_technique");

49
50
51
    Strickler = Q /(pf*L*pf^(2/3)*S^(1/2));
    print_r("Strickler");

52
endfunction