condis_sisyphe.f

Go to the documentation of this file.

!                   *************************
                    SUBROUTINE condis_sisyphe
!                   *************************
!
     &(constflow)
!
!***********************************************************************
! SISYPHE   V6P1                                   21/07/2011
!***********************************************************************
!
!brief
!
!history  BUI MINH DUC
!+        **/08/2003
!+
!+   DEVELOPED THE SUBROUTINE
!
!history  F. HUVELIN
!+        **/02/2004
!+        V5P4
!+
!
!history  B. MINH DUC; F. HUVELIN
!+
!+        V5P5
!+
!
!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  J,RIEHME (ADJOINTWARE)
!+        November 2016
!+        V7P2
!+   Replaced EXTERNAL statements to parallel functions / subroutines
!+   by the INTERFACE_PARALLEL
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| CONSTFLOW      |<->|LOGICAL, CONSTANT FLOW DISCHARGE OR NOT
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
      USE declarations_telemac
      USE declarations_sisyphe
      !
      USE declarations_special
      USE interface_parallel, ONLY : p_isum
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      LOGICAL, INTENT(INOUT) :: CONSTFLOW
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
!3/ LOCAL VARIABLES
!--------------------
!
      INTEGER          :: NZFMAX, I
      DOUBLE PRECISION :: ZFMAX
!
!=======================================================================!
!=======================================================================!
!                               PROGRAM                                 !
!=======================================================================!
!=======================================================================!
!
      nzfmax = 0
!
      IF(constflow) THEN
        CALL os('X=X+Y   ', x=ecpl, y=e)
!
        DO i=1,npoin
          zfmax = abs(ecpl%R(i)) - crit_cfd*hcpl%R(i)
          IF (zfmax.GT.1.d-8) nzfmax=nzfmax+1
        ENDDO
!
        IF (ncsize.GT.1) THEN
          nzfmax=p_isum(nzfmax)
          CALL parcom(ecpl,2,mesh)
        ENDIF
!
        IF (nzfmax.GE.1) constflow = .false.
      ENDIF
!
      IF(.NOT.constflow) THEN
        CALL os('X=C     ', x=ecpl, c=0.d0)
        CALL os('X=Y     ', x=hcpl, y=hn)
!
        IF (ncsize.GT.1) THEN
          CALL parcom(ecpl,2,mesh)
          CALL parcom(hcpl,2,mesh)
        ENDIF
      ENDIF
!
!=======================================================================!
!=======================================================================!
!
      RETURN
      END