remous_rect_trapezes_pentefaible.spec.ts

import { ParamsSectionRectang, cSnRectang } from "../src/section/section_rectang";
import { CourbeRemousParams, MethodeResolution, CourbeRemous } from "../src/remous";
import { precDigits, precDist, equalEpsilon, compareArray, compareObject, compareLog } from "./nubtest";
import { round } from "../src/base";
import { cLog } from "../src/util/log";
import { Message, MessageCode } from "../src/util/message";

/*
  Le code de modification des lignes fluviale et torrentielle a été modifié, on enlève un point de plus
  ligne d'eau complète : if (iSens * (rXCC - rX) < 0) remplacé par if (iSens * (rXCC - rX) <= 0)
  lign d'eau partielle : if (iSens * (rXCN - xRst) > 0)  remplacé par if (iSens * (rXCN - xRst) >= 0)
  Les données de validation ont été modifiées en conséquence
*/

describe('Class Remous / section rectangulaire :', () => {
    describe('méthode trapèzes :', () => {
        it("faible pente, ressaut avant l'amont", () => {
            let prms = new ParamsSectionRectang(undefined, // tirant d'eau
                2.5, // largeur de fond
                40, //  Ks=Strickler
                2, // Q=Débit
                0.001, // If=pente du fond
                precDist, // précision
                1 // YB=hauteur de berge
            );

            let sect = new cSnRectang(prms);

            let prem = new CourbeRemousParams(sect, 0.3, // Yamont = tirant amont
                0.403, // Yaval = tirant aval
                100,  // Long= Longueur du bief
                5,  // Dx=Pas d'espace
                MethodeResolution.Trapezes
            );

            let rem = new CourbeRemous(prem);

            let res = rem.calculRemous(undefined);

            // données de validation : version PHP (Oct 2017) méthode trapèzes

            let f = { 100.000: 0.403, 95.000: 0.524, 90.000: 0.558, 85.000: 0.584, 80.000: 0.604, 75.000: 0.621, 70.000: 0.637, 65.000: 0.65, 60.000: 0.662, 55.000: 0.673, 50.000: 0.684, 45.000: 0.693, 40.000: 0.701, 35.000: 0.709, 30.000: 0.717, 25.000: 0.725, 20.000: 0.731, 15.000: 0.738, 10.000: 0.744, 5.000: 0.75, 0.000: 0.755 };
            compareObject("Yfluvial", res["flu"], f, 0.03);

            let t = {};
            compareObject("Ytorrentiel", res["tor"], t, 0.03);

            let x = [0.000, 5.000, 10.000, 15.000, 20.000, 25.000, 30.000, 35.000, 40.000, 45.000, 50.000, 55.000, 60.000, 65.000, 70.000, 75.000, 80.000, 85.000, 90.000, 95.000, 100.000];
            compareArray("abscisses", res["trX"], x);

            let expLog = new cLog();
            let m = new Message(MessageCode.ERROR_REMOUS_LARGEUR_BERGE);
            m.extraVar["B"] = 2.5;
            expLog.add(m);

            m = new Message(MessageCode.ERROR_REMOUS_H_CRITIQUE);
            m.extraVar["Yc"] = 0.403;
            expLog.add(m);

            m = new Message(MessageCode.ERROR_REMOUS_H_NORMALE);
            m.extraVar["Yn"] = 0.953;
            expLog.add(m);

            expLog.add(new Message(MessageCode.ERROR_REMOUS_CALCUL_FLUVIAL));

            expLog.add(new Message(MessageCode.ERROR_REMOUS_CALCUL_TORRENTIEL));

            m = new Message(MessageCode.ERROR_REMOUS_ARRET_CRITIQUE)
            m.extraVar["x"] = 5;
            expLog.add(m);