dampvref_sum6.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "lockon.inc"
#include "lockoff.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	dampvref_sum6 (igrnod, v, a, ms, weight, nd, nnod, igr, suml6, ngrnod, numnod, ndamp)

Function/Subroutine Documentation

dampvref_sum6()

subroutine dampvref_sum6	(	type (group_), dimension(ngrnod), intent(in)	igrnod,
		dimension(3,numnod), intent(in)	v,
		dimension(3,numnod), intent(in)	a,
		dimension(numnod), intent(in)	ms,
		integer, dimension(numnod), intent(in)	weight,
		integer, intent(in)	nd,
		integer, intent(in)	nnod,
		integer, intent(in)	igr,
		double precision, dimension(7,6,ndamp), intent(inout)	suml6,
		integer, intent(in)	ngrnod,
		integer, intent(in)	numnod,
		integer, intent(in)	ndamp )

Definition at line 32 of file dampvref_sum6.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),ND,NNOD,IGR,NGRNOD,NUMNOD,NDAMP
      my_real, INTENT(IN) :: v(3,numnod),a(3,numnod),ms(numnod)
      DOUBLE PRECISION, INTENT(INOUT) :: SUML6(7,6,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,K,N
      INTEGER INDEX(NNOD),NINDEX
      my_real :: f1(nnod),f2(nnod),f3(nnod),f4(nnod),f5(nnod),
     .           f6(nnod),f7(nnod)
      DOUBLE PRECISION :: SUML6_L(7,6)
C-----------------------------------------------
C
C======================================================================|
C     Parithon computation of reference velocity for /DAMP/VREF 
C======================================================================|
C
      suml6_l(1:7,1:6) = zero
C      
      nindex = 0
      DO n=1,nnod
        i=igrnod(igr)%ENTITY(n)
        nindex = nindex+1
        index(nindex) = i
        f1(nindex) = ms(i)*v(1,i)*weight(i)
        f2(nindex) = ms(i)*v(2,i)*weight(i)
        f3(nindex) = ms(i)*v(3,i)*weight(i)
        f4(nindex) = ms(i)*a(1,i)*weight(i)
        f5(nindex) = ms(i)*a(2,i)*weight(i)
        f6(nindex) = ms(i)*a(3,i)*weight(i)
        f7(nindex) = ms(i)*weight(i)
      ENDDO
C 
      CALL sum_6_float(1, nindex, f1, suml6_l(1,1), 7)
      CALL sum_6_float(1, nindex, f2, suml6_l(2,1), 7)
      CALL sum_6_float(1, nindex, f3, suml6_l(3,1), 7)
      CALL sum_6_float(1, nindex, f4, suml6_l(4,1), 7)
      CALL sum_6_float(1, nindex, f5, suml6_l(5,1), 7)
      CALL sum_6_float(1, nindex, f6, suml6_l(6,1), 7)
      CALL sum_6_float(1, nindex, f7, suml6_l(7,1), 7)
C
#include "lockon.inc"
      DO k = 1, 6
        suml6(1,k,nd) = suml6(1,k,nd)+suml6_l(1,k)
        suml6(2,k,nd) = suml6(2,k,nd)+suml6_l(2,k)
        suml6(3,k,nd) = suml6(3,k,nd)+suml6_l(3,k)
        suml6(4,k,nd) = suml6(4,k,nd)+suml6_l(4,k)
        suml6(5,k,nd) = suml6(5,k,nd)+suml6_l(5,k)
        suml6(6,k,nd) = suml6(6,k,nd)+suml6_l(6,k)
        suml6(7,k,nd) = suml6(7,k,nd)+suml6_l(7,k)
      END DO
#include "lockoff.inc"
C
      RETURN
C