suderi3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	suderi3 (nel, vol, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8)

Function/Subroutine Documentation

suderi3()

subroutine suderi3	(	integer	nel,
			vol,
			x1,
			x2,
			x3,
			x4,
			x5,
			x6,
			x7,
			x8,
			y1,
			y2,
			y3,
			y4,
			y5,
			y6,
			y7,
			y8,
			z1,
			z2,
			z3,
			z4,
			z5,
			z6,
			z7,
			z8 )

Definition at line 28 of file suderi3.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
      my_real
     .   vol(*), 
     .   x1(*), x2(*), x3(*), x4(*), x5(*), x6(*),
     .   x7(*), x8(*), y1(*), y2(*), y3(*), y4(*), y5(*), y6(*), y7(*),
     .   y8(*), z1(*), z2(*), z3(*), z4(*), z5(*), z6(*), z7(*), z8(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J
      
      double precision
     .   xd1(mvsiz), xd2(mvsiz), xd3(mvsiz), xd4(mvsiz),
     .   xd5(mvsiz), xd6(mvsiz), xd7(mvsiz), xd8(mvsiz),
     .   yd1(mvsiz), yd2(mvsiz), yd3(mvsiz), yd4(mvsiz),
     .   yd5(mvsiz), yd6(mvsiz), yd7(mvsiz), yd8(mvsiz),
     .   zd1(mvsiz), zd2(mvsiz), zd3(mvsiz), zd4(mvsiz),
     .   zd5(mvsiz), zd6(mvsiz), zd7(mvsiz), zd8(mvsiz)
C     REAL
      my_real
     .   jac1(mvsiz), jac2(mvsiz), jac3(mvsiz), jac4(mvsiz), jac5(mvsiz), jac6(mvsiz),
     .   jac7(mvsiz), jac8(mvsiz) , jac9(mvsiz),
     .   a17(mvsiz) , a28(mvsiz) ,
     .   b17(mvsiz) , b28(mvsiz) ,
     .   c17(mvsiz) , c28(mvsiz) 
      my_real
     .   x17(mvsiz), x28(mvsiz), x35(mvsiz), x46(mvsiz),
     .   y17(mvsiz), y28(mvsiz), y35(mvsiz), y46(mvsiz),
     .   z17(mvsiz), z28(mvsiz), z35(mvsiz), z46(mvsiz),
     .   jac_59_68(mvsiz), jac_67_49(mvsiz), jac_48_57(mvsiz)
C=======================================================================
      DO i=1,nel
        xd1(i)=x1(i)
        xd2(i)=x2(i)
        xd3(i)=x3(i)
        xd4(i)=x4(i)
        xd5(i)=x5(i)
        xd6(i)=x6(i)
        xd7(i)=x7(i)
        xd8(i)=x8(i)
        yd1(i)=y1(i)
        yd2(i)=y2(i)
        yd3(i)=y3(i)
        yd4(i)=y4(i)
        yd5(i)=y5(i)
        yd6(i)=y6(i)
        yd7(i)=y7(i)
        yd8(i)=y8(i)
        zd1(i)=z1(i)
        zd2(i)=z2(i)
        zd3(i)=z3(i)
        zd4(i)=z4(i)
        zd5(i)=z5(i)
        zd6(i)=z6(i)
        zd7(i)=z7(i)
        zd8(i)=z8(i)
      ENDDO 
 
      DO i=1,nel
        x17(i)=xd7(i)-xd1(i)
        x28(i)=xd8(i)-xd2(i)
        x35(i)=xd5(i)-xd3(i)
        x46(i)=xd6(i)-xd4(i)
        y17(i)=yd7(i)-yd1(i)
        y28(i)=yd8(i)-yd2(i)
        y35(i)=yd5(i)-yd3(i)
        y46(i)=yd6(i)-yd4(i)
        z17(i)=zd7(i)-zd1(i)
        z28(i)=zd8(i)-zd2(i)
        z35(i)=zd5(i)-zd3(i)
        z46(i)=zd6(i)-zd4(i)
      ENDDO 
C
      DO i=1,nel               
        jac1(i)=x17(i)+x28(i)-x35(i)-x46(i)
        jac2(i)=y17(i)+y28(i)-y35(i)-y46(i)
        jac3(i)=z17(i)+z28(i)-z35(i)-z46(i)
        a17(i)=x17(i)+x46(i)
        a28(i)=x28(i)+x35(i)
        b17(i)=y17(i)+y46(i)
        b28(i)=y28(i)+y35(i)
        c17(i)=z17(i)+z46(i)
        c28(i)=z28(i)+z35(i)
      ENDDO 
         
      DO i=1,nel
        jac4(i)=a17(i)+a28(i)
        jac5(i)=b17(i)+b28(i)
        jac6(i)=c17(i)+c28(i)
        jac7(i)=a17(i)-a28(i)
        jac8(i)=b17(i)-b28(i)
        jac9(i)=c17(i)-c28(i)
      ENDDO 
C
      DO i=1,nel
        jac_59_68(i)=jac5(i)*jac9(i)-jac6(i)*jac8(i)
        jac_67_49(i)=jac6(i)*jac7(i)-jac4(i)*jac9(i)
        jac_48_57(i)=jac4(i)*jac8(i)-jac5(i)*jac7(i)
      ENDDO
C
      DO i=1,nel
        vol(i)=one_over_64*(jac1(i)*jac_59_68(i)+jac2(i)*jac_67_49(i)+jac3(i)*jac_48_57(i))
      ENDDO     
C-----------
      RETURN