nfacette.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	nfacette (xfas, nx, ny, nz)

Function/Subroutine Documentation

nfacette()

subroutine nfacette	(	dimension(3,4), intent(in)	xfas,
		intent(inout)	nx,
		intent(inout)	ny,
		intent(inout)	nz )

Definition at line 29 of file nfacette.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, intent(in) :: xfas(3,4)
      my_real, intent(inout) :: nx,ny,nz
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real rx,ry,rz,sx,sy,sz,sum
C-----------------------------------------------
       rx=xfas(1,3) - xfas(1,1)
       sx=xfas(1,4) - xfas(1,2)
       ry=xfas(2,3) - xfas(2,1)
       sy=xfas(2,4) - xfas(2,2)
       rz=xfas(3,3) - xfas(3,1)
       sz=xfas(3,4) - xfas(3,2)
C---------NORMAL------------
       nx = ry * sz - rz * sy 
       ny = rz * sx - rx * sz 
       nz = rx * sy - ry * sx 
       sum = sqrt(nx*nx + ny*ny + nz*nz)
C
       sum = one/max(sum,em30)
       nx = nx * sum 
       ny = ny * sum 
       nz = nz * sum 
C---
      RETURN