MAKE_WGT_SPX_NODES.f90

Go to the documentation of this file.

!    Copyright (C) 2012 The SPEED FOUNDATION
!    Author: Ilario Mazzieri
!
!    This file is part of SPEED.
!
!    SPEED is free software; you can redistribute it and/or modify it
!    under the terms of the GNU Affero General Public License as
!    published by the Free Software Foundation, either version 3 of the
!    License, or (at your option) any later version.
!
!    SPEED is distributed in the hope that it will be useful, but
!    WITHOUT ANY WARRANTY; without even the implied warranty of
!    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
!    Affero General Public License for more details.
!
!    You should have received a copy of the GNU Affero General Public License
!    along with SPEED.  If not, see <http://www.gnu.org/licenses/>.
 
 
      subroutine make_wgt_spx_nodes(nb_node,spx_con_nnz,spx_con,node_wgt,nb_nz)
      
      
      implicit none
      
      integer*4 :: nb_node,spx_con_nnz,nb_nz
      integer*4 :: i,j,ie,nb_elem,nn3
      
      integer*4, dimension(0:spx_con_nnz) :: spx_con
      integer*4, dimension(nb_node) :: node_wgt
      
      nb_nz = nb_node +1
      node_wgt = 0
      
      nb_elem = spx_con(0) -1
      do ie = 1, nb_elem
         nn3 = spx_con(ie) - spx_con(ie -1) -1
         do i = 1,nn3
            j = spx_con(spx_con(ie -1) +i)
            
            node_wgt(j) = node_wgt(j) +1
            nb_nz = nb_nz +1
         enddo
      enddo
      
      return
      
      end subroutine make_wgt_spx_nodes