!                    ***************************
                     SUBROUTINE USER_CONDIN_TRAC
!                    ***************************
!
     &(CS,CS0,NSICLA)
!
!***********************************************************************
! TELEMAC2D   V7P3
!***********************************************************************
!
!brief    USER INITIALISES THE PHYSICAL PARAMETERS TRAC
!
!history  J-M HERVOUET (LNHE)
!+        30/08/2007
!+        V6P0
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE BIEF
      USE DECLARATIONS_TELEMAC
      USE DECLARATIONS_TELEMAC2D
      USE TPXO
      USE OKADA
!
      USE DECLARATIONS_SPECIAL
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN)           :: NSICLA
      DOUBLE PRECISION,INTENT(IN)   :: CS0(1)
      TYPE(BIEF_OBJ), INTENT(INOUT) :: CS
      INTEGER ITRAC,I
!
!-----------------------------------------------------------------------
!
! MODIFIED BY QINGHUI ZHANG 15 JULY 2013
! PURPOSE IS TO PRESCRIBE A SPACE-DEPENDENT INITIAL TRACER VALUES
! THE TRACER VALUES FOR LEFT BOUNDARY NODES ARE 30
! CENCENTRATION FOR THE REST PART OF LEFT BOUNDARY ELEMENTS ARE INTERPOLATED
! FOR REST PART OF THE DOMAIN, CONCENTRATION IS 0
! THE MESH ELEMENT LENGTH ALONG CHANNEL BANK IS 40 METERS
! SO HERE IF THE X COORDINATE OF NODE LESS THAN 5 METERS, INITIAL VALUE: 30
		WRITE(*,*)NTRAC
		READ(*,*)
      IF(NTRAC.GT.0) THEN
        DO ITRAC=1,NTRAC
          CALL OS('X=0     ', X=CS%ADR(ITRAC)%P)
          DO I=1,NPOIN
!            IF(X(I).LT.5.D0) THEN
              CS%ADR(ITRAC)%P%R(I) = 0.002
!            ENDIF
          ENDDO
        ENDDO
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END