bedload_hunz_meyer.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\bedload_hunz_meyer.f
00002 !
00079                      SUBROUTINE BEDLOAD_HUNZ_MEYER
00080 !                    *****************************
00081 !
00082      &  (TOB, MU, ACLADM, UNLADM, NPOIN, DENS, XMVE, GRAV, DM, AC,
00083      &    TETAP, AHUNZI, ACP, HIDING, QSC)
00084 !
00085 !***********************************************************************
00086 ! SISYPHE   V6P1                                   21/07/2011
00087 !***********************************************************************
00088 !
00089 !
00090 !
00091 !
00092 !
00093 !
00094 !
00095 !
00096 !
00097 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00098 !| AC             |<->| CRITICAL SHIELDS PARAMETER
00099 !| ACLADM         |-->| MEAN DIAMETER OF SEDIMENT
00100 !| ACP            |<->| MODIFIED SHIELDS PARAMETER
00101 !| AHUNZI         |<->| COEFFICIENT OF HUNZIKER FORMULA
00102 !| DENS           |-->| RELATIVE DENSITY
00103 !| DM             |-->| SEDIMENT GRAIN DIAMETER
00104 !| GRAV           |-->| ACCELERATION OF GRAVITY
00105 !| HIDING         |-->| HIDING FACTOR CORRECTION
00106 !| MU             |<->| CORRECTION FACTOR FOR BED ROUGHNESS
00107 !| NPOIN          |-->| NUMBER OF POINTS
00108 !| QSC            |<->| BED LOAD TRANSPORT
00109 !| TETAP          |<->| DIMENSIONLESS BED SHEAR STRESS
00110 !| TOB            |<->| BED SHEAR STRESS (TOTAL FRICTION)
00111 !| UNLADM         |-->| MEAN DIAMETER OF ACTIVE STRATUM LAYER
00112 !| XMVE           |-->| FLUID DENSITY
00113 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00114 !
00115       USE INTERFACE_SISYPHE,
00116      &          EX_BEDLOAD_HUNZ_MEYER => BEDLOAD_HUNZ_MEYER
00117       USE BIEF
00118       IMPLICIT NONE
00119       INTEGER LNG,LU
00120       COMMON/INFO/LNG,LU
00121 !
00122       ! 2/ GLOBAL VARIABLES
00123       ! -------------------
00124       TYPE(BIEF_OBJ),   INTENT(IN)    :: TOB, MU, ACLADM, UNLADM
00125       INTEGER,          INTENT(IN)    :: NPOIN
00126       DOUBLE PRECISION, INTENT(IN)    :: DENS, XMVE, GRAV, DM, AC
00127       TYPE(BIEF_OBJ),   INTENT(INOUT) :: TETAP, AHUNZI ! WORK ARRAY T1, T2
00128       TYPE(BIEF_OBJ),   INTENT(INOUT) :: ACP           ! WORK ARRAY T3
00129       TYPE(BIEF_OBJ),   INTENT(INOUT)   :: HIDING, QSC
00130 !
00131       ! 3/ LOCAL VARIABLES
00132       ! ------------------
00133       INTEGER          :: I
00134       DOUBLE PRECISION :: C1, C2
00135 !======================================================================!
00136 !======================================================================!
00137 !                               PROGRAM                                !
00138 !======================================================================!
00139 !======================================================================!
00140       ! ************************************* !
00141       ! I - ADIMENSIONAL SKIN STRESS          !
00142       ! ************************************* !
00143       C1 = 1.D0/(DENS*XMVE*GRAV*DM)
00144       C2 = 5.D0*SQRT(GRAV*DENS*DM**3)
00145       CALL OS('X=CYZ   ', X=TETAP, Y=TOB, Z=MU, C=C1)
00146 !
00147 !     CHANGED BY JMH ON 28/10/2009 AFTER MODIFICATIONS BY
00148 !     REBEKKA KOPMANN TRANSMITTED BY JACEK JANKOWSKI
00149 !     CALL OS('X=+(Y,C)', X=TETAP , Y=TETAP, C= 1.D-06 )
00150       CALL OS('X=+(Y,C)', X=TETAP , Y=TETAP, C= 1.D-02 )
00151 !
00152       CALL OS('X=Y**C  ', X=AHUNZI, Y=TETAP, C=-1.5D0  )
00153       CALL OS('X=CX    ', X=AHUNZI, C= 0.011D0)
00154       CALL OS('X=X+C   ', X=AHUNZI, C=-0.3D0  )
00155 !
00156 ! RK COMMENT:
00157 ! BEWARE: AHUNZI CAN BECOME SO LARGE THAT THE HIDING FACTOR BECOMES
00158 ! INFINITE; HUNZIKER HIMSELF SUGGESTS THAT IT BE CAPPED TO 2.3.
00159 ! THIS INITIALLY ADOPTS A LIMIT OF APPROXIMATELY 10.
00160 ! (WHICH IS APPARENT IN TETAP BEING SET TO VALUES .GE. 0.01)
00161 !
00162 !     REMARK BY JMH: I WOULD STRONGLY RECOMMEND A SINGLE LOOP
00163 !                    WITH THE WHOLE FORMULA, INSTEAD OF PILING
00164 !                    UP CALLS TO OS
00165 !
00166       DO I = 1, NPOIN
00167         HIDING%R(I) = (DM/ACLADM%R(I))**(-AHUNZI%R(I))
00168       ENDDO
00169       ! ************************************************* !
00170       ! IV - CORRECTS THE ADIMENSIONAL CRITICAL STRESS    !
00171       ! ************************************************* !
00172       CALL OS('X=Y/Z   ', X=ACP, Y=UNLADM, Z=ACLADM)
00173       CALL OS('X=Y**C  ', X=ACP, Y=ACP   , C=0.33D0)
00174       CALL OS('X=CX    ', X=ACP, C=AC)
00175       ! ********************* !
00176       ! V - TRANSPORT RATE    !
00177       ! ********************* !
00178       CALL OS('X=Y-Z   ', X=QSC, Y=TETAP , Z=ACP )
00179       CALL OS('X=+(Y,C)', X=QSC, Y=QSC   , C=0.D0)
00180       CALL OS('X=XY    ', X=QSC, Y=HIDING)
00181       CALL OS('X=Y**C  ', X=QSC, Y=QSC   , C=1.5D0)
00182 
00183       CALL OS('X=CX    ', X=QSC, C=C2)
00184 !======================================================================!
00185 !======================================================================!
00186       RETURN
00187       END

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