w_icfield.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	w_icfield (icfield, lcfield, nodlocal, proc, llcfield_l, len_ia, p, numnod_l)

Function/Subroutine Documentation

w_icfield()

subroutine w_icfield	(	integer, dimension(sizfield,*)	icfield,
		integer, dimension(*)	lcfield,
		integer, dimension(*)	nodlocal,
		integer	proc,
		integer	llcfield_l,
		integer	len_ia,
		integer	p,
		integer	numnod_l )

Definition at line 30 of file w_icfield.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"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER PROC, LEN_IA, LLCFIELD_L,P,NUMNOD_L,
     .        ICFIELD(SIZFIELD,*),LCFIELD(*),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(SIZFIELD,NLOADC),
     .        LGRTMP(LLCFIELD_L)
C
      iadl = 1
      DO i = 1, nloadc
        nl = icfield(1,i)
        iad= icfield(4,i)
        DO j = 1, sizfield
          igrtmp(j,i) = icfield(j,i)
        ENDDO
        nll = 0
        DO j = 1, nl
          nn = lcfield(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
          ELSEIF(n == 0)THEN
            lgrtmp(iadl+nll) = 0
            nll = nll + 1
          ENDIF
        ENDDO
        igrtmp(1,i) = nll
        igrtmp(4,i) = iadl
        iadl = iadl+nll
      ENDDO
C
      CALL write_i_c(igrtmp,sizfield*nloadc)
      CALL write_i_c(lgrtmp,llcfield_l)
      len_ia = len_ia + sizfield*nloadc + iadl-1
C
      RETURN