cg1112.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\cg1112.f
00002 !
00047                      SUBROUTINE CG1112
00048 !                    *****************
00049 !
00050      &(X,DIM1,DIM2,IKLE,NELEM,NELMAX)
00051 !
00052 !***********************************************************************
00053 ! BIEF   V6P1                                   21/08/2010
00054 !***********************************************************************
00055 !
00056 !
00057 !
00058 !
00059 !
00060 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00061 !| DIM1           |-->| FIRST DIMENSION OF X
00062 !| DIM2           |-->| SECOND DIMENSION OF X
00063 !| IKLE           |-->| CONNECTIVITY TABLE
00064 !| NELEM          |-->| NUMBER OF ELEMENTS IN THE MESH
00065 !| NELMAX         |-->| FIRST DIMENSION OF IKLE AND W.
00066 !| X              |<--| VECTOR TO BE MODIFIED.
00067 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00068 !
00069       IMPLICIT NONE
00070       INTEGER LNG,LU
00071       COMMON/INFO/LNG,LU
00072 !
00073 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00074 !
00075       INTEGER         , INTENT(IN)    :: NELEM,NELMAX,DIM1,DIM2
00076       DOUBLE PRECISION, INTENT(INOUT) :: X(DIM1,DIM2)
00077       INTEGER         , INTENT(IN)    :: IKLE(NELMAX,4)
00078 !
00079 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00080 !
00081       INTEGER IELEM,IDIM1
00082 !
00083       DOUBLE PRECISION TIERS
00084 !
00085 !-----------------------------------------------------------------------
00086 !
00087       TIERS = 1.D0/3.D0
00088 !
00089 !-----------------------------------------------------------------------
00090 !
00091 !VOCL LOOP,NOVREC
00092 !DIR$ IVDEP
00093       DO IDIM1  = 1 , DIM2
00094       DO IELEM = 1 , NELEM
00095 !
00096         X(IKLE(IELEM,4),IDIM1) = TIERS * ( X(IKLE(IELEM,1),IDIM1)
00097      &                                  + X(IKLE(IELEM,2),IDIM1)
00098      &                                  + X(IKLE(IELEM,3),IDIM1) )
00099 !
00100       ENDDO ! IELEM
00101       ENDDO ! IDIM1
00102 !
00103 !-----------------------------------------------------------------------
00104 !
00105       RETURN
00106       END

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