w_frbe3.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_frbe3 (irbe3, lrbe3, proc, llrbe3_l, frbe3, frbe3ad, len_ia, nrbe3pen_l)

Function/Subroutine Documentation

w_frbe3()

subroutine w_frbe3	(	integer, dimension(nrbe3l,*)	irbe3,
		integer, dimension(*)	lrbe3,
		integer	proc,
		integer	llrbe3_l,
			frbe3,
			frbe3ad,
		integer	len_ia,
		integer, intent(in)	nrbe3pen_l )

Definition at line 31 of file w_frbe3.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      "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,
     .        IRBE3(NRBE3L,*), LRBE3(*)
      INTEGER, INTENT(IN) :: NRBE3PEN_L
      my_real
     .        frbe3(6,*),frbe3ad(*)
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, NN, M, N,IAD,IADS
      my_real
     .        fbtmp1(6,llrbe3_l),fbtmp2(2*llrbe3_l),fbtmp3(10*nrbe3pen_l)
C
      nsnl = 0
      iads=slrbe3/2
      DO i = 1, nrbe3
       k = irbe3(1,i)
       nsn = irbe3(5,i)
       DO j = 1, nsn
          n = lrbe3(k+j)
          IF (nlocal(n,proc+1)==1)THEN
            nsnl = nsnl + 1
            DO nn =1,6
             fbtmp1(nn,nsnl) = frbe3(nn,k+j)
            ENDDO
            fbtmp2(nsnl) = frbe3ad(k+j)
            fbtmp2(nsnl+llrbe3_l) = frbe3ad(k+j+iads)
          ENDIF
       ENDDO
      ENDDO
C
      CALL write_db(fbtmp1,llrbe3_l*6)
      CALL write_db(fbtmp2,llrbe3_l*2)
      IF (nrbe3pen_l>0) THEN 
        fbtmp3 =zero
        CALL write_db(fbtmp3,10*nrbe3pen_l)
      END IF
      len_ia = len_ia + llrbe3_l*8 + nrbe3pen_l*10
C
      RETURN