suspension_bijker.f

Go to the documentation of this file.

!                   ****************************
                    SUBROUTINE suspension_bijker
!                   ****************************
!
     &(taup,npoin,charr,qsc,zref,zero,cstaeq,xmve)
!
!***********************************************************************
! SISYPHE   V6P1                                   21/07/2011
!***********************************************************************
!
!brief    COMPUTES THE REFERENCE CONCENTRATION AT Z= 2*D50
!+                USING ZYSERMAN AND FREDSOE FORMULATION (1994).
!
!history  C. VILLARET
!+        14/04/2004
!+        V5P5
!+
!
!history  F. HUVELIN
!+        04/01/2005
!+        V5P6
!+
!
!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
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| CHARR          |-->| BEDLOAD
!| CSTAEQ         |<->| EQUILIBRIUM CONCENTRATION
!| NPOIN          |-->| NUMBER OF POINTS
!| QSC            |-->| BED LOAD TRANSPPORT RATE
!| TAUP           |-->| SKIN FROCTION
!| XMVE           |-->| FLUID DENSITY
!| ZERO           |-->| ZERO
!| ZREF           |-->| REFERENCE ELEVATION
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE interface_sisyphe,ex_suspension_bijker => suspension_bijker
      USE bief
      USE declarations_special
      IMPLICIT NONE
!
      ! 2/ GLOBAL VARIABLES
      ! -------------------
      TYPE(bief_obj),   INTENT(IN)    :: TAUP,QSC
      type(bief_obj),  INTENT(IN)    :: zref
      INTEGER,          INTENT(IN)    :: NPOIN
      LOGICAL,          INTENT(IN)    :: CHARR
      DOUBLE PRECISION, INTENT(IN)    :: ZERO,XMVE
!
      TYPE(bief_obj),   INTENT(INOUT)   ::  CSTAEQ
!
      ! 3/ LOCAL VARIABLES
      ! ------------------
      INTEGER                     :: I
      DOUBLE PRECISION            :: USTARP
!
!     MAXIMUM CONCENTRATION CORRESPONDING TO DENSE PACKING
!
      DOUBLE PRECISION, PARAMETER :: CMAX = 0.6d0
!
!======================================================================!
!======================================================================!
!                               PROGRAM                                !
!======================================================================!
!======================================================================!
!
      IF(.NOT.charr) THEN
        WRITE(lu,*) 'SUSPENSION_BIJKER ERROR ON CHARR'
        CALL plante(1)
        stop
      ENDIF
!
      DO i=1,npoin
!
        IF(taup%R(i).LE.zero) THEN
          cstaeq%R(i) = 0.d0
        ELSE
          ustarp=sqrt(taup%R(i)/xmve)
          cstaeq%R(i) = qsc%R(i)/(6.34d0*ustarp*zref%R(i))
          cstaeq%R(i) = min(cstaeq%R(i),cmax)
        ENDIF
!
      ENDDO
!
!======================================================================!
!======================================================================!
!
      RETURN
      END SUBROUTINE suspension_bijker