cordep.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\stbtel\cordep.f
00002 !
00031                         SUBROUTINE CORDEP
00032 !                       *****************
00033 !
00034      &(IKLE,LGVEC)
00035 !
00036 !***********************************************************************
00037 ! PROGICIEL: STBTEL V5.2          28/08/89  J-M HERVOUET (LNH) 3071 8018
00038 !***********************************************************************
00039 !
00040 ! FONCTION : CORRECTION DES DEPENDANCES ARRIERES
00041 !
00042 !-----------------------------------------------------------------------
00043 !                             ARGUMENTS
00044 ! .________________.____.______________________________________________.
00045 ! |      NOM       |MODE|                   ROLE                       |
00046 ! |________________|____|______________________________________________|
00047 ! |  IKLE          |<-->| NUMERO GLOBAUX DES NOEUDS DE CHAQUE ELEMENT  |
00048 ! |________________|____|______________________________________________
00049 ! | COMMON:        |    |
00050 ! |  GEO:          |    |
00051 ! |    MESH        | -->| TYPE DES ELEMENTS DU MAILLAGE
00052 ! |    NDP         | -->| NOMBRE DE NOEUDS PAR ELEMENTS
00053 ! |    NPOIN       | -->| NOMBRE TOTAL DE NOEUDS DU MAILLAGE
00054 ! |    NELEM       | -->| NOMBRE TOTAL D'ELEMTS DU MAILLAGE
00055 ! |    NPMAX       | -->| DIMENSION EFFECTIVE DES TABLEAUX X ET Y
00056 ! |                |    | (NPMAX = NPOIN + 0.1*NELEM)
00057 ! |    NELMAX      | -->| DIMENSION EFFECTIVE DES TABLEAUX CONCERNANT
00058 ! |                |    | LES ELEMENTS (NELMAX = NELEM + 0.2*NELEM)
00059 ! |________________|____|______________________________________________|
00060 ! MODE : -->(DONNEE NON MODIFIEE), <--(RESULTAT), <-->(DONNEE MODIFIEE)
00061 !-----------------------------------------------------------------------
00062 ! APPELE PAR : STBTEL
00063 ! APPEL DE : ECHELE
00064 !***********************************************************************
00065 !
00066       IMPLICIT NONE
00067       INTEGER LNG,LU
00068       COMMON/INFO/LNG,LU
00069 !
00070       INTEGER IELEM , NELEM , IEL1 , IEL2
00071       INTEGER I1 , I2 , I3 , J1 , J2 , J3
00072       INTEGER LGVEC , MESH , NDP , K
00073       INTEGER NPMAX , NPOIN , NELMAX
00074       INTEGER IKLE(NELMAX,4)
00075 !
00076       LOGICAL DEP
00077 !
00078       COMMON/GEO/ MESH , NDP , NPOIN , NELEM , NPMAX , NELMAX
00079 !
00080 !=======================================================================
00081 !
00082       DO IELEM = 1,NELEM
00083         IEL2 = IELEM
00084 25      CONTINUE
00085         DEP = .FALSE.
00086         I1 = IKLE(IELEM,1)
00087         I2 = IKLE(IELEM,2)
00088         I3 = IKLE(IELEM,3)
00089         DO K = 2,LGVEC
00090           IEL1 = MOD(NELEM+IELEM-K,NELEM) + 1
00091           J1 = IKLE(IEL1,1)
00092           J2 = IKLE(IEL1,2)
00093           J3 = IKLE(IEL1,3)
00094           IF (I1.EQ.J1.OR.I1.EQ.J2.OR.I1.EQ.J3.OR.
00095      &        I2.EQ.J1.OR.I2.EQ.J2.OR.I2.EQ.J3.OR.
00096      &        I3.EQ.J1.OR.I3.EQ.J2.OR.I3.EQ.J3) DEP=.TRUE.
00097         ENDDO
00098         IF (DEP) THEN
00099           IEL2 = MOD(IEL2,NELEM) + 1
00100           IF (IEL2.EQ.IELEM) GOTO 40
00101           CALL ECHELE(IKLE,IELEM,IEL2)
00102           GOTO 25
00103         ENDIF
00104       ENDDO ! IELEM
00105 !
00106 !=======================================================================
00107 !
00108 40    CONTINUE
00109 !
00110       RETURN
00111       END

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0