s10volj.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	s10volj (w, x1b, x2b, x3b, x4a, x5b, x6b, x7b, x8b, x9b, x10b, x8a, x9a, x10a, y1b, y2b, y3b, y4a, y5b, y6b, y7b, y8b, y9b, y10b, y8a, y9a, y10a, z1b, z2b, z3b, z4a, z5b, z6b, z7b, z8b, z9b, z10b, z8a, z9a, z10a, voldp, nel)

Function/Subroutine Documentation

s10volj()

subroutine s10volj	(		w,
		double precision, dimension(mvsiz)	x1b,
		double precision, dimension(mvsiz)	x2b,
		double precision, dimension(mvsiz)	x3b,
		double precision, dimension(mvsiz)	x4a,
		double precision, dimension(mvsiz)	x5b,
		double precision, dimension(mvsiz)	x6b,
		double precision, dimension(mvsiz)	x7b,
		double precision, dimension(mvsiz)	x8b,
		double precision, dimension(mvsiz)	x9b,
		double precision, dimension(mvsiz)	x10b,
		double precision, dimension(mvsiz)	x8a,
		double precision, dimension(mvsiz)	x9a,
		double precision, dimension(mvsiz)	x10a,
		double precision, dimension(mvsiz)	y1b,
		double precision, dimension(mvsiz)	y2b,
		double precision, dimension(mvsiz)	y3b,
		double precision, dimension(mvsiz)	y4a,
		double precision, dimension(mvsiz)	y5b,
		double precision, dimension(mvsiz)	y6b,
		double precision, dimension(mvsiz)	y7b,
		double precision, dimension(mvsiz)	y8b,
		double precision, dimension(mvsiz)	y9b,
		double precision, dimension(mvsiz)	y10b,
		double precision, dimension(mvsiz)	y8a,
		double precision, dimension(mvsiz)	y9a,
		double precision, dimension(mvsiz)	y10a,
		double precision, dimension(mvsiz)	z1b,
		double precision, dimension(mvsiz)	z2b,
		double precision, dimension(mvsiz)	z3b,
		double precision, dimension(mvsiz)	z4a,
		double precision, dimension(mvsiz)	z5b,
		double precision, dimension(mvsiz)	z6b,
		double precision, dimension(mvsiz)	z7b,
		double precision, dimension(mvsiz)	z8b,
		double precision, dimension(mvsiz)	z9b,
		double precision, dimension(mvsiz)	z10b,
		double precision, dimension(mvsiz)	z8a,
		double precision, dimension(mvsiz)	z9a,
		double precision, dimension(mvsiz)	z10a,
		double precision, dimension(mvsiz)	voldp,
		integer	nel )

Definition at line 28 of file s10volj.F.

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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NEL
C     REAL
      double precision
     .   x1b(mvsiz),x2b(mvsiz),x3b(mvsiz),x4a(mvsiz),x5b(mvsiz),
     .   x6b(mvsiz),x7b(mvsiz),x8b(mvsiz),x9b(mvsiz),x10b(mvsiz),
     .   x8a(mvsiz),x9a(mvsiz),x10a(mvsiz),
     .   y1b(mvsiz),y2b(mvsiz),y3b(mvsiz),y4a(mvsiz),y5b(mvsiz),
     .   y6b(mvsiz),y7b(mvsiz),y8b(mvsiz),y9b(mvsiz),y10b(mvsiz),
     .   y8a(mvsiz),y9a(mvsiz),y10a(mvsiz),
     .   z1b(mvsiz),z2b(mvsiz),z3b(mvsiz),z4a(mvsiz),z5b(mvsiz),
     .   z6b(mvsiz),z7b(mvsiz),z8b(mvsiz),z9b(mvsiz),z10b(mvsiz),
     .   z8a(mvsiz),z9a(mvsiz),z10a(mvsiz),voldp(mvsiz)
      my_real
     .   w
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      double precision
     .   a4m1,b4m1,aa
      double precision
     .   d,a4mb4,a4,b4,bb,ab,b2,
     .   drdx, drdy, drdz,det,
     .   dxdr,dxds,dxdt,dydr,dyds,dydt,dzdr,dzds,dzdt
C-----------------------------------------------
C
      DO i=1,nel
       aa = x5b(i) + x6b(i) + x7b(i)
     .    - x4a(i) - x8b(i) - x9b(i) - x10b(i)
       dxdr = x1b(i) +x8a(i) - x6b(i) + aa
       dxds = x2b(i) +x9a(i) - x7b(i) + aa
       dxdt = x3b(i) +x10a(i)- x5b(i) + aa
C
       aa = y5b(i) + y6b(i) + y7b(i)
     .    - y4a(i) - y8b(i) - y9b(i) - y10b(i)
       dydr = y1b(i) +y8a(i) - y6b(i) + aa
       dyds = y2b(i) +y9a(i) - y7b(i) + aa
       dydt = y3b(i) +y10a(i)- y5b(i) + aa
C
       aa = z5b(i) + z6b(i) + z7b(i)
     .    - z4a(i) - z8b(i) - z9b(i) - z10b(i)
       dzdr = z1b(i) +z8a(i) - z6b(i) + aa
       dzds = z2b(i) +z9a(i) - z7b(i) + aa
       dzdt = z3b(i) +z10a(i)- z5b(i) + aa
C
       drdx =dyds*dzdt-dzds*dydt
       drdy =dzds*dxdt-dxds*dzdt
       drdz =dxds*dydt-dyds*dxdt
C
       det = dxdr * drdx + dydr * drdy + dzdr * drdz
       voldp(i) = w * det * one_over_6
C
      ENDDO
C
      RETURN