ranim33.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "scr14_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	ranim33 (jft, jlt, forx, fory, forz, xmom, ymom, zmom, xl, tani, h3d_data)

Function/Subroutine Documentation

ranim33()

subroutine ranim33	(	integer	jft,
		integer	jlt,
			forx,
			fory,
			forz,
			xmom,
			ymom,
			zmom,
		double precision, dimension(mvsiz,3)	xl,
			tani,
		type(h3d_database)	h3d_data )

Definition at line 30 of file ranim33.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE h3d_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT,JLT
      my_real
     .   forx(*), fory(*), forz(*), xmom(*),
     .   ymom(*), zmom(*),tani(15,*)
      DOUBLE PRECISION XL(MVSIZ,3)
      TYPE(H3D_DATABASE) :: H3D_DATA
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "scr14_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real m1,m2,m3
C=======================================================================
C
      IF(anim_ft(1)/=0 .OR. h3d_data%UND_FORC /= 0)THEN
       DO i=jft,jlt
        tani(1,i)= forx(i)
        tani(2,i)= fory(i)
        tani(3,i)= forz(i)
        m1 = half*(xl(i,2)*forz(i)-xl(i,3)*fory(i))
        m2 = half*(xl(i,3)*forx(i)-xl(i,1)*forz(i))
        m3 = half*(xl(i,1)*fory(i)-xl(i,2)*forx(i))
        tani(4,i)= xmom(i)  + m1
        tani(5,i)= ymom(i)  + m2
        tani(6,i)= zmom(i)  + m3
        tani(7,i)= -xmom(i) + m1
        tani(8,i)= -ymom(i) + m2
        tani(9,i)= -zmom(i) + m3
       ENDDO
      ENDIF
C
      RETURN