w_ibvel.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_ibvel (ibvel, lbvel, nodlocal, proc, llbvel_l, len_ia)

Function/Subroutine Documentation

w_ibvel()

subroutine w_ibvel	(	integer, dimension(nbvelp,*)	ibvel,
		integer, dimension(*)	lbvel,
		integer, dimension(*)	nodlocal,
		integer	proc,
		integer	llbvel_l,
		integer	len_ia )

Definition at line 31 of file w_ibvel.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, LLBVEL_L, LEN_IA,
     .        IBVEL(NBVELP,*), LBVEL(*), 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(NBVELP,NIBVEL), LBVTMP(LLBVEL_L)
C
      k = 0
      kl = 0
      DO i = 1, nibvel
        nsn = ibvel(3,i)
        DO j = 1, nbvelp
          ibvtmp(j,i) = ibvel(j,i)
        ENDDO
        nsnl = 0
        DO j = 1, nsn
          n = lbvel(k+j)
          IF (nlocal(n,proc+1)==1)THEN
            nsnl = nsnl + 1
            lbvtmp(kl+nsnl) = nodlocal(n)
          ENDIF
        ENDDO
        ibvtmp(3,i) = nsnl
        m = ibvel(4,i)
Coptimisation pour ne pas traiter les rb mous sans main et sans noeud secnd
        IF(nlocal(m,proc+1)==0.AND.nsnl==0)THEN
          ibvtmp(4,i) = -m
        ELSE
          ibvtmp(4,i) = nodlocal(m)
        ENDIF
        kl = kl+nsnl
        k = k + nsn
      ENDDO
C
      CALL write_i_c(ibvtmp,nbvelp*nibvel)
      CALL write_i_c(lbvtmp,llbvel_l)
      len_ia = len_ia + nbvelp*nibvel + llbvel_l
C
      RETURN