MAKE_NH_Enhanced.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_nh_enhanced()
 
        use speed_par
 
        implicit none
 
      include 'SPEED.MPI'
 
      integer*4 :: count
      integer*4, dimension(:), allocatable :: node_nhe_flag, NN_src_ind_loc
 
 
        if (mpi_id.eq.0) write(*,'(A)')
      if (mpi_id.eq.0) write(*,'(A)')'---------------Setup Not-Honoring Enhanced ---------------' 
 
 
      allocate(node_nhe_flag(nnod_loc))
      call make_nh_enhanced_initialise(nnod_loc, nmat_nhe, val_nhe, &
                             nmat, tag_mat, nelem_loc, con_nnz_loc, con_spx_loc, &
                             xx_spx_loc, yy_spx_loc, zz_spx_loc, &
                             count, &
                             node_nhe_flag, mpi_id, mpi_comm, mpi_file)
 
 
      allocate(nn_src_ind_loc(count))
      call make_nh_enhanced_nnsearch(nnod_loc, count, mpi_id, mpi_np, mpi_comm, &
                                    mpi_file, nn_src_ind_loc)
 
 
      allocate(lambda_nhe(nnod_loc), mu_nhe(nnod_loc), rho_nhe(nnod_loc))
      allocate(qs_nhe_el(nelem_loc), qp_nhe_el(nelem_loc)) ! gamma_nhe_el(nelem_loc))
 
      call make_nh_enhanced_assign_prop(nnod_loc, nmat, prop_mat, sdeg_mat, &
                                           nelem_loc, con_nnz_loc, con_spx_loc, &
                             node_nhe_flag, count, nn_src_ind_loc, qs, qp, &
                             lambda_nhe, mu_nhe, rho_nhe, qs_nhe_el, qp_nhe_el, &
                             mpi_id, mpi_comm)
 
 
      deallocate(node_nhe_flag, nn_src_ind_loc)
 
      if (mpi_id.eq.0) write(*,'(A)')'Completed.' 
 
     end subroutine MAKE_NH_Enhanced