volint.F File Reference

#include "implicit_f.inc"
#include "scr08_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	volint (vol)

Function/Subroutine Documentation

volint()

subroutine volint

(

intent(out)

vol

)

Definition at line 37 of file volint.F.

C-----------------------------------------------
C   COMPUTE HEXAHEDRON VOLUME
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "scr08_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real, INTENT(OUT) :: vol
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real x17, x28, x35, x46,
     .        y17, y28, y35, y46,
     .        z17, z28, z35, z46
      my_real a17, a28, b17, b28, c17, c28 
      my_real jac1, jac2, jac3, jac4, jac5, jac6, jac7, jac8, jac9 
      my_real jac_59_68, jac_67_49, jac_48_57 
C-------------------------------------------------------
      x17 = xc(7)-xc(1)
      x28 = xc(8)-xc(2)
      x35 = xc(5)-xc(3)
      x46 = xc(6)-xc(4)
C
      y17 = yc(7)-yc(1)
      y28 = yc(8)-yc(2)
      y35 = yc(5)-yc(3)
      y46 = yc(6)-yc(4)
C
      z17 = zc(7)-zc(1)
      z28 = zc(8)-zc(2)
      z35 = zc(5)-zc(3)
      z46 = zc(6)-zc(4)
C
      jac4 = x17 + x28 - x35 - x46 
      jac5 = y17 + y28 - y35 - y46 
      jac6 = z17 + z28 - z35 - z46 
C
      a17 = x17 + x46 
      a28 = x28 + x35 
      b17 = y17 + y46 
      b28 = y28 + y35 
      c17 = z17 + z46 
      c28 = z28 + z35 
C
      jac7 = a17 + a28 
      jac8 = b17 + b28 
      jac9 = c17 + c28 
      jac1 = a17 - a28 
      jac2 = b17 - b28 
      jac3 = c17 - c28 
C
      jac_59_68 = jac5 * jac9 - jac6 * jac8 
      jac_67_49 = jac6 * jac7 - jac4 * jac9 
      jac_48_57 = jac4 * jac8 - jac5 * jac7 
C
      vol =  zep015625 * (jac1 * jac_59_68 + jac2 * jac_67_49 + jac3 * jac_48_57)
 
      RETURN