reaction_forces_th.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "com08_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	reaction_forces_th (nodft, nodlt, a, ar, ms, in, fthreac, iflag, nodreac)

Function/Subroutine Documentation

reaction_forces_th()

subroutine reaction_forces_th	(	integer, intent(in)	nodft,
		integer, intent(in)	nodlt,
		dimension(3,numnod), intent(in)	a,
		dimension(3,numnod), intent(in)	ar,
		dimension(numnod), intent(in)	ms,
		dimension(iroddl*numnod), intent(in)	in,
		dimension(6,*), intent(inout)	fthreac,
		integer, intent(in)	iflag,
		integer, dimension(*), intent(in)	nodreac )

Definition at line 28 of file reaction_forces_th.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      "com01_c.inc"
#include      "com04_c.inc"
#include      "com08_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER,INTENT(IN) :: NODFT, NODLT, IFLAG
      INTEGER,INTENT(IN) :: NODREAC(*)
      my_real,INTENT(IN) :: a(3,numnod) , ar(3,numnod), ms(numnod), in(iroddl*numnod)
      my_real,INTENT(INOUT) :: fthreac(6,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  N
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------     
      
      !---------------------------------------------------------!
      ! ADDITIONAL CONTRIBUTION OF REACTION FORCES for /TH/NODE !
      !   First call with iflag=1, second one with iflag=1 :    !
      !      FADD = m(A-A~)                                     !
      !---------------------------------------------------------!
      DO n=nodft,nodlt
        IF (nodreac(n) > 0) THEN
          fthreac(1,nodreac(n)) = fthreac(1,nodreac(n)) + iflag * ms(n)*a(1,n)*dt12
          fthreac(2,nodreac(n)) = fthreac(2,nodreac(n)) + iflag * ms(n)*a(2,n)*dt12
          fthreac(3,nodreac(n)) = fthreac(3,nodreac(n)) + iflag * ms(n)*a(3,n)*dt12
        ENDIF ! IF (NODREAC(N) > 0)
      ENDDO
 
      IF (iroddl/=0) THEN
        DO n=nodft,nodlt
          IF (nodreac(n) > 0) THEN
            fthreac(4,nodreac(n)) = fthreac(4,nodreac(n)) + iflag * in(n)*ar(1,n)*dt12
            fthreac(5,nodreac(n)) = fthreac(5,nodreac(n)) + iflag * in(n)*ar(2,n)*dt12
            fthreac(6,nodreac(n)) = fthreac(6,nodreac(n)) + iflag * in(n)*ar(3,n)*dt12
          ENDIF ! IF (NODREAC(N) > 0)
        ENDDO
      END IF
 
C-----------------------------------------------
      RETURN