condim_sisyphe.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\condim_sisyphe.f
00002 !
00084                      SUBROUTINE CONDIM_SISYPHE
00085 !                    *************************
00086 !
00087      & (U      , V       , QU    , QV   , H    , ZF , Z ,
00088      &  ESOMT  , THETAWR ,  Q    , HWR  , TWR  ,
00089      &  X      , Y       , NPOIN , AT   , PMAREE)
00090 !
00091 !***********************************************************************
00092 ! SISYPHE   V6P1                                   21/07/2011
00093 !***********************************************************************
00094 !
00095 !
00096 !
00097 !
00098 !
00099 !
00100 !
00101 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00102 !| AT             |-->| TIME
00103 !| ESOMT          |<->| TOTAL BED EVOLUTION
00104 !| H              |<->| WATER DEPTH
00105 !| HWR            |<->| WAVE HEIGHT
00106 !| NPOIN          |-->| NUMBER OF 2D POINTS
00107 !| PMAREE         |-->| TIDAL PERIOD
00108 !| Q              |<->| FLOW RATE
00109 !| THETAWR        |<->| WAVE DIRECTION (WRT OY AXIS)
00110 !| TWR            |<->| WAVE PERIOD
00111 !| X,Y            |-->| COORDINATES
00112 !| Z              |<->| FREE SURFACE
00113 !| ZF             |<->| BED ELEVATION
00114 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00115 !
00116       USE BIEF
00117       USE DECLARATIONS_SISYPHE, ONLY: HW,TW,THETAW
00118 !
00119       IMPLICIT NONE
00120 !
00121       INTEGER LNG,LU
00122       COMMON/INFO/LNG,LU
00123 !
00124 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00125 !
00126       INTEGER, INTENT(IN)::NPOIN
00127 !
00128       DOUBLE PRECISION, INTENT(IN):: X(NPOIN),Y(NPOIN)
00129       DOUBLE PRECISION, INTENT(IN):: AT , PMAREE
00130 ! SEDIMENT
00131       DOUBLE PRECISION, INTENT(INOUT) ::  ZF(NPOIN)
00132       DOUBLE PRECISION, INTENT (INOUT)::  ESOMT(NPOIN)
00133 ! HYDRODYNAMICS
00134       DOUBLE PRECISION, INTENT(INOUT):: Z(NPOIN) , H(NPOIN)
00135       DOUBLE PRECISION, INTENT(INOUT):: U(NPOIN) , V(NPOIN)
00136       DOUBLE PRECISION, INTENT (INOUT)::QU(NPOIN), QV(NPOIN), Q(NPOIN)
00137 ! WAVES
00138       DOUBLE PRECISION, INTENT (INOUT):: HWR(NPOIN) , TWR(NPOIN)
00139       DOUBLE PRECISION, INTENT (INOUT):: THETAWR(NPOIN)
00140 !
00141 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00142 !
00143 !      INTEGER I
00144 !-----------------------------------------------------------------------
00145 !
00146 !     ------------------------
00147 !     THE USER SHOULD BE AWARE
00148 !     ++++++++++++++++++++++++
00149 !
00150 !     SUBROUTINE CONDIM_SISYPHE IS CALLED AT EACH TIME STEP
00151 !     IN ORDER TO IMPOSE A VARIABLE FORCING
00152 !     (TIDAL CURRENT, FOR EXAMPLE)
00153 !
00154 !     IT IS NOT SUFFICIENT TO PRESCRIBE THE FLOW RATE
00155 !     THE MAIN VARIABLES ARE NOW THE 2D FLOW VELOCITY FIELD
00156 !     AND THE FLOW DEPTH
00157 !
00158 !-----------------------------------------------------------------------
00159 !
00160 !     WAVES, EXAMPLE WITH NO WAVES:
00161 !
00162 !     AMPLITUDE = 0
00163 !     CALL OS('X=0     ',X=HW)
00164 !     PERIOD = 1 S
00165 !     CALL OS('X=C     ',X=TW,C=1.D0)
00166 !     ANGLE = 0
00167 !     CALL OS('X=0     ',X=THETAW)
00168 !
00169 !     AFTER SETTING HWR, TWR AND THETAWR, PLEASE ADD:
00170 !
00171 !     HW%TYPR    ='Q'
00172 !     TW%TYPR    ='Q'
00173 !     THETAW%TYPR='Q'
00174 !
00175 !     TO ENABLE THE CONTROL OF WAVE DATA
00176 !
00177 !-----------------------------------------------------------------------
00178 !
00179       RETURN
00180       END SUBROUTINE CONDIM_SISYPHE

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0