fonvas.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac3d\fonvas.f
00002 !
00093                      SUBROUTINE FONVAS
00094 !                    *****************
00095 !
00096      &(IVIDE   , EPAI   , CONC   , TREST , TEMP   , HDEP   ,
00097      & FLUDP   , FLUDPT , FLUER  , ZF    , TA     , WC     ,
00098      & TRA01   , TRA02  , TRA03  , NPOIN2, NPOIN3 , NPFMAX ,
00099      & NCOUCH  , NPF    , LT     , DT    , DTC    , GRAV   ,
00100      & RHOS    , CFMAX  , TASSE  , ITASS , ZF_S   , ESOMT  ,
00101      & VOLU2D  , MASDEP , SETDEP , ZR    , TS     , FLUDPTC,
00102      & FLUDPTNC, FLUERC , FLUERNC, MIXTE , FLUDPC , FLUDPNC,
00103      & PVSCO   , PVSNCO , CFDEP  , EPAICO, EPAINCO)
00104 !
00105 !***********************************************************************
00106 ! TELEMAC3D   V7P0                                   21/08/2010
00107 !***********************************************************************
00108 !
00109 !
00110 !
00111 !
00112 !
00113 !
00114 !
00115 !
00116 !
00117 !
00118 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00119 !| CFDEP          |-->| CONCENTRATION OF MUD DEPOSIT (G/L)
00120 !| CFMAX          |<->| CONCENTRATION OF CONSOLIDATED MUD (G/L)
00121 !| CONC           |-->| MUD BED LAYER CONCENTRATION
00122 !|                |   | (MULTILAYER MODEL)
00123 !| DT             |-->| HYDRAULIC TIME STEP
00124 !| DTC            |-->| TIME STEP FOR CONSOLIDATION PHENOMENON
00125 !| EPAI           |<->| THICKNESS OF SOLID FRACTION OF THE BED LAYER
00126 !|                |   | (EPAI=DZ/(1+IVIDE), DZ BED LAYER THICKNESS
00127 !| EPAICO         |<->| THICKNESS OF COHESIVE SUB-LAYER
00128 !| EPAINCO        |<->| THICKNESS OF NON-COHESIVE SUB-LAYER
00129 !| ESOMT          |<->| CUMULATED BED EVOLUTION
00130 !| FLUDP          |<->| DEPOSITION FLUX
00131 !| FLUDPC         |<->| DEPOSITION FLUX FOR COHESIVE SEDIMENT IN 2D
00132 !| FLUDPNC        |<->| DEPOSITION FLUX FOR NON-COHESIVE SEDIMENT IN 2D
00133 !| FLUDPT         |<--| IMPLICIT DEPOSITION FLUX
00134 !| FLUDPTC        |<--| IMPLICIT DEPOSITION FLUX FOR COHESIVE SEDIMENT
00135 !| FLUDPTNC       |<--| IMPLICIT DEPOSITION FLUX FOR NON-COHESIVE SEDIMENT
00136 !| FLUER          |<--| EROSION FLUX FOR POINTS IN 2D
00137 !| FLUERC         |<--| EROSION FLUX FOR COHESIVE SEDIMENT IN 2D
00138 !| FLUERNC        |<--| EROSION FLUX FOR NON-COHESIVE SEDIMENT IN 2D
00139 !| GRAV           |-->| GRAVITY ACCELERATION
00140 !| HDEP           |<->| THICKNESS OF FRESH DEPOSIT (FLUID MUD LAYER)
00141 !| IVIDE          |<->| VOID RATIO
00142 !|                |   | (GIBSON MODEL ONLY)
00143 !| ITASS          |-->| INDEX OF MODEL CHOICE
00144 !| LT             |-->| CURRENT TIME STEP NUMBER
00145 !| MASDEP         |<->| DEPOSITED MASS
00146 !| MIXTE          |-->| LOGICAL, MIXED SEDIMENTS OR NOT
00147 !| NCOUCH         |-->| NUMBER OF LAYERS WITHIN THE BED
00148 !|                |   | (GIBSON MULTILAYER SETTLING MODEL)
00149 !| NPF            |-->| NUMBER OF POINTS OF THE BOTTOM ON ONE VERTICAL
00150 !| NPFMAX         |-->| MAXIMUM NUMBER OF HORIZONTAL PLANES
00151 !|                |   | DISCRETIZATION OF MUD BED (GIBSON MODEL)
00152 !| NPOIN2         |-->| NUMBER OF POINTS  (2D MESH)
00153 !| NPOIN3         |-->| NUMBER OF POINTS  (3D MESH)
00154 !| PVSCO          |<->| PERCENTAGE OF MUD
00155 !| PVSNCO         |<->| PERCENTAGE OF SAND
00156 !| RHOS           |-->| SEDIMENT DENSITY
00157 !| SETDEP         |-->| CHOICE OF ADVECTION SCHEME FOR VERTICAL SETTLING
00158 !| TA             |-->| ACTIVE TRACOR
00159 !| TASSE          |-->| MULTILAYER SETTLING MODEL LOGICAL
00160 !| TEMP           |<->| TIME COUNTER FOR CONSOLIDATION MODEL
00161 !|                |   | (MULTILAYER MODEL)
00162 !| TS             |-->| SAND CONCENTRATION
00163 !| TRA01          |<->| WORK ARRAY
00164 !| TRA02          |<->| WORK ARRAY
00165 !| TRA03          |<->| WORK ARRAY
00166 !| TREST          |<->| CONSOLIDATION TIME SCALE
00167 !|                |   | (ONLY FOR MULTILAYER MODEL)
00168 !| VOLU2D         |-->|  INTEGRAL OF TEST FUNCTIONS IN 2D (SURFACE OF ELEMENTS)
00169 !| WC             |-->| SETTLING VELOCITY
00170 !| ZF             |<->| BOTTOM ELEVATION
00171 !| ZF_S           |<->| BED EVOLUTION
00172 !| ZR             |-->| RIGID BED LEVEL
00173 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00174 !
00175       USE BIEF
00176       USE INTERFACE_TELEMAC3D, EX_FONVAS => FONVAS
00177       USE DECLARATIONS_TELEMAC3D, ONLY : IPBOT,OPTBAN,NPLAN,CGEL
00178 !
00179       IMPLICIT NONE
00180 !
00181       INTEGER LNG,LU
00182       COMMON/INFO/LNG,LU
00183 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00184 !
00185       INTEGER, INTENT(IN) ::  LT,NPOIN2,NPOIN3
00186       INTEGER, INTENT(IN) ::  NCOUCH,NPFMAX
00187 !
00188       DOUBLE PRECISION, INTENT(INOUT) :: IVIDE(NPOIN2,NCOUCH+1)
00189       DOUBLE PRECISION, INTENT(INOUT) :: TREST(NCOUCH)
00190       DOUBLE PRECISION, INTENT(INOUT) :: TEMP(NCOUCH,NPOIN2)
00191       DOUBLE PRECISION, INTENT(INOUT) :: HDEP(NPOIN2)
00192 !
00193       DOUBLE PRECISION, INTENT(INOUT) :: EPAI(NPOIN2,NCOUCH)
00194       DOUBLE PRECISION, INTENT(INOUT) :: EPAICO(NPOIN2), EPAINCO(NPOIN2)
00195 !
00196       DOUBLE PRECISION, INTENT(IN)    :: CONC(NPOIN2,NCOUCH),ZR(NPOIN2)
00197       DOUBLE PRECISION, INTENT(INOUT) :: ZF(NPOIN2)
00198       DOUBLE PRECISION, INTENT(IN)    :: TA(NPOIN3),WC(NPOIN3)
00199       DOUBLE PRECISION, INTENT(IN)    :: TS(NPOIN3)
00200       DOUBLE PRECISION, INTENT(INOUT) :: TRA01(NPFMAX,6)
00201       DOUBLE PRECISION, INTENT(INOUT) :: TRA02(NPFMAX),TRA03(NPFMAX)
00202       DOUBLE PRECISION, INTENT(IN)    :: FLUDPT(NPOIN2),FLUER(NPOIN2)
00203       DOUBLE PRECISION, INTENT(IN)    :: FLUDPTC(NPOIN2),FLUERC(NPOIN2)
00204       DOUBLE PRECISION, INTENT(IN)    :: FLUDPTNC(NPOIN2)
00205       DOUBLE PRECISION, INTENT(IN)    :: FLUERNC(NPOIN2)
00206       DOUBLE PRECISION, INTENT(INOUT) :: FLUDP(NPOIN2),ZF_S(NPOIN2)
00207       DOUBLE PRECISION, INTENT(INOUT) :: FLUDPC(NPOIN2), FLUDPNC(NPOIN2)
00208       DOUBLE PRECISION, INTENT(INOUT) :: PVSCO(NPOIN2), PVSNCO(NPOIN2)
00209       DOUBLE PRECISION, INTENT(INOUT) :: ESOMT(NPOIN2)
00210 !
00211       DOUBLE PRECISION, INTENT(IN)    :: DT,RHOS,GRAV,DTC, CFDEP
00212       DOUBLE PRECISION, INTENT(INOUT) :: CFMAX,MASDEP
00213 !
00214       INTEGER, INTENT(INOUT) ::  NPF(NPOIN2)
00215 !
00216       LOGICAL, INTENT(IN) :: TASSE, MIXTE
00217       INTEGER, INTENT(IN) :: SETDEP
00218       INTEGER, INTENT(IN) :: ITASS
00219 !
00220       TYPE(BIEF_OBJ), INTENT(IN) :: VOLU2D
00221 !
00222 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00223 !
00224       DOUBLE PRECISION C, TOTMASS, QERODE, QS, DELTAF, FLUX, SEDBED
00225       DOUBLE PRECISION MASTMP
00226       DOUBLE PRECISION TOTMASSC, QERODEC, QSC, DELTAFC, FLUXC
00227       DOUBLE PRECISION TOTMASSNC, QERODENC, QSNC, DELTAFNC, FLUXNC
00228       INTEGER IPOIN, IC
00229       DOUBLE PRECISION P_DSUM
00230       EXTERNAL         P_DSUM
00231 !
00232 !=======================================================================
00233 ! FIRST STEP
00234 !     COMPUTES THE DEPOSITED QUANTITY (IN MATERIAL COORDINATES)
00235 !     CFDEP=CONC(NCOUCH)
00236 !     COMPUTES THE ERODED QUANTITY
00237 !=======================================================================
00238 !
00239 ! Multi layer model
00240 ! +++++++++++++++++++
00241 ! deposition in the first top layer
00242 ! calculate the layers thicknesses and deposited thicknes:  HDEP = sum ( EPAI)
00243 !
00244       FLUX=0.D0
00245 
00246       IF(MIXTE) THEN
00247 
00248         IF(SETDEP.NE.1) THEN
00249           IF(OPTBAN.EQ.1) THEN
00250             DO IPOIN=1,NPOIN2
00251               IF(IPBOT%I(IPOIN).NE.NPLAN-1) THEN
00252 !         COMPUTES FIRST THE DEPOSIT FLUX OF COHESIVE SEDIMENTS
00253                 FLUDPC(IPOIN)=FLUDPTC(IPOIN)*
00254      &          TA(IPBOT%I(IPOIN)*NPOIN2+IPOIN)
00255                 FLUDPC(IPOIN)=MAX(FLUDPC(IPOIN),0.D0)
00256 !         THEN COMPUTES THE DEPOSIT FLUX OF NON COHESIVE SEDIMENTS
00257                 FLUDPNC(IPOIN)=FLUDPTNC(IPOIN)*
00258      &          TS(IPBOT%I(IPOIN)*NPOIN2+IPOIN)
00259                 FLUDPNC(IPOIN)=MAX(FLUDPNC(IPOIN),0.D0)
00260 !         THE GLOBAL DEPOSITION FLUX IS THE SUMM OF BOTH C & NC
00261                 FLUDP(IPOIN)=FLUDPC(IPOIN)+FLUDPNC(IPOIN)
00262                 FLUDP(IPOIN)=MAX(FLUDP(IPOIN),0.D0)
00263               ELSE
00264                 FLUDPC(IPOIN) = 0.D0
00265                 FLUDPNC(IPOIN)= 0.D0
00266                 FLUDP(IPOIN)  = 0.D0
00267               ENDIF
00268             ENDDO
00269           ELSE
00270             DO IPOIN=1,NPOIN2
00271 !         COMPUTES FIRST THE DEPOSIT FLUX OF COHESIVE SEDIMENTS
00272               FLUDPC(IPOIN)=FLUDPTC(IPOIN)*TA(IPOIN)
00273               FLUDPC(IPOIN)=MAX(FLUDPC(IPOIN),0.D0)
00274 !         THEN COMPUTES THE DEPOSIT FLUX OF NON COHESIVE SEDIMENTS
00275               FLUDPNC(IPOIN)=FLUDPTNC(IPOIN)*TS(IPOIN)
00276               FLUDPNC(IPOIN)=MAX(FLUDPNC(IPOIN),0.D0)
00277 !         THE GLOBAL DEPOSITION FLUX IS THE SUMM OF BOTH C & NC
00278               FLUDP(IPOIN)=FLUDPC(IPOIN)+FLUDPNC(IPOIN)
00279               FLUDP(IPOIN)=MAX(FLUDP(IPOIN),0.D0)
00280             ENDDO
00281           ENDIF
00282         ENDIF
00283 !
00284 !         BED EVOLUTION
00285 !
00286         DO IPOIN=1,NPOIN2
00287 
00288 !         COMPUTES QERODE FOR COHESIVE SEDIMENTS
00289 
00290           DELTAFC  = FLUDPC(IPOIN)-FLUERC(IPOIN)
00291           FLUXC    = FLUXC+DELTAFC*VOLU2D%R(IPOIN)
00292           TOTMASSC = 0.D0
00293           QERODEC  = FLUERC(IPOIN)*DT
00294 
00295 !         COMPUTES QERODE FOR NON-COHESIVE SEDIMENT
00296 
00297           DELTAFNC  = FLUDPNC(IPOIN)-FLUERNC(IPOIN)
00298           FLUXNC    = FLUXNC+DELTAFNC*VOLU2D%R(IPOIN)
00299           TOTMASSNC = 0.D0
00300           QERODENC  = FLUERNC(IPOIN)*DT
00301 
00302 !         COMPUTES QERODE FOR ALL THE SEDIMENTS
00303 
00304           DELTAF  = FLUDP(IPOIN)-FLUER(IPOIN)
00305           FLUX    = FLUX+DELTAF*VOLU2D%R(IPOIN)
00306           TOTMASS = 0.D0
00307           QERODE  = FLUER(IPOIN)*DT
00308 !
00309 
00310           QSC  = CONC(IPOIN,1)*EPAICO(IPOIN)
00311           QSNC = CFDEP*EPAINCO(IPOIN)
00312           QS   = QSC + QSNC
00313 !
00314           TOTMASSC  = TOTMASSC + QSC
00315           TOTMASSNC = TOTMASSNC + QSNC
00316           TOTMASS   = TOTMASS + QS
00317 
00318 !         check if we have eroded enough entire layers
00319 
00320           IF(TOTMASSC.LT.QERODEC) THEN
00321             EPAICO(IPOIN) = 0.D0
00322           ELSE
00323             QSC           = TOTMASSC - QERODEC
00324             EPAICO(IPOIN) = QSC/MAX(CONC(IPOIN,1),1.D-10)
00325           ENDIF
00326 !
00327           IF(TOTMASSNC.LT.QERODENC) THEN
00328             EPAINCO(IPOIN) = 0.D0
00329           ELSE
00330             QSNC           = TOTMASSNC - QERODENC
00331             EPAINCO(IPOIN) = QSNC/CFDEP
00332           ENDIF
00333 !
00334 !         DEPOSITION IN THE TOP LAYER
00335 !
00336         EPAICO(IPOIN)  = EPAICO(IPOIN)+FLUDPC(IPOIN)*
00337      &                DT/MAX(CONC(IPOIN,1),1.D-10)
00338         EPAINCO(IPOIN) = EPAINCO(IPOIN)+FLUDPNC(IPOIN)*
00339      &                DT/MAX(CFDEP,1.D-10)
00340         EPAI(IPOIN,1)  = EPAICO(IPOIN) + EPAINCO(IPOIN)
00341 !
00342 !         UPDATES PERCENTAGES OF EACH CLASSE
00343 !
00344         IF(EPAI(IPOIN, 1).GT.0.D0) THEN
00345           PVSCO(IPOIN)  = EPAICO(IPOIN)/EPAI(IPOIN,1)
00346           PVSNCO(IPOIN) = 1.D0-PVSCO(IPOIN)
00347         ELSE
00348           PVSCO(IPOIN)  = 0.D0
00349           PVSNCO(IPOIN) = 0.D0
00350         ENDIF
00351 !
00352 !         COMPUTING THE NEW SEDIMENT BED THICKNESS
00353 !
00354         SEDBED = EPAI(IPOIN,1)
00355 !
00356 !         EVOLUTION OBTAINED FROM OLD AND NEW SEDIMENT HEIGHT
00357 !
00358         ZF_S(IPOIN) = SEDBED-HDEP(IPOIN)
00359 !
00360 !         SEDIMENT HEIGHT UPDATED
00361 !
00362         HDEP(IPOIN) = SEDBED
00363 !
00364       ENDDO
00365 !
00366 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++
00367 !
00368 !     UPDATE THE CUMULATED BED EVOLUTION : ESOMT
00369 !     BOTTOM ELEVATION : ZF
00370 !
00371 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++
00372 !
00373         CALL OV( 'X=Y+Z   ',ESOMT,ESOMT,ZF_S,C,NPOIN2)
00374         CALL OV( 'X=Y+Z   ' , ZF   ,   ZR, HDEP, C, NPOIN2)
00375 
00376 
00377 !***********************************************************
00378 ! PREVIOUS CASE NON MIXED
00379 !************************************************************
00380       ELSE
00381 
00382 !     EXPLICIT SCHEME (SETDEP=1) FLUDP COMPUTED IN SET_DIF
00383 !
00384 !     OTHER SCHEMES : FLUDP BUILT HERE
00385 !
00386       IF(SETDEP.NE.1) THEN
00387         IF(OPTBAN.EQ.1) THEN
00388           DO IPOIN=1,NPOIN2
00389 !           correction for tidal flats: take the first point above crushed planes
00390 !           IPBOT =0  :  no tidal flats  IPBOT = NPLAN-1 : dry element
00391             IF(IPBOT%I(IPOIN).NE.NPLAN-1) THEN
00392               FLUDP(IPOIN)=FLUDPT(IPOIN)*
00393      &        TA(IPBOT%I(IPOIN)*NPOIN2+IPOIN)
00394               FLUDP(IPOIN)=MAX(FLUDP(IPOIN),0.D0)
00395             ELSE
00396               FLUDP(IPOIN)=0.D0
00397             ENDIF
00398           ENDDO
00399         ELSE
00400           DO IPOIN=1,NPOIN2
00401             FLUDP(IPOIN)=FLUDPT(IPOIN)*TA(IPOIN)
00402 !           FLUDP MUST BE POSITIVE, EVEN IF TA<0 DUE TO TRUNCATION ERRORS
00403 !           PROBLEM SEEN WITH TA=-1.D-87 !!!!!
00404             FLUDP(IPOIN)=MAX(FLUDP(IPOIN),0.D0)
00405           ENDDO
00406         ENDIF
00407       ENDIF
00408 !
00409 !     BED EVOLUTION
00410 !
00411       DO IPOIN=1,NPOIN2
00412 !
00413         DELTAF=FLUDP(IPOIN)-FLUER(IPOIN)
00414         FLUX=FLUX+DELTAF*VOLU2D%R(IPOIN)
00415         TOTMASS=0.D0
00416         QERODE = FLUER(IPOIN)*DT
00417 !
00418         DO IC=1,NCOUCH
00419 !
00420           QS = CONC(IPOIN,IC)*EPAI(IPOIN,IC)
00421 !
00422           TOTMASS = TOTMASS + QS
00423 !         check if we have eroded enough entire layers
00424           IF(TOTMASS.LT.QERODE) THEN
00425             EPAI(IPOIN,IC) = 0.D0
00426           ELSE
00427 !           we have got to the correct layer.
00428 !           How much of it do we need to erode?
00429             QS = TOTMASS - QERODE
00430 !           calculate new thickness
00431             EPAI(IPOIN,IC) = QS/MAX(CONC(IPOIN,IC),1.D-10)
00432 !           jump out of layer loop
00433             EXIT
00434           ENDIF
00435 !
00436         ENDDO
00437 !
00438 !       Then Deposition in Top layer
00439 !
00440         EPAI(IPOIN,1)=EPAI(IPOIN,1)+FLUDP(IPOIN)*
00441      &             DT/MAX(CONC(IPOIN,1),1.D-10)
00442 !
00443 !       COMPUTING THE NEW SEDIMENT BED THICKNESS
00444 !
00445         SEDBED= 0.D0
00446         DO IC=1,NCOUCH
00447           SEDBED=SEDBED+EPAI(IPOIN,IC)
00448         ENDDO
00449 !
00450 !       EVOLUTION OBTAINED FROM OLD AND NEW SEDIMENT HEIGHT
00451 !
00452         ZF_S(IPOIN)=SEDBED-HDEP(IPOIN)
00453 !
00454 !       SEDIMENT HEIGHT UPDATED
00455 !
00456         HDEP(IPOIN) = SEDBED
00457 !
00458       ENDDO
00459 !
00460 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++
00461 !
00462 !     UPDATE THE CUMULATED BED EVOLUTION : ESOMT
00463 !     BOTTOM ELEVATION : ZF
00464 !
00465 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++
00466 !
00467       CALL OV( 'X=Y+Z   ',ESOMT,ESOMT,ZF_S,C,NPOIN2)
00468 !
00469 !     NOTE JMH: THIS WAY OF WRITING THE NEW ZF ENSURES
00470 !               THAT ZF ABOVE ZR EVEN WITH TRUNCATION
00471 !               ERRORS, IF HDEP >0, THIS IS IMPORTANT
00472 !
00473 !     CALL OV( 'X=Y+Z   ' , ZF   ,   ZF, ZF_S, C, NPOIN2)
00474       CALL OV( 'X=Y+Z   ' , ZF   ,   ZR, HDEP, C, NPOIN2)
00475 !
00476       ENDIF
00477 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++
00478 ! TASSEMENT HERE
00479 ! ---> add the other models here
00480 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++
00481 !
00482 !      IF (TASS) THEN
00483 !
00484 !     -----CHECKS THE TIMESTEP TO TAKE CONSOLIDATION-----
00485 !     -----              INTO ACCOUNT               -----
00486 !
00487 !         IF(ITASS.EQ.2)
00488 !
00489 !         IF (MOD(LT*DT,DTC).LT.1.D-8) THEN
00490 !
00491 !     -----MANAGES THE DEPOSITED QUANTITY : CREATES-----
00492 !     -----     NEW LAYERS NEAR THE MUDDY BED      -----
00493 !
00494 !            CALL GESTDP( IVIDE  , EPAI   , HDEP    ,
00495 !     &                  NPOIN2 , NPFMAX , NPF     ,
00496 !     &                  EPAI0  , CONC(1)  , RHOS    )
00497 !
00498 !     -----CONSOLIDATES THE MUDDY BED-----
00499 !     -----    (GIBSON EQUATION)     -----
00500 !
00501 !           CALL TASSEM( IVIDE , EPAI ,
00502 !     &                  NPOIN2, NPFMAX, NPF  ,
00503 !     &                  GRAV  , RHOS  , DTC  , CFMAX ,
00504 !     &                  TRA01 , TRA02 ,TRA03 )
00505 !
00506 !         ENDIF
00507 !
00508 !     -----UPDATES THE BOTTOM ELEVATION-----
00509 !
00510 !         CALL ACTUZF(IVIDE,EPAI,ZF,NPOIN2,NPFMAX,NPF)
00511 !
00512 ! Will have to be rewritten
00513 !          CALL GIBSON(ZF,NPOIN2,DT,ELAY,
00514 !     &               T3,T2,LT,XMVS,XMVE,GRAV,NOMBLAY,
00515 !     &               ES,CONC_VASE,CONC,IVIDE,MS_VASE%R,XWC(1),
00516 !     &                 TRA01,TRA02,TRA03,CGEL,COEF_N,CFMAX)
00517 !      ELSEIF(ITASS.EQ.1) THEN
00518 !
00519 !     -----MODELS CONSOLIDATION (SIMPLE)-----
00520 !
00521 !         IF (MOD(LT*DT,DTC).LT.1.D-8)
00522 !     &    CALL TASSEC( CONC   , EPAI , TREST , TEMP , DTC ,
00523 !     &                 NPOIN2 , NCOUCH )
00524 !
00525 !
00526 !
00527 !     -----UPDATES THE BOTTOM ELEVATION-----
00528 !
00529 !          DO IPOIN = 1 , NPOIN2
00530 !            HDEP(IPOIN)=0.D0
00531 !            DO IC = 1 , NCOUCH
00532 !              HDEP(IPOIN)= HDEP(IPOIN) + EPAI(IPOIN,IC)
00533 !            ENDDO
00534 !            ZF(IPOIN) = ZR(IPOIN) + HDEP(IPOIN)
00535 !          ENDDO
00536 !
00537 !      ENDIF
00538 !      ENDIF
00539 !
00540 !!+++++++++++++++++++++++++++++++++++++++++++++++++++++++
00541 ! COMPUTES HERE DEPOSITED MASS
00542 !
00543 !     TOTAL DEPOSITED MASS --> MASDEP
00544 !
00545 !###>TBE - fixed bug...add up the deposition flux for all partitions (MASDEP)
00546 !       MASDEP = MASDEP + FLUX*DT
00547         MASTMP = FLUX*DT
00548         IF(NCSIZE.GT.1) MASTMP=P_DSUM(MASTMP)
00549         MASDEP = MASDEP + MASTMP
00550 !###<TBE
00551 !
00552 !=======================================================================
00553 !
00554       RETURN
00555       END

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