qvolu2.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "mvsiz_p.inc"
#include "com01_c.inc"
#include "com06_c.inc"
#include "scr17_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	qvolu2 (off, aire, volu, ngl, y1, y2, y3, y4, z1, z2, z3, z4, y234, y124, nel, jmult, jcvt)

Function/Subroutine Documentation

qvolu2()

subroutine qvolu2	(	dimension(nel), intent(inout)	off,
		dimension(nel), intent(inout)	aire,
		dimension(nel), intent(inout)	volu,
		integer, dimension(nel), intent(in)	ngl,
		dimension(nel), intent(in)	y1,
		dimension(nel), intent(in)	y2,
		dimension(nel), intent(in)	y3,
		dimension(nel), intent(in)	y4,
		dimension(nel), intent(in)	z1,
		dimension(nel), intent(in)	z2,
		dimension(nel), intent(in)	z3,
		dimension(nel), intent(in)	z4,
		dimension(nel), intent(in)	y234,
		dimension(nel), intent(in)	y124,
		integer, intent(in)	nel,
		integer, intent(in)	jmult,
		integer, intent(in)	jcvt )

Definition at line 37 of file qvolu2.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER, INTENT(IN) :: JMULT
      INTEGER, INTENT(IN) :: JCVT
      INTEGER, INTENT(IN) :: NGL(NEL)
      my_real, INTENT(IN) ::  y1(nel),y2(nel),y3(nel),y4(nel),z1(nel),z2(nel),z3(nel),z4(nel),y234(nel),y124(nel)
      my_real, INTENT(INOUT) :: off(nel),aire(nel),volu(nel)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com06_c.inc"
#include      "scr17_c.inc"
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-----------------------------------------------
C  S o u r c e   L i n e s
C-----------------------------------------------
      DO i=1,nel
        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
      ENDDO
C-----------------------------------------------
      IF(n2d == 1) THEN
        IF(jcvt == 0)THEN
          DO i=1,nel
            volu(i)=((y2(i)+y3(i)+y4(i))*a1(i)+(y1(i)+y2(i)+y4(i))*a2(i))*one_over_6
          ENDDO
        ELSE
          DO i=1,nel
           volu(i)=(y234(i)*a1(i)+y124(i)*a2(i))*one_over_6
          ENDDO
        END IF
      ELSE
        DO i=1,nel
          volu(i)=aire(i)
        ENDDO
      ENDIF
C-----------------------------------------------
      IF(jmult == 0)THEN
        DO i=1,nel
          IF(off(i) == zero)THEN
            volu(i)=one
            aire(i)=one
          ELSE
            IF(volu(i) <= volmin)THEN
              volu(i)=one
              off(i)=zero
              CALL ancmsg(msgid=142,anmode=aninfo,i1=ngl(i))
              idel7nok = 1
            ENDIF
          ENDIF
        ENDDO
      ENDIF
C-----------------------------------------------
      DO i=1,nel
        IF(off(i) >= one .AND. volu(i) <= zero) THEN
          CALL ancmsg(msgid=142,anmode=aninfo,i1=ngl(i))
          CALL arret(2)
        ENDIF
      ENDDO
C-----------------------------------------------
      RETURN