multi_i18_force_poff.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"
#include "com04_c.inc"
#include "com01_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	multi_i18_force_poff (dt, jlt, ix1, ix2, ix3, ix4, nsvg, h1, h2, h3, h4, stif, fx1, fy1, fz1, fx2, fy2, fz2, fx3, fy3, fz3, fx4, fy4, fz4, fxi, fyi, fzi, a, stifn, nin, jtask, multi_fvm, x, ixs, v, elbuf_tab, igroups, iparg, msi)

Function/Subroutine Documentation

multi_i18_force_poff()

subroutine multi_i18_force_poff	(		dt,
		integer	jlt,
		integer, dimension(mvsiz)	ix1,
		integer, dimension(mvsiz)	ix2,
		integer, dimension(mvsiz)	ix3,
		integer, dimension(mvsiz)	ix4,
		integer, dimension(mvsiz)	nsvg,
			h1,
			h2,
			h3,
			h4,
			stif,
			fx1,
			fy1,
			fz1,
			fx2,
			fy2,
			fz2,
			fx3,
			fy3,
			fz3,
			fx4,
			fy4,
			fz4,
			fxi,
			fyi,
			fzi,
			a,
			stifn,
		integer	nin,
		integer	jtask,
		type(multi_fvm_struct), intent(inout)	multi_fvm,
			x,
		integer, dimension(nixs,*)	ixs,
			v,
		type (elbuf_struct_), dimension(ngroup)	elbuf_tab,
		integer, dimension(numels)	igroups,
		integer, dimension(nparg,*)	iparg,
		intent(in)	msi )

Definition at line 32 of file multi_i18_force_poff.F.

!$COMMENT
!       MULTI_I18_FORCE_POFF description
!       accumulation of force for local and remote nodes
!       
!       MULTI_I18_FORCE_POFF organization :
!       - secondary nodes:
!            * if NSV > 0 --> local node (phantom node id = NSV - NUMNOD)
!                             accumulation in forc_int array
!            * if NSV < 0 --> remote node
!                             accumulation in AFI array
!$ENDCOMMENT
 
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE multi_fvm_mod
      USE tri7box
      USE elbufdef_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   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "param_c.inc"
#include      "com04_c.inc"
#include      "com01_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP) :: ELBUF_TAB
      INTEGER JLT, NIN,INTTH,JTASK,IGROUPS(NUMELS),IPARG(NPARG,*),
     .        IX1(MVSIZ),IX2(MVSIZ),IX3(MVSIZ),IX4(MVSIZ),NSVG(MVSIZ),IXS(NIXS,*)
      my_real, DIMENSION(2*MVSIZ), INTENT(in) :: msi
      my_real
     .    dt, h1(mvsiz),h2(mvsiz),h3(mvsiz),h4(mvsiz),stif(mvsiz),
     .    fx1(mvsiz),fy1(mvsiz),fz1(mvsiz),
     .    fx2(mvsiz),fy2(mvsiz),fz2(mvsiz),
     .    fx3(mvsiz),fy3(mvsiz),fz3(mvsiz),
     .    fx4(mvsiz),fy4(mvsiz),fz4(mvsiz),
     .    fxi(mvsiz),fyi(mvsiz),fzi(mvsiz),
     .    a(3,*), x(3,*), stifn(*), v(3,*)
      TYPE(MULTI_FVM_STRUCT), INTENT(INOUT) :: MULTI_FVM
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J1, IG, ISHIFT, NODFI, NG, ILOC, NEL, NFT
      INTEGER :: SHIFT_FORCE_INT
      my_real mass
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
c     Nodal forces on the interface
      DO i=1,jlt
         j1=ix1(i)
         a(1,j1)=a(1,j1)+fx1(i)
         a(2,j1)=a(2,j1)+fy1(i)
         a(3,j1)=a(3,j1)+fz1(i)
         stifn(j1) = stifn(j1) + stif(i)*abs(h1(i))
C     
         j1=ix2(i)
         a(1,j1)=a(1,j1)+fx2(i)
         a(2,j1)=a(2,j1)+fy2(i)
         a(3,j1)=a(3,j1)+fz2(i)
         stifn(j1) = stifn(j1) + stif(i)*abs(h2(i))
C     
         j1=ix3(i)
         a(1,j1)=a(1,j1)+fx3(i)
         a(2,j1)=a(2,j1)+fy3(i)
         a(3,j1)=a(3,j1)+fz3(i)
         stifn(j1) = stifn(j1) + stif(i)*abs(h3(i))
C     
         j1=ix4(i)
         a(1,j1)=a(1,j1)+fx4(i)
         a(2,j1)=a(2,j1)+fy4(i)
         a(3,j1)=a(3,j1)+fz4(i)
         stifn(j1) = stifn(j1) + stif(i)*abs(h4(i))
      ENDDO
C
      nodfi = nlskyfi(nin)
      ishift = nodfi*(jtask-1)
      DO i=1,jlt
         ig=nsvg(i)
         ! --------------------
         ! local node
         IF (ig > 0) THEN
            ig = ig - numnod
            ng = igroups(ig)
            nft=iparg(3,ng)
            nel=iparg(2,ng)
            iloc=ig-nft     
 
            shift_force_int = (jtask-1)*numels
            multi_fvm%FORCE_INT(1,ig+shift_force_int) = multi_fvm%FORCE_INT(1,ig+shift_force_int) - dt*fxi(i)
            multi_fvm%FORCE_INT(2,ig+shift_force_int) = multi_fvm%FORCE_INT(2,ig+shift_force_int) - dt*fyi(i)
            multi_fvm%FORCE_INT(3,ig+shift_force_int) = multi_fvm%FORCE_INT(3,ig+shift_force_int) - dt*fzi(i)
         ! --------------------
         ! remote node
         ELSE
            ig = - ig
            afi(nin)%P(1,ig+ishift) = afi(nin)%P(1,ig+ishift) - dt * fxi(i) 
            afi(nin)%P(2,ig+ishift) = afi(nin)%P(2,ig+ishift) - dt * fyi(i) 
            afi(nin)%P(3,ig+ishift) = afi(nin)%P(3,ig+ishift) - dt * fzi(i) 
         ENDIF
         ! --------------------
      ENDDO