suspension_vanrijn.f

Go to the documentation of this file.

!                     *****************************
                      SUBROUTINE suspension_vanrijn
!                     *****************************
!
     &(fdm,taup,npoin,grav,xmve,xmvs,vce,zero,ac,cstaeq,zref)
!
!***********************************************************************
! SISYPHE   V6P1                                   21/07/2011
!***********************************************************************
! NEW SUBROUTINE C. VILLARET N. HUYBRECHTS
!
!history  C.VILLARET (EDF-LNHE), P.TASSI (EDF-LNHE)
!+        19/07/2011
!+        V6P1
!+   Name of variables
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| AC             |<->| CRITICAL SHIELDS PARAMETER
!| ACLADM         |-->| MEAN DIAMETER OF SEDIMENT
!| GRAV           |-->| ACCELERATION OF GRAVITY
!| NPOIN          |-->| NUMBER OF POINTS
!| VCE            |-->| FLOW VISCOSITY
!| XMVE           |-->| FLUID DENSITY
!| XMVS           |-->| SEDIMENT DENSITY
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
      USE declarations_special
      IMPLICIT NONE
!
      ! 2/ GLOBAL VARIABLES
      ! -------------------
      TYPE(bief_obj),   INTENT(IN)    ::  TAUP,ZREF
      INTEGER,          INTENT(IN)    :: NPOIN
      DOUBLE PRECISION, INTENT(IN)    ::  GRAV,  XMVE, XMVS,VCE
      DOUBLE PRECISION, INTENT(IN)    :: ZERO,AC,FDM
      TYPE(bief_obj),   INTENT(INOUT) :: CSTAEQ
!
      ! 3/ LOCAL VARIABLES
      ! ------------------
      INTEGER            :: I
      DOUBLE PRECISION   ::  TAUC,AUX,DSTAR,DENS
!
!======================================================================!
!======================================================================!
!                               PROGRAMME                              !
!======================================================================!
!======================================================================!
!
      ! ******************************** !
      !    I - CRITICAL SHIELD PARAMETER !
      ! ******************************** !
!
!
      DO i=1,npoin
!
! ****************** !
! II - SKIN FRICTION !
! ****************** !
!
        tauc = ac * grav*(xmvs-xmve)*fdm
        dens  = (xmvs - xmve )/ xmve
        dstar = fdm*(grav*dens/vce**2)**(1.d0/3.d0)
!
! ***************** !
! IV - EROSION FLUX !
! ***************** !
! Concentration increased by AVA because it is assumed
! that it is computed only with one class of sediment
!
        IF(dstar.LE.zero) THEN
          WRITE(lu,*) 'ERROR SUSPENSION_VANRIJN'
          CALL plante(1)
          stop
        ENDIF
        aux=(taup%R(i)-tauc)/tauc
        IF(aux.GT.zero) THEN
          cstaeq%R(i)=0.015*fdm*sqrt(aux**3)/(zref%R(i)*dstar**0.3d0)
        ELSE
          cstaeq%R(i) = 0.d0
        ENDIF
!
      ENDDO
!
!======================================================================!
!======================================================================!
!======================================================================!
!
      RETURN
      END