The TELEMAC-MASCARET system  trunk
zref_gaia.f
Go to the documentation of this file.
1 ! ********************
2  SUBROUTINE zref_gaia
3 ! ********************
4 !
5 !***********************************************************************
6 ! GAIA
7 !***********************************************************************
8 !
10 !
11 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
12 !
13  USE bief
16  IMPLICIT NONE
17 !
18 !-----------------------------------------------------------------------
24 !
25 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
26 ! 1. COMPUTES THE REFERENCE ELEVATION --> ZREF
27 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
28 !
29 ! THREE OPTIONS : ICQ=1: FREDSOE REFERENCE CONC. ZREF = 2.D50
30 ! ICQ=2: BIJKER METHOD ZREF = MAX(KSP,KS)
31 ! ICQ=3: VAN RIJN ZREF= 0.5 KS
32 ! ICQ=4: SOULSBY/VAN RIJN ZREF= 0.5 KS
33 !
34 ! Each relation proposed to estimate the near-bed concentration
35 ! under equilibrium conditions has an associated reference height
36 ! called zref. For further details see Sedimentation Engineering
37 ! (edited by Garcia, M.)
38 !
39  IF(icq.EQ.1) THEN
40  CALL os('X=Y ', x=zref, y=ksp)
41  ELSEIF(icq.EQ.2) THEN
42  CALL os('X=Y ', x=zref, y=ksr)
43  ELSEIF(icq.EQ.3) THEN
44  CALL os('X=CY ', x=zref, y=ks,c=0.5d0)
45  ELSEIF(icq.EQ.4) THEN
46  CALL os('X=CY ', x=zref, y=ks,c=0.5d0)
47  ELSEIF(icq.EQ.0) THEN
48 ! CALL OS('X=CY ', X=ZREF, Y=KS,C=0.5D0)
49 ! The value of ZREF must be provided by the user according to the
50 ! near-bed concentration formula
51  WRITE(lu,201)
52 201 FORMAT(1x,'SUSPENSION_MAIN_GAIA:',/,1x,
53  & 'DEFINE ZREF FOR USER FORMULA')
54  CALL plante(1)
55  stop
56  ELSE
57  WRITE(lu,203) icq
58 203 FORMAT(1x,'SUSPENSION_MAIN_GAIA:',/,1x,
59  & 'REFERENCE CONCENTRATION FORMULA',/,1x,
60  & 'UNEXPECTED VALUE:',1i6)
61  CALL plante(1)
62  stop
63  ENDIF
64 !
65  RETURN
66  END
integer icq
Reference concentration formula.
double precision, dimension(:), pointer x
2d coordinates of the mesh
type(bief_obj), target ks
Total bed roughness.
double precision, dimension(:), pointer y
type(bief_obj), target zref
Reference elevation.
subroutine os(OP, X, Y, Z, C, IOPT, INFINI, ZERO)
Definition: os.f:7
type(bief_obj), target ksr
Ripple bed roughness.
type(bief_obj), target ksp
Bed skin roughness.
subroutine zref_gaia
Definition: zref_gaia.f:4
Definition: bief.f:3