q4deri2.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	q4deri2 (vol, ksi, eta, wi, y12, y34, y13, y24, y14, y23, z12, z34, z13, z24, z14, z23, y1, y2, y3, y4, yavg, ihbe, ngl)

Function/Subroutine Documentation

q4deri2()

subroutine q4deri2	(		vol,
			ksi,
			eta,
			wi,
			y12,
			y34,
			y13,
			y24,
			y14,
			y23,
			z12,
			z34,
			z13,
			z24,
			z14,
			z23,
			y1,
			y2,
			y3,
			y4,
			yavg,
		integer	ihbe,
		integer, dimension(*)	ngl )

Definition at line 32 of file q4deri2.F.

      USE message_mod
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      "vect01_c.inc"
#include      "com01_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IHBE, NGL(*)
      my_real
     +   vol(*),ksi,eta,wi,
     +   y12(*),y34(*),y13(*),y24(*),y14(*),y23(*),
     +   z12(*),z34(*),z13(*),z24(*),z14(*),z23(*),
     +   y1(*),y2(*),y3(*),y4(*),yavg(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     +    det(mvsiz)
      my_real
     +    qn1,qn2,qn3,qn4,yavg1
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
      DO i=lft,llt
C       8.*|J|
        det(i) = (y34(i)*z12(i)-y12(i)*z34(i))*ksi +
     +           (y23(i)*z14(i)-y14(i)*z23(i))*eta +
     +           y13(i)*z24(i)-y24(i)*z13(i)
C       W*|J|
        vol(i) = one_over_8*det(i)*wi
      ENDDO
C
      IF(n2d==1) THEN
        IF(ihbe==17) THEN
          DO i=lft,llt
            vol(i) = yavg(i)*vol(i)
          ENDDO
c        ELSE
c          QN1 = (ONE - ETA)*(ONE - KSI)
c          QN2 = (ONE - ETA)*(KSI + ONE)
c          QN3 = (ETA + ONE)*(KSI + ONE)
c          QN4 = (ETA + ONE)*(ONE - KSI)
c          DO I=LFT,LLT
c            YAVG1 = FOURTH*(QN1*Y1(I)+QN2*Y2(I)+QN3*Y3(I)+QN4*Y4(I))
c            VOL(I) = YAVG1*VOL(I)
c          ENDDO
        ENDIF
      ENDIF
C
      DO i=lft,llt
        IF(det(i)<=zero) THEN
C         WRITE(IOUT,1000) NGL(I)
C         CALL ARRET(2)
          CALL ancmsg(msgid=228,
     .                msgtype=msgerror,
     .                anmode=aninfo,
     .                i1=ngl(i))
        ENDIF
      ENDDO
C
      RETURN