MAKE_WGT_GRID_NODES.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	make_wgt_grid_nodes (nb_node, nb_elem, con_matrix, node_wgt
	Computes multeplicity for mesh nodes.

Function/Subroutine Documentation

make_wgt_grid_nodes()

subroutine make_wgt_grid_nodes	(	integer*4	nb_node,
		integer*4	nb_elem,
		integer*4, dimension(nb_elem,9)	con_matrix,
		integer*4, dimension(nb_node)	node_wgt
	)

Computes multeplicity for mesh nodes.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	nb_node	grid nodes
[in]	nb_elem	number of elements (hexes)
[in]	con_matrix	con_matrix(i,1) -> element, con_matrix(i,2),...,con_matrix(i,9) -> grid nodes of element con_matrix(i,1)
[out]	node_wgt	node_wgt(i) matrix containing the multiplicity of the i-th grid node (e.g. if the node i is shared by 4 elements node_wgt(i) = 4
[out]	nb_nz	number of grid nodes + number of grid nodes including repetitions + 1

Definition at line 32 of file MAKE_WGT_GRID_NODES.f90.

      
      implicit none
      
      integer*4 :: nb_node,nb_elem,nb_nz
      integer*4 :: i,j,ie
      
      integer*4, dimension(nb_node) :: node_wgt
      integer*4, dimension(nb_elem,9) :: con_matrix
 
      
      nb_nz = nb_node + 1
      node_wgt = 0
      
      do ie = 1,nb_elem
         do i = 1,8
            j = con_matrix(ie,i+1)         
            node_wgt(j) = node_wgt(j) + 1
            nb_nz = nb_nz + 1
         enddo
      enddo
      
      return