GET_NODE_FROM_FACE.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	get_node_from_face (nb_nodes, nb_nz, con_spc, nb_load, l
	Computes total number of boundary nodes and second derivative on a given point x.

Function/Subroutine Documentation

get_node_from_face()

subroutine get_node_from_face	(	integer*4	nb_nodes,
		integer*4	nb_nz,
		integer*4, dimension(0:nb_nz)	con_spc,
		integer*4	nb_load,
			l
	)

Computes total number of boundary nodes and second derivative on a given point x.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	nb_nodes	number of nodes
[in]	nb_nz	length of con_spc
[in]	con_spc	spectral connectivity vector
[in]	nb_load	load number for boundary condition
[in]	lab_bc	label for boundary condition
[in]	local_n_num	= local node numbering (dummy)
[in]	node_indicator	node_indicator(i) = 1 if the node belongs to the boundary
[out]	nb_nodes_face_total	number of nodes lying on the boundary

Definition at line 33 of file GET_NODE_FROM_FACE.f90.

                                    nb_nodes_face_total, node_indicator, local_n_num)
      
      implicit none
      
      integer*4 :: nb_nodes,nb_nz,nb_load
      integer*4 :: nb_nodes_face_total
      integer*4 :: i,j,ie, nb_el, nb_nodes_face, cek_el
      
      integer*4, dimension(0:nb_nz) :: con_spc
      integer*4, dimension(nb_load) :: lab_bc
      integer*4, dimension(nb_nodes) :: node_indicator, local_n_num
 
      
      node_indicator = 0
      nb_nodes_face_total = 0
      nb_el = con_spc(0) - 1
      
      if (nb_nz .gt. 0) then
         
         do ie = 1,nb_el
            nb_nodes_face = con_spc(ie) - con_spc(ie -1) -1
            
            cek_el = 0
            do j = 1,nb_load
               if (con_spc(con_spc(ie -1)) .eq. lab_bc(j)) cek_el = 1
            enddo
            
            if (cek_el .ne. 0) then          
               do i = 1, nb_nodes_face                      
                  node_indicator(con_spc(con_spc(ie -1) +i)) = 1
               enddo
            endif
         enddo
         
         do i = 1,nb_nodes
            if (node_indicator(i) .ne. 0) nb_nodes_face_total = nb_nodes_face_total + 1
         enddo
      endif
      
      return
      

Referenced by make_extint_forces().

Here is the caller graph for this function: