damping_vref.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	damping_vref (igrnod, v, a, ms, dampr, weight, v_ref, a_ref, ndamp_vrel, id_damp_vrel, fr_damp_vrel, ngrnod, numnod, ndamp, nrdamp, nspmd, iparit, ispmd)

Function/Subroutine Documentation

damping_vref()

subroutine damping_vref	(	type (group_), dimension(ngrnod), intent(in)	igrnod,
		dimension(3,numnod), intent(in)	v,
		dimension(3,numnod), intent(in)	a,
		dimension(numnod), intent(in)	ms,
		dimension(nrdamp,ndamp), intent(in)	dampr,
		integer, dimension(numnod), intent(in)	weight,
		dimension(3,ndamp), intent(out)	v_ref,
		dimension(3,ndamp), intent(out)	a_ref,
		integer, intent(in)	ndamp_vrel,
		integer, dimension(ndamp_vrel), intent(in)	id_damp_vrel,
		integer, dimension(nspmd+2,ndamp_vrel), intent(in)	fr_damp_vrel,
		integer, intent(in)	ngrnod,
		integer, intent(in)	numnod,
		integer, intent(in)	ndamp,
		integer, intent(in)	nrdamp,
		integer, intent(in)	nspmd,
		integer, intent(in)	iparit,
		integer, intent(in)	ispmd )

Definition at line 34 of file damping_vref.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"
#include      "comlock.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: WEIGHT(NUMNOD),NDAMP_VREL,ID_DAMP_VREL(NDAMP_VREL)
      INTEGER, INTENT(IN) :: FR_DAMP_VREL(NSPMD+2,NDAMP_VREL)
      INTEGER, INTENT(IN) :: NGRNOD,NUMNOD,NDAMP,NRDAMP,NSPMD,IPARIT,ISPMD
      my_real, INTENT(IN) :: a(3,numnod),v(3,numnod),ms(numnod)
      my_real, INTENT(IN) :: dampr(nrdamp,ndamp)
      my_real, INTENT(OUT) :: a_ref(3,ndamp),v_ref(3,ndamp)
C-----------------------------------------------
      TYPE (GROUP_),DIMENSION(NGRNOD), INTENT(IN) :: IGRNOD
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,N,ND,IGR,NNOD,ID
      my_real :: inv_mas 
      DOUBLE PRECISION SUML6(7,6,NDAMP_VREL),SUML(7,NDAMP_VREL)  
C-----------------------------------------------
C
C======================================================================|
C     Computation of reference velocity for /DAMP/VREF with no rbody
C======================================================================|
 
      IF (iparit == 1) THEN
C        
C--     Parith/ON assembling
C   
        DO nd = 1, ndamp_vrel
          id = id_damp_vrel(nd)
          igr   = nint(dampr(2,id))
          suml6(1:7,1:6,nd) = zero
          nnod = igrnod(igr)%NENTITY
          CALL dampvref_sum6(igrnod,v,a,ms,weight,
     .                       nd,nnod,igr,suml6,ngrnod,
     .                       numnod,ndamp)
        ENDDO
C
        DO nd = 1, ndamp_vrel
          id = id_damp_vrel(nd)
          IF ((nspmd>1).AND.(fr_damp_vrel(ispmd+1,nd)==1)) THEN
            CALL spmd_exch_fr6(fr_damp_vrel(1,nd),suml6(1,1,nd),7*6)
          ENDIF   
C
          inv_mas=one/max(em20,(suml6(7,1,nd)+suml6(7,2,nd)+suml6(7,3,nd)+
     .                          suml6(7,4,nd)+suml6(7,5,nd)+suml6(7,6,nd)))         
C          
          v_ref(1,id)=inv_mas*(suml6(1,1,nd)+suml6(1,2,nd)+suml6(1,3,nd)+
     .                         suml6(1,4,nd)+suml6(1,5,nd)+suml6(1,6,nd))
          v_ref(2,id)=inv_mas*(suml6(2,1,nd)+suml6(2,2,nd)+suml6(2,3,nd)+
     .                         suml6(2,4,nd)+suml6(2,5,nd)+suml6(2,6,nd))
          v_ref(3,id)=inv_mas*(suml6(3,1,nd)+suml6(3,2,nd)+suml6(3,3,nd)+
     .                         suml6(3,4,nd)+suml6(3,5,nd)+suml6(3,6,nd))
          a_ref(1,id)=inv_mas*(suml6(4,1,nd)+suml6(4,2,nd)+suml6(4,3,nd)+
     .                         suml6(4,4,nd)+suml6(4,5,nd)+suml6(4,6,nd))
          a_ref(2,id)=inv_mas*(suml6(5,1,nd)+suml6(5,2,nd)+suml6(5,3,nd)+
     .                         suml6(5,4,nd)+suml6(5,5,nd)+suml6(5,6,nd))
          a_ref(3,id)=inv_mas*(suml6(6,1,nd)+suml6(6,2,nd)+suml6(6,3,nd)+
     .                         suml6(6,4,nd)+suml6(6,5,nd)+suml6(6,6,nd))
        ENDDO          
C        
      ELSE
C        
C--     Parith/OFF assembling        
C
        DO nd = 1, ndamp_vrel
          id = id_damp_vrel(nd)
          suml(1:7,nd) = zero
          igr   = nint(dampr(2,id))
          DO n=1,igrnod(igr)%NENTITY
            i=igrnod(igr)%ENTITY(n)
            suml(1,nd) = suml(1,nd) + ms(i)*v(1,i)*weight(i)
            suml(2,nd) = suml(2,nd) + ms(i)*v(2,i)*weight(i)
            suml(3,nd) = suml(3,nd) + ms(i)*v(3,i)*weight(i)
            suml(4,nd) = suml(4,nd) + ms(i)*a(1,i)*weight(i)
            suml(5,nd) = suml(5,nd) + ms(i)*a(2,i)*weight(i)
            suml(6,nd) = suml(6,nd) + ms(i)*a(3,i)*weight(i)
            suml(7,nd) = suml(7,nd) + ms(i)*weight(i)
          ENDDO 
        ENDDO
C
        DO nd = 1, ndamp_vrel
          id = id_damp_vrel(nd)
          IF ((nspmd>1).AND.(fr_damp_vrel(ispmd+1,nd)==1)) THEN           
            CALL spmd_fr_poff(fr_damp_vrel(1,nd),suml(1,nd),7)
          ENDIF  
          inv_mas = one/max(em20,suml(7,nd))
          v_ref(1,id) = suml(1,nd)*inv_mas
          v_ref(2,id) = suml(2,nd)*inv_mas
          v_ref(3,id) = suml(3,nd)*inv_mas
          a_ref(1,id) = suml(4,nd)*inv_mas
          a_ref(2,id) = suml(5,nd)*inv_mas
          a_ref(3,id) = suml(6,nd)*inv_mas 
        ENDDO         
C       
      ENDIF  
C
      RETURN
C