From 97aa6891a05ddea41266309fb5ed37423f1b1ca2 Mon Sep 17 00:00:00 2001
From: Johan Fourdrinoy <johan.fourdrinoy@edf.fr>
Date: Thu, 27 Mar 2025 15:00:07 +0100
Subject: [PATCH] [fix][mascaret] Bazin's friction law corrected
---
NEWS.txt | 1 +
sources/mascaret/Modules/m_debitance_s.f90 | 2 +-
sources/mascaret/ModulesDeriv/m_debitance_s_b.f90 | 2 +-
sources/mascaret/ModulesDeriv/m_debitance_s_d.f90 | 5 ++---
4 files changed, 5 insertions(+), 5 deletions(-)
diff --git a/NEWS.txt b/NEWS.txt
index 9148be08e8..b446a1b2ae 100644
--- a/NEWS.txt
+++ b/NEWS.txt
@@ -1,5 +1,6 @@
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
diff --git a/sources/mascaret/Modules/m_debitance_s.f90 b/sources/mascaret/Modules/m_debitance_s.f90
index 1b36db02a5..343c94eba9 100644
--- a/sources/mascaret/Modules/m_debitance_s.f90
+++ b/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
diff --git a/sources/mascaret/ModulesDeriv/m_debitance_s_b.f90 b/sources/mascaret/ModulesDeriv/m_debitance_s_b.f90
index 612b5929ef..25f1a819ae 100644
--- a/sources/mascaret/ModulesDeriv/m_debitance_s_b.f90
+++ b/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)
diff --git a/sources/mascaret/ModulesDeriv/m_debitance_s_d.f90 b/sources/mascaret/ModulesDeriv/m_debitance_s_d.f90
index ca6812b029..9c1371bb54 100644
--- a/sources/mascaret/ModulesDeriv/m_debitance_s_d.f90
+++ b/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
--
GitLab