decv21.f

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

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