The TELEMAC-MASCARET system  trunk
bedload_vanrijn_gaia.f
Go to the documentation of this file.
1 ! *******************************
2  SUBROUTINE bedload_vanrijn_gaia
3 ! *******************************
4 !
5  &(tetap, npoin, dcla, dens, grav, dstar, ac, qsc, xmvs)
6 !
7 !***********************************************************************
8 ! GAIA V6P1 21/07/2011
9 !***********************************************************************
10 !
12 !
13 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
23 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
24 !
25  USE interface_gaia,ex_bedload_vanrijn => bedload_vanrijn_gaia
26  USE bief
28  IMPLICIT NONE
29 !
30  TYPE(bief_obj), INTENT(IN) :: TETAP
31  INTEGER, INTENT(IN) :: NPOIN
32  DOUBLE PRECISION, INTENT(IN) :: DCLA, DENS, GRAV, DSTAR, AC,XMVS
33  TYPE(bief_obj), INTENT(INOUT) :: QSC
34 !
35  INTEGER :: I
36  DOUBLE PRECISION :: C1, C2, T
37 !======================================================================!
38 !======================================================================!
39 ! PROGRAM !
40 !======================================================================!
41 !======================================================================!
42  c1 = dens * grav * dcla
43  c2 = 0.053d0 * sqrt(dcla**3*dens*grav) * dstar**(-0.3d0)
44  DO i = 1, npoin
45 ! ******************************
46 ! I - TRANSPORT STAGE PARAMETER
47 ! ******************************
48  IF(tetap%R(i) .LE. ac) THEN
49  t = 0.d0
50  ELSE
51  t = (tetap%R(i)-ac)/max(ac,1.d-06)
52  ENDIF
53 !
54 ! *****************************
55 ! II - BEDLOAD TRANSPORT RATE
56 ! *****************************
57  qsc%R(i) = c2 * t**2.1d0
58  ENDDO
59 !======================================================================!
60 ! SOLID DISCHARGE IS TRANSFORMED IN [kg/(m*s)]
61 !
62  CALL os('X=CX ', x=qsc, c=xmvs)
63 !======================================================================!
64  RETURN
65  END SUBROUTINE bedload_vanrijn_gaia
subroutine bedload_vanrijn_gaia(TETAP, NPOIN, DCLA, DENS, GRAV, DSTAR, AC, QSC, XMVS)
subroutine os(OP, X, Y, Z, C, IOPT, INFINI, ZERO)
Definition: os.f:7
Definition: bief.f:3