q4deric2.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	q4deric2 (y1, y2, y3, y4, z1, z2, z3, z4, y12, y34, y13, y24, y14, y23, z12, z34, z13, z24, z14, z23, pyc1, pyc2, pzc1, pzc2, aire, volu, yavg, rx, ry, rz, sx, sy, sz, nel, jhbe)

Function/Subroutine Documentation

q4deric2()

subroutine q4deric2	(		y1,
			y2,
			y3,
			y4,
			z1,
			z2,
			z3,
			z4,
			y12,
			y34,
			y13,
			y24,
			y14,
			y23,
			z12,
			z34,
			z13,
			z24,
			z14,
			z23,
			pyc1,
			pyc2,
			pzc1,
			pzc2,
			aire,
			volu,
			yavg,
			rx,
			ry,
			rz,
			sx,
			sy,
			sz,
		integer, intent(in)	nel,
		integer, intent(in)	jhbe )

Definition at line 29 of file q4deric2.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-----------------------------------------------
#include      "com01_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER, INTENT(IN) :: JHBE
C     REAL
      my_real
     .    y1(*),y2(*),y3(*),y4(*),z1(*),z2(*),z3(*),z4(*),
     .    y12(*),y34(*),y13(*),y24(*),y14(*),y23(*),
     .    z12(*),z34(*),z13(*),z24(*),z14(*),z23(*),
     .    pyc1(*),pyc2(*),pzc1(*),pzc2(*),aire(*),volu(*),
     .   rx(mvsiz) , ry(mvsiz) , rz(mvsiz) ,
     .   sx(mvsiz) , sy(mvsiz) , sz(mvsiz) ,yavg(*)
C-----------------------------------------------
c FUNCTION: 
c ARGUMENTS:  (I: input, O: output, IO: input & output, W: workspace)
c TYPE NAME                FUNCTION
c  I   Y1(*)~Z4(*)        - NODAL COORDINATES
c  O   Y12(*)~Z23(*)      - DIFFERENCE: Yi-Yj, Zi-Zj
c  O   PYC1(*)~PZC2(*)    - SHAPE DERIVATIVES AT ELEMENT CENTER
c  I   AIRE(*)            - ELEMENT AREA
c  I/O VOL(*)            - ELEMENT VOL
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
      my_real
     .    y,z,v0
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
      DO i=1,nel
        y12(i) = y1(i) - y2(i)
        y34(i) = y3(i) - y4(i)
        y13(i) = y1(i) - y3(i)
        y24(i) = y2(i) - y4(i)
        y14(i) = y1(i) - y4(i)
        y23(i) = y2(i) - y3(i)
        z12(i) = z1(i) - z2(i)
        z34(i) = z3(i) - z4(i)
        z13(i) = z1(i) - z3(i)
        z24(i) = z2(i) - z4(i)
        z14(i) = z1(i) - z4(i)
        z23(i) = z2(i) - z3(i)
      END DO
C
      DO i=1,nel
C       dNi/dY, dNi/dZ, A=4.*|J|
        pyc1(i) = half*(z2(i)-z4(i))/aire(i)
        pyc2(i) = half*(z3(i)-z1(i))/aire(i)
        pzc1(i) = half*(y4(i)-y2(i))/aire(i)
        pzc2(i) = half*(y1(i)-y3(i))/aire(i)
      ENDDO
C     COMPUTE VOLUME/(2*PI) FOR AXISYMMETRIC CASE 
      IF(n2d==1.AND.jhbe==17) THEN
          DO i=1,nel
            volu(i) = fourth*yavg(i)*aire(i)
          ENDDO
      ENDIF
      DO i=1,nel
        rx(i) = one
        ry(i) = zero
        rz(i) = zero
        sx(i) = zero
        sy(i) =half*(y2(i)+y3(i)-y1(i)-y4(i))
        sz(i) =half*(z2(i)+z3(i)-z1(i)-z4(i))
      END DO
C
      RETURN