[fix][mascaret] Bazin's friction law corrected

NEWS.txt

 Latest changes
 ==============
+MASCARET: Fix Bazin's friction Law
 
 KHIONE: Addition of a new formula to compute under ice cover friction
 A new default value of 0.02 m1/3/s is set for under cover friction

sources/mascaret/Modules/m_debitance_s.f90

@@ -104,7 +104,7 @@ module M_DEBITANCE_S
       !----------------------------------------------------
       case(LOI_FROTTEMENT_BAZIN)
 
-         chezy = 45._DOUBLE / ( 1._DOUBLE + CF1 / RH1 )
+         chezy = 87._DOUBLE / ( 1._DOUBLE + CF1 / RH1**W12 )
 
    end select
 

sources/mascaret/ModulesDeriv/m_debitance_s_b.f90

@@ -90,7 +90,7 @@ CONTAINS
 !----------------------------------------------------
 ! BAZIN   Chezy(i) en fonction de RH et mb (fixe)
 !----------------------------------------------------
-        chezy = 45._DOUBLE/(1._DOUBLE+cf1/rh1)
+        chezy = 87._DOUBLE  / ( 1._DOUBLE + CF1 / RH1**W12 )
         CALL PUSHCONTROL3B(4)
       CASE DEFAULT
         CALL PUSHCONTROL3B(0)

sources/mascaret/ModulesDeriv/m_debitance_s_d.f90

@@ -89,9 +89,8 @@ MODULE M_DEBITANCE_S_D
 !----------------------------------------------------
 ! BAZIN   Chezy(i) en fonction de RH et mb (fixe)
 !----------------------------------------------------
-        chezyd = -(45._DOUBLE*(cf1d*rh1-cf1*rh1d)/rh1**2/(1._DOUBLE+cf1/&
-&         rh1)**2)
-        chezy = 45._DOUBLE/(1._DOUBLE+cf1/rh1)
+          chezyd = 43.5_DOUBLE*(cf1*rh1d-2*cf1d*rh1)/(2._DOUBLE*rh1**W12*(cf1+rh1**W12)**2)
+          chezy  = 87._DOUBLE/(1._DOUBLE+cf1/rh1**w12)
       CASE DEFAULT
         chezyd = 0.D0
       END SELECT