The TELEMAC-MASCARET system  trunk
eqcae_bc_gaia.f
Go to the documentation of this file.
1 ! ************************
2  SUBROUTINE eqcae_bc_gaia
3 ! ************************
4 !
5  &(litbor,tbor,tn,j,kent)
6 !
7 !***********************************************************************
8 ! GAIA
9 !***********************************************************************
10 !
12 !
13 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
20 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
21 !
22  USE bief
25  IMPLICIT NONE
26 !
27 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
28 !
29  INTEGER, INTENT(IN) :: LITBOR(nptfr),J,KENT
30  DOUBLE PRECISION, INTENT(INOUT) :: TBOR(nptfr),TN(npoin)
31 !
32 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
33 !
34  INTEGER I,K
35 !
36 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
37 !
38  DO k = 1, nptfr
39  IF(litbor(k).EQ.kent) THEN
40  i = mesh%NBOR%I(k)
41  IF(.NOT.sedco(j)) THEN
42  tbor(k) = cstaeq%ADR(j)%P%R(i)/csratio%R(i)
43  ELSE
44  tbor(k) = fluer%ADR(j)%P%R(i)/xwc(j)
45  ENDIF
46 ! THIS IS THE CONDITION TO HAVE NO EVOLUTION
47 ! CS%R(I) MAY BE DIFFERENT FROM CBOR%R(K) IF UNSTEADY FLOW
48 ! OR IF DIRFLU.EQ.2 (CASE OF PRIORITY TO FLUXES)
49  fluer%ADR(j)%P%R(i) = fludpt%ADR(j)%P%R(i)*tn(i)
50  ENDIF
51  ENDDO
52 !
53  RETURN
54  END
integer, pointer nptfr
Number of boundary points.
type(bief_obj), target fludpt
Deposition flux for implicitation.
double precision, dimension(nsiclm), target xwc
Settling velocities.
subroutine eqcae_bc_gaia(LITBOR, TBOR, TN, J, KENT)
Definition: eqcae_bc_gaia.f:7
logical, dimension(nsiclm) sedco
Cohesive sediments (for each class)
type(bief_obj), target csratio
Ratio between bottom concentration and average concentration.
type(bief_obj), target cstaeq
Sediment equilibrium concentration.
type(bief_mesh), target mesh
Mesh structure.
type(bief_obj), target fluer
Erosion flux.
integer, pointer npoin
Number of 2d points in the mesh.
Definition: bief.f:3