The TELEMAC-MASCARET system  trunk
user_bedload_qb.f
Go to the documentation of this file.
1 ! **************************
2  SUBROUTINE user_bedload_qb
3 ! **************************
4  & (hn, u2d, v2d, thetac, houle, hw, tw, thetaw,
5  & tob,tobw,tobcw_mean,tobcw_max, dcla, dens, grav, dstar, ac,
6  & xmve, xmvs, tetap, mu, npoin, qsc, qss, cstaeq)
7 !
8 !***********************************************************************
9 ! GAIA
10 !***********************************************************************
11 !
14 !
17 !
18 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
44 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
45 !
46  USE interface_gaia, ex_user_bedload_qb => user_bedload_qb
47  USE bief
49  IMPLICIT NONE
50 !
51 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
52 !
53  TYPE(bief_obj), INTENT(IN) :: HN,U2D,V2D,THETAC
54  TYPE(bief_obj), INTENT(IN) :: HW, TW, THETAW
55  TYPE(bief_obj), INTENT(IN) :: TOB,TOBW,TOBCW_MEAN,TOBCW_MAX
56  DOUBLE PRECISION, INTENT(IN) :: DCLA, DENS, GRAV, DSTAR, AC
57  DOUBLE PRECISION, INTENT(IN) :: XMVE, XMVS
58  TYPE(bief_obj), INTENT(IN) :: TETAP, MU
59  TYPE(bief_obj), INTENT(IN) :: CSTAEQ
60  INTEGER, INTENT(IN) :: NPOIN
61  LOGICAL, INTENT(IN) :: HOULE
62  TYPE(bief_obj), INTENT(INOUT) :: QSC, QSS
63 !
64 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
65 !
66  INTEGER :: I
67 ! DOUBLE PRECISION :: C1, C2, T
68 !
69 !======================================================================!
70 !======================================================================!
71 ! PROGRAM !
72 !======================================================================!
73 !======================================================================!
74 !
75 ! EXAMPLE BY VAN RIJN
76 !
77 ! C1 = DENS * GRAV * DCLA
78 ! C2 = 0.053D0 * SQRT(DCLA**3*DENS*GRAV) * DSTAR**(-0.3D0)
79 !
80  DO i = 1, npoin
81 !
82 ! TRANSPORT STAGE PARAMETER
83 !
84 ! IF(TETAP%R(I) .LE. AC) THEN
85 ! T = 0.D0
86 ! ELSE
87 ! T = (TETAP%R(I)-AC)/MAX(AC,1.D-06)
88 ! ENDIF
89 !
90 ! BEDLOAD TRANSPORT RATE
91 !
92  qsc%R(i) = 0.d0 ! C2 * T**2.1D0
93  qss%R(i) = 0.d0
94 !
95  ENDDO
96 !
97 ! FOLLOWING LINES NEED TO BE COMMENTED OUT
98 !
99  WRITE(lu,53)
100 
101 53 FORMAT(/,1x,'GAIA IS STOPPED : ',/
102  & ,1x,' SAND TRANSPORT MUST BE CALCULATED IN USER_BEDLOAD_QB')
103  CALL plante(1)
104  stop
105 !
106 !-----------------------------------------------------------------------
107 !
108  RETURN
109  END
subroutine user_bedload_qb(HN, U2D, V2D, THETAC, HOULE, HW, TW, THETAW, TOB, TOBW, TOBCW_MEAN, TOBCW_MAX, DCLA, DENS, GRAV, DSTAR, AC, XMVE, XMVS, TETAP, MU, NPOIN, QSC, QSS, CSTAEQ)
Definition: bief.f:3