The TELEMAC-MASCARET system  trunk
gaia_suspension_erode.f
Go to the documentation of this file.
1 ! ********************************
2  SUBROUTINE gaia_suspension_erode
3 ! ********************************
4 !
5 !***********************************************************************
6 ! GAIA
7 !***********************************************************************
8 !
10 !
11 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
12 !
13  USE bief
16  IMPLICIT NONE
17 !
18 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
19 ! COMPUTES EROSION FLUX : FLUER
20 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
21 ! SKIN FRICTION TAUP
22 !
23  CALL os('X=CYZ ', x=taup, y=tobcw_mean, z=mu, c=1.d0)
24  CALL os('X=+(Y,C)', x=taup, y=taup, c=zero)
25 !
26  IF (bed_model.EQ.1.OR.bed_model.EQ.2)THEN
27  IF(debug.GT.0) WRITE(lu,*) 'BED1_SUSPENSION_ERODE'
29  IF(debug.GT.0) WRITE(lu,*) 'END_BED1_SUSPENSION_ERODE'
30  ELSE
31  WRITE(lu,*)'ONLY BED_MODEL 1 or 2 for tis time'
32  stop
33  ENDIF
34 !
35  RETURN
36  END
double precision, dimension(:), pointer x
2d coordinates of the mesh
integer bed_model
Bed model (3 choices, cf. dico)
double precision, dimension(:), pointer y
subroutine gaia_suspension_erode
double precision zero
Parameter used for clipping variables or testing values against zero.
subroutine bed1_suspension_erode
type(bief_obj), target z
Free surface elevation.
type(bief_obj), target mu
Skin friction correction factor for bed roughness: Ratio between shear stress due skin friction and t...
type(bief_obj), target taup
Shear stress modified by skin friction.
integer debug
Debugger.
subroutine os(OP, X, Y, Z, C, IOPT, INFINI, ZERO)
Definition: os.f:7
type(bief_obj), target tobcw_mean
Mean of total current + wave shear stress.
Definition: bief.f:3