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

286 Annexe E. Algorithmique et codes C des modèles GR4H et Hayami* LTRO [ 0 ] = L (m) Longueur du t r o n ç o n** Moussa , R . ( 1 9 9 6 ) . ” A n a l y t i c a l Hayami s o l u t i o n f o r the* d i f f u s i v e wave f l o o d r o u t i n g problem with l a t e r a l i n f l o w . ”* H y d r o l o g i c a l P r o c e s s e s 1 0 ( 9 ) : 1209 −1227.** I n s p i r e d from FORTRAN code w r i t t e n by Claude Michel*************************************************/tel-00392240, version 1 - 5 Jun 2009/* ************** Hayami k e r n e l s u b r o u t i n e s *********************I n t e g r a t i o n o f hayami k e r n e l :1 . Each time s t e p i s d i v i d e d i n TSdiv2 . Within each s u b d i v i s i o n , F u n c t i o n i s e v a l u a t e d at 5 p o i n t s3 . I n t e r g r a l i s c a l c u l a t e d by a v e r a g i n g t h e s e 5 e v a l u a t i o n sTH, Z = Hayami model p a r a m e t e r sj = time s t e p*/double KERHAY INTEG UH( double *PAS , i n t TSdiv , double TH, double Z , i n t j ){i n t l , p ;double u [ 5 ] , v [ 5 ] , x , x l ,CONVOL,SCONVOL, t , ordUH , absUH ;// D e f i n i t i o n o f i n t e g r a t i o n c o n s t a n t s u−>a b s c i s s a e v−>weightu [0]= −0.90618; u [1]= −0.53847; u [ 2 ] = 0 ; u [ 3 ] = 0 . 5 3 8 4 7 ; u [ 4 ] = 0 . 9 0 6 1 8 ;v [ 0 ] = 0 . 2 3 6 9 3 ; v [ 1 ] = 0 . 4 7 8 6 3 ; v [ 2 ] = 0 . 5 6 8 8 8 ; v [ 3 ] = 0 . 4 7 8 6 3 ; v [ 4 ] = 0 . 2 3 6 9 3 ;// V e r i f i c a t i o ni f (TH>MAXTH){TH=MAXTH; }i f ( j 0i f (ZMINZ){i f (TH>0){