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

280 Annexe E. Algorithmique et codes C des modèles GR4H et Hayamir e t u r n f ;}double INISTATE ( i n t nval , double *mPmE, double XS){i n t i ;double s t a r t , end ;// I t e r a t i o n to f i n d the s o l u t i o n o f the e q u a t i o n p r o v i d e d// by Le Moine (2008)s t a r t =0; end =1; i =0;w h i l e ( end−s t a r t >1e−3 && i 0){s t a r t +=(end−s t a r t ) / 2 ;}e l s e {end−=(end−s t a r t ) / 2 ; }tel-00392240, version 1 - 5 Jun 2009i ++;}r e t u r n s t a r t ;}/* ************** Core s u b r o u t i n e ****************************** */v o i d c GR4H ( i n t nval ,double * Xinp ,double * Pinp , double * Einp , double * mPmE,double * Qcalc ,double * Scalc , double * Rcalc , double * ECHcalc , double * ERcalc ){i n t i , l , nH , nHini , nLAC , problemUH ;double P , E ,Q, S i n i , R i n i ;double V[ XSIZE ] , X[ XSIZE ] ;double EN, ES , PS , PR,WS, S2 , PERC, ECH, TP, R2 ,QR,QD,PERCFACTOR, Suh ;double ech1 , ech2 , Expos ;/*−−−− Exponent o f the UH −−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−*/Expos =1.25;/* −−− P e r c o l a t i o n f a c t o r f o r the h o u r l y time−s t e p −−−−−− */PERCFACTOR=4;