12.07.2015 Views

Quels apports hydrologiques pour les modèles hydrauliques? Vers ...

Quels apports hydrologiques pour les modèles hydrauliques? Vers ...

Quels apports hydrologiques pour les modèles hydrauliques? Vers ...

SHOW MORE
SHOW LESS
  • No tags were found...

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Annexe E. Algorithmique et codes C des modè<strong>les</strong> GR4H et Hayami 277où i est le numéro du pas de temps suivant le pas de temps présent. i varie entre 0 et n − 1.La clef de répartition est obtenue par l’équation suivante :UH(i) = SSH(i + 1) − SSH(i)La pluie répartie P UH <strong>pour</strong> le pas de temps présent est obtenue à partir des pluies P R des pasde temps précédents :n−1∑P UH = P R (k − i) UH(i + 1)i=0avec k le numéro du pas de temps présent.7. Soit V R (mm) le niveau de remplissage du réservoir de routage au début du pas de temps. Leterme d’échange est alors calculé par l’équation suivante :( ) 7/2 VREch = IGFRtel-00392240, version 1 - 5 Jun 20098. Ce terme ainsi qu’une partie de la pluie répartie P UH permettent de mettre à jour le contenudu réservoir de routage :V + RSi la quantité V + R est négative, on fixe V + R = 0.= V R + Ech + 0.9P UHLa contribution du réservoir de routage au débit final, notée Q R , est calculée de la manièresuivante :Q R = V + R⎛⎝1 −[1 +( V+ ) 4] ⎞ −1/4R ⎠RLe niveau du réservoir de routage est ensuite mis à jour <strong>pour</strong> le pas de temps suivant :V R = V + R − Q R9. Les 10% restants de la pluie répartie P UH sont routés vers l’exutoire en prélevant directementle terme d’échange :Q D = 0.1P UH + EchSi la quantité Q D est négative, on fixe Q D = 0.10. Le débit final est obtenu en sommant <strong>les</strong> deux contributions Q R et Q D :Q = Q R + Q D

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!