qvoli2.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	qvoli2 (volu, ixq, ngl, aire, y1, y2, y3, y4, z1, z2, z3, z4)

Function/Subroutine Documentation

qvoli2()

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

Definition at line 36 of file qvoli2.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(*)
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
        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)=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