w_elmx.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"
#include "scr23_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	w_elmx (kxx, numelx_l, isumnx_l, ixx, cep, proc, nodlocal, len_ia)

Function/Subroutine Documentation

w_elmx()

subroutine w_elmx	(	integer, dimension(nixx,*)	kxx,
		integer	numelx_l,
		integer	isumnx_l,
		integer, dimension(*)	ixx,
		integer, dimension(*)	cep,
		integer	proc,
		integer, dimension(*)	nodlocal,
		integer	len_ia )

Definition at line 30 of file w_elmx.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
#include      "scr23_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NUMELX_L, PROC, ISUMNX_L,LEN_IA,
     .        KXX(NIXX,*), CEP(*), NODLOCAL(*), IXX(*)
C-----------------------------------------------
C   L o c a l  V a r i a b l e s
C-----------------------------------------------
      INTEGER I, NN, IAD, IAD_L, J, NB_L, NNOD, IE_L, OFF
      INTEGER, DIMENSION(:), ALLOCATABLE :: IXX_L
      INTEGER, DIMENSION(:,:), ALLOCATABLE :: KXX_L
C-----------------------------------------------
C
! ----------------------------
!     allocate 1d array
      ALLOCATE( ixx_l(isumnx_l) ) 
!     allocate 2d array
      ALLOCATE( kxx_l(nixx,numelx_l) )
! ----------------------------
      off = numelc+numeltg+numels+numelq+numelt+numelp+numelr
      ie_l = 0
      nb_l = 0
C
      DO i = 1, numelx
        IF(cep(off+i)==proc) THEN
          nnod  = kxx(3,i)
          iad   = kxx(4,i)
C
          ie_l = ie_l + 1
          DO j = 1, nixx
            kxx_l(j,ie_l) = kxx(j,i)
          ENDDO
          iad_l = nb_l+1
          kxx_l(4,ie_l) = iad_l
          nb_l = nb_l + nnod
          DO j = 1, nnod
            ixx_l(iad_l+j-1) = nodlocal(ixx(iad+j-1))
          ENDDO
        ENDIF
      ENDDO
C
      CALL write_i_c(kxx_l,numelx_l*nixx)
      CALL write_i_c(ixx_l,isumnx_l)
      len_ia = len_ia + numelx_l*nixx + isumnx_l
C
! ----------------------------
!     deallocate 1d array
      DEALLOCATE( ixx_l ) 
!     deallocate 2d array
      DEALLOCATE( kxx_l )
! ----------------------------
      RETURN