w_igrav.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	w_igrav (igrav, lgrav, nodlocal, proc, llgrav_l, len_ia, numnod_l)

Function/Subroutine Documentation

w_igrav()

subroutine w_igrav	(	integer, dimension(nigrv,*)	igrav,
		integer, dimension(*)	lgrav,
		integer, dimension(*)	nodlocal,
		integer	proc,
		integer	llgrav_l,
		integer	len_ia,
		integer	numnod_l )

Definition at line 30 of file w_igrav.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------  
 
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      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER PROC, LEN_IA, LLGRAV_L,NUMNOD_L,
     .        IGRAV(NIGRV,*), LGRAV(*), NODLOCAL(*)
C-----------------------------------------------
C   F u n c t i o n
C-----------------------------------------------
!      INTEGER  NLOCAL
!      EXTERNAL NLOCAL
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J, IAD, IADL, NL, N, NLL, NN, IGRTMP(NIGRV,NGRAV)
      INTEGER, DIMENSION(:), ALLOCATABLE :: LGRTMP ! LLGRAV_L
C
      ALLOCATE(lgrtmp(llgrav_l))
      iadl = 1
      DO i = 1, ngrav
        nl = igrav(1,i)
        iad= igrav(4,i)
        DO j = 1, nigrv
          igrtmp(j,i) = igrav(j,i)
        ENDDO
        nll = 0
        DO j = 1, nl
          nn = lgrav(iad+j-1)
          n = abs(nn)
          IF (nodlocal(n)/=0.AND.nodlocal(n)<=numnod_l)THEN
            IF(nn>0) THEN
              lgrtmp(iadl+nll) = nodlocal(n)
            ELSE
              lgrtmp(iadl+nll) = -nodlocal(n)
            ENDIF
            nll = nll + 1
          ENDIF
        ENDDO
        igrtmp(1,i) = nll
        igrtmp(4,i) = iadl
        iadl = iadl+nll
      ENDDO
C
      CALL write_i_c(igrtmp,nigrv*ngrav)
      CALL write_i_c(lgrtmp,iadl-1)
      len_ia = len_ia + nigrv*ngrav + iadl-1
 
      DEALLOCATE(lgrtmp)
C
      RETURN