bedload_engel.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\bedload_engel.f
00002 !
00062                      SUBROUTINE BEDLOAD_ENGEL
00063 !                    ************************
00064 !
00065      &(TOB,CF,DENS,GRAV,DM,XMVE,TETA,QSC)
00066 !
00067 !***********************************************************************
00068 ! SISYPHE   V6P2                                   21/07/2011
00069 !***********************************************************************
00070 !
00071 !
00072 !
00073 !
00074 !
00075 !
00076 !
00077 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00078 !| CF             |-->| QUADRATIC FRICTION COEFFICIENT
00079 !| DENS           |-->| RELATIVE DENSITY
00080 !| DM             |-->| SEDIMENT GRAIN DIAMETER
00081 !| GRAV           |-->| ACCELERATION OF GRAVITY
00082 !| QSC            |<->| BED LOAD TRANSPORT
00083 !| TETA           |<->| DIMENSIONLESS BED SHEAR STRESS
00084 !| TOB            |<->| BED SHEAR STRESS (TOTAL FRICTION)
00085 !| XMVE           |-->| FLUID DENSITY
00086 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00087 !
00088       USE INTERFACE_SISYPHE, EX_BEDLOAD_ENGEL => BEDLOAD_ENGEL
00089       USE BIEF
00090       IMPLICIT NONE
00091       INTEGER LNG,LU
00092       COMMON/INFO/LNG,LU
00093 !
00094 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00095 !
00096       TYPE(BIEF_OBJ),   INTENT(IN)    :: TOB, CF
00097       DOUBLE PRECISION, INTENT(IN)    :: DENS, GRAV, DM, XMVE
00098       TYPE(BIEF_OBJ),   INTENT(INOUT) :: TETA ! WORK ARRAY T1
00099       TYPE(BIEF_OBJ),   INTENT(INOUT) :: QSC
00100 !
00101 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00102 !
00103       INTEGER I
00104       DOUBLE PRECISION :: CENGEL, C1
00105 !
00106 !======================================================================!
00107 !                               PROGRAM                                !
00108 !======================================================================!
00109 !
00110       C1 = 1.D0/(DENS*XMVE*GRAV*DM)
00111       CENGEL = 0.1D0*SQRT(DENS*GRAV*DM**3)
00112 !     CALL OS('X=CY    ', X=TETA, Y=TOB , C=C1)
00113 !     CALL OS('X=Y**C  ', X=TETA, Y=TETA, C=5.D0/2.D0)
00114 !
00115 !     CALL OS('X=+(Y,C)', X=QSC , Y=CF  , C=1.D-06)
00116 !     CALL OS('X=1/Y   ', X=QSC , Y=QSC)
00117 !     CALL OS('X=CXY   ', X=QSC , Y=TETA, C=CENGEL)
00118       DO I=1,QSC%DIM1
00119 !       TOTAL NON DIMENSIONAL STRESS = SQRT(C1*TOB%R(I))**5
00120 !       BEDLOAD TRANSPORT
00121         QSC%R(I)=CENGEL*SQRT(C1*TOB%R(I))**5/MAX(CF%R(I),1.D-6)
00122       ENDDO
00123 !
00124 !======================================================================
00125 !
00126       RETURN
00127       END

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