norma1.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	norma1 (n1, n2, n3, area, xx1, xx2, xx3)
subroutine	norma1d (n1, n2, n3, area, xx1, xx2, xx3)
subroutine	norma4n (n1, n2, n3, area, irect, x)
subroutine	norma4nx (n1, n2, n3, area, irect, x, i1, i2, shf)

Function/Subroutine Documentation

norma1()

subroutine norma1	(		n1,
			n2,
			n3,
			area,
		dimension(4), intent(in)	xx1,
		dimension(4), intent(in)	xx2,
		dimension(4), intent(in)	xx3 )

Definition at line 37 of file norma1.F.

C
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
     .   n1, n2, n3, area
      my_real, INTENT(IN) :: xx1(4),xx2(4),xx3(4)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real
     .   xx13, yy13, zz13, xx24, yy24, zz24
C
      xx13 =xx1(3)-xx1(1)
      yy13 =xx2(3)-xx2(1)
      zz13 =xx3(3)-xx3(1)
      xx24 =xx1(4)-xx1(2)
      yy24 =xx2(4)-xx2(2)
      zz24 =xx3(4)-xx3(2)
C
      n1=yy13*zz24-zz13*yy24
      n2=zz13*xx24-xx13*zz24
      n3=xx13*yy24-yy13*xx24
      area= max(em20,sqrt(n1*n1+n2*n2+n3*n3))
      n1=n1/area
      n2=n2/area
      n3=n3/area
      area=half*area
      RETURN

norma1d()

subroutine norma1d	(		n1,
			n2,
			n3,
			area,
		dimension(4), intent(in)	xx1,
		dimension(4), intent(in)	xx2,
		dimension(4), intent(in)	xx3 )

Definition at line 80 of file norma1.F.

C
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
     .   n1, n2, n3, area
      my_real, INTENT(IN) :: xx1(4),xx2(4),xx3(4)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER J
      my_real
     .   xx13, yy13, zz13, xx24, yy24, zz24
C
      xx13 =xx1(3)-xx1(1)
      yy13 =xx2(3)-xx2(1)
      zz13 =xx3(3)-xx3(1)
      xx24 =xx1(4)-xx1(2)
      yy24 =xx2(4)-xx2(2)
      zz24 =xx3(4)-xx3(2)
C
      n1=yy13*zz24-zz13*yy24
      n2=zz13*xx24-xx13*zz24
      n3=xx13*yy24-yy13*xx24
      area= max(em20,sqrt(n1*n1+n2*n2+n3*n3))
      n1=n1/area
      n2=n2/area
      n3=n3/area
C
      area=half*area
      RETURN

norma4n()

subroutine norma4n	(		n1,
			n2,
			n3,
			area,
		integer, dimension(4)	irect,
			x )

Definition at line 131 of file norma1.F.

C
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 IRECT(4)
      my_real
     .   n1, n2, n3, area, x(3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER J,NN
C     REAL
      my_real
     .   x1(4),x2(4),x3(4),xx13, yy13, zz13, xx24, yy24, zz24
C
       DO j=1,4
        nn=irect(j)
        x1(j)=x(1,nn)
        x2(j)=x(2,nn)
        x3(j)=x(3,nn)
       END DO 
C       
      xx13 =x1(3)-x1(1)
      yy13 =x2(3)-x2(1)
      zz13 =x3(3)-x3(1)
      xx24 =x1(4)-x1(2)
      yy24 =x2(4)-x2(2)
      zz24 =x3(4)-x3(2)
C
      n1=yy13*zz24-zz13*yy24
      n2=zz13*xx24-xx13*zz24
      n3=xx13*yy24-yy13*xx24
      area= max(em20,sqrt(n1*n1+n2*n2+n3*n3))
      n1=n1/area
      n2=n2/area
      n3=n3/area
C
      area=half*area
C      
      RETURN

norma4nx()

subroutine norma4nx	(		n1,
			n2,
			n3,
			area,
		integer, dimension(4)	irect,
			x,
		integer	i1,
		integer	i2,
			shf )

Definition at line 182 of file norma1.F.

C
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 IRECT(4),I1,I2
      my_real
     .   n1, n2, n3, area, x(3,*),shf(3)
C-----------------------------------------------
c FUNCTION: calculate normale with perturbation
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER J,NN
C     REAL
      my_real
     .   x1(4),x2(4),x3(4),xx13, yy13, zz13, xx24, yy24, zz24
C
       DO j=1,4
        nn=irect(j)
        IF (nn == i1 .OR. nn == i2) THEN
         x1(j)=x(1,nn)
         x2(j)=x(2,nn)
         x3(j)=x(3,nn)
        ELSE
         x1(j)=x(1,nn)+shf(1)
         x2(j)=x(2,nn)+shf(2)
         x3(j)=x(3,nn)+shf(3)
        END IF        
       END DO 
C       
      xx13 =x1(3)-x1(1)
      yy13 =x2(3)-x2(1)
      zz13 =x3(3)-x3(1)
      xx24 =x1(4)-x1(2)
      yy24 =x2(4)-x2(2)
      zz24 =x3(4)-x3(2)
C
      n1=yy13*zz24-zz13*yy24
      n2=zz13*xx24-xx13*zz24
      n3=xx13*yy24-yy13*xx24
      area= max(em20,sqrt(n1*n1+n2*n2+n3*n3))
      n1=n1/area
      n2=n2/area
      n3=n3/area
C
      area=half*area
C      
      RETURN