w_dampvrel.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	w_dampvrel (dampr, igrnod, idamp_vrel_l, ndamp_vrel_l, len_ia, ngrnod, ndamp, nrdamp, nspmd)

Function/Subroutine Documentation

w_dampvrel()

subroutine w_dampvrel	(	dimension(nrdamp,ndamp), intent(in)	dampr,
		type (group_), dimension(ngrnod), intent(in)	igrnod,
		integer, dimension(ndamp), intent(in)	idamp_vrel_l,
		integer, intent(in)	ndamp_vrel_l,
		integer, intent(inout)	len_ia,
		integer, intent(in)	ngrnod,
		integer, intent(in)	ndamp,
		integer, intent(in)	nrdamp,
		integer, intent(in)	nspmd )

Definition at line 31 of file w_dampvrel.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE groupdef_mod     
C---------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NDAMP_VREL_L,NGRNOD,NDAMP,NRDAMP,NSPMD
      INTEGER, INTENT(IN) :: IDAMP_VREL_L(NDAMP)
      INTEGER, INTENT(INOUT) :: LEN_IA
      my_real, INTENT(IN) :: dampr(nrdamp,ndamp)
      TYPE (GROUP_),DIMENSION(NGRNOD),INTENT(IN) :: IGRNOD
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
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,IGR,PROC,PMAIN,CPT
      INTEGER FR_DAMP_VREL(NSPMD+2,NDAMP_VREL_L),ID_DAMP_VREL(NDAMP_VREL_L)
C----------------------------------------------- 
C
      pmain = 1
      cpt = 0
      DO i=1,ndamp
        IF (idamp_vrel_l(i) > 0) THEN
          cpt = cpt+1  
          id_damp_vrel(cpt) = i    
        ENDIF      
      ENDDO        
C
      fr_damp_vrel(1:nspmd+2,1:ndamp_vrel_l) = 0
      DO i=1,ndamp_vrel_l
        igr   = nint(dampr(2,id_damp_vrel(i)))
C--     Tag of procs having nodes in damping                 
        DO proc = 1,nspmd
          DO j = 1,igrnod(igr)%NENTITY
            IF (nlocal(igrnod(igr)%ENTITY(j),proc)==1) THEN
              fr_damp_vrel(proc,i) = 1   
            ENDIF        
          ENDDO 
        ENDDO
C--     Find main proc for damping (smallest proc id)          
        DO proc = 1,nspmd
          IF (fr_damp_vrel(proc,i) == 1) THEN
            pmain =  proc
            exit
          ENDIF   
        ENDDO
        fr_damp_vrel(nspmd+1,i) = 0 ! Not used for now
        fr_damp_vrel(nspmd+2,i) = pmain
      ENDDO
C
      CALL write_i_c(id_damp_vrel,ndamp_vrel_l)
      len_ia = len_ia + ndamp_vrel_l
C      
      CALL write_i_c(fr_damp_vrel,ndamp_vrel_l*(nspmd+2)) 
      len_ia = len_ia + ndamp_vrel_l*(nspmd+2)
C                
! --------------------------------------
      RETURN