s4derito3.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	s4derito3 (off, det, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, rx, ry, rz, sx, sy, sz, tx, ty, tz, jac_i, nel)

Function/Subroutine Documentation

s4derito3()

subroutine s4derito3	(		off,
			det,
		double precision, dimension(*)	x1,
		double precision, dimension(*)	x2,
		double precision, dimension(*)	x3,
		double precision, dimension(*)	x4,
		double precision, dimension(*)	y1,
		double precision, dimension(*)	y2,
		double precision, dimension(*)	y3,
		double precision, dimension(*)	y4,
		double precision, dimension(*)	z1,
		double precision, dimension(*)	z2,
		double precision, dimension(*)	z3,
		double precision, dimension(*)	z4,
			px1,
			px2,
			px3,
			px4,
			py1,
			py2,
			py3,
			py4,
			pz1,
			pz2,
			pz3,
			pz4,
			rx,
			ry,
			rz,
			sx,
			sy,
			sz,
			tx,
			ty,
			tz,
			jac_i,
		integer, intent(in)	nel )

Definition at line 28 of file s4derito3.F.

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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      my_real
     .   off(*),det(*),
     .   px1(*), px2(*), px3(*), px4(*),  
     .   py1(*), py2(*), py3(*), py4(*),  
     .   pz1(*), pz2(*), pz3(*), pz4(*), 
     .   rx(*), ry(*), rz(*), sx(*), sy(*), sz(*),tx(*), ty(*), tz(*),
     .   jac_i(10,*)
      double precision
     .   x1(*), x2(*), x3(*), x4(*),
     .   y1(*), y2(*), y3(*), y4(*),
     .   z1(*), z2(*), z3(*), z4(*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  I
C     REAL
C                                                                     12
      my_real
     .   x41, y41, z41, x42, y42, z42, x43, y43, z43
      my_real
     .   ji1, ji2, ji3,  
     .   ji4, ji5, ji6,  
     .   ji7, ji8, ji9
C-----------------------------------------------
      DO i=1,nel
       x43 = x4(i) - x3(i)
       y43 = y4(i) - y3(i)
       z43 = z4(i) - z3(i)
       x41 = x4(i) - x1(i)
       y41 = y4(i) - y1(i)
       z41 = z4(i) - z1(i)
       x42 = x4(i) - x2(i)
       y42 = y4(i) - y2(i)
       z42 = z4(i) - z2(i)
C
       rx(i) =  -x41
       ry(i) =  -y41
       rz(i) =  -z41
       sx(i) =  -x42
       sy(i) =  -y42
       sz(i) =  -z42
       tx(i) =  -x43
       ty(i) =  -y43
       tz(i) =  -z43
C
      ENDDO
C
      DO i=1,nel
       ji1 = jac_i(1,i)
       ji2 = jac_i(2,i)
       ji3 = jac_i(3,i)
       ji4 = jac_i(4,i)
       ji5 = jac_i(5,i)
       ji6 = jac_i(6,i)
       ji7 = jac_i(7,i)
       ji8 = jac_i(8,i)
       ji9 = jac_i(9,i)
C
        px1(i)= ji1
        py1(i)= ji4
        pz1(i)= ji7
        px2(i)= ji2
        py2(i)= ji5
        pz2(i)= ji8
        px3(i)= ji3
        py3(i)= ji6
        pz3(i)= ji9
        px4(i)=-(ji1+ji2+ji3)
        py4(i)=-(ji4+ji5+ji6)
        pz4(i)=-(ji7+ji8+ji9)
      ENDDO
C
      RETURN
C