admnorm.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	admnorm4 (n, ixc, x)
subroutine	admnorm3 (n, ixtg, x)

Function/Subroutine Documentation

admnorm3()

subroutine admnorm3	(	integer	n,
		integer, dimension(nixtg,*)	ixtg,
			x )

Definition at line 109 of file admnorm.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE remesh_mod
      use element_mod , only : nixtg
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER N, IXTG(NIXTG,*)
      my_real
     .        x(3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER N1,N2,N3
      my_real
     .       x12,y12,z12,x13,y13,z13,nx,ny,nz,aaa
C-----------------------------------------------
      n1=ixtg(2,n)
      n2=ixtg(3,n)
      n3=ixtg(4,n)
 
      tagnod(n1)=1
      tagnod(n2)=1
      tagnod(n3)=1
 
      x13 = x(1,n3) - x(1,n1)
      y13 = x(2,n3) - x(2,n1)
      z13 = x(3,n3) - x(3,n1)
 
      x12 = x(1,n2) - x(1,n1)
      y12 = x(2,n2) - x(2,n1)
      z12 = x(3,n2) - x(3,n1)
 
      nx = y12*z13 - z12*y13
      ny = z12*x13 - x12*z13
      nz = x12*y13 - y12*x13
 
      aaa=one/max(em30,sqrt(nx*nx+ny*ny+nz*nz))
      nx = nx * aaa
      ny = ny * aaa
      nz = nz * aaa
 
      nodnorm(1,n1)=nodnorm(1,n1)+nx
      nodnorm(2,n1)=nodnorm(2,n1)+ny
      nodnorm(3,n1)=nodnorm(3,n1)+nz
C  
      nodnorm(1,n2)=nodnorm(1,n2)+nx
      nodnorm(2,n2)=nodnorm(2,n2)+ny
      nodnorm(3,n2)=nodnorm(3,n2)+nz
C  
      nodnorm(1,n3)=nodnorm(1,n3)+nx
      nodnorm(2,n3)=nodnorm(2,n3)+ny
      nodnorm(3,n3)=nodnorm(3,n3)+nz
 
      RETURN

admnorm4()

subroutine admnorm4	(	integer	n,
		integer, dimension(nixc,*)	ixc,
			x )

Definition at line 31 of file admnorm.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE remesh_mod
      use element_mod , only : nixc
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER N, IXC(NIXC,*)
      my_real
     .        x(3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER SON,LEVEL,IP,NLEV,N1,N2,N3,N4
      my_real
     .        x13,y13,z13,x24,y24,z24,nx,ny,nz,aaa
C-----------------------------------------------
      n1=ixc(2,n)
      n2=ixc(3,n)
      n3=ixc(4,n)
      n4=ixc(5,n)
 
      tagnod(n1)=1
      tagnod(n2)=1
      tagnod(n3)=1
      tagnod(n4)=1
 
      x13 = x(1,n3) - x(1,n1)
      y13 = x(2,n3) - x(2,n1)
      z13 = x(3,n3) - x(3,n1)
 
      x24 = x(1,n4) - x(1,n2)
      y24 = x(2,n4) - x(2,n2)
      z24 = x(3,n4) - x(3,n2)
 
      nx = y13*z24 - z13*y24
      ny = z13*x24 - x13*z24
      nz = x13*y24 - y13*x24
 
      aaa=one/max(em30,sqrt(nx*nx+ny*ny+nz*nz))
      nx = nx * aaa
      ny = ny * aaa
      nz = nz * aaa
 
      nodnorm(1,n1)=nodnorm(1,n1)+nx
      nodnorm(2,n1)=nodnorm(2,n1)+ny
      nodnorm(3,n1)=nodnorm(3,n1)+nz
C  
      nodnorm(1,n2)=nodnorm(1,n2)+nx
      nodnorm(2,n2)=nodnorm(2,n2)+ny
      nodnorm(3,n2)=nodnorm(3,n2)+nz
C  
      nodnorm(1,n3)=nodnorm(1,n3)+nx
      nodnorm(2,n3)=nodnorm(2,n3)+ny
      nodnorm(3,n3)=nodnorm(3,n3)+nz
C  
      nodnorm(1,n4)=nodnorm(1,n4)+nx
      nodnorm(2,n4)=nodnorm(2,n4)+ny
      nodnorm(3,n4)=nodnorm(3,n4)+nz
 
      RETURN