qsform.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\qsform.f
00002 !
00062                      SUBROUTINE QSFORM
00063 !                    *****************
00064 !
00065      &(U2D, V2D, TOB, HN, XMVE, TETAP, MU, NPOIN, DM,
00066      & DENS, GRAV, DSTAR, AC, QSC, QSS)
00067 !
00068 !***********************************************************************
00069 ! SISYPHE   V6P2                                   21/07/2011
00070 !***********************************************************************
00071 !
00072 !
00073 !
00074 !
00075 !
00076 !
00077 !
00078 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00079 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00080 !
00081       USE INTERFACE_SISYPHE, EX_QSFORM => QSFORM
00082 !     USE DECLARATIONS_SISYPHE
00083       USE BIEF
00084       IMPLICIT NONE
00085       INTEGER LNG,LU
00086       COMMON/INFO/LNG,LU
00087 !
00088 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00089 !
00090       TYPE(BIEF_OBJ),   INTENT(IN)    :: U2D,V2D,TOB,HN,TETAP,MU
00091       TYPE(BIEF_OBJ),   INTENT(INOUT) :: QSC, QSS
00092       INTEGER,          INTENT(IN)    :: NPOIN
00093       DOUBLE PRECISION, INTENT(IN)    :: XMVE, DM, DENS, GRAV, DSTAR, AC
00094 !
00095 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00096 !
00097 !     INTEGER          :: I
00098 !     DOUBLE PRECISION :: C1, C2, T
00099 !
00100 !======================================================================!
00101 !======================================================================!
00102 !                               PROGRAM                                !
00103 !======================================================================!
00104 !======================================================================!
00105 !
00106 !     EXAMPLE BY VAN RIJN
00107 !
00108 !     C1 = DENS * GRAV * DM
00109 !     C2 = 0.053D0 * SQRT(DM**3*DENS*GRAV) * DSTAR**(-0.3D0)
00110 !
00111 !     DO I = 1, NPOIN
00112 !
00113 !       TRANSPORT STAGE PARAMETER
00114 !
00115 !       IF(TETAP%R(I) .LE. AC) THEN
00116 !         T = 0.D0
00117 !       ELSE
00118 !         T = (TETAP%R(I)-AC)/MAX(AC,1.D-06)
00119 !       ENDIF
00120 !
00121 !       BEDLOAD TRANSPORT RATE
00122 !
00123 !       QSC%R(I) = C2 * T**2.1D0
00124 !       QSS%R(I) = 0.D0
00125 !
00126 !     ENDDO
00127 !
00128 !  FOLLOWING LINES NEED TO BE COMMENTED OUT
00129 !
00130       IF(LNG.EQ.1) WRITE(LU,52)
00131       IF(LNG.EQ.2) WRITE(LU,53)
00132 !
00133 52    FORMAT(/,1X,' STOP :',/
00134      &     ,1X,' LE TAUX DE TRANSPORT DOIT ETRE CALCULE DANS QSFORM')
00135 53    FORMAT(/,1X,'SISYPHE IS STOPPED : ',/
00136      &      ,1X,' SAND TRANSPORT MUST BE CALCULATED IN QSFORM')
00137       CALL PLANTE(1)
00138       STOP
00139 !
00140 !-----------------------------------------------------------------------
00141 !
00142       RETURN
00143       END

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