q3.f

Go to the documentation of this file.

!                   ****************************
                    DOUBLE PRECISION FUNCTION q3
!                   ****************************
!
     &( i , time , entet )
!
!***********************************************************************
! TELEMAC3D   V7P1
!***********************************************************************
!
!brief    PRESCRIBES THE DISCHARGE FOR FLOW IMPOSED
!+                LIQUID BOUNDARIES.
!
!history  J-M HERVOUET (LNHE)
!+        08/04/09
!+        V6P0
!+   Original version.
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        13/07/2010
!+        V6P0
!+   Translation of French comments within the FORTRAN sources into
!+   English comments
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        21/08/2010
!+        V6P0
!+   Creation of DOXYGEN tags for automated documentation and
!+   cross-referencing of the FORTRAN sources
!
!history  C. COULET (ARTELIA GROUP)
!+        08/11/2011
!+        V6P2
!+   Modification of FCT size due to modification of TRACER numbering
!
!history  J-M HERVOUET (LNHE)
!+        28/09/2015
!+        V7P1
!+   Removing hardcoded array OK.
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| ENTET          |-->| IF YES, LISTING PRINTOUTS ALLOWED
!| I              |-->| NUMBER OF THE LIQUID BOUNDARY
!| TIME           |-->| TIME OF TIME STEP
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE interface_telemac3d, ex_q3 => q3
      USE declarations_special
      USE bief
      USE declarations_telemac
      USE declarations_telemac3d
!
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER          , INTENT(IN) :: I
      DOUBLE PRECISION , INTENT(IN) :: TIME
      LOGICAL          , INTENT(IN) :: ENTET
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      CHARACTER(LEN=9) :: FCT
!
!-----------------------------------------------------------------------
!
!     IF LIQUID BOUNDARY FILE EXISTS, ATTEMPTS TO FIND
!     THE VALUE IN IT. IF YES, OKQ3 REMAINS TO .TRUE. FOR NEXT CALLS
!                      IF  NO, OKQ3 IS SET  TO .FALSE.
!
      IF(okq3(i).AND.t3d_files(t3dimp)%NAME(1:1).NE.' ') THEN
!
!       FCT WILL BE Q(1), Q(2), ETC, Q(99), DEPENDING ON I
        fct='Q(       '
        IF(i.LT.10) THEN
          WRITE(fct(3:3),fmt='(I1)') i
          fct(4:4)=')'
        ELSEIF(i.LT.100) THEN
          WRITE(fct(3:4),fmt='(I2)') i
          fct(5:5)=')'
        ELSE
          WRITE(lu,*) 'Q3 NOT PROGRAMMED FOR MORE
     &                 THAN 99 BOUNDARIES'
          CALL plante(1)
          stop
        ENDIF
        CALL read_fic_frliq(q3,fct,time,t3d_files(t3dimp)%LU,
     &                      entet,okq3(i))
!
      ENDIF
!
      IF(.NOT.okq3(i).OR.t3d_files(t3dimp)%NAME(1:1).EQ.' ') THEN
!
!     Q IS TAKEN FROM THE STEERING FILE, BUT MAY BE CHANGED
!
        q3 = debimp(i)

        ! USER UPDATE
        CALL user_q3(q3, i, time, entet)
!
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END