i2rep.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	i2rep (x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, nir)

Function/Subroutine Documentation

i2rep()

subroutine i2rep	(		x1,
			x2,
			x3,
			x4,
			y1,
			y2,
			y3,
			y4,
			z1,
			z2,
			z3,
			z4,
			e1x,
			e1y,
			e1z,
			e2x,
			e2y,
			e2z,
			e3x,
			e3y,
			e3z,
		integer	nir )

Definition at line 41 of file i2rep.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-----------------------------------------------
        my_real x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, e1x, e1y, e1z, e2x,e2y, e2z, e3x, e3y, e3z
        INTEGER NIR
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        my_real :: sum_,s1,s2
C-----------------------------------------------
 
        IF(nir==4)THEN
          e1x = (x2 + x3 - x1 - x4)
          e1y = (y2 + y3 - y1 - y4)
          e1z = (z2 + z3 - z1 - z4)
          e2x = (x3 + x4 - x1 - x2)
          e2y = (y3 + y4 - y1 - y2)
          e2z = (z3 + z4 - z1 - z2)
        ELSE
          e1x = (x2 + x3 - two * x1)
          e1y = (y2 + y3 - two * y1)
          e1z = (z2 + z3 - two * z1)
          e2x = (two * x3 - x1 - x2)
          e2y = (two * y3 - y1 - y2)
          e2z = (two * z3 - z1 - z2)
        END IF
        e3x = e1y*e2z - e1z*e2y
        e3y = e1z*e2x - e1x*e2z
        e3z = e1x*e2y - e1y*e2x
C
        sum_ = e3x*e3x+e3y*e3y+e3z*e3z
        sum_ = one / max(sqrt(sum_),em20)
        e3x = e3x * sum_
        e3y = e3y * sum_
        e3z = e3z * sum_
C
        IF (nir == 4) THEN
          s1  = e1x*e1x+e1y*e1y+e1z*e1z
          s2  = e2x*e2x+e2y*e2y+e2z*e2z
          sum_ = sqrt(s1/s2)
          e1x = e1x + (e2y*e3z-e2z*e3y)*sum_
          e1y = e1y + (e2z*e3x-e2x*e3z)*sum_
          e1z = e1z + (e2x*e3y-e2y*e3x)*sum_
        ENDIF
C
        sum_ = e1x*e1x+e1y*e1y+e1z*e1z
        sum_ = one / max(sqrt(sum_),em20)
        e1x = e1x * sum_
        e1y = e1y * sum_
        e1z = e1z * sum_
C
        e2x = e3y * e1z - e3z * e1y
        e2y = e3z * e1x - e3x * e1z
        e2z = e3x * e1y - e3y * e1x
C---
        RETURN