decv21.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE decv21
!                   *****************
!
     &(teta,sl,zf,ikle,nelem,nelmax)
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    IDENTIFIES TIDAL FLATS.
!+
!+            DRYING ELEMENT : TETA = 0,
!+
!+            NORMAL ELEMENT : TETA = 1.
!+
!+            THE CRITERION FOR DRYING ELEMENTS IS THAT OF
!+                J.-M. JANIN : BOTTOM ELEVATION OF A POINT IN AN
!+                ELEMENT BEING HIGHER THAN THE FREE SURFACE
!+                ELEVATION OF ANOTHER.
!
!history  J-M HERVOUET (LNH)
!+        09/12/94
!+        V5P1
!+
!
!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
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| IKLE           |-->| CONNECTIVITY TABLE
!| NELEM          |-->| NUMBER OF ELEMENTS IN THE MESH
!| NELMAX         |-->| MAXIMUM NUMBER OF ELEMENTS IN THE MESH
!| SL             |-->| FREE SURFACE
!| TETA           |<--| STATES IF AN ELEMENT IS DRY OR NOT
!| ZF             |-->| BOTTOM ELEVATION
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_special
      USE bief, ex_decv21 => decv21
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER         , INTENT(IN)  :: NELEM,NELMAX
      INTEGER         , INTENT(IN)  :: IKLE(nelmax,*)
      DOUBLE PRECISION, INTENT(OUT) :: TETA(nelem)
      DOUBLE PRECISION, INTENT(IN)  :: SL(*),ZF(*)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER IELEM
!
      DOUBLE PRECISION SL1,SL2,SL3,SL4,ZF1,ZF2,ZF3,ZF4
!
      INTRINSIC max,min
!
!-----------------------------------------------------------------------
!
      CALL ov('X=C     ', x=teta, c=1.d0, dim1=nelem)
!
!-----------------------------------------------------------------------
!
      DO ielem = 1 , nelem
!
        sl1 = sl(ikle(ielem,1))
        sl2 = sl(ikle(ielem,2))
        sl3 = sl(ikle(ielem,3))
        sl4 = sl(ikle(ielem,4))
!
        zf1 = zf(ikle(ielem,1))
        zf2 = zf(ikle(ielem,2))
        zf3 = zf(ikle(ielem,3))
        zf4 = zf(ikle(ielem,4))
!
        IF(max(zf1,zf2,zf3,zf4).GT.min(sl1,sl2,sl3,sl4)) THEN
          teta(ielem) = 0.d0
        ENDIF
!
      ENDDO ! IELEM
!
!-----------------------------------------------------------------------
!
      RETURN
      END