gestio.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac2d\gestio.f
00002 !
00054                      SUBROUTINE GESTIO
00055 !                    *****************
00056 !
00057      &(U,V,C,T,AK,EP,UTILD,VTILD,CTILD,TTILD,AKTILD,EPTILD,
00058      & TRAC,PROPA,CONVV,ITURB,IETAPE)
00059 !
00060 !***********************************************************************
00061 ! TELEMAC2D   V6P1                                   21/08/2010
00062 !***********************************************************************
00063 !
00064 !
00065 !
00066 !
00067 !
00068 !
00069 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00070 !| AK             |<->| TURBULENT KINETIC ENERGY
00071 !| AKTILD         |-->| TURBULENT KINETIC BEFORE CURRENT STEP
00072 !| C              |<->| CELERITY
00073 !| CONVV          |-->| LOGICAL ARRAY SAYING IF A VARIABLE IS ADVECTED OR NOT
00074 !| CTILD          |-->| CELERITY BEFORE CURRENT STEP
00075 !| EP             |<->| TURBULENT ENERGY DISSIPASSION
00076 !| EPTILD         |-->| TURBULENT ENERGY DISSIPASSION BEFORE CURRENT STEP
00077 !| IETAPE         |-->| FRACTIONAL STEP NUMBER
00078 !| ITURB          |-->| TURBULENCE MODEL 1: LAMINAR, CONSTANT COEFFICIENT
00079 !|                |   |                  2: MIXING LENGTH
00080 !|                |   |                  3: K-EPSILON
00081 !| PROPA          |-->| IF PROPA=.FALSE. : NO PROPAGATION STEP.
00082 !| T              |<->| BLOCK OF TRACERS
00083 !| TRAC           |-->| LOGICAL, YES IF THERE ARE TRACERS
00084 !| TTILD          |-->| BLOCK OF TRACERS BEFORE CURRENT STEP
00085 !| U              |<->| X-COMPONENT OF VELOCITY
00086 !| V              |<->| Y-COMPONENT OF VELOCITY
00087 !| UTILD          |-->| X-COMPONENT OF VELOCITY BEFORE CURRENT STEP
00088 !| VTILD          |-->| Y-COMPONENT OF VELOCITY BEFORE CURRENT STEP
00089 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00090 !
00091       USE BIEF
00092 !
00093       IMPLICIT NONE
00094       INTEGER LNG,LU
00095       COMMON/INFO/LNG,LU
00096 !
00097 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00098 !
00099       INTEGER, INTENT(IN)           :: ITURB,IETAPE
00100       LOGICAL, INTENT(IN)           :: TRAC,CONVV(4),PROPA
00101       TYPE(BIEF_OBJ), INTENT(IN)    :: T,AK,EP
00102       TYPE(BIEF_OBJ), INTENT(INOUT) :: U,V,C
00103       TYPE(BIEF_OBJ), INTENT(INOUT) :: UTILD,VTILD,CTILD,TTILD
00104       TYPE(BIEF_OBJ), INTENT(INOUT) :: AKTILD,EPTILD
00105 !
00106 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00107 !
00108 !----------------------------------------------------------------------
00109 !    ARRAYS DEPENDING ON THE SELECTED EQUATIONS
00110 !-----------------------------------------------------------------------
00111 !
00112 !    ADVECTION
00113 !
00114       IF(IETAPE.EQ.3) THEN
00115 !
00116         IF(.NOT.CONVV(1)) THEN
00117           CALL OS( 'X=Y     ' , X=UTILD , Y=U )
00118           CALL OS( 'X=Y     ' , X=VTILD , Y=V )
00119         ENDIF
00120         IF(.NOT.CONVV(2)) THEN
00121           CALL OS( 'X=Y     ' , X=CTILD , Y=C )
00122         ENDIF
00123         IF(TRAC.AND.(.NOT.CONVV(3))) THEN
00124           CALL OS( 'X=Y     ' , X=TTILD , Y=T )
00125         ENDIF
00126         IF(ITURB.EQ.3.AND.(.NOT.CONVV(4))) THEN
00127           CALL OS( 'X=Y     ' , X=AKTILD , Y=AK )
00128           CALL OS( 'X=Y     ' , X=EPTILD , Y=EP )
00129         ENDIF
00130 !
00131       ENDIF
00132 !
00133 !-----------------------------------------------------------------------
00134 !
00135 !    PROPAGATION
00136 !
00137       IF(IETAPE.EQ.6) THEN
00138 !
00139             IF(.NOT.PROPA) THEN
00140 !
00141                    CALL OS( 'X=Y     ' , X=U , Y=UTILD )
00142                    CALL OS( 'X=Y     ' , X=V , Y=VTILD )
00143                    CALL OS( 'X=Y     ' , X=C , Y=CTILD )
00144 !
00145             ENDIF
00146 !
00147       ENDIF
00148 !
00149 !-----------------------------------------------------------------------
00150 !
00151       RETURN
00152       END

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