cdk6defo3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	cdk6defo3 (jft, jlt, vlx, vly, vlz, vz4, vz5, vz6, px2, py2, px3, py3, pb1, pb2, pb3, vdef)

Function/Subroutine Documentation

cdk6defo3()

subroutine cdk6defo3	(	integer	jft,
		integer	jlt,
			vlx,
			vly,
			vlz,
			vz4,
			vz5,
			vz6,
			px2,
			py2,
			px3,
			py3,
			pb1,
			pb2,
			pb3,
			vdef )

Definition at line 29 of file cdk6defo3.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 JFT, JLT
      my_real
     .   px2(*),py2(*),px3(*),py3(*),pb1(mvsiz,3,3),pb2(mvsiz,3,3),pb3(mvsiz,3,6),
     .   vlx(mvsiz,2),vly(mvsiz,2),vlz(mvsiz,6),vz4(mvsiz,3),vz5(mvsiz,3),vz6(mvsiz,3),vdef(mvsiz,8)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real a5,a6,a4
C-----------------------------------------------
C
      DO i=jft,jlt
       vdef(i,1)=px2(i)*vlx(i,1)+px3(i)*vlx(i,2)
       vdef(i,2)=py2(i)*vly(i,1)+py3(i)*vly(i,2)
       vdef(i,3)=px2(i)*vly(i,1)+px3(i)*vly(i,2)
     .          +py2(i)*vlx(i,1)+py3(i)*vlx(i,2)
       vdef(i,4)=zero
       vdef(i,5)=zero
      ENDDO
C
      DO i=jft,jlt
       a5 = pb3(i,1,4)*vz5(i,1)+pb3(i,1,5)*vz5(i,2)
     .     +pb3(i,1,6)*vz5(i,3)
       a6 = pb3(i,2,4)*vz6(i,1)+pb3(i,2,5)*vz6(i,2)
     .     +pb3(i,2,6)*vz6(i,3)
       a4 = pb3(i,3,4)*vz4(i,1)+pb3(i,3,5)*vz4(i,2)
     .     +pb3(i,3,6)*vz4(i,3)
       vdef(i,6)=pb1(i,1,1)*vlz(i,1)+pb1(i,1,2)*vlz(i,2)
     .          +pb1(i,1,3)*vlz(i,3)+pb2(i,1,1)*vlz(i,4)
     .          +pb2(i,1,2)*vlz(i,5)+pb2(i,1,3)*vlz(i,6)
     .          +pb3(i,1,1)*a5+pb3(i,1,2)*a6+pb3(i,1,3)*a4
       vdef(i,7)=pb1(i,2,1)*vlz(i,1)+pb1(i,2,2)*vlz(i,2)
     .          +pb1(i,2,3)*vlz(i,3)+pb2(i,2,1)*vlz(i,4)
     .          +pb2(i,2,2)*vlz(i,5)+pb2(i,2,3)*vlz(i,6)
     .          +pb3(i,2,1)*a5+pb3(i,2,2)*a6+pb3(i,2,3)*a4
       vdef(i,8)=pb1(i,3,1)*vlz(i,1)+pb1(i,3,2)*vlz(i,2)
     .          +pb1(i,3,3)*vlz(i,3)+pb2(i,3,1)*vlz(i,4)
     .          +pb2(i,3,2)*vlz(i,5)+pb2(i,3,3)*vlz(i,6)
     .          +pb3(i,3,1)*a5+pb3(i,3,2)*a6+pb3(i,3,3)*a4
      ENDDO
C
      RETURN