suspension_fredsoe.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\suspension_fredsoe.f
00002 !
00073                      SUBROUTINE SUSPENSION_FREDSOE
00074 !                    *****************************
00075 !
00076      &(DM,TAUP,NPOIN,GRAV,XMVE,XMVS,ZERO,AC,CSTAEQ)
00077 !
00078 !***********************************************************************
00079 ! SISYPHE   V6P1                                   21/07/2011
00080 !***********************************************************************
00081 !
00082 !
00083 !
00084 !
00085 !
00086 !
00087 !
00088 !
00089 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00090 !| AC             |<->| CRITICAL SHIELDS PARAMETER
00091 !| DM             |-->| SEDIMENT GRAIN DIAMETER
00092 !| CSTAEQ         |<->| EQUILIBRIUM CONCENTRATION
00093 !| GRAV           |-->| ACCELERATION OF GRAVITY
00094 !| NPOIN          |-->| NUMBER OF POINTS
00095 !| TAUP           |-->| CRITICAL SHEAR STRESS
00096 !| XMVE           |-->| FLUID DENSITY
00097 !| XMVS           |-->| WATER DENSITY
00098 !| ZERO           |-->| ZERO
00099 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00100 !
00101       USE INTERFACE_SISYPHE,EX_SUSPENSION_FREDSOE => SUSPENSION_FREDSOE
00102       USE BIEF
00103       IMPLICIT NONE
00104       INTEGER LNG,LU
00105       COMMON/INFO/LNG,LU
00106 !
00107 !
00108       ! 2/ GLOBAL VARIABLES
00109       ! -------------------
00110       TYPE(BIEF_OBJ),   INTENT(IN)    ::  TAUP
00111       INTEGER,          INTENT(IN)    :: NPOIN
00112       DOUBLE PRECISION, INTENT(IN)    :: GRAV, XMVE, XMVS
00113       DOUBLE PRECISION, INTENT(IN)    :: ZERO,AC,DM
00114       TYPE(BIEF_OBJ),   INTENT(INOUT) :: CSTAEQ
00115 
00116 
00117       ! 3/ LOCAL VARIABLES
00118       ! ------------------
00119       INTEGER            :: I
00120       DOUBLE PRECISION   ::  TETAP,AUX
00121 !
00122       DOUBLE PRECISION   :: CMAX
00123 !
00124 !     MAXIMUM CONCENTRATION CORRESPONDING TO DENSE PACKING
00125 !
00126       DATA CMAX/0.6D0/
00127       INTRINSIC MAX
00128 !
00129 !======================================================================!
00130 !======================================================================!
00131 !                               PROGRAM                                !
00132 !======================================================================!
00133 !======================================================================!
00134 !
00135       ! ******************************** !
00136       !    I - CRITICAL SHIELDS PARAMETER!
00137       ! ******************************** !
00138 !
00139       DO I=1,NPOIN
00140 !
00141         ! ****************** !
00142         ! II - SKIN FRICTION !
00143         ! ****************** !
00144 !
00145         TETAP = TAUP%R(I) / (GRAV*(XMVS-XMVE)*DM)
00146 !
00147         ! ***************** !
00148         ! IV - EROSION FLUX ! (_IMP_)
00149         ! ***************** !
00150         ! CONCENTRATION INCREASED BY AVA BECAUSE IT IS COMPUTED
00151         ! ONLY WITH ONE CLASS OF SEDIMENT (ASSUMPTION)
00152 !
00153         IF(TETAP.GT.AC) THEN
00154           AUX=(TETAP-AC)**1.75D0
00155           CSTAEQ%R(I) = 0.331D0*AUX/(1.D0+0.72D0*AUX)
00156           CSTAEQ%R(I) = MIN(CSTAEQ%R(I),CMAX)
00157         ELSE
00158           CSTAEQ%R(I) = 0.D0
00159         ENDIF
00160 !
00161       ENDDO
00162 !
00163 !======================================================================!
00164 !======================================================================!
00165 !
00166       RETURN
00167       END

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0