cg1113.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\cg1113.f
00002 !
00047                      SUBROUTINE CG1113
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,6)
00078 !
00079 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00080 !
00081       INTEGER IELEM,IDIM
00082 !
00083 !-----------------------------------------------------------------------
00084 !
00085 !VOCL LOOP,NOVREC
00086 !DIR$ IVDEP
00087       DO IDIM  = 1 , DIM2
00088       DO IELEM = 1 , NELEM
00089 !
00090         X(IKLE(IELEM,4),IDIM) = 0.5D0 * ( X(IKLE(IELEM,1),IDIM)
00091      &                                  + X(IKLE(IELEM,2),IDIM) )
00092         X(IKLE(IELEM,5),IDIM) = 0.5D0 * ( X(IKLE(IELEM,2),IDIM)
00093      &                                  + X(IKLE(IELEM,3),IDIM) )
00094         X(IKLE(IELEM,6),IDIM) = 0.5D0 * ( X(IKLE(IELEM,3),IDIM)
00095      &                                  + X(IKLE(IELEM,1),IDIM) )
00096 !
00097       ENDDO
00098       ENDDO
00099 !
00100 !-----------------------------------------------------------------------
00101 !
00102       RETURN
00103       END

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