CHECK_NORMAL.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	check_normal (nor_x, nor_y, nor_z, material, element, face, nel_glo, normalxyz, mat_el_fac, yon)
	Verifies if 2 different faces are opposite or not.

Function/Subroutine Documentation

check_normal()

subroutine check_normal	(	real*8	nor_x,
		real*8	nor_y,
		real*8	nor_z,
		integer*4	material,
		integer*4	element,
		integer*4	face,
		integer*4	nel_glo,
		real*8, dimension(nel_glo,3)	normalxyz,
		integer*4, dimension(nel_glo,3)	mat_el_fac,
		integer*4	yon
	)

Verifies if 2 different faces are opposite or not.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	nor_x	x-comp of the normal vector for the current element
[in]	nor_y	y-comp of the normal vector for the current element
[in]	nor_z	z-comp of the normal vector for the current element
[in]	material	block id for the neigh. element
[in]	element	element id for the neigh. element
[in]	face	face id for the neigh. element
[in]	nel_glo	number of global element dg
[in]	nomalxyz	nomals for dg elements
[in]	mat_el_fac	id for dg elements
[out]	yon	1 if the normals are opposite 0 otherwise

Definition at line 35 of file CHECK_NORMAL.f90.

      
      implicit none
      
      integer*4 :: i
      integer*4 :: material, element, face, yon, nel_glo
      integer*4, dimension(nel_glo,3) :: mat_el_fac
      
      real*8 :: nor_x, nor_y, nor_z, res                        
      real*8, dimension(nel_glo,3) :: normalxyz
                                 
      yon = 0              
      
      do i = 1, nel_glo
      
         if(material .eq. mat_el_fac(i,1) .and. &
            element .eq. mat_el_fac(i,2)  .and. &
            face .eq. mat_el_fac(i,3) ) then
      
            res = dsqrt(dabs(nor_x + normalxyz(i,1))**2 + dabs(nor_y + normalxyz(i,2))**2 + dabs(nor_z + normalxyz(i,3))**2)
            
            if(res .le. 5.d-2) yon = 1
            
         endif
      
      enddo                     
                    
                    
                    

Referenced by write_file_dgfs().

Here is the caller graph for this function: