The TELEMAC-MASCARET system  trunk
settling_vel.f
Go to the documentation of this file.
1 ! ***********************
2  SUBROUTINE settling_vel
3 ! ***********************
4 !***********************************************************************
5 ! GAIA
6 !***********************************************************************
7 !
13 !
14 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
15 !
17  IMPLICIT NONE
18 !
19 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
20  INTEGER :: ICLA
21  DOUBLE PRECISION :: DENS
22 !
23 !
24 ! SETTLING VELOCITY
25 !
26  DO icla = 1, nsicla
27  IF(xwc0(icla).LT.0) THEN
28 ! SETTLING VELOCITY IS NOT GIVEN IN THE PARAMETER FILE AND INITIALISED HERE
29 ! IF T3D: SETTLING VELOCITY IS POTENTIALY MODIFIED IN 3D AND THEN GIVEN TO GAIA
30 ! IF T2D: ?
31  dens = (xmvs0(icla) - xmve) / xmve
32  IF (dcla(icla).LT.1.d-4) THEN
33  xwc0(icla) = dens * dcla(icla) *
34  & dcla(icla) * grav / ( 18.d0 * vce )
35  ELSEIF (dcla(icla).LT.1d-3) THEN
36  xwc0(icla) = 10.d0 * vce / dcla(icla) *
37  & (sqrt( 1.d0 + 0.01d0* dens * grav *
38  & dcla(icla)**3.d0 / (vce*vce) ) -1.d0 )
39  ELSE
40  xwc0(icla) = 1.1d0 * sqrt( dens * grav * dcla(icla) )
41  ENDIF
42  ENDIF
43  xwc(icla)=xwc0(icla)
44  ENDDO
45 !
46  RETURN
47  END
subroutine settling_vel
Definition: settling_vel.f:4
double precision, dimension(nsiclm), target xwc
Settling velocities.
double precision xmve
Water density (from steering file of T2D or T3D)
double precision vce
Water viscosity: it is defined here because the viscosity set in TELEMAC2D or TELEMAC3D may not b the...
double precision, dimension(nsiclm) xmvs0
Sand density.
double precision, dimension(nsiclm), target dcla
Sediment diameter for each class It is only relevant for non-cohesive sediments. For the bedload...
integer, target nsicla
Number of sediment classes of bed material (less than NISCLM)
double precision grav
Gravity acceleration.
double precision, dimension(nsiclm) xwc0
Initial settling velocities.