bedload_hiding_factor.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\bedload_hiding_factor.f
00002 !
00075                      SUBROUTINE BEDLOAD_HIDING_FACTOR
00076 !                    ********************************
00077 !
00078      &(ACLADM, HIDFAC, NPOIN, HIDI, DM, KARIM_HOLLY_YANG, HIDING)
00079 !
00080 !***********************************************************************
00081 ! SISYPHE   V6P1                                   21/07/2011
00082 !***********************************************************************
00083 !
00084 !
00085 !
00086 !
00087 !
00088 !
00089 !
00090 !
00091 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00092 !| ACLADM         |-->| MEAN DIAMETER OF SEDIMENT
00093 !| DM             |-->| SEDIMENT GRAIN DIAMETER
00094 !| HIDFAC         |-->| HIDING FACTOR FORMULAS
00095 !| HIDI           |-->| HIDING FACTOR FOR PARTICULAR SIZE CLASS (HIDFAC =0)
00096 !| HIDING         |-->| HIDING FACTOR CORRECTION
00097 !| NPOIN          |-->| NUMBER OF POINTS
00098 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00099 !
00100       USE INTERFACE_SISYPHE,
00101      &    EX_BEDLOAD_HIDING_FACTOR => BEDLOAD_HIDING_FACTOR
00102       USE BIEF
00103       IMPLICIT NONE
00104       INTEGER LNG,LU
00105       COMMON/INFO/LNG,LU
00106       !
00107       !
00108       ! 2/ GLOBAL VARIABLES
00109       ! -------------------
00110       !
00111       TYPE(BIEF_OBJ),   INTENT(IN)    :: ACLADM
00112       INTEGER,          INTENT(IN)    :: HIDFAC, NPOIN
00113       DOUBLE PRECISION, INTENT(IN)    :: HIDI, DM, KARIM_HOLLY_YANG
00114       TYPE(BIEF_OBJ),   INTENT(INOUT) :: HIDING
00115       !
00116       !
00117       ! 3/ LOCAL VARIABLES
00118       ! ------------------
00119       !
00120       INTEGER          :: J
00121       DOUBLE PRECISION :: C1, C2
00122 !
00123 !======================================================================!
00124 !======================================================================!
00125 !                               PROGRAM                                !
00126 !======================================================================!
00127 !======================================================================!
00128 !
00129 ! *************************** !
00130 ! IA - CONSTANT HIDING FACTOR !
00131 ! *************************** !
00132 !
00133       IF (HIDFAC == 0) THEN
00134 !
00135         CALL OS('X=C     ', X=HIDING, C=HIDI)
00136 !
00137 ! ************************** !
00138 ! IB - EGIAZAROFF FORMULATION !
00139 ! ************************** !
00140 !
00141       ELSEIF (HIDFAC == 1) THEN
00142 !
00143         C1 = LOG10(19.D0)
00144         C2 = 19.D0*DM
00145         DO J = 1, NPOIN
00146           HIDING%R(J) = (C1/LOG10(C2/ACLADM%R(J)))**2
00147         ENDDO
00148 !
00149 ! ********************************** !
00150 ! IC - ASHIDA AND MICHIUE FORMULATION !
00151 ! ********************************** !
00152 !
00153       ELSEIF (HIDFAC == 2) THEN
00154 !
00155         C1 = LOG10(19.D0)
00156         C2 = 19.D0*DM
00157         DO J = 1, NPOIN
00158 !
00159           IF(DM/ACLADM%R(J) >= 0.4D0) THEN
00160             HIDING%R(J) = (C1 / LOG10(C2/ACLADM%R(J)) )**2
00161           ELSE
00162             HIDING%R(J) = 0.85D0*(ACLADM%R(J)/DM)
00163           ENDIF
00164 !
00165         ENDDO
00166 !
00167 ! ************************************* !
00168 ! IE - KARIM, HOLLY AND YANG FORMULATION !
00169 ! ************************************* !
00170 !
00171       ELSEIF (HIDFAC == 4) THEN
00172 !
00173         CALL OS('X=1/Y   ', X=HIDING, Y=ACLADM)
00174         CALL OS('X=CX    ', X=HIDING, C=DM)
00175         CALL OS('X=Y**C  ', X=HIDING, Y=HIDING, C=KARIM_HOLLY_YANG)
00176 !
00177       ELSE
00178 !
00179         IF(LNG.EQ.1) THEN
00180           WRITE(LU,*) 'FORMULE DE MASQUAGE INCONNUE : ',HIDFAC
00181         ENDIF
00182         IF(LNG.EQ.2) THEN
00183           WRITE(LU,*) 'UNKNOWN HIDING FACTOR FORMULA: ',HIDFAC
00184         ENDIF
00185         CALL PLANTE(1)
00186         STOP
00187 !
00188       ENDIF
00189 !
00190 !======================================================================!
00191 !======================================================================!
00192 !
00193       RETURN
00194       END SUBROUTINE BEDLOAD_HIDING_FACTOR

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