corrsl.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE corrsl
!                   *****************
!
     &(newsl,oldsl,zf,mesh)
!
!***********************************************************************
! BIEF   V7P0                                   21/08/2010
!***********************************************************************
!
!brief    CORRECTS THE FREE SURFACE COMPUTATION BY ELEMENTS
!+                TO TAKE ACCOUNT OF THE TIDAL FLATS.
!
!history  J-M JANIN (LNH)    ; J-M HERVOUET (LNH)
!+        27/11/92
!+        V5P1
!+
!
!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  J-M HERVOUET (EDF LAB, LNHE)
!+        12/05/2014
!+        V7P0
!+   Discontinuous elements better treated: new types 15, 16 and 17 for
!+   discontinuous linear, quasi-bubble, and quadratic, rather than
!+   using component DIMDISC=11, 12 or 13.
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| MESH           |-->| MESH STRUCTURE
!| NEWSL          |<->| MODIFIED FREE SURFACE, GIVEN PER ELEMENT
!| OLDSL          |-->| ORIGINAL FREE SURFACE
!| ZF             |-->| BATHYMETRY
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_corrsl => corrsl
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      TYPE(bief_obj) , INTENT(INOUT) :: NEWSL
      TYPE(bief_obj) , INTENT(IN)    :: OLDSL,ZF
      TYPE(bief_mesh), INTENT(IN)    :: MESH
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER NELEM,NELMAX,IELM
!
!-----------------------------------------------------------------------
!
      nelem = mesh%NELEM
      nelmax= mesh%NELMAX
!
      ielm=oldsl%ELM
!
!-----------------------------------------------------------------------
!
      IF(ielm.EQ.11) THEN
!
        CALL crsl11(newsl%R,oldsl%R,zf%R,mesh%IKLE%I,nelem,nelmax)
!
!-----------------------------------------------------------------------
!
      ELSEIF(ielm.EQ.12) THEN
!
        CALL crsl12(newsl%R,oldsl%R,zf%R,mesh%IKLE%I,nelem,nelmax)
!
!-----------------------------------------------------------------------
!
      ELSE
!
        WRITE(lu,11) ielm
11      FORMAT(1x,'CORRSL: UNKNOWN DISCRETIZATION :',i6)
        CALL plante(1)
        stop
!
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END