vitchu_sisyphe.f

Go to the documentation of this file.

!                   *************************
                    SUBROUTINE vitchu_sisyphe
!                   *************************
!
     & ( ws , dens , dm , grav , vce )
!
!***********************************************************************
! SISYPHE   V6P1                                   21/07/2011
!***********************************************************************
!
!brief    COMPUTES THE FALL VELOCITY.
!
!history
!+        20/05/96
!+        V5P1
!+
!
!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
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| DENS           |-->| POIDS DEJAUGE
!| DM             |-->| MEAN DIAMETER OF THE SEDIMENT
!| GRAV           |-->| ACCELERATION OF GRAVITY
!| VCE            |-->| FLOW VISCOSITY
!| WS             |-->| SETTLING/FALL VELOCITY
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      DOUBLE PRECISION, INTENT(IN)    :: DENS,  DM,  GRAV, VCE
      DOUBLE PRECISION, INTENT(INOUT) :: WS
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
!
! FALL VELOCITY
! ================
!
      IF (dm.LT.1.d-4) THEN
        ws = dens * dm * dm * grav / ( 18.d0 * vce )
      ELSEIF (dm.LT.1d-3) THEN
        ws = 10.d0 * vce / dm * (sqrt( 1.d0 + 0.01d0* dens * grav *
     &       dm**3.d0 / (vce*vce) ) -1.d0 )
      ELSE
        ws = 1.1d0 * sqrt( dens * grav * dm )
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END SUBROUTINE vitchu_sisyphe