bedload_engel_cc.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\bedload_engel_cc.f
00002 !
00067                      SUBROUTINE BEDLOAD_ENGEL_CC
00068 !                    ***************************
00069 !
00070      &(TETAP,CF,NPOIN,GRAV,DM,DENS,TETA,QSC)
00071 !
00072 !***********************************************************************
00073 ! SISYPHE   V6P1                                   21/07/2011
00074 !***********************************************************************
00075 !
00076 !
00077 !
00078 !
00079 !
00080 !
00081 !
00082 !
00083 !
00084 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00085 !| CF             |-->| QUADRATIC FRICTION COEFFICIENT
00086 !| DENS           |-->| RELATIVE DENSITY
00087 !| DM             |-->| SEDIMENT GRAIN DIAMETER
00088 !| GRAV           |-->| ACCELERATION OF GRAVITY
00089 !| NPOIN          |-->| NUMBER OF POINTS
00090 !| QSC            |<->| BED LOAD TRANSPORT
00091 !| TETA           |<->| DIMENSIONLESS BED SHEAR STRESS
00092 !| TETAP          |-->| ADIMENSIONAL SKIN FRICTION
00093 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00094 !
00095       USE INTERFACE_SISYPHE,EX_BEDLOAD_ENGEL_CC => BEDLOAD_ENGEL_CC
00096       USE BIEF
00097 !
00098       IMPLICIT NONE
00099       INTEGER LNG,LU
00100       COMMON/INFO/LNG,LU
00101 !
00102 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00103 !
00104       TYPE(BIEF_OBJ),   INTENT(IN)    :: TETAP,CF
00105       INTEGER,          INTENT(IN)    :: NPOIN
00106       DOUBLE PRECISION, INTENT(IN)    :: GRAV, DM, DENS
00107       TYPE(BIEF_OBJ),   INTENT(INOUT) :: TETA! WORK ARRAY T1
00108       TYPE(BIEF_OBJ),   INTENT(INOUT) :: QSC
00109 !
00110 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00111 !
00112       INTEGER          :: I
00113       DOUBLE PRECISION :: CENGEL
00114 !
00115       INTRINSIC SQRT
00116 !
00117 !======================================================================!
00118 !======================================================================!
00119 !                               PROGRAM                                !
00120 !======================================================================!
00121 !======================================================================!
00122 !
00123 !     ADIMENSIONAL SKIN STRESS: TETAP
00124 !
00125 !     ADIMENSIONAL TOTAL STRESS
00126 !
00127       DO I = 1, NPOIN
00128         IF(TETAP%R(I) <= 0.06D0) THEN
00129           TETA%R(I) = 0.D0
00130         ELSEIF(TETAP%R(I) <  0.384D0) THEN
00131           TETA%R(I) = SQRT( 2.5D0 * (TETAP%R(I) - 0.06D0))
00132         ELSEIF(TETAP%R(I) <  1.080D0) THEN
00133           TETA%R(I) = 1.066D0 * TETAP%R(I)**0.176D0
00134         ELSE
00135           TETA%R(I) = TETAP%R(I)
00136         ENDIF
00137       ENDDO
00138 !
00139 !     BEDLOAD TRANSPORT
00140 !
00141       CENGEL = 0.1D0*SQRT(DENS*GRAV*DM**3)
00142       DO I=1,NPOIN
00143         QSC%R(I)=CENGEL*SQRT(TETA%R(I)**5)/MAX(CF%R(I),1.D-6)
00144       ENDDO
00145 !
00146 !-----------------------------------------------------------------------
00147 !
00148       RETURN
00149       END

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