bedload_einst.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\bedload_einst.f
00002 !
00062                      SUBROUTINE BEDLOAD_EINST
00063 !                    ************************
00064 !
00065      &(TETAP, NPOIN, DENS, GRAV, DM, DSTAR, QSC)
00066 !
00067 !***********************************************************************
00068 ! SISYPHE   V6P1                                   21/07/2011
00069 !***********************************************************************
00070 !
00071 !
00072 !
00073 !
00074 !
00075 !
00076 !
00077 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00078 !| DENS           |-->| RELATIVE SENSITY OF SEDIMENT
00079 !| DM             |-->| SEDIMENT GRAIN DIAMETER
00080 !| DSTAR          |-->| NON-DIMENSIONAL DIAMETER
00081 !| GRAV           |-->| ACCELERATION OF GRAVITY
00082 !| NPOIN          |-->| NUMBER OF POINTS
00083 !| QSC            |<->| BED LOAD TRANSPORT
00084 !| TETAP          |-->| ADIMENSIONAL SKIN FRICTION
00085 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00086 !
00087       USE INTERFACE_SISYPHE, EX_BEDLOAD_EINST => BEDLOAD_EINST
00088       USE BIEF
00089       IMPLICIT NONE
00090       INTEGER LNG,LU
00091       COMMON/INFO/LNG,LU
00092 !
00093 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00094 !
00095       TYPE(BIEF_OBJ),   INTENT(IN)    :: TETAP
00096       INTEGER,          INTENT(IN)    :: NPOIN
00097       DOUBLE PRECISION, INTENT(IN)    :: DENS, GRAV, DM, DSTAR
00098       TYPE(BIEF_OBJ),   INTENT(INOUT)   :: QSC
00099 !
00100 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00101 !
00102       INTEGER          :: I
00103       DOUBLE PRECISION :: CEINST
00104 !
00105 !======================================================================!
00106 !                               PROGRAM                                !
00107 !======================================================================!
00108 !
00109 !     BEDLOAD TRANSPORT
00110 !
00111       CEINST = 36.D0/(DSTAR**3)
00112       CEINST = SQRT(2.D0/3.D0+CEINST) -  SQRT(CEINST)
00113       CEINST = CEINST * SQRT(DENS*GRAV*(DM**3))
00114       DO I = 1, NPOIN
00115         IF(TETAP%R(I) < 2.5D-3) THEN
00116           QSC%R(I) = 0.D0
00117         ELSE IF (TETAP%R(I) < 0.2D0) THEN
00118           QSC%R(I) = 2.15D0* CEINST * EXP(-0.391D0/TETAP%R(I))
00119         ELSE
00120           QSC%R(I) = 40.D0 * CEINST * (TETAP%R(I)**3.D0)
00121         ENDIF
00122       ENDDO
00123 !
00124 !=======================================================================
00125 !
00126       RETURN
00127       END

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