cgl.f

Go to the documentation of this file.

!                   *****************************
                    DOUBLE PRECISION FUNCTION cgl
!                   *****************************
!
     &(i,at)
!
!***********************************************************************
! TELEMAC2D   V7P1
!***********************************************************************
!
!brief    PRESCRIBES THE FREE SURFACE ELEVATION FOR LEVEL IMPOSED
!+                LIQUID BOUNDARIES.
!
!history  J-M HERVOUET (LNHE)
!+        17/08/1994
!+        V6P0
!+
!
!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
!
!history  C. COULET (ARTELIA GROUP)
!+        08/11/2011
!+        V6P2
!+   Modification size FCT due to modification of TRACER numbering
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| I              |-->| NUMBER OF LIQUID BOUNDARY
!| N              |-->| GLOBAL NUMBER OF POINT
!|                |   | IN PARALLEL NUMBER IN THE ORIGINAL MESH
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
      USE declarations_sisyphe
      USE declarations_special
!
!      USE INTERFACE_SISYPHE, EX_CGL => CGL
!
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN) :: I
      DOUBLE PRECISION, INTENT(IN):: AT
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      CHARACTER(LEN=9) FCT
!
!-----------------------------------------------------------------------
!
!     IF THE LIQUID BOUNDARY FILE EXISTS, ATTEMPTS TO FIND
!     THE VALUE IN IT. IF YES, OK REMAINS TO .TRUE. FOR NEXT CALLS
!                      IF  NO, OK IS SET  TO .FALSE.
!
      IF(okcgl(i).AND.sis_files(sisliq)%NAME(1:1).NE.' ') THEN
!
!       FCT WILL BE CGL(1), CGL(2), ETC, CGL(9), DEPENDING ON I
        fct='CG(      '
        IF(i.LT.10) THEN
          WRITE(fct(4:4),fmt='(I1)') i
          fct(5:5)=')'
        ELSEIF(i.LT.100) THEN
          WRITE(fct(4:5),fmt='(I2)') i
          fct(6:6)=')'
        ELSE
          WRITE(lu,*)'I=',i
          WRITE(lu,*) 'CGL NOT PROGRAMMED FOR MORE THAN 99 BOUNDARIES'
          CALL plante(1)
          stop
        ENDIF
!
        CALL read_fic_conc(cgl,fct,at,sis_files(sisliq)%LU,
     &                     okcgl(i))
!
      ENDIF
!
      IF(.NOT.okcgl(i).OR.sis_files(sisliq)%NAME(1:1).EQ.' ') THEN
!       USER DEFINED
        CALL user_cgl(cgl, i, at)
        WRITE(lu,101) i
101     FORMAT(1x,/,1x,'CG: MORE PRESCRIBED ELEVATIONS ARE REQUIRED'
     &           ,/,1x,'     IN THE PARAMETER FILE'
     &           ,/,1x,'     AT LEAST ',1i6,' MUST BE GIVEN')
        CALL plante(1)
        stop
      ENDIF
!
!
!-----------------------------------------------------------------------
!
      RETURN
      END