GET_ELEM_FROM_SURF.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	get_elem_from_surf (nnz_bc, cs_bc, n1, n2, n3, n4, ie_surf)
	Find element index from verteces number.

Function/Subroutine Documentation

get_elem_from_surf()

subroutine get_elem_from_surf	(	integer*4	nnz_bc,
		integer*4, dimension(0:nnz_bc)	cs_bc,
		integer*4	n1,
		integer*4	n2,
		integer*4	n3,
		integer*4	n4,
		integer*4	ie_surf
	)

Find element index from verteces number.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	nnz_bc	length of cs_bc
[in]	cs_bc	spectral connectivity vector
[in]	n1	node index of 1 vertex of the element (quad)
[in]	n2	node index of 1 vertex of the element (quad)
[in]	n3	node index of 1 vertex of the element (quad)
[in]	n4	node index of 1 vertex of the element (quad)
[out]	ie_surf	index of the element found

Definition at line 32 of file GET_ELEM_FROM_SURF.f90.

     
 
     implicit none
 
     integer*4 :: i, ne, nel, find
     integer*4 :: nnz_bc, n1, n2, n3, n4, ie_surf
     
     integer*4, dimension(:), allocatable :: vet
     integer*4, dimension(0:nnz_bc) :: cs_bc
     
 
     
     ne = cs_bc(0) - 1 
     
     do i = 1, ne
        nel = cs_bc(i) - cs_bc(i-1)  - 1 
        allocate(vet(nel))
        
        vet = cs_bc(cs_bc(i-1) + 1 : cs_bc(i) - 1 )
        
        if (find(nel, vet, n1) .eq. 1 .and. &
            find(nel, vet, n2) .eq. 1 .and. &
            find(nel, vet, n3) .eq. 1 .and. &
            find(nel, vet, n4) .eq. 1 ) then
           
            ie_surf = i
            return
        endif    
        
        deallocate(vet)   
     enddo
     
     ie_surf = 0
     return
     
     

Referenced by make_extint_forces().

Here is the caller graph for this function: