The TELEMAC-MASCARET system  trunk
gaia_bedload_update.f
Go to the documentation of this file.
1 ! ******************************
2  SUBROUTINE gaia_bedload_update
3 ! ******************************
4 !
5 !***********************************************************************
6 ! GAIA
7 !***********************************************************************
8 !
10 !
11 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
12 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
13  USE interface_gaia
14  USE bief
19 !
20  IMPLICIT NONE
21 !
22 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
23 !
24 !
25 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
26 !
27  INTEGER I
28 !
29 !-----------------------------------------------------------------------
30 !
31 !
32 ! RECONSTITUTES THE BEDLOAD AND/OR SUSPENSION DATA
33 ! -----------------------------------------------------
34 !
35  IF(debug.GT.0) WRITE(lu,*) 'QS_RESULT'
36 !
37  DO i = 1, nsicla
38  CALL os('X=Y ',x=t1,y=qsclxc%ADR(i)%P)
39  CALL os('X=Y ',x=t2,y=qsclyc%ADR(i)%P)
40  CALL os('X=N(Y,Z)', x=qscl%ADR(i)%P,y=t1,z=t2)
41  IF(i.EQ.1) THEN
42  CALL os('X=Y ', x=qsx, y=t1)
43  CALL os('X=Y ', x=qsy, y=t2)
44  ELSE
45  CALL os('X=X+Y ', x=qsx, y=t1)
46  CALL os('X=X+Y ', x=qsy, y=t2)
47  ENDIF
48  ENDDO
49  CALL os('X=N(Y,Z)', x=qs, y=qsx, z=qsy)
50  IF(debug.GT.0) WRITE(lu,*) 'END_QS_RESULT'
51 !
52 !-----------------------------------------------------------------------
53 !
54  RETURN
55  END
type(bief_obj), target qs
Total solid discharge (bedload+suspension)
subroutine gaia_bedload_update
double precision, dimension(:), pointer x
2d coordinates of the mesh
type(bief_obj), target qscl
Total transport rate for a sediment class : bedload+suspended load.
type(bief_obj), target qsx
Solid discharge, along x and y.
double precision, dimension(:), pointer y
integer, target nsicla
Number of sediment classes of bed material (less than NISCLM)
type(bief_obj), pointer t2
type(bief_obj), target z
Free surface elevation.
type(bief_obj), target qsclyc
type(bief_obj), pointer t1
Aliases for work vectors in tb.
type(bief_obj), target qsy
type(bief_obj), target qsclxc
Bedload transport rate in the x and y direction for a sediment class [kg*(m-1*s-1)].
integer debug
Debugger.
subroutine os(OP, X, Y, Z, C, IOPT, INFINI, ZERO)
Definition: os.f:7
Definition: bief.f:3