q4voli2.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	q4voli2 (volu, ixq, ngl, aire, y1, y2, y3, y4, z1, z2, z3, z4, y234, y124)

Function/Subroutine Documentation

q4voli2()

subroutine q4voli2	(		volu,
		integer, dimension(nixq,*)	ixq,
		integer, dimension(*)	ngl,
			aire,
			y1,
			y2,
			y3,
			y4,
			z1,
			z2,
			z3,
			z4,
			y234,
			y124 )

Definition at line 32 of file q4voli2.F.

      USE message_mod
      use element_mod , only : nixq
C
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"
#include      "vect01_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IXQ(NIXQ,*), NGL(*)
      my_real
     .   volu(*), aire(*), 
     .   y1(*), y2(*), y3(*), y4(*),
     .   z1(*), z2(*), z3(*), z4(*),y234(*),y124(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   a1(mvsiz), a2(mvsiz)
C----------------------------------------------------------------------
      DO i=lft,llt
        a1(i) =y2(i)*(z3(i)-z4(i))+
     .         y3(i)*(z4(i)-z2(i))+y4(i)*(z2(i)-z3(i))
        a2(i) =y2(i)*(z4(i)-z1(i))+
     .         y4(i)*(z1(i)-z2(i))+y1(i)*(z2(i)-z4(i))
        aire(i)=(a1(i)+a2(i))* half
        IF (aire(i)==zero) THEN
           CALL ancmsg(msgid=227,
     .                 msgtype=msgerror,
     .                 anmode=aninfo_blind_1,
     .                 i1=ixq(nixq,i))
        ELSEIF(aire(i)<zero.AND.n2d==1) THEN
                     CALL ancmsg(msgid=1228,
     .                 msgtype=msgerror,
     .                 anmode=aninfo_blind_1,
     .                 i1=ngl(i),
     .                 prmod=msg_cumu)
         aire(i) = -aire(i)
        ENDIF
      END DO
 
      IF(n2d==1) THEN
       IF(jcvt==0)THEN
        DO i=lft,llt
         volu(i)=((y2(i)+y3(i)+y4(i))*a1(i)
     .                 +(y1(i)+y2(i)+y4(i))*a2(i))*one_over_6
        END DO
       ELSE
        DO i=lft,llt
         volu(i)=(y234(i)*a1(i)+y124(i)*a2(i))*one_over_6
        END DO
       END IF
      ELSE
        DO i=lft,llt
          volu(i)=aire(i)
          END DO
      ENDIF
 
      DO i=lft,llt
        IF(volu(i)<=0.0)THEN
          CALL ancmsg(msgid=228,
     .                msgtype=msgerror,
     .                anmode=aninfo,
     .                i1=ngl(i))
            EXIT
        END IF
      END DO
 
      RETURN