decvrt.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE decvrt
!                   *****************
!
     &(teta,sl,zf,mesh)
!
!***********************************************************************
! 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
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| MESH           |-->| MESH STRUCTURE
!| SL             |-->| FREE SURFACE
!| TETA           |<--| STATES IF AN ELEMENT IS DRY OR NOT
!| ZF             |-->| BOTTOM ELEVATION
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_decvrt => decvrt
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      TYPE(bief_obj) , INTENT(INOUT) :: TETA
      TYPE(bief_obj) , INTENT(IN)    :: SL,ZF
      TYPE(bief_mesh), INTENT(IN)    :: MESH
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER NELEM,NELMAX,IELMS,IELMZ
!
!-----------------------------------------------------------------------
!
      ielms=sl%ELM
      ielmz=zf%ELM
!
!-----------------------------------------------------------------------
!
!  DEPLOYS THE MESH STRUCTURE
!
      nelem = mesh%NELEM
      nelmax= mesh%NELMAX
!
!-----------------------------------------------------------------------
!
      IF(ielms.EQ.11.AND.ielmz.EQ.11) THEN
!
        CALL decv11(teta%R,sl%R,zf%R,mesh%IKLE%I,nelem,nelmax)
!
      ELSEIF((ielms.EQ.21.AND.ielmz.EQ.21).OR.
     &       (ielms.EQ.12.AND.ielmz.EQ.12)      ) THEN
!
        CALL decv21(teta%R,sl%R,zf%R,mesh%IKLE%I,nelem,nelmax)
!
      ELSE
!
        WRITE(lu,11) ielms,ielmz
11      FORMAT(1x,'DECVRT : DISCRETIZATION NOT IMPLEMENTED:',i6,' ',i6)
        CALL plante(1)
        stop
!
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END