condis_sisyphe.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\condis_sisyphe.f
00002 !
00080                      SUBROUTINE CONDIS_SISYPHE
00081 !                    *************************
00082 !
00083      &(CONSTFLOW)
00084 !
00085 !***********************************************************************
00086 ! SISYPHE   V6P1                                   21/07/2011
00087 !***********************************************************************
00088 !
00089 !
00090 !
00091 !
00092 !
00093 !
00094 !
00095 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00096 !| CONSTFLOW      |<->|LOGICAL, CONSTANT FLOW DISCHARGE OR NOT
00097 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00098 !
00099       USE BIEF
00100       USE DECLARATIONS_TELEMAC
00101       USE DECLARATIONS_SISYPHE
00102       !
00103       IMPLICIT NONE
00104       INTEGER LNG,LU
00105       COMMON/INFO/LNG,LU
00106 !
00107 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00108 !
00109       LOGICAL, INTENT(INOUT) :: CONSTFLOW
00110 !
00111 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00112 !
00113       INTEGER, EXTERNAL      :: P_ISUM
00114 !
00115 !3/ LOCAL VARIABLES
00116 !--------------------
00117 !
00118       INTEGER          :: NZFMAX, I
00119       DOUBLE PRECISION :: ZFMAX, C
00120 !
00121 !=======================================================================!
00122 !=======================================================================!
00123 !                               PROGRAM                                 !
00124 !=======================================================================!
00125 !=======================================================================!
00126 !
00127       NZFMAX = 0
00128 !
00129       IF(CONSTFLOW) THEN
00130         CALL OS('X=X+Y   ', ECPL, E, S, C)
00131 !
00132         DO I=1,NPOIN
00133           ZFMAX = ABS(ECPL%R(I)) - CRIT_CFD*HCPL%R(I)
00134           IF (ZFMAX.GT.1.D-8) NZFMAX=NZFMAX+1
00135         ENDDO
00136 !
00137         IF (NCSIZE.GT.1) THEN
00138           NZFMAX=P_ISUM(NZFMAX)
00139           CALL PARCOM(ECPL,2,MESH)
00140         ENDIF
00141 !
00142         IF (NZFMAX.GE.1) CONSTFLOW = .FALSE.
00143       ENDIF
00144 !
00145       IF(.NOT.CONSTFLOW) THEN
00146         CALL OS('X=C     ', ECPL,  S, S, 0.D0)
00147         CALL OS('X=Y     ', HCPL, HN, S,    C)
00148 !
00149         IF (NCSIZE.GT.1) THEN
00150           CALL PARCOM(ECPL,2,MESH)
00151           CALL PARCOM(HCPL,2,MESH)
00152         ENDIF
00153       ENDIF
00154 !
00155 !=======================================================================!
00156 !=======================================================================!
00157 !
00158       RETURN
00159       END

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