s6cortho3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	s6cortho3 (x1, x2, x3, x4, x5, x6, y1, y2, y3, y4, y5, y6, z1, z2, z3, z4, z5, z6, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, rx, ry, rz, sx, sy, sz, tx, ty, tz, nel)

Function/Subroutine Documentation

s6cortho3()

subroutine s6cortho3	(	dimension(mvsiz)	x1,
		dimension(mvsiz)	x2,
		dimension(mvsiz)	x3,
		dimension(mvsiz)	x4,
		dimension(mvsiz)	x5,
		dimension(mvsiz)	x6,
		dimension(mvsiz)	y1,
		dimension(mvsiz)	y2,
		dimension(mvsiz)	y3,
		dimension(mvsiz)	y4,
		dimension(mvsiz)	y5,
		dimension(mvsiz)	y6,
		dimension(mvsiz)	z1,
		dimension(mvsiz)	z2,
		dimension(mvsiz)	z3,
		dimension(mvsiz)	z4,
		dimension(mvsiz)	z5,
		dimension(mvsiz)	z6,
		dimension(mvsiz)	e1x,
		dimension(mvsiz)	e2x,
		dimension(mvsiz)	e3x,
		dimension(mvsiz)	e1y,
		dimension(mvsiz)	e2y,
		dimension(mvsiz)	e3y,
		dimension(mvsiz)	e1z,
		dimension(mvsiz)	e2z,
		dimension(mvsiz)	e3z,
		dimension(mvsiz)	rx,
		dimension(mvsiz)	ry,
		dimension(mvsiz)	rz,
		dimension(mvsiz)	sx,
		dimension(mvsiz)	sy,
		dimension(mvsiz)	sz,
		dimension(mvsiz)	tx,
		dimension(mvsiz)	ty,
		dimension(mvsiz)	tz,
		integer, intent(in)	nel )

Definition at line 29 of file s6cortho3.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, INTENT(IN) :: NEL
C     REAL
      my_real, INTENT(IN),  DIMENSION(MVSIZ) ::
     .   x1, x2, x3, x4, x5, x6, 
     .   y1, y2, y3, y4, y5, y6,   
     .   z1, z2, z3, z4, z5, z6
      my_real, INTENT(OUT),  DIMENSION(MVSIZ) ::
     .   rx, ry, rz, sx, sy, sz, tx, ty, tz,
     .   e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z
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
      my_real
     .   x14 ,y14 , z14
      my_real
     .   det,c1,c2 
C-----------------------------------------------
      DO i=1,nel
        x14=x1(i)+x4(i)
        y14=y1(i)+y4(i)
        z14=z1(i)+z4(i)
C--------ksi---------
        tx(i)=x2(i)+x5(i)-x14
        ty(i)=y2(i)+y5(i)-y14
        tz(i)=z2(i)+z5(i)-z14
C--------eta---------
        rx(i)=x3(i)+x6(i)-x14
        ry(i)=y3(i)+y6(i)-y14
        rz(i)=z3(i)+z6(i)-z14
       ENDDO
      DO i=1,nel
        sx(i)= (x4(i)+x5(i)+x6(i)-x1(i)-x2(i)-x3(i))*third
        sy(i)= (y4(i)+y5(i)+y6(i)-y1(i)-y2(i)-y3(i))*third
        sz(i)= (z4(i)+z5(i)+z6(i)-z1(i)-z2(i)-z3(i))*third
      ENDDO
c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
      DO 100 i=1,nel
C
       e3x(i) = ty(i) * rz(i) - tz(i) * ry(i) 
       e3y(i) = tz(i) * rx(i) - tx(i) * rz(i) 
       e3z(i) = tx(i) * ry(i) - ty(i) * rx(i) 
C
       det = sqrt(e3x(i)*e3x(i) + e3y(i)*e3y(i) + e3z(i)*e3z(i))
       IF ( det/=zero) det = one / det
       e3x(i) = e3x(i) * det
       e3y(i) = e3y(i) * det
       e3z(i) = e3z(i) * det
C
       c1=sqrt(tx(i)*tx(i)+ty(i)*ty(i)+tz(i)*tz(i))
       c2=sqrt(rx(i)*rx(i)+ry(i)*ry(i)+rz(i)*rz(i))
       e1x(i)=tx(i)*c2 +(ry(i) * e3z(i) - rz(i) * e3y(i))*c1  
       e1y(i)=ty(i)*c2 +(rz(i) * e3x(i) - rx(i) * e3z(i))*c1  
       e1z(i)=tz(i)*c2 +(rx(i) * e3y(i) - ry(i) * e3x(i))*c1
       det = sqrt(e1x(i)*e1x(i) + e1y(i)*e1y(i) + e1z(i)*e1z(i))
       IF ( det/=zero) det = one / det
       e1x(i) = e1x(i)*det
       e1y(i) = e1y(i)*det
       e1z(i) = e1z(i)*det
C
       e2x(i) = e3y(i) * e1z(i) - e3z(i) * e1y(i)
       e2y(i) = e3z(i) * e1x(i) - e3x(i) * e1z(i)
       e2z(i) = e3x(i) * e1y(i) - e3y(i) * e1x(i)
 100  CONTINUE
c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
      RETURN