s8zked3.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	s8zked3 (lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, dd, kij, is)

Function/Subroutine Documentation

s8zked3()

subroutine s8zked3	(	integer	lft,
		integer	llt,
			pxi,
			pyi,
			pzi,
			pxj,
			pyj,
			pzj,
			dd,
			kij,
		integer	is )

Definition at line 32 of file s8zked3.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 LFT,LLT,IS
C     REAL
      my_real
     .   pxi(*), pxj(*),pyi(*), pyj(*), pzi(*), pzj(*),  
     .   dd(3,3,*),kij(3,3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
C-----------------------------------------------
      IF (is==1) THEN
       DO i=lft,llt
        kij(1,1,i)=kij(1,1,i)+pxi(i)*dd(1,1,i)*pxj(i)
        kij(2,2,i)=kij(2,2,i)+pyi(i)*dd(2,2,i)*pyj(i)
        kij(3,3,i)=kij(3,3,i)+pzi(i)*dd(3,3,i)*pzj(i)
        kij(1,2,i)=kij(1,2,i)+pxi(i)*dd(1,2,i)*pyj(i)
        kij(1,3,i)=kij(1,3,i)+pxi(i)*dd(1,3,i)*pzj(i)
        kij(2,3,i)=kij(2,3,i)+pyi(i)*dd(2,3,i)*pzj(i)
       ENDDO
      ELSE
       DO i=lft,llt
        kij(1,1,i)=kij(1,1,i)+pxi(i)*dd(1,1,i)*pxj(i)
        kij(2,2,i)=kij(2,2,i)+pyi(i)*dd(2,2,i)*pyj(i)
        kij(3,3,i)=kij(3,3,i)+pzi(i)*dd(3,3,i)*pzj(i)
        kij(1,2,i)=kij(1,2,i)+pxi(i)*dd(1,2,i)*pyj(i)
        kij(2,1,i)=kij(2,1,i)+pxj(i)*dd(2,1,i)*pyi(i)
        kij(1,3,i)=kij(1,3,i)+pxi(i)*dd(1,3,i)*pzj(i)
        kij(3,1,i)=kij(3,1,i)+pxj(i)*dd(3,1,i)*pzi(i)
        kij(2,3,i)=kij(2,3,i)+pyi(i)*dd(2,3,i)*pzj(i)
        kij(3,2,i)=kij(3,2,i)+pyj(i)*dd(3,2,i)*pzi(i)
       ENDDO
      ENDIF 
      RETURN