mean_node_norm2.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	mean_node_norm2 (in, nod_normal, nx, ny, nz)

Function/Subroutine Documentation

mean_node_norm2()

subroutine mean_node_norm2	(	integer, dimension(4), intent(in)	in,
		dimension(3,*), intent(inout)	nod_normal,
		intent(in)	nx,
		intent(in)	ny,
		intent(in)	nz )

Definition at line 28 of file mean_node_norm2.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,INTENT(IN) :: IN(4)
      my_real,INTENT(IN) :: nx,ny,nz
      my_real,INTENT(INOUT) :: nod_normal(3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER r,NBNOD
      my_real aaa, x0,y0,z0
      my_real x01,y01,z01, x02,y02,z02, x03,y03,z03, x04,y04,z04
      my_real xn1,yn1,zn1, xn2,yn2,zn2, xn3,yn3,zn3, xn4,yn4,zn4
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
      IF(in(4) /= in(3))THEN
        nod_normal(1,in(1))=nod_normal(1,in(1))+nx
        nod_normal(2,in(1))=nod_normal(2,in(1))+ny
        nod_normal(3,in(1))=nod_normal(3,in(1))+nz
 
        nod_normal(1,in(2))=nod_normal(1,in(2))+nx
        nod_normal(2,in(2))=nod_normal(2,in(2))+ny
        nod_normal(3,in(2))=nod_normal(3,in(2))+nz
 
        nod_normal(1,in(3))=nod_normal(1,in(3))+nx
        nod_normal(2,in(3))=nod_normal(2,in(3))+ny
        nod_normal(3,in(3))=nod_normal(3,in(3))+nz
 
        nod_normal(1,in(4))=nod_normal(1,in(4))+nx
        nod_normal(2,in(4))=nod_normal(2,in(4))+ny
        nod_normal(3,in(4))=nod_normal(3,in(4))+nz
      ELSE
        nod_normal(1,in(1))=nod_normal(1,in(1))+nx
        nod_normal(2,in(1))=nod_normal(2,in(1))+ny
        nod_normal(3,in(1))=nod_normal(3,in(1))+nz
 
        nod_normal(1,in(2))=nod_normal(1,in(2))+nx
        nod_normal(2,in(2))=nod_normal(2,in(2))+ny
        nod_normal(3,in(2))=nod_normal(3,in(2))+nz
 
        nod_normal(1,in(3))=nod_normal(1,in(3))+nx
        nod_normal(2,in(3))=nod_normal(2,in(3))+ny
        nod_normal(3,in(3))=nod_normal(3,in(3))+nz
      ENDIF
C-----------------------------------------------
      RETURN