sfint43.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	sfint43 (ipg, npg, nel, hh, areapg, soft, fx1, fx2, fx3, fx4, fx5, fx6, fx7, fx8, fy1, fy2, fy3, fy4, fy5, fy6, fy7, fy8, fz1, fz2, fz3, fz4, fz5, fz6, fz7, fz8, signzz, signyz, signzx)

Function/Subroutine Documentation

sfint43()

subroutine sfint43	(	integer	ipg,
		integer	npg,
		integer	nel,
		dimension(npg,npg)	hh,
		dimension(nel)	areapg,
		dimension(nel)	soft,
		dimension(nel)	fx1,
		dimension(nel)	fx2,
		dimension(nel)	fx3,
		dimension(nel)	fx4,
		dimension(nel)	fx5,
		dimension(nel)	fx6,
		dimension(nel)	fx7,
		dimension(nel)	fx8,
		dimension(nel)	fy1,
		dimension(nel)	fy2,
		dimension(nel)	fy3,
		dimension(nel)	fy4,
		dimension(nel)	fy5,
		dimension(nel)	fy6,
		dimension(nel)	fy7,
		dimension(nel)	fy8,
		dimension(nel)	fz1,
		dimension(nel)	fz2,
		dimension(nel)	fz3,
		dimension(nel)	fz4,
		dimension(nel)	fz5,
		dimension(nel)	fz6,
		dimension(nel)	fz7,
		dimension(nel)	fz8,
		dimension(nel)	signzz,
		dimension(nel)	signyz,
		dimension(nel)	signzx )

Definition at line 28 of file sfint43.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IPG,NEL,NPG
      my_real, DIMENSION(NEL) :: areapg,signzz,signyz,signzx,soft,
     .  fx1,fx2,fx3,fx4,fx5,fx6,fx7,fx8,fy1,fy2,fy3,fy4,fy5,fy6,fy7,fy8,         
     .  fz1,fz2,fz3,fz4,fz5,fz6,fz7,fz8
      my_real, DIMENSION(NPG,NPG) :: hh 
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real fzz,fyz,fzx,rx,ry
      my_real, DIMENSION(NEL) :: dfx1,dfx2,dfx3,dfx4,dfx5,dfx6,dfx7,dfx8,
     .                           dfy1,dfy2,dfy3,dfy4,dfy5,dfy6,dfy7,dfy8,
     .                           dfz1,dfz2,dfz3,dfz4,dfz5,dfz6,dfz7,dfz8 
C=======================================================================
      DO i=1,nel
        fzz = signzz(i)*areapg(i)*soft(i)
        fyz = signyz(i)*areapg(i)*soft(i)
        fzx = signzx(i)*areapg(i)*soft(i)
C       inf
        dfx1(i) = fzx*hh(1,ipg)
        dfx2(i) = fzx*hh(2,ipg)
        dfx3(i) = fzx*hh(3,ipg)
        dfx4(i) = fzx*hh(4,ipg)
        dfy1(i) = fyz*hh(1,ipg)
        dfy2(i) = fyz*hh(2,ipg)
        dfy3(i) = fyz*hh(3,ipg)
        dfy4(i) = fyz*hh(4,ipg)
        dfz1(i) = fzz*hh(1,ipg)
        dfz2(i) = fzz*hh(2,ipg)
        dfz3(i) = fzz*hh(3,ipg)
        dfz4(i) = fzz*hh(4,ipg)
C       sup
        dfx5(i) =-fzx*hh(1,ipg)
        dfx6(i) =-fzx*hh(2,ipg)
        dfx7(i) =-fzx*hh(3,ipg)
        dfx8(i) =-fzx*hh(4,ipg)
        dfy5(i) =-fyz*hh(1,ipg)
        dfy6(i) =-fyz*hh(2,ipg)
        dfy7(i) =-fyz*hh(3,ipg)
        dfy8(i) =-fyz*hh(4,ipg)
        dfz5(i) =-fzz*hh(1,ipg)
        dfz6(i) =-fzz*hh(2,ipg)
        dfz7(i) =-fzz*hh(3,ipg)
        dfz8(i) =-fzz*hh(4,ipg)
      ENDDO
C-----------------------------
      DO i=1,nel
C       inf
        fx1(i) = fx1(i) + dfx1(i)
        fx2(i) = fx2(i) + dfx2(i)
        fx3(i) = fx3(i) + dfx3(i)
        fx4(i) = fx4(i) + dfx4(i)
        fy1(i) = fy1(i) + dfy1(i)
        fy2(i) = fy2(i) + dfy2(i)
        fy3(i) = fy3(i) + dfy3(i)
        fy4(i) = fy4(i) + dfy4(i)
        fz1(i) = fz1(i) + dfz1(i)
        fz2(i) = fz2(i) + dfz2(i)
        fz3(i) = fz3(i) + dfz3(i)
        fz4(i) = fz4(i) + dfz4(i)
C       sup             
        fx5(i) = fx5(i) + dfx5(i)
        fx6(i) = fx6(i) + dfx6(i)
        fx7(i) = fx7(i) + dfx7(i)
        fx8(i) = fx8(i) + dfx8(i)
        fy5(i) = fy5(i) + dfy5(i)
        fy6(i) = fy6(i) + dfy6(i)
        fy7(i) = fy7(i) + dfy7(i)
        fy8(i) = fy8(i) + dfy8(i)
        fz5(i) = fz5(i) + dfz5(i)
        fz6(i) = fz6(i) + dfz6(i)
        fz7(i) = fz7(i) + dfz7(i)
        fz8(i) = fz8(i) + dfz8(i)
      ENDDO
C-----------
      RETURN