w_irbe3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	w_irbe3 (irbe3, lrbe3, nodlocal, proc, llrbe3_l, itabrbe3m, nbddrbe3m, numnod_l, len_ia)

Function/Subroutine Documentation

w_irbe3()

subroutine w_irbe3	(	integer, dimension(nrbe3l,*)	irbe3,
		integer, dimension(*)	lrbe3,
		integer, dimension(*)	nodlocal,
		integer	proc,
		integer	llrbe3_l,
		integer, dimension(*)	itabrbe3m,
		integer	nbddrbe3m,
		integer	numnod_l,
		integer	len_ia )

Definition at line 31 of file w_irbe3.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"
#include      "tabsiz_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER PROC, LLRBE3_L, LEN_IA,ITABRBE3M(*),NBDDRBE3M,NUMNOD_L,
     .        IRBE3(NRBE3L,*), LRBE3(*), 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, K, KL, NSNL, NSN, ILP, NLL, NN, M, N,
     .        IBVTMP(NRBE3L,NRBE3), LBVTMP(2*LLRBE3_L),IAD
C
      kl=0
      nsnl = 0
      DO i = 1, nrbe3
       k = irbe3(1,i)
       m = irbe3(3,i)
       nsn = irbe3(5,i)
       iad = nsnl
       ilp = 0
       DO j = 1, nsn
          n = lrbe3(k+j)
          IF (nlocal(n,proc+1)==1)THEN
            ilp = ilp + 1
            nsnl = nsnl + 1
            lbvtmp(nsnl) = nodlocal(n)
            lbvtmp(nsnl+llrbe3_l) = lrbe3(k+j+slrbe3/2)
          ENDIF
       ENDDO
       IF (nlocal(m,proc+1)==1.OR.ilp>0)THEN
        kl = kl + 1
        DO j = 1, nrbe3l
          ibvtmp(j,kl) = irbe3(j,i)
        ENDDO
        ibvtmp(1,kl) = iad
        ibvtmp(5,kl) = nsnl-iad
        IF (nlocal(m,proc+1)==1.AND.ibvtmp(5,kl)>0)THEN
         ibvtmp(3,kl) = nodlocal(m)
        ELSE
         ibvtmp(3,kl) = 0
        ENDIF
C SPMD issue global RBE2 Number
        ibvtmp(10,kl) = i
       ENDIF
      ENDDO
C
      CALL write_i_c(ibvtmp,nrbe3l*kl)
      CALL write_i_c(lbvtmp,llrbe3_l*2)
      len_ia = len_ia + nrbe3l*kl + llrbe3_l*2
C
      IF(nbddrbe3m>0) THEN
        DO i = 1, numnod_l
          itabrbe3m(i) = 0
        ENDDO
          ilp=0
          DO k = 1, nsnl
            n = lbvtmp(k)
             IF(itabrbe3m(n)==0) THEN
              ilp = ilp + 1
              itabrbe3m(n) = ilp
            ENDIF
          END DO
      ENDIF
C
      RETURN