build_ikle_ext.f

Go to the documentation of this file.

!                   *************************
                    SUBROUTINE build_ikle_ext
!                   *************************
!
     &(ikle_ext,nelmax,ikle,nelem)
!
!***********************************************************************
! BIEF   V6P3                                   05/10/2012
!***********************************************************************
!
!brief    Copy of IKLE2 into a larger array with a larger number of
!+        elements. Values of 1 put in extended area.
!         This is used in post_interp (module streamline)
!
!history  J-M HERVOUET (EDF-LNHE)
!+        05/10/2012
!+        V6P3
!+   First version
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| IKLE_EXT       |<->| EXTENDED IKLE
!| NELMAX         |-->| FIRST DIMENSION OF EXTENDED IKLE
!| IKLE           |-->| ORIGINAL IKLE
!| NELEM          |-->| FIRST DIMENSION OF ORIGINAL IKLE
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN)    :: NELMAX,NELEM
      INTEGER, INTENT(IN)    :: IKLE(nelem,3)
      INTEGER, INTENT(INOUT) :: IKLE_EXT(nelmax,3)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER IELEM
!
!-----------------------------------------------------------------------
!
      DO ielem=1,nelem
        ikle_ext(ielem,1)=ikle(ielem,1)
        ikle_ext(ielem,2)=ikle(ielem,2)
        ikle_ext(ielem,3)=ikle(ielem,3)
      ENDDO
!
      IF(nelmax.GT.nelem) THEN
        DO ielem=nelem+1,nelmax
          ikle_ext(ielem,1)=1
          ikle_ext(ielem,2)=1
          ikle_ext(ielem,3)=1
        ENDDO
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END