36 material, element, face, &
37 nel_glo, normalxyz, mat_el_fac, yon)
42 integer*4 :: material, element, face, yon, nel_glo
43 integer*4,
dimension(nel_glo,3) :: mat_el_fac
45 real*8 :: nor_x, nor_y, nor_z, res
46 real*8,
dimension(nel_glo,3) :: normalxyz
52 if(material .eq. mat_el_fac(i,1) .and. &
53 element .eq. mat_el_fac(i,2) .and. &
54 face .eq. mat_el_fac(i,3) )
then
56 res = dsqrt(dabs(nor_x + normalxyz(i,1))**2 + dabs(nor_y + normalxyz
58 if(res .le. 5.d-2) yon = 1
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.