bedload_engel.f

Go to the documentation of this file.

!                   ************************
                    SUBROUTINE bedload_engel
!                   ************************
!
     &(tob,cf,dens,grav,dm,xmve,qsc)
!
!***********************************************************************
! SISYPHE   V7P3                                   10/01/2018
!***********************************************************************
!
!brief    ENGELUND-HANSEN BEDLOAD TRANSPORT FORMULATION.
!
!history  E. PELTIER; C. LENORMANT; J.-M. HERVOUET
!+        11/09/1995
!+        V5P1
!+
!
!history  C.VILLARET
!+        **/10/2003
!+        V5P4
!+
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        13/07/2010
!+        V6P0
!+   Translation of French comments within the FORTRAN sources into
!+   English comments
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        21/08/2010
!+        V6P0
!+   Creation of DOXYGEN tags for automated documentation and
!+   cross-referencing of the FORTRAN sources
!
!history  C.VILLARET (EDF-LNHE), P.TASSI (EDF-LNHE)
!+        19/07/2011
!+        V6P1
!+  Name of variables
!+
!history  P.TASSI (EDF-LNHE)
!+        10/01/2018
!+        V7P3
!+  Correction of the coefficient of the Engelund and Hansen formula
!+ (thanks to Alexander Breugem)
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| CF             |-->| QUADRATIC FRICTION COEFFICIENT
!| DENS           |-->| RELATIVE DENSITY
!| DM             |-->| SEDIMENT GRAIN DIAMETER
!| GRAV           |-->| ACCELERATION OF GRAVITY
!| QSC            |<->| BED LOAD TRANSPORT
!| TOB            |<->| BED SHEAR STRESS (TOTAL FRICTION)
!| XMVE           |-->| FLUID DENSITY
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE interface_sisyphe, ex_bedload_engel => bedload_engel
      USE bief
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      TYPE(bief_obj),   INTENT(IN)    :: TOB, CF
      DOUBLE PRECISION, INTENT(IN)    :: DENS, GRAV, DM, XMVE
      TYPE(bief_obj),   INTENT(INOUT) :: QSC
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER I
      DOUBLE PRECISION :: CENGEL, C1
!
!======================================================================!
!                               PROGRAM                                !
!======================================================================!
!
      c1 = 1.d0/(dens*xmve*grav*dm)
      cengel = 0.05d0*sqrt(dens*grav*dm**3)
!     CALL OS('X=CY    ', X=TETA, Y=TOB , C=C1)
!     CALL OS('X=Y**C  ', X=TETA, Y=TETA, C=5.D0/2.D0)
!
!     CALL OS('X=+(Y,C)', X=QSC , Y=CF  , C=1.D-06)
!     CALL OS('X=1/Y   ', X=QSC , Y=QSC)
!     CALL OS('X=CXY   ', X=QSC , Y=TETA, C=CENGEL)
      DO i=1,qsc%DIM1
!       TOTAL NON DIMENSIONAL STRESS = SQRT(C1*TOB%R(I))**5
!       BEDLOAD TRANSPORT
        qsc%R(i)=cengel*sqrt(c1*tob%R(i))**5/max(cf%R(i),1.d-6)
      ENDDO
!
!======================================================================
!
      RETURN
      END