fxbvini.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"
#include "fxbcom.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	fxbvini (fxbipm, fxbvit, fxbrpm, v, vr)

Function/Subroutine Documentation

fxbvini()

subroutine fxbvini	(	integer, dimension(nbipm,nfxbody), intent(in)	fxbipm,
		dimension(lenvar), intent(inout)	fxbvit,
		dimension(lenrpm), intent(in)	fxbrpm,
		dimension(3,numnod), intent(in)	v,
		dimension(3,numnod), intent(in)	vr )

Definition at line 29 of file fxbvini.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      "com04_c.inc"
#include      "fxbcom.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER, INTENT(IN) :: FXBIPM(NBIPM,NFXBODY)
        my_real, INTENT(IN) :: fxbrpm(lenrpm),v(3,numnod),vr(3,numnod)
        my_real, INTENT(INOUT) :: fxbvit(lenvar)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER NFX,I,AVAR,ARPM,IDMAST
        my_real p11,p12,p13,p21,p22,p23,p31,p32,p33,w1,w2,w3,wl1,wl2,wl3
 
        DO nfx=1,nfxbody
          idmast=fxbipm(2,nfx)
          avar=fxbipm(13,nfx)
          arpm=fxbipm(14,nfx)
          DO i=1,3
            fxbvit(avar+9+i-1)=v(i,idmast)
          END DO
          p11=fxbrpm(arpm+1)
          p12=fxbrpm(arpm+2)
          p13=fxbrpm(arpm+3)
          p21=fxbrpm(arpm+4)
          p22=fxbrpm(arpm+5)
          p23=fxbrpm(arpm+6)
          p31=fxbrpm(arpm+7)
          p32=fxbrpm(arpm+8)
          p33=fxbrpm(arpm+9)
          w1=vr(1,idmast)
          w2=vr(2,idmast)
          w3=vr(3,idmast)
          wl1=p11*w1+p21*w2+p31*w3
          wl2=p12*w1+p22*w2+p32*w3
          wl3=p13*w1+p23*w2+p33*w3
          fxbvit(avar)  =v(1,idmast)+p12*wl3-p13*wl2
          fxbvit(avar+1)=v(2,idmast)+p22*wl3-p23*wl2
          fxbvit(avar+2)=v(3,idmast)+p32*wl3-p33*wl2
          fxbvit(avar+3)=v(1,idmast)-p11*wl3+p13*wl1
          fxbvit(avar+4)=v(2,idmast)-p21*wl3+p23*wl1
          fxbvit(avar+5)=v(3,idmast)-p31*wl3+p33*wl1
          fxbvit(avar+6)=v(1,idmast)+p11*wl2-p12*wl1
          fxbvit(avar+7)=v(2,idmast)+p21*wl2-p22*wl1
          fxbvit(avar+8)=v(3,idmast)+p31*wl2-p32*wl1
        ENDDO
C
        RETURN