The TELEMAC-MASCARET system  trunk
bedload_direction_gaia.f
Go to the documentation of this file.
1 ! *********************************
2  SUBROUTINE bedload_direction_gaia
3 ! *********************************
4 !
5  & (u2d, v2d, npoin, pi, thetac)
6 !
7 !***********************************************************************
8 ! GAIA
9 !***********************************************************************
10 !
12 !
13 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
19 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
20 !
21  USE interface_gaia,
22  & ex_bedload_direction => bedload_direction_gaia
23  USE bief
25  IMPLICIT NONE
26 !
27 ! GLOBAL VARIABLES
28 ! -------------------
29  TYPE(bief_obj), INTENT(IN) :: U2D, V2D
30  INTEGER, INTENT(IN) :: NPOIN
31  DOUBLE PRECISION, INTENT(IN) :: PI
32  TYPE(bief_obj), INTENT(INOUT) :: THETAC
33 !
34 ! LOCAL VARIABLES
35 ! ------------------
36  INTEGER :: I
37  DOUBLE PRECISION, PARAMETER :: LOCAL_ZERO = 1.d-6
38 !
39 !======================================================================!
40 !======================================================================!
41 !
42  DO i = 1, npoin
43  IF (abs(u2d%R(i)) <= local_zero) THEN
44  IF (v2d%R(i) < = local_zero) THEN
45  thetac%R(i) = -pi*0.5d0
46  ELSE
47  thetac%R(i) = pi*0.5d0
48  ENDIF
49  ELSE
50  thetac%R(i) = atan(v2d%R(i) / u2d%R(i))
51  IF (u2d%R(i) < 0.d0) THEN
52  thetac%R(i) = pi + thetac%R(i)
53  ENDIF
54  ENDIF
55  END DO
56 !======================================================================!
57 !======================================================================!
58  RETURN
59  END SUBROUTINE bedload_direction_gaia
subroutine bedload_direction_gaia(U2D, V2D, NPOIN, PI, THETAC)
Definition: bief.f:3