massed.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac3d\massed.f
00002 !
00085                      SUBROUTINE MASSED
00086 !                    *****************
00087 !
00088      &(MASBED,EPAI,CONC,HDEP,TRA02,NPOIN2,
00089      & NPFMAX,NCOUCH,NPF,TASSE,GIBSON,RHOS,VOLU2D,
00090      & CFDEP,EPAICO,EPAINCO,MIXTE)
00091 !
00092 !***********************************************************************
00093 ! TELEMAC3D   V7P0                                   21/08/2010
00094 !***********************************************************************
00095 !
00096 !
00097 !
00098 !
00099 !
00100 !
00101 !
00102 !
00103 !
00104 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00105 !| CFDEP          |-->| CONCENTRATION OF MUD DEPOSIT (G/L)
00106 !| CONC           |-->| CONCENTRATION OF MUD BED LAYER
00107 !|                |   | (MULTILAYER MODEL)
00108 !| EPAI           |-->| THICKNESS OF SOLID BED LAYER
00109 !|                |   | (EPAI=DZ/(1+IVIDE), DZ TOTAL BED THICKNESS)
00110 !| EPAICO         |-->| THICKNESS OF COHESIVE SUB-LAYER
00111 !| EPAINCO        |-->| THICKNESS OF NON-COHESIVE SUB-LAYER
00112 !| GIBSON         |-->| LOGICAL FOR GIBSON MODEL
00113 !| HDEP           |-->| THICKNESS OF FRESH DEPOSIT (FLUID MUD LAYER)
00114 !| MASBED         |<->| MASS OF SEDIMENT BED
00115 !| MIXTE          |-->| LOGICAL, MIXED SEDIMENTS OR NOT
00116 !| NCOUCH         |-->| NUMBER OF LAYERS WITHIN THE BED
00117 !|                |   | (GIBSON MODEL)
00118 !| NPF            |-->| NUMBER OF POINTS OF THE BOTTOM ON ONE VERTICAL
00119 !| NPFMAX         |-->| MAXIMUM NUMBER OF HORIZONTAL PLANES THAT
00120 !|                |   | DISCRETISE MUD BED
00121 !| NPOIN2         |-->| NUMBER OF POINTS IN 2D
00122 !| RHOS           |-->| DENSITY OF SEDIMENT
00123 !| TASSE          |-->| MULTILAYER SETTLING MODEL LOGICAL
00124 !| TRA02          |<->| WORK ARRAY
00125 !| VOLU2D         |-->| INTEGRAL OF TEST FUNCTIONS IN 2D
00126 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00127 !
00128       USE BIEF
00129       USE INTERFACE_TELEMAC3D, EX_MASSED => MASSED
00130       IMPLICIT NONE
00131       INTEGER LNG,LU
00132       COMMON/INFO/LNG,LU
00133 !
00134 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00135 !
00136       INTEGER, INTENT(IN)             :: NPFMAX,NPOIN2,NCOUCH
00137       DOUBLE PRECISION, INTENT(IN)    :: CFDEP
00138       INTEGER, INTENT(IN)             :: NPF(NPOIN2)
00139       DOUBLE PRECISION, INTENT(INOUT) :: MASBED
00140       DOUBLE PRECISION, INTENT(IN)    :: EPAI(NPOIN2,NCOUCH)
00141       DOUBLE PRECISION, INTENT(IN)    :: EPAICO(NPOIN2), EPAINCO(NPOIN2)
00142       DOUBLE PRECISION, INTENT(IN)    :: VOLU2D(NPOIN2)
00143       DOUBLE PRECISION, INTENT(IN)    :: HDEP(NPOIN2)
00144       DOUBLE PRECISION, INTENT(IN)    :: CONC(NPOIN2,NCOUCH)
00145       DOUBLE PRECISION, INTENT(INOUT) :: TRA02(NPOIN2)
00146       DOUBLE PRECISION, INTENT(IN)    :: RHOS
00147       LOGICAL, INTENT(IN)             :: TASSE, GIBSON, MIXTE
00148 !
00149 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00150 !
00151       INTEGER IPOIN,IPF
00152       DOUBLE PRECISION MASSE6
00153       DOUBLE PRECISION P_DSUM
00154       EXTERNAL         P_DSUM
00155 !
00156 !=======================================================================
00157 !
00158 ! MASS OF MUDDY DEPOSITS ON THE RIGID BED (MASSE6)
00159 !
00160 !=======================================================================
00161 !
00162       IF(MIXTE) THEN
00163 
00164         DO IPOIN=1,NPOIN2
00165           TRA02(IPOIN)=0.D0
00166           TRA02(IPOIN)=TRA02(IPOIN) + CONC(IPOIN,1)*EPAICO(IPOIN)
00167      &      + CFDEP*EPAINCO(IPOIN)
00168         ENDDO
00169 
00170       ELSE
00171 
00172         DO IPOIN=1,NPOIN2
00173           TRA02(IPOIN)=0.D0
00174           DO IPF=1,NCOUCH
00175             TRA02(IPOIN)=TRA02(IPOIN)+CONC(IPOIN,IPF)*EPAI(IPOIN,IPF)
00176           ENDDO
00177         ENDDO
00178 !
00179       ENDIF
00180 
00181       MASSE6=0.D0
00182       DO IPOIN=1,NPOIN2
00183         MASSE6=MASSE6+VOLU2D(IPOIN)*TRA02(IPOIN)
00184       ENDDO
00185       IF(NCSIZE.GT.1) MASSE6=P_DSUM(MASSE6)
00186       MASBED = MASSE6
00187 !
00188 !-----------------------------------------------------------------------
00189 !
00190       RETURN
00191       END

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