The TELEMAC-MASCARET system  trunk
bedload_vanrijn.f
Go to the documentation of this file.
1 ! **************************
2  SUBROUTINE bedload_vanrijn
3 ! **************************
4 !
5  &(tetap,npoin, dm, dens, grav, dstar, ac, qsc)
6 !
7 !***********************************************************************
8 ! SISYPHE V6P1 21/07/2011
9 !***********************************************************************
10 !
11 !brief VAN RIJN BEDLOAD TRANSPORT FORMULATION.
12 !
13 !history BUI MINH DUC
14 !+ **/10/2001
15 !+ V5P2
16 !+
17 !
18 !history C. VILLARET
19 !+ **/**/2004
20 !+ V5P4
21 !+
22 !
23 !history N.DURAND (HRW), S.E.BOURBAN (HRW)
24 !+ 13/07/2010
25 !+ V6P0
26 !+ Translation of French comments within the FORTRAN sources into
27 !+ English comments
28 !
29 !history N.DURAND (HRW), S.E.BOURBAN (HRW)
30 !+ 21/08/2010
31 !+ V6P0
32 !+ Creation of DOXYGEN tags for automated documentation and
33 !+ cross-referencing of the FORTRAN sources
34 !
35 !history C.VILLARET (EDF-LNHE), P.TASSI (EDF-LNHE)
36 !+ 19/07/2011
37 !+ V6P1
38 !+ Name of variables
39 !+
40 !
41 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
42 !| AC |<->| CRITICAL SHIELDS PARAMETER
43 !| DENS |-->| RELATIVE DENSITY
44 !| DM |-->| SEDIMENT GRAIN DIAMETER
45 !| DSTAR |-->| NON-DIMENSIONAL DIAMETER
46 !| GRAV |-->| ACCELERATION OF GRAVITY
47 !| NPOIN |-->| NUMBER OF POINTS
48 !| QSC |<->| BEDLOAD TRANSPORT RATE
49 !| TETAP |-->| ADIMENSIONAL SKIN FRICTION
50 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
51 !
52  USE interface_sisyphe,ex_bedload_vanrijn => bedload_vanrijn
53  USE bief
55  IMPLICIT NONE
56  ! 2/ GLOBAL VARIABLES
57  ! -------------------
58  TYPE(bief_obj), INTENT(IN) :: TETAP
59  INTEGER, INTENT(IN) :: NPOIN
60  DOUBLE PRECISION, INTENT(IN) :: DM, DENS, GRAV, DSTAR, AC
61  TYPE(bief_obj), INTENT(INOUT) :: QSC
62  ! 3/ LOCAL VARIABLES
63  ! ------------------
64  INTEGER :: I
65  DOUBLE PRECISION :: C1, C2, T
66 !======================================================================!
67 !======================================================================!
68 ! PROGRAM !
69 !======================================================================!
70 !======================================================================!
71  c1 = dens * grav * dm
72  c2 = 0.053d0 * sqrt(dm**3*dens*grav) * dstar**(-0.3d0)
73  DO i = 1, npoin
74  ! ****************************** !
75  ! I - TRANSPORT STAGE PARAMETER !
76  ! ****************************** !
77  IF(tetap%R(i) .LE. ac) THEN
78  t = 0.d0
79  ELSE
80  t = (tetap%R(i)-ac)/max(ac,1.d-06)
81  ENDIF
82 !
83  ! ***************************** !
84  ! II - BEDLOAD TRANSPORT RATE !
85  ! ***************************** !
86  qsc%R(i) = c2 * t**2.1d0
87  ENDDO
88 !======================================================================!
89 !======================================================================!
90  RETURN
91  END SUBROUTINE bedload_vanrijn
subroutine bedload_vanrijn(TETAP, NPOIN, DM, DENS, GRAV, DSTAR, AC, QSC)
Definition: bief.f:3