diffcl.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE diffcl
!                   *****************
!
     &(litbor,ttild,tbor,nbor,iconv,npoin,nptfr)
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    INITIALISES THE VALUE OF THE TRACER FOR BOUNDARY
!+                CONDITIONS OF TYPE DIRICHLET, IN THE DIFFUSION STEP.
!
!history  J-M HERVOUET (LNHE)
!+        09/10/2009
!+        V6P0
!+   MOVED FROM TELEMAC-2D TO ALLOW CALL BY SISYPHE
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        13/07/2010
!+        V6P0
!+   Translation of French comments within the FORTRAN sources into
!+   English comments
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        21/08/2010
!+        V6P0
!+   Creation of DOXYGEN tags for automated documentation and
!+   cross-referencing of the FORTRAN sources
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| ICONV          |-->| THE CONVENTION FOR THE ADVECTION SCHEME
!|                |   | E.G. 1 FOR CHARACTERISTICS
!| LITBOR         |-->| TYPES OF BOUNDARY CONDITIONS FOR THE TRACER
!| NBOR           |-->| GLOBAL NUMBERS OF BOUNDARY POINTS
!| NPOIN          |-->| NUMBER OF POINTS
!| NPTFR          |-->| NUMBER OF BOUNDARY POINTS
!| TBOR           |<->| PRESCRIBED BOUNDARY CONDITIONS ON T.
!| TTILD          |-->| TRACER AFTER ADVECTION.
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_telemac
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN)             :: NPOIN,NPTFR,ICONV
      INTEGER, INTENT(IN)             :: NBOR(nptfr)
      INTEGER, INTENT(IN)             :: LITBOR(nptfr)
      DOUBLE PRECISION, INTENT(IN)    :: TTILD(npoin)
      DOUBLE PRECISION, INTENT(INOUT) :: TBOR(nptfr)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER K
!
!----------------------------------------------------------------------
!
!  INITIALISES TBOR (SEE CONSTRUCTION OF LIMTRA IN DIFFIN)
!
      IF(iconv.EQ.adv_car) THEN
!
      DO k=1,nptfr
!
!  IMPOSES THE RESULT OF ADVECTION AT FREE EXITS
!                      WITH THE METHOD OF CHARACTERISTICS
!
        IF(litbor(k).EQ.ksort) tbor(k) = ttild(nbor(k))
!
      ENDDO ! K
!
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END