The TELEMAC-MASCARET system  trunk
bedload_engel_gaia.f
Go to the documentation of this file.
1 ! *****************************
2  SUBROUTINE bedload_engel_gaia
3 ! *****************************
4 !
5  &(tob,cf,dens,grav,dcla,xmve,teta,qsc,xmvs)
6 !
7 !***********************************************************************
8 ! GAIA V7P3 10/01/2018
9 !***********************************************************************
10 !
12 !
13 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
23 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
24 !
25  USE interface_gaia, ex_bedload_engel => bedload_engel_gaia
26  USE bief
28  IMPLICIT NONE
29 !
30 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
31 !
32  TYPE(bief_obj), INTENT(IN) :: TOB, CF
33  DOUBLE PRECISION, INTENT(IN) :: DENS, GRAV, DCLA, XMVE, XMVS
34  TYPE(bief_obj), INTENT(INOUT) :: TETA ! WORK ARRAY T1
35  TYPE(bief_obj), INTENT(INOUT) :: QSC
36 !
37 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
38 !
39  INTEGER I
40  DOUBLE PRECISION :: CENGEL, C1
41 !
42 !======================================================================!
43 ! PROGRAM !
44 !======================================================================!
45 !
46  c1 = 1.d0/(dens*xmve*grav*dcla)
47  cengel = 0.05d0*sqrt(dens*grav*dcla**3)
48 ! CALL OS('X=CY ', X=TETA, Y=TOB , C=C1)
49 ! CALL OS('X=Y**C ', X=TETA, Y=TETA, C=5.D0/2.D0)
50 !
51 ! CALL OS('X=+(Y,C)', X=QSC , Y=CF , C=1.D-06)
52 ! CALL OS('X=1/Y ', X=QSC , Y=QSC)
53 ! CALL OS('X=CXY ', X=QSC , Y=TETA, C=CENGEL)
54  DO i=1,qsc%DIM1
55 ! TOTAL NON DIMENSIONAL STRESS = SQRT(C1*TOB%R(I))**5
56 ! BEDLOAD TRANSPORT
57  qsc%R(i)=cengel*sqrt(c1*tob%R(i))**5/max(cf%R(i),1.d-6)
58  ENDDO
59 !
60 ! SOLID DISCHARGE IS TRANSFORMED IN [kg/(m*s)]
61 !
62  CALL os('X=CX ', x=qsc, c=xmvs)
63 !======================================================================
64 !
65  RETURN
66  END
subroutine bedload_engel_gaia(TOB, CF, DENS, GRAV, DCLA, XMVE, TETA, QSC, XMVS)
subroutine os(OP, X, Y, Z, C, IOPT, INFINI, ZERO)
Definition: os.f:7
Definition: bief.f:3