s8zkehij.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	s8zkehij (lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, kij, fxx, fyy, fzz, fxy, fyz, fxz)

Function/Subroutine Documentation

s8zkehij()

subroutine s8zkehij	(	integer	lft,
		integer	llt,
			pxi,
			pyi,
			pzi,
			pxj,
			pyj,
			pzj,
			kij,
			fxx,
			fyy,
			fzz,
			fxy,
			fyz,
			fxz )

Definition at line 28 of file s8zkehij.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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER LFT,LLT
C     REAL
      my_real
     .   pxi(*), pxj(*),pyi(*), pyj(*), pzi(*), pzj(*),  
     .   fxx(*),fyy(*),fzz(*),fxy(*),fxz(*),fyz(*),kij(3,3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J
C     REAL
      my_real
     .   hij(mvsiz)
C-----------------------------------------------
      DO i=lft,llt
        hij(i)=fxx(i)*pxi(i)*pxj(i)+fyy(i)*pyi(i)*pyj(i)+
     .         fzz(i)*pzi(i)*pzj(i)+
     .         fxy(i)*(pxi(i)*pyj(i)+pyi(i)*pxj(i))+
     .         fyz(i)*(pyi(i)*pzj(i)+pzi(i)*pyj(i))+
     .         fxz(i)*(pxi(i)*pzj(i)+pzi(i)*pxj(i))
      ENDDO
       DO j = 1,3
        DO i=lft,llt
         kij(j,j,i)=kij(j,j,i)+hij(i)
        ENDDO
       ENDDO
C
      RETURN