bedload_vanrijn.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\bedload_vanrijn.f
00002 !
00062                      SUBROUTINE BEDLOAD_VANRIJN !
00063 !                    ****************************
00064 !
00065      &  (TETAP,MU, NPOIN, DM, DENS, GRAV, DSTAR, AC, QSC)
00066 !
00067 !***********************************************************************
00068 ! SISYPHE   V6P1                                   21/07/2011
00069 !***********************************************************************
00070 !
00071 !
00072 !
00073 !
00074 !
00075 !
00076 !
00077 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00078 !| AC             |<->| CRITICAL SHIELDS PARAMETER
00079 !| DENS           |-->| RELATIVE DENSITY
00080 !| DM             |-->| SEDIMENT GRAIN DIAMETER
00081 !| DSTAR          |-->| NON-DIMENSIONAL DIAMETER
00082 !| GRAV           |-->| ACCELERATION OF GRAVITY
00083 !| MU             |<->| CORRECTION FACTOR FOR BED ROUGHNESS
00084 !| NPOIN          |-->| NUMBER OF POINTS
00085 !| QSC            |<->| BEDLOAD TRANSPORT RATE
00086 !| TETAP          |-->| ADIMENSIONAL SKIN FRICTION
00087 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00088 !
00089       USE INTERFACE_SISYPHE,EX_BEDLOAD_VANRIJN => BEDLOAD_VANRIJN
00090       USE BIEF
00091       IMPLICIT NONE
00092       INTEGER LNG,LU
00093       COMMON/INFO/LNG,LU
00094       ! 2/ GLOBAL VARIABLES
00095       ! -------------------
00096       TYPE(BIEF_OBJ),   INTENT(IN)  :: TETAP,MU
00097       INTEGER,          INTENT(IN)  :: NPOIN
00098       DOUBLE PRECISION, INTENT(IN)  :: DM, DENS, GRAV, DSTAR, AC
00099       TYPE(BIEF_OBJ),   INTENT(INOUT) :: QSC
00100       ! 3/ LOCAL VARIABLES
00101       ! ------------------
00102       INTEGER          :: I
00103       DOUBLE PRECISION :: C1, C2, T
00104 !======================================================================!
00105 !======================================================================!
00106 !                               PROGRAM                                !
00107 !======================================================================!
00108 !======================================================================!
00109       C1 = DENS * GRAV * DM
00110       C2 = 0.053D0 * SQRT(DM**3*DENS*GRAV) * DSTAR**(-0.3D0)
00111       DO I = 1, NPOIN
00112         ! ****************************** !
00113         ! I - TRANSPORT STAGE PARAMETER  !
00114         ! ****************************** !
00115         IF(TETAP%R(I) .LE. AC) THEN
00116           T = 0.D0
00117         ELSE
00118           T = (TETAP%R(I)-AC)/MAX(AC,1.D-06)
00119         ENDIF
00120 !
00121         ! ***************************** !
00122         ! II - BEDLOAD TRANSPORT RATE   !
00123         ! ***************************** !
00124         QSC%R(I) = C2 * T**2.1D0
00125       ENDDO
00126 !======================================================================!
00127 !======================================================================!
00128       RETURN
00129       END SUBROUTINE BEDLOAD_VANRIJN

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