The TELEMAC-MASCARET system  trunk
suspension_bijker_gaia.f
Go to the documentation of this file.
1 ! *********************************
2  SUBROUTINE suspension_bijker_gaia
3 ! *********************************
4 !
5  &(taup,npoin,charr,qsc,zref,zero,cstaeq,xmve,ratio_toce)
6 !
7 !***********************************************************************
8 ! GAIA
9 !***********************************************************************
10 !
13 !
14 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
27 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
28 !
29  USE interface_gaia,ex_suspension_bijker => suspension_bijker_gaia
30  USE bief
32  IMPLICIT NONE
33 !
34 ! GLOBAL VARIABLES
35 ! -------------------
36  TYPE(bief_obj), INTENT(IN) :: TAUP,QSC
37  type(bief_obj), INTENT(IN) :: zref
38  INTEGER, INTENT(IN) :: NPOIN
39  LOGICAL, INTENT(IN) :: CHARR
40  DOUBLE PRECISION, INTENT(IN) :: ZERO,XMVE
41 !
42  TYPE(bief_obj), INTENT(INOUT) :: CSTAEQ
43  TYPE(bief_obj), INTENT(IN) :: RATIO_TOCE
44 !
45 ! LOCAL VARIABLES
46 ! ------------------
47  INTEGER :: I
48  DOUBLE PRECISION :: USTARP
49 !
50 ! MAXIMUM CONCENTRATION CORRESPONDING TO DENSE PACKING
51 !
52  DOUBLE PRECISION, PARAMETER :: CMAX = 0.6d0
53 !
54 !======================================================================!
55 !======================================================================!
56 ! PROGRAM !
57 !======================================================================!
58 !======================================================================!
59 !
60  IF(.NOT.charr) THEN
61  WRITE(lu,*) 'SUSPENSION_BIJKER_GAIA ERROR ON CHARR'
62  CALL plante(1)
63  stop
64  ENDIF
65 !
66  DO i=1,npoin
67 ! CSTAEQ IS FUNCTION OF QB, AND QB TAKE IN ACCOUNT THE RATIO OF MUD
68 ! RATIO_TOCE IS NOT USE HERE
69  IF(taup%R(i).LE.zero) THEN
70  cstaeq%R(i) = 0.d0
71  ELSE
72  ustarp=sqrt(taup%R(i)/xmve)
73  cstaeq%R(i) = qsc%R(i)/(6.34d0*ustarp*zref%R(i))
74  cstaeq%R(i) = min(cstaeq%R(i),cmax)
75  ENDIF
76 !------------ CSTAEQ%R is in kg/m3 because QSC is already in kg/m3
77 !
78  ENDDO
79 !
80 !======================================================================!
81 !======================================================================!
82 !
83  RETURN
84  END SUBROUTINE suspension_bijker_gaia
subroutine suspension_bijker_gaia(TAUP, NPOIN, CHARR, QSC, ZREF, ZERO, CSTAEQ, XMVE, RATIO_TOCE)
Definition: bief.f:3