decv11.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\decv11.f
00002 !
00060                      SUBROUTINE DECV11
00061 !                    *****************
00062 !
00063      &(TETA,SL,ZF,IKLE,NELEM,NELMAX)
00064 !
00065 !***********************************************************************
00066 ! BIEF   V6P1                                   21/08/2010
00067 !***********************************************************************
00068 !
00069 !
00070 !
00071 !
00072 !
00073 !
00074 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00075 !| IKLE           |-->| CONNECTIVITY TABLE
00076 !| NELEM          |-->| NUMBER OF ELEMENTS IN THE MESH
00077 !| NELMAX         |-->| MAXIMUM NUMBER OF ELEMENTS IN THE MESH
00078 !| SL             |-->| FREE SURFACE
00079 !| TETA           |<--| STATES IF AN ELEMENT IS DRY OR NOT
00080 !| ZF             |-->| BOTTOM ELEVATION
00081 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00082 !
00083       IMPLICIT NONE
00084       INTEGER LNG,LU
00085       COMMON/INFO/LNG,LU
00086 !
00087 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00088 !
00089       INTEGER         , INTENT(IN)  :: NELEM,NELMAX
00090       INTEGER         , INTENT(IN)  :: IKLE(NELMAX,*)
00091       DOUBLE PRECISION, INTENT(OUT) :: TETA(NELEM)
00092       DOUBLE PRECISION, INTENT(IN)  :: SL(*),ZF(*)
00093 !
00094 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00095 !
00096       INTEGER IELEM
00097 !
00098       DOUBLE PRECISION SL1,SL2,SL3,ZF1,ZF2,ZF3
00099 !
00100       INTRINSIC MAX,MIN
00101 !
00102 !-----------------------------------------------------------------------
00103 !
00104       DO IELEM = 1 , NELEM
00105 !
00106         SL1 = SL(IKLE(IELEM,1))
00107         SL2 = SL(IKLE(IELEM,2))
00108         SL3 = SL(IKLE(IELEM,3))
00109 !
00110         ZF1 = ZF(IKLE(IELEM,1))
00111         ZF2 = ZF(IKLE(IELEM,2))
00112         ZF3 = ZF(IKLE(IELEM,3))
00113 !
00114         IF(MAX(ZF1,ZF2,ZF3).GT.MIN(SL1,SL2,SL3)) THEN
00115           TETA(IELEM) = 0.D0
00116         ELSE
00117           TETA(IELEM) = 1.D0
00118         ENDIF
00119 !
00120       ENDDO
00121 !
00122 !-----------------------------------------------------------------------
00123 !
00124       RETURN
00125       END

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