s8zfintp3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	s8zfintp3 (px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, vol, pp, idege, nel)

Function/Subroutine Documentation

s8zfintp3()

subroutine s8zfintp3	(		px1,
			px2,
			px3,
			px4,
			py1,
			py2,
			py3,
			py4,
			pz1,
			pz2,
			pz3,
			pz4,
			f11,
			f21,
			f31,
			f12,
			f22,
			f32,
			f13,
			f23,
			f33,
			f14,
			f24,
			f34,
			f15,
			f25,
			f35,
			f16,
			f26,
			f36,
			f17,
			f27,
			f37,
			f18,
			f28,
			f38,
			vol,
			pp,
		integer, dimension(*)	idege,
		integer, intent(in)	nel )

Definition at line 31 of file s8zfintp3.F.

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-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER IDEGE(*)
C     REAL
      my_real
     .   px1(*), px2(*), px3(*), px4(*),  
     .   py1(*), py2(*), py3(*), py4(*),  
     .   pz1(*), pz2(*), pz3(*), pz4(*),  
     .   f11(*),f21(*),f31(*),f12(*),f22(*),f32(*),
     .   f13(*),f23(*),f33(*),f14(*),f24(*),f34(*),
     .   f15(*),f25(*),f35(*),f16(*),f26(*),f36(*),
     .   f17(*),f27(*),f37(*),f18(*),f28(*),f38(*),
     .   vol(*),pp(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J
C     REAL
      my_real
     .   sp,sx,sy,sz
C-----------------------------------------------
       DO i=1,nel
        IF (idege(i)>10) cycle
        sp =pp(i)*vol(i)
        sx =sp*px1(i)
        sy =sp*py1(i)
        sz =sp*pz1(i)
        f11(i)=f11(i)-sx
        f17(i)=f17(i)+sx
        f21(i)=f21(i)-sy
        f27(i)=f27(i)+sy
        f31(i)=f31(i)-sz
        f37(i)=f37(i)+sz
        sx =sp*px2(i)
        sy =sp*py2(i)
        sz =sp*pz2(i)
        f12(i)=f12(i)-sx
        f18(i)=f18(i)+sx
        f22(i)=f22(i)-sy
        f28(i)=f28(i)+sy
        f32(i)=f32(i)-sz
        f38(i)=f38(i)+sz
        sx =sp*px3(i)
        sy =sp*py3(i)
        sz =sp*pz3(i)
        f13(i)=f13(i)-sx
        f15(i)=f15(i)+sx
        f23(i)=f23(i)-sy
        f25(i)=f25(i)+sy
        f33(i)=f33(i)-sz
        f35(i)=f35(i)+sz
        sx =sp*px4(i)
        sy =sp*py4(i)
        sz =sp*pz4(i)
        f14(i)=f14(i)-sx
        f16(i)=f16(i)+sx
        f24(i)=f24(i)-sy
        f26(i)=f26(i)+sy
        f34(i)=f34(i)-sz
        f36(i)=f36(i)+sz
       ENDDO
c-----------
      RETURN