cordep.f

Go to the documentation of this file.

!                       *****************
                        SUBROUTINE cordep
!                       *****************
!
     &(ikle,lgvec)
!
!***********************************************************************
! PROGICIEL: STBTEL V5.2          28/08/89  J-M HERVOUET (LNH) 3071 8018
!***********************************************************************
!
! FONCTION : CORRECTION DES DEPENDANCES ARRIERES
!
!-----------------------------------------------------------------------
!                             ARGUMENTS
! .________________.____.______________________________________________.
! |      NOM       |MODE|                   ROLE                       |
! |________________|____|______________________________________________|
! |  IKLE          |<-->| NUMERO GLOBAUX DES NOEUDS DE CHAQUE ELEMENT  |
! |________________|____|______________________________________________
! | COMMON:        |    |
! |  GEO:          |    |
! |    MESH        | -->| TYPE DES ELEMENTS DU MAILLAGE
! |    NDP         | -->| NOMBRE DE NOEUDS PAR ELEMENTS
! |    NPOIN       | -->| NOMBRE TOTAL DE NOEUDS DU MAILLAGE
! |    NELEM       | -->| NOMBRE TOTAL D'ELEMTS DU MAILLAGE
! |    NPMAX       | -->| DIMENSION EFFECTIVE DES TABLEAUX X ET Y
! |                |    | (NPMAX = NPOIN + 0.1*NELEM)
! |    NELMAX      | -->| DIMENSION EFFECTIVE DES TABLEAUX CONCERNANT
! |                |    | LES ELEMENTS (NELMAX = NELEM + 0.2*NELEM)
! |________________|____|______________________________________________|
! MODE : -->(DONNEE NON MODIFIEE), <--(RESULTAT), <-->(DONNEE MODIFIEE)
!-----------------------------------------------------------------------
! APPELE PAR : STBTEL
! APPEL DE : ECHELE
!***********************************************************************
!
      USE declarations_special
      USE declarations_stbtel, ONLY: nelem,nelmax
      USE interface_stbtel, ex_cordep => cordep
      IMPLICIT NONE
!
      INTEGER, INTENT(INOUT) :: IKLE(nelmax,4)
      INTEGER, INTENT(IN) :: LGVEC
!
      INTEGER IELEM , IEL1 , IEL2
      INTEGER I1 , I2 , I3 , J1 , J2 , J3
      INTEGER K
!
      LOGICAL DEP
!
!
!=======================================================================
!
      DO ielem = 1,nelem
        iel2 = ielem
25      CONTINUE
        dep = .false.
        i1 = ikle(ielem,1)
        i2 = ikle(ielem,2)
        i3 = ikle(ielem,3)
        DO k = 2,lgvec
          iel1 = mod(nelem+ielem-k,nelem) + 1
          j1 = ikle(iel1,1)
          j2 = ikle(iel1,2)
          j3 = ikle(iel1,3)
          IF (i1.EQ.j1.OR.i1.EQ.j2.OR.i1.EQ.j3.OR.
     &        i2.EQ.j1.OR.i2.EQ.j2.OR.i2.EQ.j3.OR.
     &        i3.EQ.j1.OR.i3.EQ.j2.OR.i3.EQ.j3) dep=.true.
        ENDDO
        IF (dep) THEN
          iel2 = mod(iel2,nelem) + 1
          IF (iel2.EQ.ielem) GOTO 40
          CALL echele(ikle,ielem,iel2)
          GOTO 25
        ENDIF
      ENDDO ! IELEM
!
!=======================================================================
!
40    CONTINUE
!
      RETURN
      END