deparr.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\stbtel\deparr.f
00002 !
00037                         SUBROUTINE DEPARR
00038 !                       *****************
00039 !
00040      &(IKLE,NDEPAR,LGVEC)
00041 !
00042 !***********************************************************************
00043 ! STBTEL    V5P2                                   28/08/89
00044 !***********************************************************************
00045 !
00046 !
00047 !
00048 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00049 !| IKLE      |<--| GLOBAL NUMBERING OF ALL NODES BY ELEMENT
00050 !| NDEPAR    |<--| NUMBER OF BACKWARD DEPENDENCIES
00051 !| MESH      |-->| TYPE OF ELEMENTS IN THE MESH
00052 !| NDP       |-->| NUMBER OF NODE PER ELEMENT
00053 !| NPOIN     |-->| TOTAL NUMBER OF NODES IN THE MESH
00054 !| NELEM     |-->| TOTAL NUMBER OF ELEMENTS IN THE MESH
00055 !| NPMAX     |-->| EFFECTIVE DIMENSION OF ARRAYS X AND Y
00056 !|           |   | (NPMAX = NPOIN + 0.1*NELEM)
00057 !| NELMAX    |-->| EFFECTIVE DIMENSION OF THE ARRAYS DEALING WITH
00058 !|           |   | THE ELEMENTS (NELMAX = NELEM + 0.2*NELEM)
00059 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00060 !
00061       IMPLICIT NONE
00062       INTEGER LNG,LU
00063       COMMON/INFO/LNG,LU
00064 !
00065       INTEGER LGVEC , MESH , NDP , NDEPAR , IELEM , NELEM , K
00066       INTEGER NPMAX , NPOIN , NELMAX
00067       INTEGER I1 , I2 , I3 , J1 , J2 , J3 , IEL1
00068       INTEGER IKLE(NELMAX,4)
00069 !
00070       COMMON/GEO/ MESH , NDP , NPOIN , NELEM , NPMAX , NELMAX
00071 !
00072 !=======================================================================
00073 !
00074       NDEPAR = 0
00075       DO IELEM = 1,NELEM
00076         I1 = IKLE(IELEM,1)
00077         I2 = IKLE(IELEM,2)
00078         I3 = IKLE(IELEM,3)
00079         DO K = 2,LGVEC
00080           IEL1 = MOD(NELEM+IELEM-K,NELEM) + 1
00081           J1 = IKLE(IEL1,1)
00082           J2 = IKLE(IEL1,2)
00083           J3 = IKLE(IEL1,3)
00084           IF (I1.EQ.J1.OR.I1.EQ.J2.OR.I1.EQ.J3.OR.
00085      &        I2.EQ.J1.OR.I2.EQ.J2.OR.I2.EQ.J3.OR.
00086      &        I3.EQ.J1.OR.I3.EQ.J2.OR.I3.EQ.J3) NDEPAR = NDEPAR + 1
00087         ENDDO
00088       ENDDO
00089 !
00090 !=======================================================================
00091 !
00092       RETURN
00093       END

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