loinoy.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac2d\loinoy.f
00002 !
00048                      SUBROUTINE LOINOY
00049 !                    *****************
00050 !
00051      &(YAM,YAV,YS,PHI,DEB,G)
00052 !
00053 !***********************************************************************
00054 ! TELEMAC2D   V6P1                                   21/08/2010
00055 !***********************************************************************
00056 !
00057 !
00058 !
00059 !
00060 !
00061 !
00062 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00063 !| DEB            |<--| DISCHARGE OF WEIR
00064 !| G              |-->| GRAVITY.
00065 !| PHI            |-->| DISCHARGE COEFFICIENT OF WEIR.
00066 !| YAM            |-->| UPSTREAM ELEVATION
00067 !| YAV            |-->| DOWNSTREAM ELEVATION
00068 !| YS             |-->| ELEVATION OF WEIR
00069 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00070 !
00071       IMPLICIT NONE
00072       INTEGER LNG,LU
00073       COMMON/INFO/LNG,LU
00074 !
00075 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00076 !
00077       DOUBLE PRECISION, INTENT(INOUT) :: DEB
00078       DOUBLE PRECISION, INTENT(IN)    :: G,YAM,YAV,PHI,YS
00079 !
00080 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00081 !
00082       IF(YAM.LT.YS.AND.YAV.LT.YS) THEN
00083         DEB=0.D0
00084       ELSE
00085         DEB=2.598D0*PHI*SQRT(2.D0*G)*(YAV-YS)*SQRT(YAM-YAV)
00086 !           2.598D0 IS THE INVERSE OF SQRT(1/3)*2/3
00087       ENDIF
00088 !
00089 !-----------------------------------------------------------------------
00090 !
00091       RETURN
00092       END

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