sderito3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	sderito3 (off, det, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, jac_i, nel)

Function/Subroutine Documentation

sderito3()

subroutine sderito3	(		off,
			det,
			px1,
			px2,
			px3,
			px4,
			py1,
			py2,
			py3,
			py4,
			pz1,
			pz2,
			pz3,
			pz4,
			jac_i,
		integer	nel )

Definition at line 28 of file sderito3.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
     .   off(*),det(*),  
     .   px1(*), px2(*), px3(*), px4(*),  
     .   py1(*), py2(*), py3(*), py4(*),  
     .   pz1(*), pz2(*), pz3(*), pz4(*),  
     .   jac_i(10,mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   jaci1, jaci2, jaci3,
     .   jaci4, jaci5, jaci6,
     .   jaci7, jaci8, jaci9,
     .   jac12, jac45, jac78
C-----------------------------------------------
      DO i=1,nel
        jaci1=jac_i(1,i)
        jaci4=jac_i(4,i)
        jaci7=jac_i(7,i)
        jaci2=jac_i(2,i)
        jaci5=jac_i(5,i)
        jaci8=jac_i(8,i)
        jaci3=jac_i(3,i)
        jaci6=jac_i(6,i)
        jaci9=jac_i(9,i)
        det(i)=jac_i(10,i)
C      
        jac12=jaci1+jaci2
        jac45=jaci4+jaci5
        jac78=jaci7+jaci8
 
        px1(i)=-jac12-jaci3
        py1(i)=-jac45-jaci6
        pz1(i)=-jac78-jaci9
C
        px2(i)=-jac12+jaci3
        py2(i)=-jac45+jaci6
        pz2(i)=-jac78+jaci9
C
        jac12=jaci1-jaci2
        jac45=jaci4-jaci5
        jac78=jaci7-jaci8
 
        px3(i)= jac12+jaci3
        py3(i)= jac45+jaci6
        pz3(i)= jac78+jaci9
C
        px4(i)= jac12-jaci3
        py4(i)= jac45-jaci6
        pz4(i)= jac78-jaci9
 
 
      ENDDO
      RETURN
C