sztorth3.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	sztorth3 (jft, jlt, isorth, nel, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, g1x, g1y, g1z, g2x, g2y, g2z, g3x, g3y, g3z, gama)

Function/Subroutine Documentation

sztorth3()

subroutine sztorth3	(	integer	jft,
		integer	jlt,
		integer	isorth,
		integer	nel,
			e1x,
			e2x,
			e3x,
			e1y,
			e2y,
			e3y,
			e1z,
			e2z,
			e3z,
			g1x,
			g1y,
			g1z,
			g2x,
			g2y,
			g2z,
			g3x,
			g3y,
			g3z,
			gama )

Definition at line 29 of file sztorth3.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C----------------------------------------------- 
      INTEGER JFT, JLT,ISORTH,NEL
C     REAL
      my_real
     .   e1x(*),e1y(*),e1z(*),
     .   e2x(*),e2y(*),e2z(*),e3x(*),e3y(*),e3z(*),
     .   g1x(*),g1y(*),g1z(*),
     .   g2x(*),g2y(*),g2z(*),g3x(*),g3y(*),g3z(*),gama(nel,6)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   r11, r12, r13, r21, r22, r23, r31, r32, r33
C=======================================================================
      IF (isorth > 0) THEN 
       DO i=jft,jlt
         r11=gama(i,1)
         r12=gama(i,2)
         r13=gama(i,3)
         r21=gama(i,4)
         r22=gama(i,5)
         r23=gama(i,6)
         r31=r12*r23-r13*r22
         r32=r13*r21-r11*r23
         r33=r11*r22-r12*r21
         g1x(i)=r11*e1x(i)+r12*e2x(i)+r13*e3x(i)
         g1y(i)=r11*e1y(i)+r12*e2y(i)+r13*e3y(i)
         g1z(i)=r11*e1z(i)+r12*e2z(i)+r13*e3z(i)
         g2x(i)=r21*e1x(i)+r22*e2x(i)+r23*e3x(i)
         g2y(i)=r21*e1y(i)+r22*e2y(i)+r23*e3y(i)
         g2z(i)=r21*e1z(i)+r22*e2z(i)+r23*e3z(i)
         g3x(i)=r31*e1x(i)+r32*e2x(i)+r33*e3x(i)
         g3y(i)=r31*e1y(i)+r32*e2y(i)+r33*e3y(i)
         g3z(i)=r31*e1z(i)+r32*e2z(i)+r33*e3z(i)    
       ENDDO
      ELSE
       DO i=jft,jlt
         g1x(i)=e1x(i)
         g1y(i)=e1y(i)
         g1z(i)=e1z(i)
         g2x(i)=e2x(i)
         g2y(i)=e2y(i)
         g2z(i)=e2z(i)
         g3x(i)=e3x(i)
         g3y(i)=e3y(i)
         g3z(i)=e3z(i)    
       ENDDO
      ENDIF 
C-----------
      RETURN