q4kep2.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	q4kep2 (py1, py2, pz1, pz2, ay, r22, r23, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, hh, air, fac, icp, off, nel, jcvt)

Function/Subroutine Documentation

q4kep2()

subroutine q4kep2	(		py1,
			py2,
			pz1,
			pz2,
			ay,
			r22,
			r23,
			k11,
			k12,
			k13,
			k14,
			k22,
			k23,
			k24,
			k33,
			k34,
			k44,
			hh,
			air,
			fac,
		integer	icp,
			off,
		integer, intent(in)	nel,
		integer, intent(in)	jcvt )

Definition at line 30 of file q4kep2.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) :: JCVT
      INTEGER, INTENT(IN) :: NEL
      INTEGER ICP
      my_real
     .    py1(*), py2(*), pz1(*), pz2(*), ay(*), r22(*), r23(*),
     .    k11(2,2,*), k12(2,2,*), k13(2,2,*), k14(2,2,*), k22(2,2,*), 
     .    k23(2,2,*), k24(2,2,*), k33(2,2,*), k34(2,2,*), k44(2,2,*), 
     .    hh(2,*), air(*), fac(*), off(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J, IS
      my_real
     .    bh(mvsiz), 
     .    py3(mvsiz), py4(mvsiz), pz3(mvsiz), pz4(mvsiz)
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
      DO i=1,nel
        py3(i) = -py1(i)
        py4(i) = -py2(i)
        pz3(i) = -pz1(i)
        pz4(i) = -pz2(i)
      ENDDO
C
      IF (icp==1) THEN
        DO i=1,nel
          bh(i) = (hh(1,i)+two_third*hh(2,i))*air(i)
        ENDDO
      ELSE
        DO i=1,nel
          bh(i) = (hh(1,i)*fac(i)+two_third*hh(2,i)*(one-fac(i)))*air(i)
        ENDDO
      ENDIF
C
      is = 1
      CALL q4kepij2(
     1   py1,     pz1,     py1,     pz1,
     2   ay,      r22,     r23,     bh,
     3   k11,     nel,     jcvt)
      CALL q4kepij2(
     1   py2,     pz2,     py2,     pz2,
     2   ay,      r22,     r23,     bh,
     3   k22,     nel,     jcvt)
      CALL q4kepij2(
     1   py3,     pz3,     py3,     pz3,
     2   ay,      r22,     r23,     bh,
     3   k33,     nel,     jcvt)
      CALL q4kepij2(
     1   py4,     pz4,     py4,     pz4,
     2   ay,      r22,     r23,     bh,
     3   k44,     nel,     jcvt)
      is = 0
      CALL q4kepij2(
     1   py1,     pz1,     py2,     pz2,
     2   ay,      r22,     r23,     bh,
     3   k12,     nel,     jcvt)
      CALL q4kepij2(
     1   py1,     pz1,     py3,     pz3,
     2   ay,      r22,     r23,     bh,
     3   k13,     nel,     jcvt)
      CALL q4kepij2(
     1   py1,     pz1,     py4,     pz4,
     2   ay,      r22,     r23,     bh,
     3   k14,     nel,     jcvt)
      CALL q4kepij2(
     1   py2,     pz2,     py3,     pz3,
     2   ay,      r22,     r23,     bh,
     3   k23,     nel,     jcvt)
      CALL q4kepij2(
     1   py2,     pz2,     py4,     pz4,
     2   ay,      r22,     r23,     bh,
     3   k24,     nel,     jcvt)
      CALL q4kepij2(
     1   py3,     pz3,     py4,     pz4,
     2   ay,      r22,     r23,     bh,
     3   k34,     nel,     jcvt)
C
      RETURN