The TELEMAC-MASCARET system  trunk
suspension_vanrijn_gaia.f
Go to the documentation of this file.
1 ! **********************************
2  SUBROUTINE suspension_vanrijn_gaia
3 ! **********************************
4 !
5  &(dcla,taup,npoin,grav,xmve,xmvs,vce,zero,ac,cstaeq,zref,
6  & ratio_toce)
7 !
8 !***********************************************************************
9 ! GAIA
10 !***********************************************************************
11 !
13 !
14 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
27 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
28 !
29  USE bief
31  IMPLICIT NONE
32 !
33  ! 2/ GLOBAL VARIABLES
34  ! -------------------
35  TYPE(bief_obj), INTENT(IN) :: TAUP,ZREF
36  INTEGER, INTENT(IN) :: NPOIN
37  DOUBLE PRECISION, INTENT(IN) :: GRAV, XMVE, XMVS,VCE
38  DOUBLE PRECISION, INTENT(IN) :: ZERO,AC,DCLA
39  TYPE(bief_obj), INTENT(INOUT) :: CSTAEQ
40  TYPE(bief_obj), INTENT(IN) :: RATIO_TOCE
41 !
42  ! 3/ LOCAL VARIABLES
43  ! ------------------
44  INTEGER :: I
45  DOUBLE PRECISION :: TAUC,AUX,DSTAR,DENS
46 !
47 !======================================================================!
48 !======================================================================!
49 ! PROGRAMME !
50 !======================================================================!
51 !======================================================================!
52 !
53  ! ******************************** !
54  ! I - CRITICAL SHIELD PARAMETER !
55  ! ******************************** !
56 !
57 !
58  DO i=1,npoin
59 !
60 ! ****************** !
61 ! II - SKIN FRICTION !
62 ! ****************** !
63 !
64 ! TOCE IS MODIFIED IN TOCE_MIX : TOCE_MIX=TOCE_SAND*RATIO_TOCE
65 ! WE CAN APPLICATE DIRECTLY THIS RATIO ON TAUC
66 ! (ratio_toce=1 if no mud)
67 ! why recalculate TAUC?? and not use TOCE
68 ! TAUC = AC * GRAV*(XMVS-XMVE)*DCLA
69  tauc = ac * grav*(xmvs-xmve)*dcla*ratio_toce%R(i)
70  dens = (xmvs - xmve )/ xmve
71  dstar = dcla*(grav*dens/vce**2)**(1.d0/3.d0)
72 !
73 ! ***************** !
74 ! IV - EROSION FLUX !
75 ! ***************** !
76 ! Concentration increased by AVA because it is assumed
77 ! that it is computed only with one class of sediment
78 !
79  IF(dstar.LE.zero) THEN
80  WRITE(lu,*) 'ERROR SUSPENSION_VANRIJN_GAIA'
81  CALL plante(1)
82  stop
83  ENDIF
84  aux=(taup%R(i)-tauc)/tauc
85  IF(aux.GT.zero) THEN
86  cstaeq%R(i)=0.015*dcla*sqrt(aux**3)/(zref%R(i)*dstar**0.3d0)
87  ELSE
88  cstaeq%R(i) = 0.d0
89  ENDIF
90 !------------ in kg/m3
91  cstaeq%R(i) = xmvs*cstaeq%R(i)
92 !
93  ENDDO
94 !
95 !======================================================================!
96 !======================================================================!
97 !======================================================================!
98 !
99  RETURN
100  END
subroutine suspension_vanrijn_gaia(DCLA, TAUP, NPOIN, GRAV, XMVE, XMVS, VCE, ZERO, AC, CSTAEQ, ZREF, RATIO_TOCE)
Definition: bief.f:3