The TELEMAC-MASCARET system  trunk
user_forcing_gaia.f
Go to the documentation of this file.
1 ! ****************************
2  SUBROUTINE user_forcing_gaia
3 ! ****************************
4 !
5 !***********************************************************************
6 ! GAIA
7 !***********************************************************************
8 !
16 !
18 !! EXAMPLE WITH NO WAVES:! AMPLITUDE = 0! PERIOD = 1 S! ANGLE = 0
31 !
32 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
33 !
34  USE bief
36 !
38  IMPLICIT NONE
39 !
40 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
41 !
42 ! REMOVE THIS RETURN IF YOU WANT TO ENTER INTO THIS USER SUBROUTINE
43  IF (.true.) RETURN
44 !
45 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
46 !
47 !-----------------------------------------------------------------------
48 !
49 ! ------------------------
50 ! THE USER SHOULD BE AWARE
51 ! ++++++++++++++++++++++++
52 !
53 ! SUBROUTINE USER_FORCING_GAIA IS CALLED AT EACH TIME STEP
54 ! IN ORDER TO IMPOSE A VARIABLE FORCING
55 ! (TIDAL CURRENT, FOR EXAMPLE)
56 !
57 ! IT IS NOT SUFFICIENT TO PRESCRIBE THE FLOW RATE
58 ! THE MAIN VARIABLES ARE NOW THE 2D FLOW VELOCITY FIELD
59 ! AND THE FLOW DEPTH
60 !
61 !-----------------------------------------------------------------------
62 !
63 ! EXAMPLE 1: WITH NO WAVES
64 !
65  IF (.false.) THEN
66 ! AMPLITUDE = 0
67  CALL os('X=0 ',x=hw)
68 ! PERIOD = 1 S
69  CALL os('X=C ',x=tw,c=1.d0)
70 ! ANGLE = 0
71  CALL os('X=0 ',x=thetaw)
72 
73 ! AFTER SETTING HWR, TWR AND THETAWR, PLEASE ADD:
74 
75  hw%TYPR ='Q'
76  tw%TYPR ='Q'
77  thetaw%TYPR='Q'
78 
79 ! TO ENABLE THE CONTROL OF WAVE DATA
80  ENDIF
81 !
82 !-----------------------------------------------------------------------
83 !
84  RETURN
85  END SUBROUTINE user_forcing_gaia
double precision, dimension(:), pointer x
2d coordinates of the mesh
type(bief_obj), target thetaw
Wave direction (deg wrt ox axis) !!!!!some say oy axis!!!!!
type(bief_obj), target hw
Significant wave height.
subroutine user_forcing_gaia
type(bief_obj), target tw
Mean wave period.
subroutine os(OP, X, Y, Z, C, IOPT, INFINI, ZERO)
Definition: os.f:7
Definition: bief.f:3