The TELEMAC-MASCARET system  trunk
depos_fx.f
Go to the documentation of this file.
1 ! ***********************
2  SUBROUTINE depos_fx
3 ! **********************
4 !
5  &(sedp,taub,csus,taus,vitchu,npoin)
6 !
7 !***********************************************************************
8 ! TELEMAC2D V7P2
9 !***********************************************************************
10 !
11 !brief COMPUTES DEPOSITION FLUX
12 !
13 !
14 !
15 !history R. ATA (LNHE)
16 !+ 02/09/2015
17 !+ V7P1
18 !+
19 !
20 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
21 !| NPOIN |-->| TOTAL NUMBER OF MESH NODES
22 !| SEDP |<--| DEPOSITION FLUX
23 !| TAUB |-->| BED SHEAR STRESS
24 !| TAUS |-->| SEDIMENTATION CRITICAL STRESS
25 !| CSUS |-->| SUSPENDED LOAD (TRACER 1)
26 !| VITCHU |-->| SEDIMENT SETTLING VELOCITY
27 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
28 !
29  USE bief
31  USE interface_waqtel, ex_depos_fx => depos_fx
32 !
33  IMPLICIT NONE
34 !
35 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
36 !
37  INTEGER , INTENT(IN) :: npoin
38  DOUBLE PRECISION, INTENT(IN) :: taus,vitchu
39  TYPE(bief_obj) , INTENT(IN ) :: taub,csus
40  TYPE(bief_obj) , INTENT(INOUT) :: sedp
41  INTRINSIC max
42 !
43 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
44 !
45 !
46 ! LOCAL VARIABLES
47  INTEGER i
48 !
49 !
50  IF (abs(taus).LT.1.e-10)THEN
51  WRITE(lu,*) 'DEPOS_FX: CRITICAL STRESS OF RESUSPENSION '
52  WRITE(lu,*) ' TAUS VERY SMALL OR NIL - VERIFY !!!'
53  WRITE(lu,*) ' TAUS = ',taus
54  CALL plante(1)
55  stop
56  ENDIF
57 !
58 ! THIS WAY WORKS WELL FOR 2D AND FOR 3D AS WELL SINCE BED LAYER IS
59 ! FOR I=1 TO NPOIN
60  DO i=1,npoin
61  sedp%R(i)=vitchu*csus%R(i)*max(1.d0-taub%R(i)/taus,0.d0)
62  ENDDO
63 !
64  RETURN
65  END
66 !
67 !-----------------------------------------------------------------------
68 !
subroutine vitchu(S3D_WCHU, S3D_WCHU0, U, V, S3D_TURBA, S3D_TURBB, TRAV1, TRAV2, TRAV3, S, MESH3D, IELM3, NPOIN2, NPOIN3, NPLAN, NTRAC, MSK, MASKEL, UETCAR, TA, HN, S3D_FLOC, S3D_FLOC_TYPE, S3D_HINDER, S3D_HIND_TYPE, S3D_CGEL, S3D_CINI)
Definition: vitchu.f:10
Definition: bief.f:3