The TELEMAC-MASCARET system  trunk
suspension_compute_cae.f
Go to the documentation of this file.
1 ! *********************************
2  SUBROUTINE suspension_compute_cae
3 ! *********************************
4 !
5  &(taup,hn,dcla,npoin,charr,xmve,xmvs,vce,grav,
6  & zero,zref,ac,cstaeq,qsc,icq,u2d,v2d,csratio,debug,ratio_toce)
7 !
8 !***********************************************************************
9 ! GAIA
10 !***********************************************************************
11 !
13 !
14 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
34 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
35 !
36  USE interface_gaia, ex_suspension_compute_cae =>
38  USE bief
39  USE declarations_gaia, ONLY : nomblay
41  IMPLICIT NONE
42 ! GLOBAL VARIABLES
43 ! -------------------
44  type(bief_obj), INTENT(IN) :: taup,hn,zref,qsc
45  type(bief_obj), INTENT(IN) :: u2d,v2d,csratio
46  INTEGER, INTENT(IN) :: NPOIN,DEBUG,ICQ
47  LOGICAL, INTENT(IN) :: CHARR
48  DOUBLE PRECISION, INTENT(IN) :: XMVE,XMVS,GRAV,VCE
49  DOUBLE PRECISION, INTENT(IN) :: ZERO,DCLA
50  DOUBLE PRECISION, INTENT(IN) :: AC
51  type(bief_obj), INTENT(INOUT) :: cstaeq,ratio_toce
52 !======================================================================!
53 !======================================================================!
54 ! PROGRAM !
55 !======================================================================!
56 !======================================================================!
57 !
58 ! COMPUTES THE NEAR BED EQUILIBRIUM CONCENTRATION --> CSTAEQ
59 !
60  IF(icq.EQ.1) THEN
61 !
62  IF(debug > 0) WRITE(lu,*) 'SUSPENSION_FREDSOE_GAIA'
63  CALL suspension_fredsoe_gaia(dcla,taup,npoin,
64  & grav,xmve,xmvs,ac,cstaeq,ratio_toce)
65  IF(debug > 0) WRITE(lu,*) 'END SUSPENSION_FREDSOE_GAIA'
66 !
67  ELSEIF(icq.EQ.2) THEN
68  IF(debug > 0) WRITE(lu,*) 'SUSPENSION_BIJKER_GAIA'
69  IF(nomblay.GE.1)THEN
71  WRITE(lu,*)'QSC IS COMING FROM FIRST LAYER OF BEDLOAD PART
72  & IN THE CASE OF MULTILAYER OR ACTIVE LAYER,
73  & QSC HAVE BEEN TO BE
74  & COMPUTED FOR OTHERS LAYER IF THERE IS EROSION
75  & OF THE FIRST LAYER IN BED1_SUSPENSION_ERODE'
76  CALL plante(1)
77  stop
78  ENDIF
79  CALL suspension_bijker_gaia(taup,npoin,charr,qsc,zref,
80  & zero,cstaeq,xmve,ratio_toce)
81  IF(debug > 0) WRITE(lu,*) 'END SUSPENSION_BIJKER_GAIA'
82 !
83  ELSEIF(icq.EQ.3) THEN
84 !
85  IF(debug > 0) WRITE(lu,*) 'SUSPENSION_VANRIJN_GAIA'
86  CALL suspension_vanrijn_gaia(dcla,taup,npoin,grav,xmve,xmvs,
87  & vce,zero,ac,cstaeq,zref,ratio_toce)
88  IF(debug > 0) WRITE(lu,*) 'END SUSPENSION_VANRIJN_GAIA'
89 
90  ELSEIF(icq.EQ.4) THEN
91  IF(debug > 0) WRITE(lu,*) 'SUSPENSION_SANDFLOW_GAIA'
92  CALL suspension_sandflow_gaia(dcla,npoin,grav,xmve,xmvs,
93  & zero,cstaeq,hn,u2d,v2d,csratio,ratio_toce)
94  IF(debug > 0) WRITE(lu,*) 'END SUSPENSION_SANDFLOW_GAIA'
95 !
96  ELSEIF(icq.EQ.0) THEN
97  IF(debug > 0) WRITE(lu,*) 'SUSPENSION_CAE_USER_GAIA'
98  CALL user_suspension_cae_gaia(dcla,npoin,xmvs,cstaeq)
99  IF(debug > 0) WRITE(lu,*) 'END SUSPENSION_CAE_USER_GAIA'
100 !
101  ENDIF
102 !
103 !======================================================================!
104 !======================================================================!
105 !
106  RETURN
107  END
108 !
109 !#######################################################################
110 !
subroutine suspension_sandflow_gaia(DCLA, NPOIN, GRAV, XMVE, XMVS, ZERO, CSTAEQ, HN, U2D, V2D, CSRATIO, RATIO_TOCE)
integer, target nomblay
Number of bed load model layers = NUMSTRAT+1 to take the active layer into account.
subroutine suspension_compute_cae(TAUP, HN, DCLA, NPOIN, CHARR, XMVE, XMVS, VCE, GRAV, ZERO, ZREF, AC, CSTAEQ, QSC, ICQ, U2D, V2D, CSRATIO, DEBUG, RATIO_TOCE)
subroutine suspension_vanrijn_gaia(DCLA, TAUP, NPOIN, GRAV, XMVE, XMVS, VCE, ZERO, AC, CSTAEQ, ZREF, RATIO_TOCE)
subroutine user_suspension_cae_gaia(DCLA, NPOIN, XMVS, CSTAEQ)
subroutine suspension_bijker_gaia(TAUP, NPOIN, CHARR, QSC, ZREF, ZERO, CSTAEQ, XMVE, RATIO_TOCE)
subroutine suspension_fredsoe_gaia(DCLA, TAUP, NPOIN, GRAV, XMVE, XMVS, AC, CSTAEQ, RATIO_TOCE)
Definition: bief.f:3