The TELEMAC-MASCARET system  trunk
suspension_fredsoe_gaia.f
Go to the documentation of this file.
1 ! **********************************
2  SUBROUTINE suspension_fredsoe_gaia
3 ! **********************************
4 !
5  &(dcla,taup,npoin,grav,xmve,xmvs,ac,cstaeq,ratio_toce)
6 !
7 !***********************************************************************
8 ! GAIA
9 !***********************************************************************
10 !
13 !
14 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
25 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
26 !
27  USE interface_gaia,ex_suspension_fredsoe =>
29  USE bief
31  IMPLICIT NONE
32 !
33 ! GLOBAL VARIABLES
34 ! -------------------
35  TYPE(bief_obj), INTENT(IN) :: TAUP
36  INTEGER, INTENT(IN) :: NPOIN
37  DOUBLE PRECISION, INTENT(IN) :: GRAV, XMVE, XMVS
38  DOUBLE PRECISION, INTENT(IN) :: AC,DCLA
39  TYPE(bief_obj), INTENT(INOUT) :: CSTAEQ
40  TYPE(bief_obj), INTENT(IN) :: RATIO_TOCE
41 ! LOCAL VARIABLES
42 ! ------------------
43  INTEGER :: I
44  DOUBLE PRECISION :: TETAP,AUX
45 !
46 ! MAXIMUM CONCENTRATION CORRESPONDING TO DENSE PACKING
47 !
48  DOUBLE PRECISION, PARAMETER :: CMAX = 0.6d0
49 !!
50 !======================================================================!
51 !======================================================================!
52 ! PROGRAM !
53 !======================================================================!
54 !======================================================================!
55 !
56 ! ********************************
57 ! I - CRITICAL SHIELDS PARAMETER
58 ! ********************************
59 !
60  DO i=1,npoin
61 !
62 ! ****************** !
63 ! II - SKIN FRICTION !
64 ! ****************** !
65 !
66  tetap = taup%R(i) / (grav*(xmvs-xmve)*dcla)
67 !
68 ! *****************
69 ! IV - EROSION FLUX
70 ! *****************
71 ! CONCENTRATION INCREASED BY AVA BECAUSE IT IS COMPUTED
72 ! ONLY WITH ONE CLASS OF SEDIMENT (ASSUMPTION)
73 !
74 ! TOCE IS MODIFIED IN TOCE_MIX : TOCE_MIX=TOCE_SAND*RATIO_TOCE
75 ! WE CAN APPLICATE DIRECTLY THIS RATIO ON AC
76 ! because TOCE=AC*DM*GRAV*(XMVS-XMVE)
77 ! (ratio_toce=1 if no mud)
78 ! IF(TETAP.GT.(AC) THEN
79 ! AUX=(TETAP-AC)**1.75D0
80  IF(tetap.GT.(ac*ratio_toce%R(i))) THEN
81  aux=(tetap-(ac*ratio_toce%R(i)))**1.75d0
82  cstaeq%R(i) = 0.331d0*aux/(1.d0+0.72d0*aux)
83  cstaeq%R(i) = min(cstaeq%R(i),cmax)
84  ELSE
85  cstaeq%R(i) = 0.d0
86  ENDIF
87 !------------ in kg/m3
88  cstaeq%R(i) = xmvs*cstaeq%R(i)
89 !
90  ENDDO
91 !
92 !======================================================================!
93 !======================================================================!
94 !
95  RETURN
96  END
subroutine suspension_fredsoe_gaia(DCLA, TAUP, NPOIN, GRAV, XMVE, XMVS, AC, CSTAEQ, RATIO_TOCE)
Definition: bief.f:3